Skip to contents

Table for complex design cox model.

Usage

svycox.display(svycoxph.obj, decimal = 2, pcut.univariate = NULL)

Arguments

svycoxph.obj

svycoxph.object

decimal

digit, Default: 2

pcut.univariate

pcut.univariate, Default: NULL

Value

List including table, metric, caption

Details

DETAILS

See also

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')"
#>