# \donttest{
library(Hmisc)
my_summary <- summary(Species ~ ., data = iris, method = "reverse")
tidy_summary(my_summary)
# }
# \donttest{
library(rms)
options(datadist = "dd")
n <- 1000L
set.seed(731L)
age <- 50L + 12L * rnorm(n)
sex <- factor(sample(c("Male", "Female"), n,
rep = TRUE,
prob = c(.6, .4)
))
cens <- 15L * runif(n)
h <- .02 * exp(.04 * (age - 50L) + .8 * (sex == "Female"))
dt <- -log(runif(n)) / h
e <- ifelse(dt <= cens, 1L, 0L)
dt <- pmin(dt, cens)
dd <- datadist(age, sex)
S <- survival::Surv(dt, e)
f <- rms::cph(S ~ age + sex)
my_summary <- summary(f)
tidy_summary(my_summary)
# }
Run the code above in your browser using DataLab