我正在努力研究R中的coxme包.我想使用像survfit()这样的函数 - 通常用于coxph()模型的方法 - 绘制调整后的生存曲线并找出不同参数值下的中位生存期.
如果我使用没有随机效果的coxph来拟合模型,我可以执行以下操作:
library(KMsurv)
data(burn)
my.surv <- with(burn, Surv(T1, D1))
cox_nr = coxph(my.surv ~ Z1 , data = burn)
survfit(cox_nr, newdata = data.frame(Z1 =1))
这提供了生存估计.但如果我与coxme适合相同的模型:
library(coxme)
cox_r = coxme(my.surv ~ Z1 + (1|Z11), data = burn)
survfit(cox_r, newdata = data.frame(Z1 = 1))
UseMethod中的错误("幸存",公式):没有适用于"幸运"的方法应用于类"coxme"的对象
所以survfit.coxme似乎不存在,从阅读coxme包文档,我没有看到一个等价物.我试图做的事情是否存在根本错误？如果没有,我怎样才能获得这些估算？
r ×1