response: (chd1 %in% c('1','2'))
covariates: dsmoker
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= 3083
coef exp(coef) se(coef) robust se
allele[notmissing] 0.00279 1.003 0.0546 0.0683
wz[notmissing, selectedcovariates] 0.49898 1.647 0.0755 0.1012
wruaco[notmissing]FIN-EAS_03 0.73278 2.081 0.1292 0.1622
wruaco[notmissing]FIN-EAS_24 0.32040 1.378 0.1144 0.1451
wruaco[notmissing]FIN-WES_03 0.24810 1.282 0.1419 0.1761
wruaco[notmissing]FIN-WES_24 0.17375 1.190 0.1255 0.1622
wruaco[notmissing]FRA-LIL_21 -0.22011 0.802 0.2207 0.2959
wruaco[notmissing]FRA-STR_21 -0.04926 0.952 0.1949 0.2476
wruaco[notmissing]FRA-TOU_21 0.01092 1.011 0.1991 0.2596
wruaco[notmissing]SWE-NSW_02 0.47998 1.616 0.2052 0.2550
wruaco[notmissing]SWE-NSW_03 -0.20407 0.815 0.3048 0.3486
wruaco[notmissing]UNK-BEL_21 0.56897 1.766 0.1426 0.1849
z p
allele[notmissing] 0.0409 9.7e-01
wz[notmissing, selectedcovariates] 4.9300 8.2e-07
wruaco[notmissing]FIN-EAS_03 4.5179 6.2e-06
wruaco[notmissing]FIN-EAS_24 2.2080 2.7e-02
wruaco[notmissing]FIN-WES_03 1.4087 1.6e-01
wruaco[notmissing]FIN-WES_24 1.0710 2.8e-01
wruaco[notmissing]FRA-LIL_21 -0.7438 4.6e-01
wruaco[notmissing]FRA-STR_21 -0.1990 8.4e-01
wruaco[notmissing]FRA-TOU_21 0.0421 9.7e-01
wruaco[notmissing]SWE-NSW_02 1.8821 6.0e-02
wruaco[notmissing]SWE-NSW_03 -0.5853 5.6e-01
wruaco[notmissing]UNK-BEL_21 3.0766 2.1e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.003 0.997 0.877 1.15
wz[notmissing, selectedcovariates] 1.647 0.607 1.351 2.01
wruaco[notmissing]FIN-EAS_03 2.081 0.481 1.514 2.86
wruaco[notmissing]FIN-EAS_24 1.378 0.726 1.037 1.83
wruaco[notmissing]FIN-WES_03 1.282 0.780 0.907 1.81
wruaco[notmissing]FIN-WES_24 1.190 0.841 0.866 1.64
wruaco[notmissing]FRA-LIL_21 0.802 1.246 0.449 1.43
wruaco[notmissing]FRA-STR_21 0.952 1.050 0.586 1.55
wruaco[notmissing]FRA-TOU_21 1.011 0.989 0.608 1.68
wruaco[notmissing]SWE-NSW_02 1.616 0.619 0.980 2.66
wruaco[notmissing]SWE-NSW_03 0.815 1.226 0.412 1.61
wruaco[notmissing]UNK-BEL_21 1.766 0.566 1.229 2.54
Rsquare= 0.029 (max possible= 0.995 )
Likelihood ratio test= 90.3 on 12 df, p=4.3e-14
Wald test = 52.4 on 12 df, p=5.32e-07
Score (logrank) test = 91.8 on 12 df, p=2.19e-14, Robust = 49 p=2.11e-06
(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