response: (stroke4 %in% c('1','2','3','4','5','6'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: 19Q12_MG_446
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=2426 (44 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.06590 1.068 0.09663 0.10989
wz[notmissing, selectedcovariates]1 0.07606 1.079 0.04397 0.05258
wz[notmissing, selectedcovariates]2 -0.70973 0.492 0.17538 0.20107
wz[notmissing, selectedcovariates]3 -0.00202 0.998 0.01229 0.01372
wz[notmissing, selectedcovariates]4 0.01612 1.016 0.00367 0.00423
wz[notmissing, selectedcovariates]5 0.35323 1.424 0.12423 0.13974
wz[notmissing, selectedcovariates]6 0.83272 2.300 0.15595 0.19610
wruaco[notmissing]FIN-EAS_03 0.25373 1.289 0.21527 0.23635
wruaco[notmissing]FIN-EAS_24 -0.20669 0.813 0.19234 0.22292
wruaco[notmissing]FIN-WES_03 -0.11389 0.892 0.22201 0.24347
wruaco[notmissing]FIN-WES_24 -0.23200 0.793 0.20169 0.23809
wruaco[notmissing]FRA-LIL_21 -1.09109 0.336 0.51991 0.57104
wruaco[notmissing]FRA-STR_21 -1.21106 0.298 0.47750 0.52274
wruaco[notmissing]FRA-TOU_21 -0.81328 0.443 0.52918 0.57187
wruaco[notmissing]SWE-NSW_02 -0.01523 0.985 0.28339 0.33943
wruaco[notmissing]SWE-NSW_03 0.06220 1.064 0.25711 0.28902
wruaco[notmissing]UNK-BEL_21 -0.82392 0.439 0.37447 0.40065
z p
allele[notmissing] 0.5997 5.5e-01
wz[notmissing, selectedcovariates]1 1.4465 1.5e-01
wz[notmissing, selectedcovariates]2 -3.5298 4.2e-04
wz[notmissing, selectedcovariates]3 -0.1474 8.8e-01
wz[notmissing, selectedcovariates]4 3.8142 1.4e-04
wz[notmissing, selectedcovariates]5 2.5278 1.1e-02
wz[notmissing, selectedcovariates]6 4.2464 2.2e-05
wruaco[notmissing]FIN-EAS_03 1.0736 2.8e-01
wruaco[notmissing]FIN-EAS_24 -0.9272 3.5e-01
wruaco[notmissing]FIN-WES_03 -0.4678 6.4e-01
wruaco[notmissing]FIN-WES_24 -0.9744 3.3e-01
wruaco[notmissing]FRA-LIL_21 -1.9107 5.6e-02
wruaco[notmissing]FRA-STR_21 -2.3168 2.1e-02
wruaco[notmissing]FRA-TOU_21 -1.4221 1.5e-01
wruaco[notmissing]SWE-NSW_02 -0.0449 9.6e-01
wruaco[notmissing]SWE-NSW_03 0.2152 8.3e-01
wruaco[notmissing]UNK-BEL_21 -2.0565 4.0e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.068 0.936 0.861 1.325
wz[notmissing, selectedcovariates]1 1.079 0.927 0.973 1.196
wz[notmissing, selectedcovariates]2 0.492 2.033 0.332 0.729
wz[notmissing, selectedcovariates]3 0.998 1.002 0.972 1.025
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.424 0.702 1.083 1.872
wz[notmissing, selectedcovariates]6 2.300 0.435 1.566 3.377
wruaco[notmissing]FIN-EAS_03 1.289 0.776 0.811 2.048
wruaco[notmissing]FIN-EAS_24 0.813 1.230 0.525 1.259
wruaco[notmissing]FIN-WES_03 0.892 1.121 0.554 1.438
wruaco[notmissing]FIN-WES_24 0.793 1.261 0.497 1.264
wruaco[notmissing]FRA-LIL_21 0.336 2.978 0.110 1.029
wruaco[notmissing]FRA-STR_21 0.298 3.357 0.107 0.830
wruaco[notmissing]FRA-TOU_21 0.443 2.255 0.145 1.360
wruaco[notmissing]SWE-NSW_02 0.985 1.015 0.506 1.916
wruaco[notmissing]SWE-NSW_03 1.064 0.940 0.604 1.875
wruaco[notmissing]UNK-BEL_21 0.439 2.279 0.200 0.962
Rsquare= 0.046 (max possible= 0.928 )
Likelihood ratio test= 114 on 17 df, p=2.22e-16
Wald test = 86 on 17 df, p=3.3e-11
Score (logrank) test = 117 on 17 df, p=0, Robust = 74.6 p=3.42e-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