response: (chd2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MTHFD1L_MG_449
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=2358 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.00874 1.009 0.2484 0.29050
wz[notmissing, selectedcovariates]1 0.10034 1.106 0.0477 0.05260
wz[notmissing, selectedcovariates]2 -0.74192 0.476 0.2001 0.23606
wz[notmissing, selectedcovariates]3 0.00497 1.005 0.0140 0.01665
wz[notmissing, selectedcovariates]4 0.01630 1.016 0.0041 0.00453
wz[notmissing, selectedcovariates]5 0.77169 2.163 0.1278 0.15030
wz[notmissing, selectedcovariates]6 0.87719 2.404 0.1742 0.20965
wrua2[notmissing]FRA -1.28694 0.276 0.4037 0.43373
wrua2[notmissing]SWE -0.46693 0.627 0.2707 0.29284
wrua2[notmissing]UNK -0.16631 0.847 0.3225 0.35871
z p
allele[notmissing] 0.0301 9.8e-01
wz[notmissing, selectedcovariates]1 1.9077 5.6e-02
wz[notmissing, selectedcovariates]2 -3.1429 1.7e-03
wz[notmissing, selectedcovariates]3 0.2987 7.7e-01
wz[notmissing, selectedcovariates]4 3.6025 3.2e-04
wz[notmissing, selectedcovariates]5 5.1343 2.8e-07
wz[notmissing, selectedcovariates]6 4.1840 2.9e-05
wrua2[notmissing]FRA -2.9671 3.0e-03
wrua2[notmissing]SWE -1.5945 1.1e-01
wrua2[notmissing]UNK -0.4636 6.4e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.009 0.991 0.571 1.783
wz[notmissing, selectedcovariates]1 1.106 0.905 0.997 1.226
wz[notmissing, selectedcovariates]2 0.476 2.100 0.300 0.756
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.163 0.462 1.611 2.905
wz[notmissing, selectedcovariates]6 2.404 0.416 1.594 3.626
wrua2[notmissing]FRA 0.276 3.622 0.118 0.646
wrua2[notmissing]SWE 0.627 1.595 0.353 1.113
wrua2[notmissing]UNK 0.847 1.181 0.419 1.710
Rsquare= 0.052 (max possible= 0.875 )
Likelihood ratio test= 126 on 10 df, p=0
Wald test = 91.5 on 10 df, p=2.66e-15
Score (logrank) test = 126 on 10 df, p=0, Robust = 80 p=4.86e-13
(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