我准备了一个数据集来执行事件历史分析。当我计算事件的数量时,
> 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 中计算事件吗