3

イベント履歴分析を実行するためのデータ セットを準備しました。イベント数を数えてみると、

 > event_info[, .(sum(event))]
 V1
 1: 4102

走ると、

library(eha)
reliability <- aftreg(Surv(tstart, tstop, event) ~ 1, 
                dist = "weibull", data = event_info)
summary(reliability)

私は、

Call:
aftreg(formula = Surv(tstart, tstop, event) ~ 1, data = event_info, 
    dist = "weibull")

Covariate          W.mean      Coef Time-Accn  se(Coef)    Wald p
Baseline parameters:
log(scale)                   11.975 158696.930     0.086     0.000 
log(shape)                   -0.978     0.376     0.025     0.000 
Baseline life expectancy:  

Events                    1541 
Total time at risk        68581432 
Max. log. likelihood      -16890 

イベント数が event_info データ セットの数と異なる理由がわかりません。

残念ながら、データを投稿することはできません。

イベントがehaでどのようにカウントされるかについて誰かが光を当てることができますか

4

0 に答える 0