Delegates to survival::survfit() with Surv(left, right, type = "interval2").
Hand-rolled EM is left as a stub for environments without survival.
Usage
morie_survival_turnbull(left, right, max_iter = 200, tol = 1e-06)
Arguments
- left
Left-bracket times for interval-censored data
(morie_survival_turnbull).
- right
Right-bracket times for interval-censored data.
- max_iter
Iteration cap for the Turnbull NPMLE EM loop.
- tol
Convergence tolerance for the Turnbull EM.
Value
A named list with the NPMLE times, the survival
function estimates, and method = "Turnbull NPMLE".
Examples
set.seed(1)
left <- rexp(40); right <- left + runif(40, 0.1, 1)
str(morie_survival_turnbull(left, right), max.level = 1)
#> List of 3
#> $ times : num [1:19] 0.378 0.482 0.602 0.647 0.665 ...
#> $ survival: num [1:19] 0.861 0.861 0.73 0.73 0.728 ...
#> $ method : chr "Turnbull NPMLE"
