response: (stroke4 %in% c('2','4','6'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: CELSR2_MG_438
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=2429 (44 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.01259 1.013 0.09150 0.10429
wz[notmissing, selectedcovariates]1 0.05990 1.062 0.04557 0.05302
wz[notmissing, selectedcovariates]2 -0.74709 0.474 0.18107 0.20437
wz[notmissing, selectedcovariates]3 0.00805 1.008 0.01240 0.01340
wz[notmissing, selectedcovariates]4 0.01526 1.015 0.00379 0.00437
wz[notmissing, selectedcovariates]5 0.39680 1.487 0.12848 0.14326
wz[notmissing, selectedcovariates]6 0.79150 2.207 0.16075 0.19882
wruaco[notmissing]FIN-EAS_03 0.34506 1.412 0.22034 0.23873
wruaco[notmissing]FIN-EAS_24 -0.09609 0.908 0.19566 0.22605
wruaco[notmissing]FIN-WES_03 -0.00888 0.991 0.22653 0.24733
wruaco[notmissing]FIN-WES_24 -0.11765 0.889 0.20648 0.24159
wruaco[notmissing]FRA-LIL_21 -0.96613 0.381 0.52107 0.57157
wruaco[notmissing]FRA-STR_21 -1.12268 0.325 0.47907 0.52086
wruaco[notmissing]FRA-TOU_21 -0.73027 0.482 0.53148 0.57497
wruaco[notmissing]SWE-NSW_02 0.05516 1.057 0.29258 0.35105
wruaco[notmissing]SWE-NSW_03 0.09184 1.096 0.26944 0.30022
wruaco[notmissing]UNK-BEL_21 -0.71658 0.488 0.37668 0.39760
z p
allele[notmissing] 0.1207 9.0e-01
wz[notmissing, selectedcovariates]1 1.1297 2.6e-01
wz[notmissing, selectedcovariates]2 -3.6555 2.6e-04
wz[notmissing, selectedcovariates]3 0.6008 5.5e-01
wz[notmissing, selectedcovariates]4 3.4923 4.8e-04
wz[notmissing, selectedcovariates]5 2.7698 5.6e-03
wz[notmissing, selectedcovariates]6 3.9809 6.9e-05
wruaco[notmissing]FIN-EAS_03 1.4454 1.5e-01
wruaco[notmissing]FIN-EAS_24 -0.4251 6.7e-01
wruaco[notmissing]FIN-WES_03 -0.0359 9.7e-01
wruaco[notmissing]FIN-WES_24 -0.4870 6.3e-01
wruaco[notmissing]FRA-LIL_21 -1.6903 9.1e-02
wruaco[notmissing]FRA-STR_21 -2.1554 3.1e-02
wruaco[notmissing]FRA-TOU_21 -1.2701 2.0e-01
wruaco[notmissing]SWE-NSW_02 0.1571 8.8e-01
wruaco[notmissing]SWE-NSW_03 0.3059 7.6e-01
wruaco[notmissing]UNK-BEL_21 -1.8023 7.2e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.013 0.987 0.825 1.242
wz[notmissing, selectedcovariates]1 1.062 0.942 0.957 1.178
wz[notmissing, selectedcovariates]2 0.474 2.111 0.317 0.707
wz[notmissing, selectedcovariates]3 1.008 0.992 0.982 1.035
wz[notmissing, selectedcovariates]4 1.015 0.985 1.007 1.024
wz[notmissing, selectedcovariates]5 1.487 0.672 1.123 1.969
wz[notmissing, selectedcovariates]6 2.207 0.453 1.495 3.258
wruaco[notmissing]FIN-EAS_03 1.412 0.708 0.884 2.255
wruaco[notmissing]FIN-EAS_24 0.908 1.101 0.583 1.415
wruaco[notmissing]FIN-WES_03 0.991 1.009 0.610 1.609
wruaco[notmissing]FIN-WES_24 0.889 1.125 0.554 1.427
wruaco[notmissing]FRA-LIL_21 0.381 2.628 0.124 1.167
wruaco[notmissing]FRA-STR_21 0.325 3.073 0.117 0.903
wruaco[notmissing]FRA-TOU_21 0.482 2.076 0.156 1.487
wruaco[notmissing]SWE-NSW_02 1.057 0.946 0.531 2.103
wruaco[notmissing]SWE-NSW_03 1.096 0.912 0.609 1.974
wruaco[notmissing]UNK-BEL_21 0.488 2.047 0.224 1.065
Rsquare= 0.043 (max possible= 0.915 )
Likelihood ratio test= 108 on 17 df, p=3.11e-15
Wald test = 84.8 on 17 df, p=5.42e-11
Score (logrank) test = 111 on 17 df, p=8.88e-16, Robust = 73.2 p=5.95e-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