查找重叠间隔

Question

上下文：据我所知，R 缺乏一致的函数，这些函数有助于在生存/事件历史分析的上下文中进行数据准备，例如片段分割以包含时变协变量（有时称为'计算过程数据'）。

对于每个人 (id)，给出每集的开始时间 (start.cp) 和结束时间 (stop.cp)。此外，对于每个 1,2, ..., p 时变协变量 (TVC)，我们知道剧集何时开始 (tvc.start_) 以及何时开始它结束（tvc.stop_）。

在我的示例中（见下文），TVC 的数量为 2，但通常该数量可能会有所不同（从 1 到 p）。

示例：

输入数据：

  id start.cp stop.cp tvc.start1 tvc.start2 tvc.stop1 tvc.stop2
1  1        1       2          2          3         4         7
2  1        2       3          2          3         4         7
3  1        3       4          2          3         4         7
4  1        4       7          2          3         4         7
5  1        7      12          2          3         4         7

structure(list(id = c(1, 1, 1, 1, 1), start.cp = c(1, 2, 3, 4, 
7), stop.cp = c(2, 3, 4, 7, 12), tvc.start1 = c(2, 2, 2, 2, 2
), tvc.start2 = c(3, 3, 3, 3, 3), tvc.stop1 = c(4, 4, 4, 4, 4
), tvc.stop2 = c(7, 7, 7, 7, 7)), .Names = c("id", "start.cp", 
"stop.cp", "tvc.start1", "tvc.start2", "tvc.stop1", "tvc.stop2"), 
row.names = c(NA, 5L), class = "data.frame")

TVC 的名称已知，即在此示例中已知

tvc.start <- c("tvc.start1", "tvc.start2") 
tvc.stop <- c("tvc.stop1", "tvc.stop2")

预期结果：

  id start.cp stop.cp tvc.start1 tvc.start2 tvc.stop1 tvc.stop2 tvc.d1 tvc.d2
1  1        1       2          2          3         4         7      0      0
2  1        2       3          2          3         4         7      1      0
3  1        3       4          2          3         4         7      1      0
4  1        4       7          2          3         4         7      0      1
5  1        7      12          2          3         4         7      0      1

structure(list(id = c(1, 1, 1, 1, 1), start.cp = c(1, 2, 3, 4, 
7), stop.cp = c(2, 3, 4, 7, 12), tvc.start1 = c(2, 2, 2, 2, 2
), tvc.start2 = c(3, 3, 3, 3, 3), tvc.stop1 = c(4, 4, 4, 4, 4
), tvc.stop2 = c(7, 7, 7, 7, 7), tvc.d1 = c(0, 1, 1, 0, 0), tvc.d2 = c(0, 
0, 0, 1, 1)), .Names = c("id", "start.cp", "stop.cp", "tvc.start1", 
"tvc.start2", "tvc.stop1", "tvc.stop2", "tvc.d1", "tvc.d2"), row.names = c(NA, 
5L), class = "data.frame")

问题： 对于每个 TVC，我想创建一个新向量（tvc.d1、tvc.d2，请参阅示例），该向量指示给定的剧集（定义通过 start.cp和 stop.cp) 与 TVC 的间隔重叠 (=1)。假设[start.cp, stop.cp)。如何在不循环 TVC 集的情况下完成此操作，即我正在寻找矢量化解决方案。

PS：请随意更改标题...

Answer 1

我认为 Terry Therneau 可能想对您的主张提出异议，推荐的生存包中的 tcut 函数和 pyears 早在 他与 Cindy Crowson 合作撰写的关于处理时间相关协变量的技术文章。我很难理解为什么 tcv.d1 应该在间隔 2 -> 2 期间贡献曝光。 3 当其停止时间为 2 时？但对后来读者的解释是在问题的评论中。

您实际上只需要 start.cp stop.cp 向量和第一行作为输入数据。您将间隔定义向量与每个组件/个体的开始和停止向量的向量进行比较，并找到 == '1' 的间隔。我想知道数据是否真的不是以这种方式出现的，您可能不需要在设置中重复开始和停止时间。

tvec <- with(dat, c(start.cp[1], stop.cp))
dat$tvc.d1 <- 1*( findInterval(tvec,      # the "1*" converts to numeric
                               as.numeric( dat[ 1, c("tvc.start1", "tvc.stop1")]) ,  
                               all.inside=FALSE)[1:5] == 1)
dat$tvc.d2 <- 1*( findInterval(tvec, 
                               as.numeric( dat[ 1, c("tvc.start2", "tvc.stop2")]) ,  
                               all.inside=FALSE)[1:5] == 1)