从稀疏生存模型提取系数与稀疏= T

Question

我适合与时间依赖协变量和使用frailty.gaussian（）一些变量的生存模型。

一个例子是呼叫

coxph1 <- coxph(Surv(tstart,tstop,dstatus)~x+frailty.gaussian(y,sparse=T),data=data) 

其中，y是与许多水平的因子变量。我想从y中提取随机效果蓝图（或任何估计值）。有没有一个好的方法来做到这一点？请注意，我坚持使用生存包，因为它允许拟合时间相关的协变量。只要能够处理这些问题，我们会很乐意转向另一个脆弱的方案。

谢谢！

Answer 1

这里是一个重复的例子，在那里我有没有问题，提取的脆弱估计..

library(survival) 
    data(rats) 
    coxph1 <- coxph(Surv(time,status)~rx+frailty.gaussian(litter,sparse=T),data=rats) 
    hist(coxph1$frail)