response: (stroke3 %in% c('1','2'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: SMAD3_MG_452
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=2442 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.019961 1.020 0.08548 0.09867
wz[notmissing, selectedcovariates]1 0.062949 1.065 0.04477 0.05307
wz[notmissing, selectedcovariates]2 -0.753167 0.471 0.17739 0.20139
wz[notmissing, selectedcovariates]3 -0.000838 0.999 0.01245 0.01377
wz[notmissing, selectedcovariates]4 0.016206 1.016 0.00372 0.00429
wz[notmissing, selectedcovariates]5 0.368894 1.446 0.12556 0.13842
wz[notmissing, selectedcovariates]6 0.743419 2.103 0.16139 0.19956
wruaco[notmissing]FIN-EAS_03 0.264442 1.303 0.21853 0.23794
wruaco[notmissing]FIN-EAS_24 -0.199176 0.819 0.19321 0.22215
wruaco[notmissing]FIN-WES_03 -0.199970 0.819 0.23226 0.24909
wruaco[notmissing]FIN-WES_24 -0.230886 0.794 0.20652 0.24295
wruaco[notmissing]FRA-LIL_21 -1.020647 0.360 0.52152 0.57241
wruaco[notmissing]FRA-STR_21 -1.125979 0.324 0.47824 0.51901
wruaco[notmissing]FRA-TOU_21 -0.765025 0.465 0.52990 0.57134
wruaco[notmissing]SWE-NSW_02 0.013509 1.014 0.28453 0.33886
wruaco[notmissing]SWE-NSW_03 0.045361 1.046 0.26260 0.29317
wruaco[notmissing]UNK-BEL_21 -0.699144 0.497 0.36125 0.38255
z p
allele[notmissing] 0.2023 0.84000
wz[notmissing, selectedcovariates]1 1.1861 0.24000
wz[notmissing, selectedcovariates]2 -3.7398 0.00018
wz[notmissing, selectedcovariates]3 -0.0609 0.95000
wz[notmissing, selectedcovariates]4 3.7751 0.00016
wz[notmissing, selectedcovariates]5 2.6651 0.00770
wz[notmissing, selectedcovariates]6 3.7254 0.00020
wruaco[notmissing]FIN-EAS_03 1.1114 0.27000
wruaco[notmissing]FIN-EAS_24 -0.8966 0.37000
wruaco[notmissing]FIN-WES_03 -0.8028 0.42000
wruaco[notmissing]FIN-WES_24 -0.9504 0.34000
wruaco[notmissing]FRA-LIL_21 -1.7831 0.07500
wruaco[notmissing]FRA-STR_21 -2.1695 0.03000
wruaco[notmissing]FRA-TOU_21 -1.3390 0.18000
wruaco[notmissing]SWE-NSW_02 0.0399 0.97000
wruaco[notmissing]SWE-NSW_03 0.1547 0.88000
wruaco[notmissing]UNK-BEL_21 -1.8276 0.06800
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.020 0.980 0.841 1.238
wz[notmissing, selectedcovariates]1 1.065 0.939 0.960 1.182
wz[notmissing, selectedcovariates]2 0.471 2.124 0.317 0.699
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.446 0.691 1.103 1.897
wz[notmissing, selectedcovariates]6 2.103 0.475 1.422 3.110
wruaco[notmissing]FIN-EAS_03 1.303 0.768 0.817 2.077
wruaco[notmissing]FIN-EAS_24 0.819 1.220 0.530 1.266
wruaco[notmissing]FIN-WES_03 0.819 1.221 0.502 1.334
wruaco[notmissing]FIN-WES_24 0.794 1.260 0.493 1.278
wruaco[notmissing]FRA-LIL_21 0.360 2.775 0.117 1.107
wruaco[notmissing]FRA-STR_21 0.324 3.083 0.117 0.897
wruaco[notmissing]FRA-TOU_21 0.465 2.149 0.152 1.426
wruaco[notmissing]SWE-NSW_02 1.014 0.987 0.522 1.969
wruaco[notmissing]SWE-NSW_03 1.046 0.956 0.589 1.859
wruaco[notmissing]UNK-BEL_21 0.497 2.012 0.235 1.052
Rsquare= 0.042 (max possible= 0.922 )
Likelihood ratio test= 106 on 17 df, p=6.77e-15
Wald test = 82.1 on 17 df, p=1.61e-10
Score (logrank) test = 107 on 17 df, p=4.22e-15, Robust = 72.8 p=7.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