response: (stroke4 %in% c('2','4','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=2423 (42 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.07342 1.076 0.08639 0.10003
wz[notmissing, selectedcovariates]1 0.05744 1.059 0.04583 0.05347
wz[notmissing, selectedcovariates]2 -0.75955 0.468 0.18198 0.20580
wz[notmissing, selectedcovariates]3 0.00819 1.008 0.01245 0.01341
wz[notmissing, selectedcovariates]4 0.01533 1.015 0.00382 0.00442
wz[notmissing, selectedcovariates]5 0.39109 1.479 0.12882 0.14260
wz[notmissing, selectedcovariates]6 0.76479 2.149 0.16208 0.19978
wruaco[notmissing]FIN-EAS_03 0.36228 1.437 0.22081 0.23967
wruaco[notmissing]FIN-EAS_24 -0.10189 0.903 0.19557 0.22532
wruaco[notmissing]FIN-WES_03 -0.02273 0.978 0.22903 0.24893
wruaco[notmissing]FIN-WES_24 -0.14098 0.869 0.21069 0.24705
wruaco[notmissing]FRA-LIL_21 -0.97447 0.377 0.52201 0.57510
wruaco[notmissing]FRA-STR_21 -1.05873 0.347 0.47938 0.52215
wruaco[notmissing]FRA-TOU_21 -0.66691 0.513 0.53107 0.57363
wruaco[notmissing]SWE-NSW_02 0.02699 1.027 0.29248 0.34900
wruaco[notmissing]SWE-NSW_03 0.09894 1.104 0.26955 0.30025
wruaco[notmissing]UNK-BEL_21 -0.71775 0.488 0.37769 0.39743
z p
allele[notmissing] 0.7339 0.46000
wz[notmissing, selectedcovariates]1 1.0743 0.28000
wz[notmissing, selectedcovariates]2 -3.6908 0.00022
wz[notmissing, selectedcovariates]3 0.6104 0.54000
wz[notmissing, selectedcovariates]4 3.4661 0.00053
wz[notmissing, selectedcovariates]5 2.7425 0.00610
wz[notmissing, selectedcovariates]6 3.8281 0.00013
wruaco[notmissing]FIN-EAS_03 1.5116 0.13000
wruaco[notmissing]FIN-EAS_24 -0.4522 0.65000
wruaco[notmissing]FIN-WES_03 -0.0913 0.93000
wruaco[notmissing]FIN-WES_24 -0.5706 0.57000
wruaco[notmissing]FRA-LIL_21 -1.6944 0.09000
wruaco[notmissing]FRA-STR_21 -2.0276 0.04300
wruaco[notmissing]FRA-TOU_21 -1.1626 0.24000
wruaco[notmissing]SWE-NSW_02 0.0773 0.94000
wruaco[notmissing]SWE-NSW_03 0.3295 0.74000
wruaco[notmissing]UNK-BEL_21 -1.8060 0.07100
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.076 0.929 0.885 1.309
wz[notmissing, selectedcovariates]1 1.059 0.944 0.954 1.176
wz[notmissing, selectedcovariates]2 0.468 2.137 0.313 0.700
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.479 0.676 1.118 1.955
wz[notmissing, selectedcovariates]6 2.149 0.465 1.452 3.178
wruaco[notmissing]FIN-EAS_03 1.437 0.696 0.898 2.298
wruaco[notmissing]FIN-EAS_24 0.903 1.107 0.581 1.405
wruaco[notmissing]FIN-WES_03 0.978 1.023 0.600 1.592
wruaco[notmissing]FIN-WES_24 0.869 1.151 0.535 1.410
wruaco[notmissing]FRA-LIL_21 0.377 2.650 0.122 1.165
wruaco[notmissing]FRA-STR_21 0.347 2.883 0.125 0.965
wruaco[notmissing]FRA-TOU_21 0.513 1.948 0.167 1.580
wruaco[notmissing]SWE-NSW_02 1.027 0.973 0.518 2.036
wruaco[notmissing]SWE-NSW_03 1.104 0.906 0.613 1.989
wruaco[notmissing]UNK-BEL_21 0.488 2.050 0.224 1.063
Rsquare= 0.043 (max possible= 0.914 )
Likelihood ratio test= 105 on 17 df, p=8.66e-15
Wald test = 83.4 on 17 df, p=9.56e-11
Score (logrank) test = 108 on 17 df, p=2.78e-15, Robust = 71.9 p=1.01e-08
(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