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.55)" #> 3 "9.78 (1.37,69.94)" "0.023" "17.03 (1.69,171.6)" #> 4 "22.89 (3.2,163.48)" "0.002" "22.56 (2.25,226.42)" #> 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 1480.29 NA NA NA #> #> $caption #> [1] "Survey cox model on time ('time') to event ('status > 0')" #>