[英]Extracting coefficients from frailty survival model with sparse=T
我在R中使用时间依赖协变量拟合生存模型,并对某些变量使用frailty.gaussian()。
一个示例电话是
coxph1 <- coxph(Surv(tstart,tstop,dstatus)~x+frailty.gaussian(y,sparse=T),data=data)
其中y是具有多个级别的因子变量。 我想从y中提取随机效应blups(或任何估计)。 有没有办法做到这一点? 请注意,我坚持使用生存包,因为它允许拟合时间依赖的协变量。 只要可以解决这些问题,他们会很乐意转向另一个脆弱的方案。
谢谢!
这是一个可重复的例子,我可以毫无问题地提取脆弱的估计值。
library(survival)
data(rats)
coxph1 <- coxph(Surv(time,status)~rx+frailty.gaussian(litter,sparse=T),data=rats)
hist(coxph1$frail)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.