Cox-Snell residuals from a fitted morie Cox model.
Examples
set.seed(1)
df <- data.frame(t = rexp(60), e = rbinom(60, 1, 0.7),
x1 = rnorm(60), x2 = rnorm(60))
fit <- morie_survival_cox(df, "t", "e", c("x1", "x2"))
str(morie_survival_coxsnell(fit), max.level = 1)
#> num [1:60] 0.4434 1.1938 0.086 0.0616 0.2749 ...
