Table for complex design cox model.
svycox.display(svycoxph.obj, decimal = 2)
| svycoxph.obj | svycoxph.object |
|---|---|
| decimal | digit, Default: 2 |
List including table, metric, caption
DETAILS
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')" #>