Question

我在R中使用时间依赖协变量拟合生存模型，并对某些变量使用frailty.gaussian（）。

示例电话是

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

其中y是具有多个级别的因子变量。我想从y中提取随机效应blups（或任何估计）。有没有办法做到这一点？请注意，我坚持使用生存包，因为它允许拟合时间依赖的协变量。只要可以处理这些问题，他们会很乐意转向另一个脆弱的方案。

谢谢！

Answer 1

这是一个可重复的例子，我可以毫无问题地提取脆弱的估计值。

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