R 当绘制 svycoxph 对象的预测结果时，如何辨别哪条曲线是哪条

发布于07月20日

我使用调查包中的svycoxph创建了一个基于复杂调查设计的对象:

object <- svycoxph(Surv(time, status) ~ group, design=design)

然后，我按照文档here on page 88中的说明使用Forecate()来预测曲线并绘制它们:

prediction <- predict(object, type="curve", newdata = data.frame(group=c("A", "B")) )

plot(prediction[[1]],ci=FALSE,col="sienna") 
lines(prediction[[2]],ci=FALSE,col="royalblue")

我的组变量是A级和B级的因子.列表1是A，列表2是B吗？我以为情况就是这样，但结果与我的预期相反.当我在源代码窗格中打开该对象时，列表1和2之外没有任何标签(任何地方都没有标记A或B).

有没有办法让曲线发生反转？结果是否可能按字母顺序而不是按因素顺序列出？

推荐答案

这是我的解决方案，我没有因为只有一个关卡而出错.我在调查中使用了这个例子，刚刚修改了预测语句.

library(survey)
library(survival)
 
 data(pbc, package="survival")
 pbc$randomized<-with(pbc, !is.na(trt) & trt>0)
 biasmodel<-glm(randomized~age*edema,data=pbc,family=binomial)
 pbc$randprob<-fitted(biasmodel)
 if (is.null(pbc$albumin)) pbc$albumin<-pbc$alb ##pre2.9.0
 dpbc<-svydesign(id=~1, prob=~randprob, strata=~edema, data=subset(pbc,randomized))
 rpbc<-as.svrepdesign(dpbc)
 
 (model<-svycoxph(Surv(time,status>0)~log(bili)+protime+albumin,design=dpbc))
 svycoxph(Surv(time,status>0)~log(bili)+protime+albumin,design=rpbc)
 s1<-predict(model,se=TRUE, type="curve",
            newdata=data.frame(bili=c(3), protime=c(10), albumin=c(3.5)))
 s2<-predict(model,se=TRUE, type="curve",
             newdata=data.frame(bili=c(9), protime=c(10), albumin=c(3.5)))