イベント履歴分析を実行するためのデータ セットを準備しました。イベント数を数えてみると、
> 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でどのようにカウントされるかについて誰かが光を当てることができますか