response: (chd1=='2')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MIA3_MG_454
Call:
coxph(formula = y[notmissing] ~ allele[notmissing] + wz[notmissing,
selectedcovariates] + wruaco[notmissing] + strata(wsex[notmissing]) +
cluster(wdata$id[notmissing]), weights = weight[notmissing],
robust = TRUE)
n=2728 (41 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.20835 0.812 0.07990 0.09963
wz[notmissing, selectedcovariates]1 0.23322 1.263 0.02884 0.05282
wz[notmissing, selectedcovariates]2 -0.99349 0.370 0.13548 0.16914
wz[notmissing, selectedcovariates]3 0.00931 1.009 0.00945 0.01127
wz[notmissing, selectedcovariates]4 0.01402 1.014 0.00269 0.00348
wz[notmissing, selectedcovariates]5 0.45214 1.572 0.09015 0.12443
wz[notmissing, selectedcovariates]6 0.77402 2.168 0.11919 0.17519
wruaco[notmissing]FIN-EAS_03 1.16362 3.202 0.14927 0.18835
wruaco[notmissing]FIN-EAS_24 0.71861 2.052 0.13887 0.17378
wruaco[notmissing]FIN-WES_03 0.50193 1.652 0.16899 0.20569
wruaco[notmissing]FIN-WES_24 0.45724 1.580 0.15627 0.22089
wruaco[notmissing]FRA-LIL_21 0.49941 1.648 0.23339 0.35199
wruaco[notmissing]FRA-STR_21 0.42473 1.529 0.21986 0.29530
wruaco[notmissing]FRA-TOU_21 0.91191 2.489 0.22530 0.29856
wruaco[notmissing]SWE-NSW_02 0.72888 2.073 0.23932 0.31578
wruaco[notmissing]SWE-NSW_03 0.21038 1.234 0.32785 0.39338
wruaco[notmissing]UNK-BEL_21 0.92583 2.524 0.16690 0.21888
z p
allele[notmissing] -2.091 3.7e-02
wz[notmissing, selectedcovariates]1 4.415 1.0e-05
wz[notmissing, selectedcovariates]2 -5.874 4.3e-09
wz[notmissing, selectedcovariates]3 0.826 4.1e-01
wz[notmissing, selectedcovariates]4 4.032 5.5e-05
wz[notmissing, selectedcovariates]5 3.634 2.8e-04
wz[notmissing, selectedcovariates]6 4.418 1.0e-05
wruaco[notmissing]FIN-EAS_03 6.178 6.5e-10
wruaco[notmissing]FIN-EAS_24 4.135 3.5e-05
wruaco[notmissing]FIN-WES_03 2.440 1.5e-02
wruaco[notmissing]FIN-WES_24 2.070 3.8e-02
wruaco[notmissing]FRA-LIL_21 1.419 1.6e-01
wruaco[notmissing]FRA-STR_21 1.438 1.5e-01
wruaco[notmissing]FRA-TOU_21 3.054 2.3e-03
wruaco[notmissing]SWE-NSW_02 2.308 2.1e-02
wruaco[notmissing]SWE-NSW_03 0.535 5.9e-01
wruaco[notmissing]UNK-BEL_21 4.230 2.3e-05
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.812 1.232 0.668 0.987
wz[notmissing, selectedcovariates]1 1.263 0.792 1.138 1.400
wz[notmissing, selectedcovariates]2 0.370 2.701 0.266 0.516
wz[notmissing, selectedcovariates]3 1.009 0.991 0.987 1.032
wz[notmissing, selectedcovariates]4 1.014 0.986 1.007 1.021
wz[notmissing, selectedcovariates]5 1.572 0.636 1.232 2.006
wz[notmissing, selectedcovariates]6 2.168 0.461 1.538 3.057
wruaco[notmissing]FIN-EAS_03 3.202 0.312 2.213 4.631
wruaco[notmissing]FIN-EAS_24 2.052 0.487 1.459 2.884
wruaco[notmissing]FIN-WES_03 1.652 0.605 1.104 2.472
wruaco[notmissing]FIN-WES_24 1.580 0.633 1.025 2.436
wruaco[notmissing]FRA-LIL_21 1.648 0.607 0.827 3.285
wruaco[notmissing]FRA-STR_21 1.529 0.654 0.857 2.728
wruaco[notmissing]FRA-TOU_21 2.489 0.402 1.386 4.469
wruaco[notmissing]SWE-NSW_02 2.073 0.482 1.116 3.849
wruaco[notmissing]SWE-NSW_03 1.234 0.810 0.571 2.668
wruaco[notmissing]UNK-BEL_21 2.524 0.396 1.644 3.876
Rsquare= 0.111 (max possible= 0.986 )
Likelihood ratio test= 321 on 17 df, p=0
Wald test = 156 on 17 df, p=0
Score (logrank) test = 351 on 17 df, p=0, Robust = 178 p=0
(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