Wraps survival::cox.zph() (scaled Schoenfeld residuals).
Arguments
- cox_result
Object returned by
morie_survival_cox().
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_schoenfeld(fit), max.level = 1)
#> List of 3
#> $ residuals: num [1:47, 1:2] 1.855 -1.034 -0.198 -0.209 -0.662 ...
#> ..- attr(*, "dimnames")=List of 2
#> $ scaled : num [1:47, 1:2] 2.526 -1.825 -0.587 -0.61 -1.295 ...
#> ..- attr(*, "dimnames")=List of 2
#> $ zph_table:'data.frame': 3 obs. of 3 variables:
