使用'向生存数据添加时变协变量；tmerge'；在'；生存'；包裹

我试图使用生存软件包中的

tmerge

向数据集添加几个时间相关的协变量，以进行生存分析。我的意思是按照主题中的建议按顺序添加每一项，但是第一次添加的输出并没有按照我的预期工作

更具体地说，我有一个简单的data.frame，其中包含个人（组织）的ID和组织停止活动之前的天数（年龄）。第二个data.frame包含组织经历“过渡”事件之前的ID和天数。并非所有组织都经历了转换，因此并非所有组织都出现在第二个data.frame中

在对

tmerge

I的第一次调用中，我按照包使用的格式格式化第一个data.frame。在第二个示例中，我尝试添加一个变量，用于统计组织经历的过渡次数。对于大多数组织来说，结果与我预期的一样，但对于少数组织来说，结果没有意义，我也没有明显的理由解释为什么失败

frames很小，所以我将它们与下面的代码一起发布

ages <- structure(list(id = c(1L, 2L, 5L, 6L, 9L, 10L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 20L, 21L, 24L, 26L, 27L, 28L, 29L, 30L, 31L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 42L, 45L, 46L, 43L, 48L, 49L, 50L, 51L, 52L, 54L, 55L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 8L, 19L, 22L, 23L, 33L, 41L), age = c(13668, 21550, 15249, 21550, 16045, 21550, 14976, 14976, 6574, 21550, 4463, 16927, 16927, 15706, 4567, 21306, 17235, 22158, 19692, 17632, 17597, 4383, 5811, 7704, 5063, 17351, 17015, 16801, 4383, 5080, 13185, 12604, 19784, 5310, 15369, 13239, 1638, 21323, 10914, 21262, 7297, 17214, 17508, 14199, 14062, 2227, 8434, 4593, 14429, 21323, 4782, 10813, 2667, 2853, 5709, 3140, 12237, 7882, 21550, 15553, 16466, 16621, 19534, 21842)), .Names = c("id", "age"), row.names = c(NA, 64L), class = "data.frame")
ages1 <- tmerge(ages, ages, id=id, tstop=age)
transitions <- structure(list(id = c(2L, 2L, 6L, 8L, 10L, 19L, 22L, 23L, 24L, 31L, 33L, 41L, 43L, 43L, 52L, 55L, 66L), transition = structure(c(18993, 13668, 15249, 15706, 15887, 11609, 4023, 9316, 16193, 1461, 4584, 17824, 3713, 11261, 16818, 10670, 15479), class = "difftime", units = "days")), .Names = c("id", "transition"), row.names = c(3L, 4L, 7L, 8L, 11L, 20L, 25L, 27L, 28L, 35L, 38L, 47L, 49L, 51L, 59L, 61L, 73L), class = "data.frame")
newdata <- tmerge(ages1, transitions, id=id, transition=cumtdc(transition))

---

ages这个问题通过简单的id排序解决。
ages1感谢您的回答。我真是抓狂。顺便问一下，软件包的创建者在哪里提到这个？