我正在尝试使用R survival
包来拟合带有间隔检查和左截断的模型。
在流行病学分析中,开始时间可能无关紧要(填写表格本身不是事件),因此可以根据年龄将人们包括在研究中。这样,x
轴不再是time
,而是age
。由于人们仅从进入年龄开始就进行随访,因此数据为left-truncated(在进入年龄之前我们没有关于其特征的信息)。您可以找到有关here的更多信息。
我的事件由问卷收集,因此事件可以在问卷n
和n-1
之间的任何地方发生。因此,它们是interval-censored。解决此问题的常用方法是考虑在问卷n
和n-1
之间的中点时间对事件进行右删失,但是我想知道是否有更好的方法。
AFAIK,可以使用以下结构实现右删截的左截断(按年龄)生存分析:
db$event = ifelse(is.event, 0, 1)
Surv(age_at_inclusion, age_at_event_or_censoring, event, type="counting")
#here "counting" could be omitted as it is the default when time2 is used
对于纵向数据集Surv
:heart
,这与Surv(heart$start, heart$stop, heart$event)
帮助示例使用相同的逻辑,但是在横向数据集上。
但是阅读Surv()
help note看来,您可以通过使用以下结构来实现间隔检查:
db$event = ifelse(is.event, 0, 3)
Surv(time_interval_left, time_interval_right, event, type="interval")
由于它们似乎使用相同的结构,因此我该如何用左截断和右删截对事件进行建模?