response: (chd1 %in% c('1','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=3010 (51 observations deleted due to missingness)
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.1899 0.827 0.06662 0.08532 -2.226
wz[notmissing, selectedcovariates]1 0.1972 1.218 0.02532 0.04508 4.375
wz[notmissing, selectedcovariates]2 -0.9108 0.402 0.11440 0.14977 -6.081
wz[notmissing, selectedcovariates]3 0.0112 1.011 0.00791 0.01019 1.101
wz[notmissing, selectedcovariates]4 0.0142 1.014 0.00230 0.00307 4.620
wz[notmissing, selectedcovariates]5 0.5410 1.718 0.07760 0.11036 4.902
wz[notmissing, selectedcovariates]6 0.7674 2.154 0.10173 0.15645 4.905
wruaco[notmissing]FIN-EAS_03 0.8159 2.261 0.12999 0.17243 4.732
wruaco[notmissing]FIN-EAS_24 0.3216 1.379 0.11928 0.15707 2.047
wruaco[notmissing]FIN-WES_03 0.2444 1.277 0.14406 0.18495 1.321
wruaco[notmissing]FIN-WES_24 0.1735 1.189 0.13040 0.19026 0.912
wruaco[notmissing]FRA-LIL_21 -0.0135 0.987 0.22363 0.34597 -0.039
wruaco[notmissing]FRA-STR_21 0.0405 1.041 0.20104 0.27751 0.146
wruaco[notmissing]FRA-TOU_21 0.5112 1.667 0.20601 0.27819 1.838
wruaco[notmissing]SWE-NSW_02 0.4089 1.505 0.21357 0.28981 1.411
wruaco[notmissing]SWE-NSW_03 -0.1086 0.897 0.29735 0.36223 -0.300
wruaco[notmissing]UNK-BEL_21 0.5973 1.817 0.14640 0.20048 2.979
p
allele[notmissing] 2.6e-02
wz[notmissing, selectedcovariates]1 1.2e-05
wz[notmissing, selectedcovariates]2 1.2e-09
wz[notmissing, selectedcovariates]3 2.7e-01
wz[notmissing, selectedcovariates]4 3.8e-06
wz[notmissing, selectedcovariates]5 9.5e-07
wz[notmissing, selectedcovariates]6 9.3e-07
wruaco[notmissing]FIN-EAS_03 2.2e-06
wruaco[notmissing]FIN-EAS_24 4.1e-02
wruaco[notmissing]FIN-WES_03 1.9e-01
wruaco[notmissing]FIN-WES_24 3.6e-01
wruaco[notmissing]FRA-LIL_21 9.7e-01
wruaco[notmissing]FRA-STR_21 8.8e-01
wruaco[notmissing]FRA-TOU_21 6.6e-02
wruaco[notmissing]SWE-NSW_02 1.6e-01
wruaco[notmissing]SWE-NSW_03 7.6e-01
wruaco[notmissing]UNK-BEL_21 2.9e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.827 1.209 0.700 0.978
wz[notmissing, selectedcovariates]1 1.218 0.821 1.115 1.331
wz[notmissing, selectedcovariates]2 0.402 2.486 0.300 0.539
wz[notmissing, selectedcovariates]3 1.011 0.989 0.991 1.032
wz[notmissing, selectedcovariates]4 1.014 0.986 1.008 1.020
wz[notmissing, selectedcovariates]5 1.718 0.582 1.384 2.133
wz[notmissing, selectedcovariates]6 2.154 0.464 1.585 2.927
wruaco[notmissing]FIN-EAS_03 2.261 0.442 1.613 3.170
wruaco[notmissing]FIN-EAS_24 1.379 0.725 1.014 1.877
wruaco[notmissing]FIN-WES_03 1.277 0.783 0.889 1.835
wruaco[notmissing]FIN-WES_24 1.189 0.841 0.819 1.727
wruaco[notmissing]FRA-LIL_21 0.987 1.014 0.501 1.944
wruaco[notmissing]FRA-STR_21 1.041 0.960 0.604 1.794
wruaco[notmissing]FRA-TOU_21 1.667 0.600 0.967 2.876
wruaco[notmissing]SWE-NSW_02 1.505 0.664 0.853 2.656
wruaco[notmissing]SWE-NSW_03 0.897 1.115 0.441 1.825
wruaco[notmissing]UNK-BEL_21 1.817 0.550 1.227 2.692
Rsquare= 0.116 (max possible= 0.995 )
Likelihood ratio test= 372 on 17 df, p=0
Wald test = 174 on 17 df, p=0
Score (logrank) test = 401 on 17 df, p=0, Robust = 195 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