response: (stroke3 %in% c('1','2'))
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=2449 (45 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.00913 1.009 0.08889 0.10197
wz[notmissing, selectedcovariates]1 0.06122 1.063 0.04459 0.05271
wz[notmissing, selectedcovariates]2 -0.74812 0.473 0.17653 0.20008
wz[notmissing, selectedcovariates]3 -0.00171 0.998 0.01242 0.01381
wz[notmissing, selectedcovariates]4 0.01589 1.016 0.00369 0.00425
wz[notmissing, selectedcovariates]5 0.36818 1.445 0.12515 0.13890
wz[notmissing, selectedcovariates]6 0.76605 2.151 0.15994 0.19821
wruaco[notmissing]FIN-EAS_03 0.25371 1.289 0.21793 0.23680
wruaco[notmissing]FIN-EAS_24 -0.16607 0.847 0.19185 0.22157
wruaco[notmissing]FIN-WES_03 -0.17751 0.837 0.22931 0.24669
wruaco[notmissing]FIN-WES_24 -0.20375 0.816 0.20237 0.23738
wruaco[notmissing]FRA-LIL_21 -1.02792 0.358 0.52058 0.56901
wruaco[notmissing]FRA-STR_21 -1.18814 0.305 0.47792 0.51735
wruaco[notmissing]FRA-TOU_21 -0.82034 0.440 0.53027 0.57233
wruaco[notmissing]SWE-NSW_02 0.03544 1.036 0.28405 0.33894
wruaco[notmissing]SWE-NSW_03 0.04557 1.047 0.26219 0.29233
wruaco[notmissing]UNK-BEL_21 -0.70897 0.492 0.36023 0.38247
z p
allele[notmissing] 0.0896 0.93000
wz[notmissing, selectedcovariates]1 1.1615 0.25000
wz[notmissing, selectedcovariates]2 -3.7392 0.00018
wz[notmissing, selectedcovariates]3 -0.1238 0.90000
wz[notmissing, selectedcovariates]4 3.7427 0.00018
wz[notmissing, selectedcovariates]5 2.6507 0.00800
wz[notmissing, selectedcovariates]6 3.8648 0.00011
wruaco[notmissing]FIN-EAS_03 1.0714 0.28000
wruaco[notmissing]FIN-EAS_24 -0.7495 0.45000
wruaco[notmissing]FIN-WES_03 -0.7196 0.47000
wruaco[notmissing]FIN-WES_24 -0.8583 0.39000
wruaco[notmissing]FRA-LIL_21 -1.8065 0.07100
wruaco[notmissing]FRA-STR_21 -2.2966 0.02200
wruaco[notmissing]FRA-TOU_21 -1.4333 0.15000
wruaco[notmissing]SWE-NSW_02 0.1046 0.92000
wruaco[notmissing]SWE-NSW_03 0.1559 0.88000
wruaco[notmissing]UNK-BEL_21 -1.8537 0.06400
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.009 0.991 0.826 1.23
wz[notmissing, selectedcovariates]1 1.063 0.941 0.959 1.18
wz[notmissing, selectedcovariates]2 0.473 2.113 0.320 0.70
wz[notmissing, selectedcovariates]3 0.998 1.002 0.972 1.03
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.02
wz[notmissing, selectedcovariates]5 1.445 0.692 1.101 1.90
wz[notmissing, selectedcovariates]6 2.151 0.465 1.459 3.17
wruaco[notmissing]FIN-EAS_03 1.289 0.776 0.810 2.05
wruaco[notmissing]FIN-EAS_24 0.847 1.181 0.549 1.31
wruaco[notmissing]FIN-WES_03 0.837 1.194 0.516 1.36
wruaco[notmissing]FIN-WES_24 0.816 1.226 0.512 1.30
wruaco[notmissing]FRA-LIL_21 0.358 2.795 0.117 1.09
wruaco[notmissing]FRA-STR_21 0.305 3.281 0.111 0.84
wruaco[notmissing]FRA-TOU_21 0.440 2.271 0.143 1.35
wruaco[notmissing]SWE-NSW_02 1.036 0.965 0.533 2.01
wruaco[notmissing]SWE-NSW_03 1.047 0.955 0.590 1.86
wruaco[notmissing]UNK-BEL_21 0.492 2.032 0.233 1.04
Rsquare= 0.043 (max possible= 0.923 )
Likelihood ratio test= 107 on 17 df, p=4.33e-15
Wald test = 82.5 on 17 df, p=1.36e-10
Score (logrank) test = 108 on 17 df, p=2.44e-15, Robust = 73.3 p=5.81e-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