response: (stroke4 %in% c('1','2','3','4','5','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=2461 (45 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.01446 1.015 0.08785 0.10106
wz[notmissing, selectedcovariates]1 0.06656 1.069 0.04394 0.05221
wz[notmissing, selectedcovariates]2 -0.71204 0.491 0.17360 0.19770
wz[notmissing, selectedcovariates]3 -0.00246 0.998 0.01224 0.01362
wz[notmissing, selectedcovariates]4 0.01599 1.016 0.00364 0.00418
wz[notmissing, selectedcovariates]5 0.34733 1.415 0.12337 0.13866
wz[notmissing, selectedcovariates]6 0.81023 2.248 0.15554 0.19438
wruaco[notmissing]FIN-EAS_03 0.24971 1.284 0.21467 0.23546
wruaco[notmissing]FIN-EAS_24 -0.19992 0.819 0.19049 0.22112
wruaco[notmissing]FIN-WES_03 -0.10493 0.900 0.22092 0.24114
wruaco[notmissing]FIN-WES_24 -0.23706 0.789 0.20119 0.23717
wruaco[notmissing]FRA-LIL_21 -1.08857 0.337 0.51979 0.56925
wruaco[notmissing]FRA-STR_21 -1.24315 0.288 0.47689 0.51726
wruaco[notmissing]FRA-TOU_21 -0.86188 0.422 0.52934 0.57142
wruaco[notmissing]SWE-NSW_02 -0.00934 0.991 0.28299 0.33883
wruaco[notmissing]SWE-NSW_03 0.06286 1.065 0.25646 0.28874
wruaco[notmissing]UNK-BEL_21 -0.74631 0.474 0.35886 0.38112
z p
allele[notmissing] 0.1430 8.9e-01
wz[notmissing, selectedcovariates]1 1.2750 2.0e-01
wz[notmissing, selectedcovariates]2 -3.6015 3.2e-04
wz[notmissing, selectedcovariates]3 -0.1806 8.6e-01
wz[notmissing, selectedcovariates]4 3.8258 1.3e-04
wz[notmissing, selectedcovariates]5 2.5049 1.2e-02
wz[notmissing, selectedcovariates]6 4.1683 3.1e-05
wruaco[notmissing]FIN-EAS_03 1.0606 2.9e-01
wruaco[notmissing]FIN-EAS_24 -0.9041 3.7e-01
wruaco[notmissing]FIN-WES_03 -0.4352 6.6e-01
wruaco[notmissing]FIN-WES_24 -0.9996 3.2e-01
wruaco[notmissing]FRA-LIL_21 -1.9123 5.6e-02
wruaco[notmissing]FRA-STR_21 -2.4034 1.6e-02
wruaco[notmissing]FRA-TOU_21 -1.5083 1.3e-01
wruaco[notmissing]SWE-NSW_02 -0.0276 9.8e-01
wruaco[notmissing]SWE-NSW_03 0.2177 8.3e-01
wruaco[notmissing]UNK-BEL_21 -1.9582 5.0e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.015 0.986 0.832 1.237
wz[notmissing, selectedcovariates]1 1.069 0.936 0.965 1.184
wz[notmissing, selectedcovariates]2 0.491 2.038 0.333 0.723
wz[notmissing, selectedcovariates]3 0.998 1.002 0.971 1.025
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.024
wz[notmissing, selectedcovariates]5 1.415 0.707 1.078 1.857
wz[notmissing, selectedcovariates]6 2.248 0.445 1.536 3.291
wruaco[notmissing]FIN-EAS_03 1.284 0.779 0.809 2.036
wruaco[notmissing]FIN-EAS_24 0.819 1.221 0.531 1.263
wruaco[notmissing]FIN-WES_03 0.900 1.111 0.561 1.444
wruaco[notmissing]FIN-WES_24 0.789 1.268 0.496 1.256
wruaco[notmissing]FRA-LIL_21 0.337 2.970 0.110 1.028
wruaco[notmissing]FRA-STR_21 0.288 3.467 0.105 0.795
wruaco[notmissing]FRA-TOU_21 0.422 2.368 0.138 1.294
wruaco[notmissing]SWE-NSW_02 0.991 1.009 0.510 1.925
wruaco[notmissing]SWE-NSW_03 1.065 0.939 0.605 1.875
wruaco[notmissing]UNK-BEL_21 0.474 2.109 0.225 1.001
Rsquare= 0.044 (max possible= 0.928 )
Likelihood ratio test= 112 on 17 df, p=5.55e-16
Wald test = 84.8 on 17 df, p=5.28e-11
Score (logrank) test = 114 on 17 df, p=2.22e-16, Robust = 74.5 p=3.58e-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