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