问题
我有四列(ID)的数据,测试完成日期(TestDT) --一列具有兴趣测试(Test1)的数值结果,另一列有不同兴趣测试(Test2)的数值结果。
ID TestDT Test1 Test2
1 2011-03-02 23 NA
2 2011-03-14 NA 16
2 2011-03-15 NA 52
2 2011-11-22 24 NA
2 2011-12-12 NA 77
3 2009-03-02 NA 23
3 2009-05-02 15 NA
3 2011-15-02 NA 66
4 2017-07-03 NA 22期望结果
我想得到任何Test2的平均结果,这发生在3个月的时间内(即之前或之后)测试1的每个ID。
ID TestDT Test1 Test2 Av_of_test2_within_range
1 2011-03-02 23 NA 34
2 2011-11-22 24 NA 77
3 2009-05-02 15 NA 23我很难过滤在时间范围内发生的Test2结果。
尝试
我尝试从filter_time包中使用tibbletime,如下所示:
library(tibbletime)
FB <- as_tbl_time(myData, index = TestDT)
FB %>% group_by(ID) %>%filter_time(TestDT ~ TestDT+84)但是得到错误:
Error: Problem with `filter()` input `..1`.
x object 'TestDT' not found
i Input `..1` is `{ ... }`.
The error occured in group 1: 发布于 2021-03-06 13:34:12
data.table具有foverlaps函数,它在日期范围内合并两个数据集。
您需要将数据拆分为test1和test2,并执行如下操作:
library(data.table)
df <- read.table(text = "ID TestDT Test1 Test2
1 2011-03-02 23 NA
2 2011-03-14 NA 16
2 2011-03-15 NA 52
2 2011-11-22 24 NA
2 2011-12-12 NA 77
3 2009-03-02 NA 23
3 2009-05-02 15 NA
3 2011-12-02 NA 66
4 2017-07-03 NA 22", header = TRUE)
dt <- data.table(df)
dt[, TestDT := as.Date(TestDT)]
test1 <- dt[!is.na(Test1), .(ID, TestDT, Test1)]
test2 <- dt[!is.na(Test2), .(ID, TestDT, Test2)]
test1[, start.date := TestDT - 91]
test1[, end.date := TestDT + 91]
test2[, start.date := TestDT]
test2[, end.date := TestDT]
setkey(test2, ID, start.date, end.date)
res <- foverlaps(
test1,
test2,
by.x = c("ID", "start.date", "end.date"),
by.y = c("ID", "start.date", "end.date")
)https://stackoverflow.com/questions/66505421
复制相似问题