是否可以使用survival包从R中的另一个生物对象中减去一个生存对象?
我想绘制一个数字,显示一条生存曲线在哪个区间内高于/低于另一条生存曲线和多少。
答案 0 :(得分:0)
survA
和survB
作为survdiff
的一种可能解决方案 - 对象:
interval <- 0:2500
# choose a different time interval if you want
sumA <- summary(survA, time = interval)
sumB <- summary(survB, time = interval)
both <- data.frame(time = interval, A = sumA$surv, B = sumB$surv)
both$diff <- both$B - both$A
# or both$diff <- both$A - both$B
plot(x = both$time, y = both$diff, type = "line")