response: (chd1=='2')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MTHFD1L_MG_449
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=2744 (41 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.03750 0.963 0.17353 0.21667
wz[notmissing, selectedcovariates]1 0.23204 1.261 0.02892 0.05188
wz[notmissing, selectedcovariates]2 -0.96587 0.381 0.13522 0.16985
wz[notmissing, selectedcovariates]3 0.00941 1.009 0.00947 0.01134
wz[notmissing, selectedcovariates]4 0.01450 1.015 0.00268 0.00346
wz[notmissing, selectedcovariates]5 0.45665 1.579 0.08973 0.12349
wz[notmissing, selectedcovariates]6 0.78457 2.191 0.11919 0.17441
wruaco[notmissing]FIN-EAS_03 1.10250 3.012 0.15000 0.18892
wruaco[notmissing]FIN-EAS_24 0.72323 2.061 0.13897 0.17388
wruaco[notmissing]FIN-WES_03 0.54992 1.733 0.16775 0.20434
wruaco[notmissing]FIN-WES_24 0.47683 1.611 0.15611 0.21945
wruaco[notmissing]FRA-LIL_21 0.48840 1.630 0.23346 0.34802
wruaco[notmissing]FRA-STR_21 0.42259 1.526 0.21726 0.29160
wruaco[notmissing]FRA-TOU_21 0.92547 2.523 0.22513 0.29821
wruaco[notmissing]SWE-NSW_02 0.73175 2.079 0.23850 0.31491
wruaco[notmissing]SWE-NSW_03 0.22594 1.253 0.32759 0.39299
wruaco[notmissing]UNK-BEL_21 0.95707 2.604 0.16637 0.21736
z p
allele[notmissing] -0.173 8.6e-01
wz[notmissing, selectedcovariates]1 4.472 7.7e-06
wz[notmissing, selectedcovariates]2 -5.687 1.3e-08
wz[notmissing, selectedcovariates]3 0.830 4.1e-01
wz[notmissing, selectedcovariates]4 4.192 2.8e-05
wz[notmissing, selectedcovariates]5 3.698 2.2e-04
wz[notmissing, selectedcovariates]6 4.498 6.9e-06
wruaco[notmissing]FIN-EAS_03 5.836 5.4e-09
wruaco[notmissing]FIN-EAS_24 4.159 3.2e-05
wruaco[notmissing]FIN-WES_03 2.691 7.1e-03
wruaco[notmissing]FIN-WES_24 2.173 3.0e-02
wruaco[notmissing]FRA-LIL_21 1.403 1.6e-01
wruaco[notmissing]FRA-STR_21 1.449 1.5e-01
wruaco[notmissing]FRA-TOU_21 3.103 1.9e-03
wruaco[notmissing]SWE-NSW_02 2.324 2.0e-02
wruaco[notmissing]SWE-NSW_03 0.575 5.7e-01
wruaco[notmissing]UNK-BEL_21 4.403 1.1e-05
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.963 1.038 0.630 1.473
wz[notmissing, selectedcovariates]1 1.261 0.793 1.139 1.396
wz[notmissing, selectedcovariates]2 0.381 2.627 0.273 0.531
wz[notmissing, selectedcovariates]3 1.009 0.991 0.987 1.032
wz[notmissing, selectedcovariates]4 1.015 0.986 1.008 1.022
wz[notmissing, selectedcovariates]5 1.579 0.633 1.239 2.011
wz[notmissing, selectedcovariates]6 2.191 0.456 1.557 3.085
wruaco[notmissing]FIN-EAS_03 3.012 0.332 2.080 4.361
wruaco[notmissing]FIN-EAS_24 2.061 0.485 1.466 2.898
wruaco[notmissing]FIN-WES_03 1.733 0.577 1.161 2.587
wruaco[notmissing]FIN-WES_24 1.611 0.621 1.048 2.477
wruaco[notmissing]FRA-LIL_21 1.630 0.614 0.824 3.224
wruaco[notmissing]FRA-STR_21 1.526 0.655 0.862 2.702
wruaco[notmissing]FRA-TOU_21 2.523 0.396 1.406 4.526
wruaco[notmissing]SWE-NSW_02 2.079 0.481 1.121 3.853
wruaco[notmissing]SWE-NSW_03 1.253 0.798 0.580 2.708
wruaco[notmissing]UNK-BEL_21 2.604 0.384 1.701 3.987
Rsquare= 0.106 (max possible= 0.986 )
Likelihood ratio test= 309 on 17 df, p=0
Wald test = 151 on 17 df, p=0
Score (logrank) test = 338 on 17 df, p=0, Robust = 172 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