response: (chd1 %in% c('1','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=3027 (51 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.0116 0.988 0.14375 0.18988
wz[notmissing, selectedcovariates]1 0.1973 1.218 0.02536 0.04439
wz[notmissing, selectedcovariates]2 -0.8928 0.410 0.11426 0.14996
wz[notmissing, selectedcovariates]3 0.0115 1.012 0.00793 0.01020
wz[notmissing, selectedcovariates]4 0.0146 1.015 0.00229 0.00306
wz[notmissing, selectedcovariates]5 0.5447 1.724 0.07728 0.10964
wz[notmissing, selectedcovariates]6 0.7735 2.167 0.10167 0.15563
wruaco[notmissing]FIN-EAS_03 0.7601 2.139 0.13053 0.17264
wruaco[notmissing]FIN-EAS_24 0.3261 1.386 0.11933 0.15702
wruaco[notmissing]FIN-WES_03 0.2791 1.322 0.14330 0.18401
wruaco[notmissing]FIN-WES_24 0.1877 1.207 0.13030 0.18932
wruaco[notmissing]FRA-LIL_21 -0.0235 0.977 0.22372 0.34236
wruaco[notmissing]FRA-STR_21 0.0388 1.040 0.19865 0.27389
wruaco[notmissing]FRA-TOU_21 0.5254 1.691 0.20583 0.27808
wruaco[notmissing]SWE-NSW_02 0.4160 1.516 0.21278 0.28911
wruaco[notmissing]SWE-NSW_03 -0.0923 0.912 0.29712 0.36211
wruaco[notmissing]UNK-BEL_21 0.6251 1.868 0.14587 0.19920
z p
allele[notmissing] -0.0611 9.5e-01
wz[notmissing, selectedcovariates]1 4.4448 8.8e-06
wz[notmissing, selectedcovariates]2 -5.9536 2.6e-09
wz[notmissing, selectedcovariates]3 1.1311 2.6e-01
wz[notmissing, selectedcovariates]4 4.7611 1.9e-06
wz[notmissing, selectedcovariates]5 4.9686 6.7e-07
wz[notmissing, selectedcovariates]6 4.9703 6.7e-07
wruaco[notmissing]FIN-EAS_03 4.4029 1.1e-05
wruaco[notmissing]FIN-EAS_24 2.0767 3.8e-02
wruaco[notmissing]FIN-WES_03 1.5166 1.3e-01
wruaco[notmissing]FIN-WES_24 0.9916 3.2e-01
wruaco[notmissing]FRA-LIL_21 -0.0687 9.5e-01
wruaco[notmissing]FRA-STR_21 0.1418 8.9e-01
wruaco[notmissing]FRA-TOU_21 1.8893 5.9e-02
wruaco[notmissing]SWE-NSW_02 1.4391 1.5e-01
wruaco[notmissing]SWE-NSW_03 -0.2549 8.0e-01
wruaco[notmissing]UNK-BEL_21 3.1382 1.7e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.988 1.012 0.681 1.434
wz[notmissing, selectedcovariates]1 1.218 0.821 1.117 1.329
wz[notmissing, selectedcovariates]2 0.410 2.442 0.305 0.549
wz[notmissing, selectedcovariates]3 1.012 0.989 0.992 1.032
wz[notmissing, selectedcovariates]4 1.015 0.986 1.009 1.021
wz[notmissing, selectedcovariates]5 1.724 0.580 1.391 2.137
wz[notmissing, selectedcovariates]6 2.167 0.461 1.598 2.940
wruaco[notmissing]FIN-EAS_03 2.139 0.468 1.525 3.000
wruaco[notmissing]FIN-EAS_24 1.386 0.722 1.018 1.885
wruaco[notmissing]FIN-WES_03 1.322 0.756 0.922 1.896
wruaco[notmissing]FIN-WES_24 1.207 0.829 0.832 1.749
wruaco[notmissing]FRA-LIL_21 0.977 1.024 0.499 1.911
wruaco[notmissing]FRA-STR_21 1.040 0.962 0.608 1.778
wruaco[notmissing]FRA-TOU_21 1.691 0.591 0.981 2.917
wruaco[notmissing]SWE-NSW_02 1.516 0.660 0.860 2.672
wruaco[notmissing]SWE-NSW_03 0.912 1.097 0.448 1.854
wruaco[notmissing]UNK-BEL_21 1.868 0.535 1.265 2.761
Rsquare= 0.112 (max possible= 0.995 )
Likelihood ratio test= 361 on 17 df, p=0
Wald test = 171 on 17 df, p=0
Score (logrank) test = 389 on 17 df, p=0, Robust = 190 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