Nelson-Aalen cumulative-hazard estimator.

Usage

morie_survival_nelsonaalen(time, event, confidence = 0.95)

Arguments

time: Numeric vector of observation times.
event: 0/1 event indicator (1 = event observed).
confidence: Confidence level (default 0.95).

Value

list with times, cumhaz, ci_lower, ci_upper, at_risk, events, censored.

Examples

set.seed(1)
time <- rexp(60); event <- rbinom(60, 1, 0.7)
str(morie_survival_nelsonaalen(time, event), max.level = 1)
#> List of 8
#>  $ times   : num [1:47] 0.0594 0.1061 0.1398 0.1457 0.147 ...
#>  $ cumhaz  : num [1:47] 0.0169 0.0342 0.0517 0.0696 0.0878 ...
#>  $ ci_lower: num [1:47] 0 0 0 0.00138 0.01081 ...
#>  $ ci_upper: num [1:47] 0.0502 0.0816 0.1103 0.1378 0.1647 ...
#>  $ at_risk : num [1:47] 59 58 57 56 55 54 53 52 51 50 ...
#>  $ events  : num [1:47] 1 1 1 1 1 1 1 1 1 1 ...
#>  $ censored: num [1:47] 0 0 0 0 0 0 0 0 0 0 ...
#>  $ method  : chr "Nelson-Aalen"