從稀疏生存模型提取係數與稀疏= 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)