疎な生存モデルからの係数の抽出= T

Question

私は、時間依存共変量でRに生存モデルをフィッティングし、いくつかの変数にfrailty.gaussian（）を使用しています。

例コールは、Yは、多くのレベルの要因変数である

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

あります。私はyからランダムな効果のブラス（または任意の推定値）を抽出したいと思います。これを行う良い方法はありますか？時間依存の共変量を当てはめることができるので、私は生存パッケージに固執していることに注意してください。それがそれらを扱うことができる限り、別の虚弱パッケージに移動して幸せでしょうか。

ありがとうございます！ここで

Answer 1

は

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