response: (stroke3 %in% c('1','2'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: 10Q11_21_MG_458
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=2423 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.15244 1.165 0.09967 0.1088
wz[notmissing, selectedcovariates]1 0.06721 1.070 0.04449 0.0526
wz[notmissing, selectedcovariates]2 -0.77637 0.460 0.17761 0.2017
wz[notmissing, selectedcovariates]3 -0.00274 0.997 0.01243 0.0138
wz[notmissing, selectedcovariates]4 0.01656 1.017 0.00373 0.0043
wz[notmissing, selectedcovariates]5 0.37628 1.457 0.12616 0.1395
wz[notmissing, selectedcovariates]6 0.75535 2.128 0.16163 0.2008
wruaco[notmissing]FIN-EAS_03 0.26930 1.309 0.21776 0.2363
wruaco[notmissing]FIN-EAS_24 -0.19029 0.827 0.19354 0.2228
wruaco[notmissing]FIN-WES_03 -0.18590 0.830 0.23256 0.2496
wruaco[notmissing]FIN-WES_24 -0.21628 0.806 0.20619 0.2431
wruaco[notmissing]FRA-LIL_21 -1.03003 0.357 0.52093 0.5718
wruaco[notmissing]FRA-STR_21 -1.11778 0.327 0.47909 0.5229
wruaco[notmissing]FRA-TOU_21 -0.72820 0.483 0.52932 0.5718
wruaco[notmissing]SWE-NSW_02 -0.09319 0.911 0.29948 0.3581
wruaco[notmissing]SWE-NSW_03 0.08980 1.094 0.26278 0.2955
wruaco[notmissing]UNK-BEL_21 -0.70422 0.494 0.36080 0.3827
z p
allele[notmissing] 1.401 0.16000
wz[notmissing, selectedcovariates]1 1.278 0.20000
wz[notmissing, selectedcovariates]2 -3.849 0.00012
wz[notmissing, selectedcovariates]3 -0.199 0.84000
wz[notmissing, selectedcovariates]4 3.852 0.00012
wz[notmissing, selectedcovariates]5 2.697 0.00700
wz[notmissing, selectedcovariates]6 3.761 0.00017
wruaco[notmissing]FIN-EAS_03 1.140 0.25000
wruaco[notmissing]FIN-EAS_24 -0.854 0.39000
wruaco[notmissing]FIN-WES_03 -0.745 0.46000
wruaco[notmissing]FIN-WES_24 -0.890 0.37000
wruaco[notmissing]FRA-LIL_21 -1.801 0.07200
wruaco[notmissing]FRA-STR_21 -2.138 0.03300
wruaco[notmissing]FRA-TOU_21 -1.273 0.20000
wruaco[notmissing]SWE-NSW_02 -0.260 0.79000
wruaco[notmissing]SWE-NSW_03 0.304 0.76000
wruaco[notmissing]UNK-BEL_21 -1.840 0.06600
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.165 0.859 0.941 1.442
wz[notmissing, selectedcovariates]1 1.070 0.935 0.965 1.186
wz[notmissing, selectedcovariates]2 0.460 2.174 0.310 0.683
wz[notmissing, selectedcovariates]3 0.997 1.003 0.971 1.025
wz[notmissing, selectedcovariates]4 1.017 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.457 0.686 1.108 1.915
wz[notmissing, selectedcovariates]6 2.128 0.470 1.436 3.155
wruaco[notmissing]FIN-EAS_03 1.309 0.764 0.824 2.080
wruaco[notmissing]FIN-EAS_24 0.827 1.210 0.534 1.279
wruaco[notmissing]FIN-WES_03 0.830 1.204 0.509 1.354
wruaco[notmissing]FIN-WES_24 0.806 1.241 0.500 1.297
wruaco[notmissing]FRA-LIL_21 0.357 2.801 0.116 1.095
wruaco[notmissing]FRA-STR_21 0.327 3.058 0.117 0.911
wruaco[notmissing]FRA-TOU_21 0.483 2.071 0.157 1.481
wruaco[notmissing]SWE-NSW_02 0.911 1.098 0.452 1.838
wruaco[notmissing]SWE-NSW_03 1.094 0.914 0.613 1.952
wruaco[notmissing]UNK-BEL_21 0.494 2.022 0.234 1.047
Rsquare= 0.044 (max possible= 0.922 )
Likelihood ratio test= 110 on 17 df, p=1.11e-15
Wald test = 84.4 on 17 df, p=6.39e-11
Score (logrank) test = 111 on 17 df, p=7.77e-16, Robust = 75 p=2.98e-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