response: (stroke4 %in% c('1','2','3','4','5','6'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: SMAD3_MG_452
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=2453 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.028879 1.029 0.08413 0.09752
wz[notmissing, selectedcovariates]1 0.063338 1.065 0.04428 0.05257
wz[notmissing, selectedcovariates]2 -0.727658 0.483 0.17477 0.19935
wz[notmissing, selectedcovariates]3 -0.000938 0.999 0.01226 0.01355
wz[notmissing, selectedcovariates]4 0.016280 1.016 0.00368 0.00423
wz[notmissing, selectedcovariates]5 0.352231 1.422 0.12394 0.13836
wz[notmissing, selectedcovariates]6 0.764880 2.149 0.15818 0.19629
wruaco[notmissing]FIN-EAS_03 0.257988 1.294 0.21557 0.23705
wruaco[notmissing]FIN-EAS_24 -0.236508 0.789 0.19205 0.22205
wruaco[notmissing]FIN-WES_03 -0.123879 0.883 0.22357 0.24311
wruaco[notmissing]FIN-WES_24 -0.265993 0.766 0.20549 0.24276
wruaco[notmissing]FRA-LIL_21 -1.080859 0.339 0.52072 0.57206
wruaco[notmissing]FRA-STR_21 -1.178022 0.308 0.47725 0.51853
wruaco[notmissing]FRA-TOU_21 -0.808015 0.446 0.52901 0.57050
wruaco[notmissing]SWE-NSW_02 -0.029746 0.971 0.28361 0.33874
wruaco[notmissing]SWE-NSW_03 0.014311 1.014 0.26164 0.29315
wruaco[notmissing]UNK-BEL_21 -0.739811 0.477 0.35993 0.38132
z p
allele[notmissing] 0.2961 7.7e-01
wz[notmissing, selectedcovariates]1 1.2049 2.3e-01
wz[notmissing, selectedcovariates]2 -3.6502 2.6e-04
wz[notmissing, selectedcovariates]3 -0.0692 9.4e-01
wz[notmissing, selectedcovariates]4 3.8443 1.2e-04
wz[notmissing, selectedcovariates]5 2.5458 1.1e-02
wz[notmissing, selectedcovariates]6 3.8966 9.8e-05
wruaco[notmissing]FIN-EAS_03 1.0883 2.8e-01
wruaco[notmissing]FIN-EAS_24 -1.0651 2.9e-01
wruaco[notmissing]FIN-WES_03 -0.5096 6.1e-01
wruaco[notmissing]FIN-WES_24 -1.0957 2.7e-01
wruaco[notmissing]FRA-LIL_21 -1.8894 5.9e-02
wruaco[notmissing]FRA-STR_21 -2.2719 2.3e-02
wruaco[notmissing]FRA-TOU_21 -1.4163 1.6e-01
wruaco[notmissing]SWE-NSW_02 -0.0878 9.3e-01
wruaco[notmissing]SWE-NSW_03 0.0488 9.6e-01
wruaco[notmissing]UNK-BEL_21 -1.9401 5.2e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.029 0.972 0.850 1.246
wz[notmissing, selectedcovariates]1 1.065 0.939 0.961 1.181
wz[notmissing, selectedcovariates]2 0.483 2.070 0.327 0.714
wz[notmissing, selectedcovariates]3 0.999 1.001 0.973 1.026
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.422 0.703 1.084 1.865
wz[notmissing, selectedcovariates]6 2.149 0.465 1.463 3.157
wruaco[notmissing]FIN-EAS_03 1.294 0.773 0.813 2.060
wruaco[notmissing]FIN-EAS_24 0.789 1.267 0.511 1.220
wruaco[notmissing]FIN-WES_03 0.883 1.132 0.549 1.423
wruaco[notmissing]FIN-WES_24 0.766 1.305 0.476 1.233
wruaco[notmissing]FRA-LIL_21 0.339 2.947 0.111 1.041
wruaco[notmissing]FRA-STR_21 0.308 3.248 0.111 0.851
wruaco[notmissing]FRA-TOU_21 0.446 2.243 0.146 1.364
wruaco[notmissing]SWE-NSW_02 0.971 1.030 0.500 1.885
wruaco[notmissing]SWE-NSW_03 1.014 0.986 0.571 1.802
wruaco[notmissing]UNK-BEL_21 0.477 2.096 0.226 1.008
Rsquare= 0.044 (max possible= 0.926 )
Likelihood ratio test= 110 on 17 df, p=1.44e-15
Wald test = 83.3 on 17 df, p=9.72e-11
Score (logrank) test = 111 on 17 df, p=7.77e-16, Robust = 73.4 p=5.47e-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