response: (stroke4 %in% c('1','2','3','4','5','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=2454 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.250130 0.779 0.24751 0.28026
wz[notmissing, selectedcovariates]1 0.064063 1.066 0.04420 0.05241
wz[notmissing, selectedcovariates]2 -0.723928 0.485 0.17475 0.19953
wz[notmissing, selectedcovariates]3 -0.000964 0.999 0.01228 0.01361
wz[notmissing, selectedcovariates]4 0.016354 1.016 0.00368 0.00423
wz[notmissing, selectedcovariates]5 0.355331 1.427 0.12413 0.13875
wz[notmissing, selectedcovariates]6 0.760870 2.140 0.15810 0.19604
wruaco[notmissing]FIN-EAS_03 0.244595 1.277 0.21533 0.23642
wruaco[notmissing]FIN-EAS_24 -0.242191 0.785 0.19226 0.22227
wruaco[notmissing]FIN-WES_03 -0.119973 0.887 0.22411 0.24345
wruaco[notmissing]FIN-WES_24 -0.265481 0.767 0.20570 0.24315
wruaco[notmissing]FRA-LIL_21 -1.064900 0.345 0.51981 0.56948
wruaco[notmissing]FRA-STR_21 -1.181069 0.307 0.47724 0.51872
wruaco[notmissing]FRA-TOU_21 -0.821882 0.440 0.52893 0.57049
wruaco[notmissing]SWE-NSW_02 -0.028659 0.972 0.28321 0.33820
wruaco[notmissing]SWE-NSW_03 0.015511 1.016 0.26135 0.29239
wruaco[notmissing]UNK-BEL_21 -0.730152 0.482 0.35978 0.38169
z p
allele[notmissing] -0.8925 0.37000
wz[notmissing, selectedcovariates]1 1.2223 0.22000
wz[notmissing, selectedcovariates]2 -3.6282 0.00029
wz[notmissing, selectedcovariates]3 -0.0709 0.94000
wz[notmissing, selectedcovariates]4 3.8633 0.00011
wz[notmissing, selectedcovariates]5 2.5609 0.01000
wz[notmissing, selectedcovariates]6 3.8811 0.00010
wruaco[notmissing]FIN-EAS_03 1.0346 0.30000
wruaco[notmissing]FIN-EAS_24 -1.0896 0.28000
wruaco[notmissing]FIN-WES_03 -0.4928 0.62000
wruaco[notmissing]FIN-WES_24 -1.0918 0.27000
wruaco[notmissing]FRA-LIL_21 -1.8700 0.06100
wruaco[notmissing]FRA-STR_21 -2.2769 0.02300
wruaco[notmissing]FRA-TOU_21 -1.4407 0.15000
wruaco[notmissing]SWE-NSW_02 -0.0847 0.93000
wruaco[notmissing]SWE-NSW_03 0.0530 0.96000
wruaco[notmissing]UNK-BEL_21 -1.9129 0.05600
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.779 1.284 0.450 1.349
wz[notmissing, selectedcovariates]1 1.066 0.938 0.962 1.182
wz[notmissing, selectedcovariates]2 0.485 2.063 0.328 0.717
wz[notmissing, selectedcovariates]3 0.999 1.001 0.973 1.026
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.427 0.701 1.087 1.873
wz[notmissing, selectedcovariates]6 2.140 0.467 1.457 3.143
wruaco[notmissing]FIN-EAS_03 1.277 0.783 0.803 2.030
wruaco[notmissing]FIN-EAS_24 0.785 1.274 0.508 1.213
wruaco[notmissing]FIN-WES_03 0.887 1.127 0.550 1.429
wruaco[notmissing]FIN-WES_24 0.767 1.304 0.476 1.235
wruaco[notmissing]FRA-LIL_21 0.345 2.901 0.113 1.053
wruaco[notmissing]FRA-STR_21 0.307 3.258 0.111 0.848
wruaco[notmissing]FRA-TOU_21 0.440 2.275 0.144 1.345
wruaco[notmissing]SWE-NSW_02 0.972 1.029 0.501 1.886
wruaco[notmissing]SWE-NSW_03 1.016 0.985 0.573 1.801
wruaco[notmissing]UNK-BEL_21 0.482 2.075 0.228 1.018
Rsquare= 0.044 (max possible= 0.926 )
Likelihood ratio test= 110 on 17 df, p=1e-15
Wald test = 84.3 on 17 df, p=6.47e-11
Score (logrank) test = 111 on 17 df, p=6.66e-16, Robust = 74.8 p=3.17e-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