response: (stroke4 %in% c('2','4','6'))
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=2424 (42 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.17469 0.840 0.24773 0.28078
wz[notmissing, selectedcovariates]1 0.05945 1.061 0.04572 0.05326
wz[notmissing, selectedcovariates]2 -0.75656 0.469 0.18183 0.20569
wz[notmissing, selectedcovariates]3 0.00809 1.008 0.01246 0.01345
wz[notmissing, selectedcovariates]4 0.01532 1.015 0.00382 0.00442
wz[notmissing, selectedcovariates]5 0.39460 1.484 0.12891 0.14282
wz[notmissing, selectedcovariates]6 0.76263 2.144 0.16206 0.19966
wruaco[notmissing]FIN-EAS_03 0.35009 1.419 0.22040 0.23862
wruaco[notmissing]FIN-EAS_24 -0.10664 0.899 0.19574 0.22541
wruaco[notmissing]FIN-WES_03 -0.02132 0.979 0.22945 0.24913
wruaco[notmissing]FIN-WES_24 -0.14248 0.867 0.21082 0.24732
wruaco[notmissing]FRA-LIL_21 -0.94369 0.389 0.52108 0.57184
wruaco[notmissing]FRA-STR_21 -1.05931 0.347 0.47931 0.52225
wruaco[notmissing]FRA-TOU_21 -0.68639 0.503 0.53089 0.57343
wruaco[notmissing]SWE-NSW_02 0.03966 1.040 0.29233 0.34934
wruaco[notmissing]SWE-NSW_03 0.09819 1.103 0.26946 0.30017
wruaco[notmissing]UNK-BEL_21 -0.69926 0.497 0.37734 0.39756
z p
allele[notmissing] -0.6222 0.53000
wz[notmissing, selectedcovariates]1 1.1163 0.26000
wz[notmissing, selectedcovariates]2 -3.6781 0.00023
wz[notmissing, selectedcovariates]3 0.6011 0.55000
wz[notmissing, selectedcovariates]4 3.4695 0.00052
wz[notmissing, selectedcovariates]5 2.7629 0.00570
wz[notmissing, selectedcovariates]6 3.8197 0.00013
wruaco[notmissing]FIN-EAS_03 1.4671 0.14000
wruaco[notmissing]FIN-EAS_24 -0.4731 0.64000
wruaco[notmissing]FIN-WES_03 -0.0856 0.93000
wruaco[notmissing]FIN-WES_24 -0.5761 0.56000
wruaco[notmissing]FRA-LIL_21 -1.6503 0.09900
wruaco[notmissing]FRA-STR_21 -2.0284 0.04300
wruaco[notmissing]FRA-TOU_21 -1.1970 0.23000
wruaco[notmissing]SWE-NSW_02 0.1135 0.91000
wruaco[notmissing]SWE-NSW_03 0.3271 0.74000
wruaco[notmissing]UNK-BEL_21 -1.7589 0.07900
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.840 1.191 0.484 1.456
wz[notmissing, selectedcovariates]1 1.061 0.942 0.956 1.178
wz[notmissing, selectedcovariates]2 0.469 2.131 0.314 0.702
wz[notmissing, selectedcovariates]3 1.008 0.992 0.982 1.035
wz[notmissing, selectedcovariates]4 1.015 0.985 1.007 1.024
wz[notmissing, selectedcovariates]5 1.484 0.674 1.122 1.963
wz[notmissing, selectedcovariates]6 2.144 0.466 1.450 3.171
wruaco[notmissing]FIN-EAS_03 1.419 0.705 0.889 2.265
wruaco[notmissing]FIN-EAS_24 0.899 1.113 0.578 1.398
wruaco[notmissing]FIN-WES_03 0.979 1.022 0.601 1.595
wruaco[notmissing]FIN-WES_24 0.867 1.153 0.534 1.408
wruaco[notmissing]FRA-LIL_21 0.389 2.569 0.127 1.194
wruaco[notmissing]FRA-STR_21 0.347 2.884 0.125 0.965
wruaco[notmissing]FRA-TOU_21 0.503 1.987 0.164 1.549
wruaco[notmissing]SWE-NSW_02 1.040 0.961 0.525 2.063
wruaco[notmissing]SWE-NSW_03 1.103 0.906 0.613 1.987
wruaco[notmissing]UNK-BEL_21 0.497 2.012 0.228 1.083
Rsquare= 0.042 (max possible= 0.914 )
Likelihood ratio test= 105 on 17 df, p=1.04e-14
Wald test = 82.9 on 17 df, p=1.17e-10
Score (logrank) test = 107 on 17 df, p=4.11e-15, Robust = 72.2 p=9.12e-09
(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