response: (stroke4 %in% c('1','2','3','4','5','6'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MIA3_MG_454
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=2439 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -4.70e-02 0.954 0.09701 0.10916
wz[notmissing, selectedcovariates]1 7.30e-02 1.076 0.04418 0.05258
wz[notmissing, selectedcovariates]2 -7.09e-01 0.492 0.17537 0.20055
wz[notmissing, selectedcovariates]3 6.59e-05 1.000 0.01225 0.01357
wz[notmissing, selectedcovariates]4 1.63e-02 1.016 0.00369 0.00425
wz[notmissing, selectedcovariates]5 3.62e-01 1.437 0.12453 0.13996
wz[notmissing, selectedcovariates]6 7.74e-01 2.167 0.15851 0.19709
wruaco[notmissing]FIN-EAS_03 2.54e-01 1.289 0.21701 0.23725
wruaco[notmissing]FIN-EAS_24 -2.50e-01 0.779 0.19286 0.22374
wruaco[notmissing]FIN-WES_03 -7.51e-02 0.928 0.22086 0.24142
wruaco[notmissing]FIN-WES_24 -3.06e-01 0.737 0.20774 0.24449
wruaco[notmissing]FRA-LIL_21 -1.06e+00 0.346 0.51990 0.56969
wruaco[notmissing]FRA-STR_21 -1.14e+00 0.319 0.47794 0.52144
wruaco[notmissing]FRA-TOU_21 -8.07e-01 0.446 0.52932 0.57182
wruaco[notmissing]SWE-NSW_02 -4.62e-02 0.955 0.28466 0.34432
wruaco[notmissing]SWE-NSW_03 8.61e-04 1.001 0.26218 0.29466
wruaco[notmissing]UNK-BEL_21 -7.26e-01 0.484 0.35991 0.38370
z p
allele[notmissing] -0.43078 6.7e-01
wz[notmissing, selectedcovariates]1 1.38897 1.6e-01
wz[notmissing, selectedcovariates]2 -3.53485 4.1e-04
wz[notmissing, selectedcovariates]3 0.00485 1.0e+00
wz[notmissing, selectedcovariates]4 3.84357 1.2e-04
wz[notmissing, selectedcovariates]5 2.58906 9.6e-03
wz[notmissing, selectedcovariates]6 3.92480 8.7e-05
wruaco[notmissing]FIN-EAS_03 1.07163 2.8e-01
wruaco[notmissing]FIN-EAS_24 -1.11762 2.6e-01
wruaco[notmissing]FIN-WES_03 -0.31107 7.6e-01
wruaco[notmissing]FIN-WES_24 -1.24955 2.1e-01
wruaco[notmissing]FRA-LIL_21 -1.86332 6.2e-02
wruaco[notmissing]FRA-STR_21 -2.18924 2.9e-02
wruaco[notmissing]FRA-TOU_21 -1.41158 1.6e-01
wruaco[notmissing]SWE-NSW_02 -0.13421 8.9e-01
wruaco[notmissing]SWE-NSW_03 0.00292 1.0e+00
wruaco[notmissing]UNK-BEL_21 -1.89294 5.8e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.954 1.048 0.770 1.182
wz[notmissing, selectedcovariates]1 1.076 0.930 0.970 1.193
wz[notmissing, selectedcovariates]2 0.492 2.032 0.332 0.729
wz[notmissing, selectedcovariates]3 1.000 1.000 0.974 1.027
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.437 0.696 1.092 1.890
wz[notmissing, selectedcovariates]6 2.167 0.461 1.473 3.189
wruaco[notmissing]FIN-EAS_03 1.289 0.775 0.810 2.053
wruaco[notmissing]FIN-EAS_24 0.779 1.284 0.502 1.207
wruaco[notmissing]FIN-WES_03 0.928 1.078 0.578 1.489
wruaco[notmissing]FIN-WES_24 0.737 1.357 0.456 1.190
wruaco[notmissing]FRA-LIL_21 0.346 2.891 0.113 1.057
wruaco[notmissing]FRA-STR_21 0.319 3.132 0.115 0.887
wruaco[notmissing]FRA-TOU_21 0.446 2.242 0.145 1.368
wruaco[notmissing]SWE-NSW_02 0.955 1.047 0.486 1.875
wruaco[notmissing]SWE-NSW_03 1.001 0.999 0.562 1.783
wruaco[notmissing]UNK-BEL_21 0.484 2.067 0.228 1.026
Rsquare= 0.044 (max possible= 0.926 )
Likelihood ratio test= 111 on 17 df, p=1e-15
Wald test = 83.3 on 17 df, p=9.7e-11
Score (logrank) test = 112 on 17 df, p=4.44e-16, Robust = 73.3 p=5.7e-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