response: (stroke4 %in% c('2','4','6'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MIA3_MG_454
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=2409 (42 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.10416 0.901 0.10330 0.11555
wz[notmissing, selectedcovariates]1 0.06982 1.072 0.04566 0.05353
wz[notmissing, selectedcovariates]2 -0.73581 0.479 0.18236 0.20646
wz[notmissing, selectedcovariates]3 0.00919 1.009 0.01241 0.01342
wz[notmissing, selectedcovariates]4 0.01529 1.015 0.00383 0.00443
wz[notmissing, selectedcovariates]5 0.40250 1.496 0.12940 0.14425
wz[notmissing, selectedcovariates]6 0.77419 2.169 0.16259 0.20107
wruaco[notmissing]FIN-EAS_03 0.35676 1.429 0.22221 0.23943
wruaco[notmissing]FIN-EAS_24 -0.11761 0.889 0.19633 0.22692
wruaco[notmissing]FIN-WES_03 0.02378 1.024 0.22606 0.24710
wruaco[notmissing]FIN-WES_24 -0.18539 0.831 0.21291 0.24875
wruaco[notmissing]FRA-LIL_21 -0.93788 0.391 0.52115 0.57304
wruaco[notmissing]FRA-STR_21 -1.02199 0.360 0.48003 0.52496
wruaco[notmissing]FRA-TOU_21 -0.67917 0.507 0.53140 0.57513
wruaco[notmissing]SWE-NSW_02 0.01613 1.016 0.29400 0.35650
wruaco[notmissing]SWE-NSW_03 0.07912 1.082 0.27032 0.30213
wruaco[notmissing]UNK-BEL_21 -0.70115 0.496 0.37762 0.40020
z p
allele[notmissing] -0.9014 0.37000
wz[notmissing, selectedcovariates]1 1.3045 0.19000
wz[notmissing, selectedcovariates]2 -3.5639 0.00037
wz[notmissing, selectedcovariates]3 0.6848 0.49000
wz[notmissing, selectedcovariates]4 3.4497 0.00056
wz[notmissing, selectedcovariates]5 2.7903 0.00530
wz[notmissing, selectedcovariates]6 3.8504 0.00012
wruaco[notmissing]FIN-EAS_03 1.4900 0.14000
wruaco[notmissing]FIN-EAS_24 -0.5183 0.60000
wruaco[notmissing]FIN-WES_03 0.0963 0.92000
wruaco[notmissing]FIN-WES_24 -0.7453 0.46000
wruaco[notmissing]FRA-LIL_21 -1.6367 0.10000
wruaco[notmissing]FRA-STR_21 -1.9468 0.05200
wruaco[notmissing]FRA-TOU_21 -1.1809 0.24000
wruaco[notmissing]SWE-NSW_02 0.0453 0.96000
wruaco[notmissing]SWE-NSW_03 0.2619 0.79000
wruaco[notmissing]UNK-BEL_21 -1.7520 0.08000
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.901 1.110 0.718 1.130
wz[notmissing, selectedcovariates]1 1.072 0.933 0.966 1.191
wz[notmissing, selectedcovariates]2 0.479 2.087 0.320 0.718
wz[notmissing, selectedcovariates]3 1.009 0.991 0.983 1.036
wz[notmissing, selectedcovariates]4 1.015 0.985 1.007 1.024
wz[notmissing, selectedcovariates]5 1.496 0.669 1.127 1.984
wz[notmissing, selectedcovariates]6 2.169 0.461 1.462 3.216
wruaco[notmissing]FIN-EAS_03 1.429 0.700 0.894 2.284
wruaco[notmissing]FIN-EAS_24 0.889 1.125 0.570 1.387
wruaco[notmissing]FIN-WES_03 1.024 0.976 0.631 1.662
wruaco[notmissing]FIN-WES_24 0.831 1.204 0.510 1.353
wruaco[notmissing]FRA-LIL_21 0.391 2.555 0.127 1.204
wruaco[notmissing]FRA-STR_21 0.360 2.779 0.129 1.007
wruaco[notmissing]FRA-TOU_21 0.507 1.972 0.164 1.565
wruaco[notmissing]SWE-NSW_02 1.016 0.984 0.505 2.044
wruaco[notmissing]SWE-NSW_03 1.082 0.924 0.599 1.957
wruaco[notmissing]UNK-BEL_21 0.496 2.016 0.226 1.087
Rsquare= 0.043 (max possible= 0.914 )
Likelihood ratio test= 107 on 17 df, p=5.22e-15
Wald test = 83 on 17 df, p=1.13e-10
Score (logrank) test = 110 on 17 df, p=1.33e-15, Robust = 71.8 p=1.07e-08
(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