我通读了所有类似"如何计算到事件的时间/事件之间的时间"的问题。以及Surv()、coxph、comme的文档。我还需要确认我的数据设置是否正确。这表示我的数据:
library(data.table)
library(survival)
library(coxme)
snippet <- data.table(age=c(rep(74,3),rep(73,3),rep(69,3)),
sex=c(rep(0,3),rep(1,3),rep(0,3)),
year=c(0:2,0:2,0:2),
CHEM_376=c(0.036469233, 0.262701338, 0.120690616, -0.197499658, -0.231236627, -0.368920289, -0.106358638, -0.055245055, -0.141187835),
subject=c(rep(1,3),rep(2,3),rep(3,3)),
status=c(0,0,1,0,0,1,0,0,1))
我想运行cox比例风险回归，受试者ID作为随机效应(来自包装)，化学物质chem376每年测量，一些协变量如年龄和男性。下面的运行，但我不确定它是否正确:
coxph(Surv(time=year, event=status) ~ CHEM_376 + age + sex + cluster(subject), data=snippet)
包是否理解年是发生这些观察的年份，而不是状态变化的时间?我是否需要计算年数，直到每个人的状态变为1 ?我不知道该怎么做。请注意，年份不是一个正x日期时间格式。谢谢你的帮助。
当使用Surv创建Survival对象时，time参数是观察开始和事件发生时间之间的年数间隔(或天、周等)。您还可以提供time和time2参数来指示生存间隔的开始和结束时间。
参数event表示事件是否发生:1表示是(受试者死亡)，0表示否(受试者活着或被审查-从研究中删除等)。
Google "生存分析"查看相关文章