models.Rmd
library(mRclwhip)
df <- arsenal::mockstudy
df_selected <- df %>%
dplyr::select(fu.time, fu.stat, mdquality.s, ps, sex) %>%
dplyr::mutate_at(dplyr::vars(mdquality.s, ps), as.factor) %>%
dplyr::mutate(fu.time = fu.time / 365)
Hmisc::label(df_selected$sex) <- "Sex"
Hmisc::label(df_selected$mdquality.s) <- "Quality"
Hmisc::label(df_selected$ps) <- "PS"
modelsum_cox(event = fu.stat, time = fu.time, data = df_selected)
Level |
N |
HR (95% CI) |
P-value |
|
Quality |
1 |
1,499 |
0.65 (0.54, 0.79) |
<0.001 |
PS |
1 |
1,499 |
1.64 (1.45, 1.85) |
<0.001 |
PS |
2 |
1,499 |
2.68 (2.07, 3.47) |
<0.001 |
Sex |
Female |
1,499 |
1.00 (0.90, 1.12) |
0.98 |
modelsum_km(event = fu.stat, time = fu.time, times = 1:3, data = df_selected, freedom.from = "death")
|
Total |
Total Events |
1-year freedom from death |
2-year freedom from death |
3-year freedom from death |
P value |
All |
1,499 |
1,356 |
67% (65%, 70%) |
34% (31%, 36%) |
17% (15%, 19%) |
- |
Quality |
<0.001 |
|||||
0 |
124 |
119 |
47.6% (39.6%, 57.2%) |
18.9% (13.1%, 27.3%) |
10.7% (6.4%, 17.9%) |
- |
1 |
1,123 |
1,013 |
69.4% (66.8%, 72.2%) |
34.5% (31.8%, 37.4%) |
17.7% (15.6%, 20.1%) |
- |
PS |
<0.001 |
|||||
0 |
635 |
554 |
76.7% (73.44%, 80.03%) |
41.4% (37.77%, 45.47%) |
23.2% (20.07%, 26.70%) |
- |
1 |
531 |
502 |
57.2% (53.09%, 61.53%) |
25.1% (21.62%, 29.05%) |
10.5% (8.17%, 13.52%) |
- |
2 |
67 |
65 |
32.8% (23.31%, 46.25%) |
13.4% (7.31%, 24.67%) |
6.0% (2.31%, 15.44%) |
- |
Sex |
0.98 |
|||||
Female |
583 |
829 |
65.31% (61.55%, 69.29%) |
32.77% (29.17%, 36.82%) |
17.01% (14.20%, 20.38%) |
- |
Male |
916 |
527 |
68.70% (65.76%, 71.77%) |
34.37% (31.42%, 37.61%) |
17.48% (15.16%, 20.15%) |
- |