R-生存竞争风险COX模型

Question

我正在使用survival软件包进行竞争风险分析，并希望使用预测进行模拟。

i计划从[0，1]中采样一个均匀的随机x，从与x相交的生存曲线中找到终止时间；然后从[0，hazard1 危险2]中采样另一个随机y，以决定选择哪个端状态。

我可以提取生存曲线，但不知道对各个竞争状态有什么危害。我从下面的survival软件包中复制了Vignettes"竞争"的示例：

data(mgus2)
cfit1 <- coxph(Surv(etime, event=="pcm") ~ age + sex + mspike, mgus2)
etime <- with(mgus2, ifelse(pstat==0, futime, ptime))
event <- with(mgus2, ifelse(pstat==0, 2*death, 1))
event <- factor(event, 0:2, labels=c("censor", "pcm", "death"))
cfit2 <- coxph(Surv(etime, event=="death") ~ age + sex + mspike, mgus2)
cfit1 <- coxph(Surv(etime, event=="pcm") ~ age + sex + mspike, mgus2)
newdata <- expand.grid(sex=c("F", "M"), age=c(60, 80), mspike=1.2)
newdata
temp <- matrix(list(), 3,3)
dimnames(temp) <- list(from=c("Entry", "PCM", "Death"),
to =c("Entry", "PCM", "Death"))
temp[1,2] <- list(survfit(cfit1, newdata, std.err=FALSE))
temp[1,3] <- list(survfit(cfit2, newdata, std.err=FALSE))
csurv <- survfit(temp, p0 =c(1,0,0))

看起来csurv$pstate包含生存曲线和累积危害。但是我不明白他们的关系。我试图通过增加两个累积危害来退出生存，并如下所示，但结果与生存曲线不同。

exp(-(csurv$pstate[,2]+csurv$pstate[,3]))
csurv$pstate[,1]

他们有什么关系？我如何使用它们进行仿真？

Answer 1

文档说，幸存对象的预科成分是"状态中的流行率"。(也许Terry更改名称时忘了修复该文档？(由于这是概率，在任何给定时间，概率的总和应为1。因此您不会指出它们。当我将PSTATE组件的行列制成表格时，所有总和(每行3个值(均在1.0的数值误差范围内。因此，总和的指数将等于e，而不是您所期望的。我会丢弃exp操作。它已经为您"幕后"完成了。