[Survival Analysis]Survfit equivalent for coxme in R?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

I am struggling with the coxme package in R. I would like to use a function like survfit() - the way it would ordinarily be used for a coxph() model - to plot adjusted survival curves and find the median survival at different parameter values.

If I fit the model using coxph without random effects I can do the following:

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))

This provides survival estimates. But if I fit the same model with coxme:

library(coxme)cox_r = coxme(my.surv ~ Z1 + (1|Z11), data = burn)survfit(cox_r, newdata = data.frame(Z1 = 1))

Error in UseMethod("survfit", formula) : no applicable method for 'survfit' applied to an object of class "coxme"

So survfit.coxme doesn't seem to exist and from reading the coxme package documentation, I don't see an equivalent. Is there something fundamentally wrong about what I am attempting to do? If not, how can I get these estimates?

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：equivalent Analysis Survival Analysi alysis different following function adjusted provides

survival package