response: (stroke3 %in% c('1','2'))
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=2428 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.052670 0.949 0.09850 0.11083
wz[notmissing, selectedcovariates]1 0.072339 1.075 0.04468 0.05311
wz[notmissing, selectedcovariates]2 -0.732956 0.480 0.17789 0.20242
wz[notmissing, selectedcovariates]3 0.000377 1.000 0.01243 0.01378
wz[notmissing, selectedcovariates]4 0.016280 1.016 0.00373 0.00431
wz[notmissing, selectedcovariates]5 0.377986 1.459 0.12619 0.14007
wz[notmissing, selectedcovariates]6 0.746900 2.110 0.16172 0.20052
wruaco[notmissing]FIN-EAS_03 0.261094 1.298 0.22001 0.23793
wruaco[notmissing]FIN-EAS_24 -0.212139 0.809 0.19399 0.22362
wruaco[notmissing]FIN-WES_03 -0.152683 0.858 0.22914 0.24678
wruaco[notmissing]FIN-WES_24 -0.270706 0.763 0.20875 0.24457
wruaco[notmissing]FRA-LIL_21 -1.003590 0.367 0.52069 0.56994
wruaco[notmissing]FRA-STR_21 -1.090191 0.336 0.47895 0.52176
wruaco[notmissing]FRA-TOU_21 -0.763979 0.466 0.53022 0.57262
wruaco[notmissing]SWE-NSW_02 -0.001082 0.999 0.28537 0.34355
wruaco[notmissing]SWE-NSW_03 0.032449 1.033 0.26310 0.29447
wruaco[notmissing]UNK-BEL_21 -0.688478 0.502 0.36127 0.38502
z p
allele[notmissing] -0.47525 0.63000
wz[notmissing, selectedcovariates]1 1.36217 0.17000
wz[notmissing, selectedcovariates]2 -3.62102 0.00029
wz[notmissing, selectedcovariates]3 0.02732 0.98000
wz[notmissing, selectedcovariates]4 3.77993 0.00016
wz[notmissing, selectedcovariates]5 2.69865 0.00700
wz[notmissing, selectedcovariates]6 3.72483 0.00020
wruaco[notmissing]FIN-EAS_03 1.09737 0.27000
wruaco[notmissing]FIN-EAS_24 -0.94865 0.34000
wruaco[notmissing]FIN-WES_03 -0.61870 0.54000
wruaco[notmissing]FIN-WES_24 -1.10685 0.27000
wruaco[notmissing]FRA-LIL_21 -1.76087 0.07800
wruaco[notmissing]FRA-STR_21 -2.08945 0.03700
wruaco[notmissing]FRA-TOU_21 -1.33418 0.18000
wruaco[notmissing]SWE-NSW_02 -0.00315 1.00000
wruaco[notmissing]SWE-NSW_03 0.11019 0.91000
wruaco[notmissing]UNK-BEL_21 -1.78815 0.07400
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.949 1.054 0.763 1.179
wz[notmissing, selectedcovariates]1 1.075 0.930 0.969 1.193
wz[notmissing, selectedcovariates]2 0.480 2.081 0.323 0.714
wz[notmissing, selectedcovariates]3 1.000 1.000 0.974 1.028
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.459 0.685 1.109 1.920
wz[notmissing, selectedcovariates]6 2.110 0.474 1.425 3.126
wruaco[notmissing]FIN-EAS_03 1.298 0.770 0.814 2.070
wruaco[notmissing]FIN-EAS_24 0.809 1.236 0.522 1.254
wruaco[notmissing]FIN-WES_03 0.858 1.165 0.529 1.392
wruaco[notmissing]FIN-WES_24 0.763 1.311 0.472 1.232
wruaco[notmissing]FRA-LIL_21 0.367 2.728 0.120 1.120
wruaco[notmissing]FRA-STR_21 0.336 2.975 0.121 0.935
wruaco[notmissing]FRA-TOU_21 0.466 2.147 0.152 1.431
wruaco[notmissing]SWE-NSW_02 0.999 1.001 0.509 1.959
wruaco[notmissing]SWE-NSW_03 1.033 0.968 0.580 1.840
wruaco[notmissing]UNK-BEL_21 0.502 1.991 0.236 1.068
Rsquare= 0.043 (max possible= 0.922 )
Likelihood ratio test= 107 on 17 df, p=5.11e-15
Wald test = 81.8 on 17 df, p=1.8e-10
Score (logrank) test = 108 on 17 df, p=2.55e-15, Robust = 72.5 p=8.03e-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