Table for complex design cox model.
Examples
library(survival)
data(pbc)
pbc$sex <- factor(pbc$sex)
pbc$stage <- factor(pbc$stage)
pbc$randomized <- with(pbc, !is.na(trt) & trt > 0)
biasmodel <- glm(randomized ~ age * edema, data = pbc, family = binomial)
pbc$randprob <- fitted(biasmodel)
if (is.null(pbc$albumin)) pbc$albumin <- pbc$alb ## pre2.9.0
dpbc <- survey::svydesign(
id = ~1, prob = ~randprob, strata = ~edema,
data = subset(pbc, randomized)
)
model <- survey::svycoxph(Surv(time, status > 0) ~ sex + protime + albumin + stage,
design = dpbc
)
svycox.display(model)
#> Stratified Independent Sampling design (with replacement)
#> survey::svydesign(id = ~1, prob = ~randprob, strata = ~edema,
#> data = subset(pbc, randomized))
#> Stratified Independent Sampling design (with replacement)
#> survey::svydesign(id = ~1, prob = ~randprob, strata = ~edema,
#> data = subset(pbc, randomized))
#> Stratified Independent Sampling design (with replacement)
#> survey::svydesign(id = ~1, prob = ~randprob, strata = ~edema,
#> data = subset(pbc, randomized))
#> Stratified Independent Sampling design (with replacement)
#> survey::svydesign(id = ~1, prob = ~randprob, strata = ~edema,
#> data = subset(pbc, randomized))
#> Stratified Independent Sampling design (with replacement)
#> survey::svydesign(id = ~1, prob = ~randprob, strata = ~edema,
#> data = subset(pbc, randomized))
#> $table
#> crude HR(95%CI) crude P value adj. HR(95%CI)
#> sex: f vs m "0.62 (0.4,0.97)" "0.038" "0.55 (0.33,0.9)"
#> protime "1.37 (1.09,1.72)" "0.006" "1.52 (1.2,1.91)"
#> albumin "0.2 (0.14,0.29)" "< 0.001" "0.31 (0.2,0.47)"
#> stage: ref.=1 NA NA NA
#> 2 "5.67 (0.77,41.78)" "0.089" "10.94 (1.01,118.54)"
#> 3 "9.78 (1.37,69.94)" "0.023" "17.03 (1.69,171.6)"
#> 4 "22.89 (3.2,163.47)" "0.002" "22.56 (2.25,226.41)"
#> adj. P value
#> sex: f vs m "0.017"
#> protime "< 0.001"
#> albumin "< 0.001"
#> stage: ref.=1 NA
#> 2 "0.049"
#> 3 "0.016"
#> 4 "0.008"
#>
#> $metric
#> [,1] [,2] [,3] [,4]
#> <NA> NA NA NA NA
#> No. of observations 312.00 NA NA NA
#> No. of events 144.00 NA NA NA
#> AIC 1483.12 NA NA NA
#>
#> $caption
#> [1] "Survey cox model on time ('time') to event ('status > 0')"
#>