地标分析,在一定时间后获得最小事件时间值

Landmark analysis and obtaining the least event time value after a certain time

提问人:Mohamed Rahouma 提问时间:11/18/2023 最后编辑:Stu SztukowskiMohamed Rahouma 更新时间:11/18/2023 访问量:11

问:

我有一个复合结果,称为包括死亡、中风、心肌梗塞和复发,其发生时间称为 。 我想做一个具有里程碑意义的分析,我看到了这个链接,但我不想排除任何患者。我想排除 90 天之前的任何事件,并考虑之后最近的时间点以创建新的事件发生时间变量。所以我需要创建和.primaryprimarydaysprimary90primarydays90

以下是我的示例数据和代码

df<-read.table(text="primary    primarydays dth dthdays stroke  strokedays  mi  midays  recurrence  recurrencedays
0   1   0   1   0   1   0   1   0   1
1   1   0   966 0   966 1   1   0   966
0   2   0   1525    0   2   0   2   0   2
0   2   0   1407    0   2   0   2   0   2
0   2   0   16  0   2   0   2   0   2
1   2   0   689 0   689 1   2   0   689
0   2   0   499 0   2   0   2   0   2
0   2   0   2   0   2   0   2   0   2
0   2   0   38  0   2   0   2   0   2
0   2   0   2   0   2   0   2   0   2
0   2   0   12  0   2   0   2   0   2
0   2   0   534 0   2   0   2   0   2
0   2   0   540 0   2   0   2   0   2
1   3   0   1659    0   1659    1   3   0   1659
1   3   0   1430    1   1430    0   1430    1   3
1   6   0   743 1   6   0   743 1   743
1   6   0   1121    1   6   1   1121    0   1121
1   7   0   833 0   833 1   7   0   833
1   95  1   95  0   36  0   36  0   36
1   95  1   101 1   95  0   101 0   101
1   95  0   753 0   753 1   95  0   753
1   96  0   879 1   96  0   879 0   879
0   96  0   128 0   96  0   96  0   96
0   97  0   1060    0   97  0   97  0   97
0   98  0   98  0   98  0   98  0   98
0   100 0   185 0   100 0   100 0   100
1   101 0   554 0   554 0   101 0   554
1   102 0   565 0   565 1   102 0   565
1   102 0   1331    0   1331    0   1331    1   102
1   102 0   1285    1   102 0   1285    0   1285", sep="\t", header=T)

#I tried to add ">90" in the code but I am sure that it should be modified:
df %>% rowwise() %>% mutate(primary90 = pmax(dth    , stroke    ,   mi, recurrence, na.rm=T) >90)->df ;df
df %>% rowwise() %>% mutate(primary90days = pmin(dthdays    , strokedays    ,   midays, recurrencedays, na.rm=T)>90 )->df ;df


我也可以尝试 sas,但仍然在为正确的代码而苦苦挣扎。

任何建议将不胜感激。

R 生存分析 Cox-回归

评论


答: 暂无答案