response: (chd1=='2')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: CELSR2_MG_438
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 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.08356 0.920 0.06581 0.08121
wz[notmissing, selectedcovariates]1 0.23906 1.270 0.02918 0.05151
wz[notmissing, selectedcovariates]2 -0.97246 0.378 0.13592 0.17046
wz[notmissing, selectedcovariates]3 0.00929 1.009 0.00951 0.01134
wz[notmissing, selectedcovariates]4 0.01343 1.014 0.00268 0.00344
wz[notmissing, selectedcovariates]5 0.46627 1.594 0.09003 0.12390
wz[notmissing, selectedcovariates]6 0.80203 2.230 0.11917 0.17398
wruaco[notmissing]FIN-EAS_03 1.09506 2.989 0.15065 0.18987
wruaco[notmissing]FIN-EAS_24 0.72710 2.069 0.14036 0.17511
wruaco[notmissing]FIN-WES_03 0.55150 1.736 0.16773 0.20456
wruaco[notmissing]FIN-WES_24 0.46182 1.587 0.15575 0.21749
wruaco[notmissing]FRA-LIL_21 0.49209 1.636 0.23339 0.34661
wruaco[notmissing]FRA-STR_21 0.40253 1.496 0.21768 0.28917
wruaco[notmissing]FRA-TOU_21 0.89637 2.451 0.22587 0.29858
wruaco[notmissing]SWE-NSW_02 0.74842 2.114 0.23882 0.31625
wruaco[notmissing]SWE-NSW_03 0.23043 1.259 0.32775 0.39265
wruaco[notmissing]UNK-BEL_21 0.97125 2.641 0.16611 0.21680
z p
allele[notmissing] -1.029 3.0e-01
wz[notmissing, selectedcovariates]1 4.641 3.5e-06
wz[notmissing, selectedcovariates]2 -5.705 1.2e-08
wz[notmissing, selectedcovariates]3 0.819 4.1e-01
wz[notmissing, selectedcovariates]4 3.899 9.7e-05
wz[notmissing, selectedcovariates]5 3.763 1.7e-04
wz[notmissing, selectedcovariates]6 4.610 4.0e-06
wruaco[notmissing]FIN-EAS_03 5.767 8.1e-09
wruaco[notmissing]FIN-EAS_24 4.152 3.3e-05
wruaco[notmissing]FIN-WES_03 2.696 7.0e-03
wruaco[notmissing]FIN-WES_24 2.123 3.4e-02
wruaco[notmissing]FRA-LIL_21 1.420 1.6e-01
wruaco[notmissing]FRA-STR_21 1.392 1.6e-01
wruaco[notmissing]FRA-TOU_21 3.002 2.7e-03
wruaco[notmissing]SWE-NSW_02 2.367 1.8e-02
wruaco[notmissing]SWE-NSW_03 0.587 5.6e-01
wruaco[notmissing]UNK-BEL_21 4.480 7.5e-06
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.920 1.087 0.784 1.079
wz[notmissing, selectedcovariates]1 1.270 0.787 1.148 1.405
wz[notmissing, selectedcovariates]2 0.378 2.644 0.271 0.528
wz[notmissing, selectedcovariates]3 1.009 0.991 0.987 1.032
wz[notmissing, selectedcovariates]4 1.014 0.987 1.007 1.020
wz[notmissing, selectedcovariates]5 1.594 0.627 1.250 2.032
wz[notmissing, selectedcovariates]6 2.230 0.448 1.586 3.136
wruaco[notmissing]FIN-EAS_03 2.989 0.335 2.060 4.337
wruaco[notmissing]FIN-EAS_24 2.069 0.483 1.468 2.916
wruaco[notmissing]FIN-WES_03 1.736 0.576 1.163 2.592
wruaco[notmissing]FIN-WES_24 1.587 0.630 1.036 2.430
wruaco[notmissing]FRA-LIL_21 1.636 0.611 0.829 3.227
wruaco[notmissing]FRA-STR_21 1.496 0.669 0.849 2.636
wruaco[notmissing]FRA-TOU_21 2.451 0.408 1.365 4.400
wruaco[notmissing]SWE-NSW_02 2.114 0.473 1.137 3.928
wruaco[notmissing]SWE-NSW_03 1.259 0.794 0.583 2.718
wruaco[notmissing]UNK-BEL_21 2.641 0.379 1.727 4.040
Rsquare= 0.107 (max possible= 0.986 )
Likelihood ratio test= 310 on 17 df, p=0
Wald test = 156 on 17 df, p=0
Score (logrank) test = 341 on 17 df, p=0, Robust = 173 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