response: (chd2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: 10Q11_21_MG_458
Call:
coxph(formula = y[notmissing] ~ allele[notmissing] + wz[notmissing,
selectedcovariates] + wrua2[notmissing] + strata(wsex[notmissing]) +
cluster(wdata$id[notmissing]), weights = weight[notmissing],
robust = TRUE)
n=2339 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.10481 1.110 0.10919 0.12629
wz[notmissing, selectedcovariates]1 0.09910 1.104 0.04757 0.05247
wz[notmissing, selectedcovariates]2 -0.75049 0.472 0.20008 0.23602
wz[notmissing, selectedcovariates]3 0.00514 1.005 0.01398 0.01655
wz[notmissing, selectedcovariates]4 0.01597 1.016 0.00412 0.00456
wz[notmissing, selectedcovariates]5 0.75846 2.135 0.12804 0.15111
wz[notmissing, selectedcovariates]6 0.89631 2.451 0.17452 0.21061
wrua2[notmissing]FRA -1.25758 0.284 0.40422 0.43480
wrua2[notmissing]SWE -0.40053 0.670 0.27097 0.29298
wrua2[notmissing]UNK -0.14016 0.869 0.32293 0.35841
z p
allele[notmissing] 0.830 4.1e-01
wz[notmissing, selectedcovariates]1 1.889 5.9e-02
wz[notmissing, selectedcovariates]2 -3.180 1.5e-03
wz[notmissing, selectedcovariates]3 0.311 7.6e-01
wz[notmissing, selectedcovariates]4 3.503 4.6e-04
wz[notmissing, selectedcovariates]5 5.019 5.2e-07
wz[notmissing, selectedcovariates]6 4.256 2.1e-05
wrua2[notmissing]FRA -2.892 3.8e-03
wrua2[notmissing]SWE -1.367 1.7e-01
wrua2[notmissing]UNK -0.391 7.0e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.110 0.900 0.867 1.422
wz[notmissing, selectedcovariates]1 1.104 0.906 0.996 1.224
wz[notmissing, selectedcovariates]2 0.472 2.118 0.297 0.750
wz[notmissing, selectedcovariates]3 1.005 0.995 0.973 1.038
wz[notmissing, selectedcovariates]4 1.016 0.984 1.007 1.025
wz[notmissing, selectedcovariates]5 2.135 0.468 1.588 2.871
wz[notmissing, selectedcovariates]6 2.451 0.408 1.622 3.703
wrua2[notmissing]FRA 0.284 3.517 0.121 0.667
wrua2[notmissing]SWE 0.670 1.493 0.377 1.190
wrua2[notmissing]UNK 0.869 1.150 0.431 1.755
Rsquare= 0.051 (max possible= 0.876 )
Likelihood ratio test= 123 on 10 df, p=0
Wald test = 89 on 10 df, p=8.33e-15
Score (logrank) test = 124 on 10 df, p=0, Robust = 78 p=1.25e-12
(Note: the likelihood ratio and score tests assume independence of
observations within a cluster, the Wald and robust score tests do not).
Updated on 2008-06-04