response: (stroke2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: CELSR2_MG_438
Call:
coxph(formula = y[notmissing] ~ allele[notmissing] + wz[notmissing,
selectedcovariates] + wrua2[notmissing] + strata(wsex[notmissing]) +
cluster(wdata$id[notmissing]), weights = weight[notmissing],
robust = TRUE)
n=2124 (34 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.39309 1.482 0.2124 0.2370
wz[notmissing, selectedcovariates]1 0.00704 1.007 0.1003 0.1158
wz[notmissing, selectedcovariates]2 -0.58895 0.555 0.3693 0.4305
wz[notmissing, selectedcovariates]3 -0.12403 0.883 0.0327 0.0392
wz[notmissing, selectedcovariates]4 0.02868 1.029 0.0076 0.0084
wz[notmissing, selectedcovariates]5 0.32274 1.381 0.2382 0.2528
wz[notmissing, selectedcovariates]6 1.06943 2.914 0.3417 0.4007
wrua2[notmissing]FRA -0.93617 0.392 0.7789 0.7454
wrua2[notmissing]SWE 0.45487 1.576 0.3725 0.3932
wrua2[notmissing]UNK 0.47525 1.608 0.6181 0.6030
z p
allele[notmissing] 1.6585 0.09700
wz[notmissing, selectedcovariates]1 0.0608 0.95000
wz[notmissing, selectedcovariates]2 -1.3681 0.17000
wz[notmissing, selectedcovariates]3 -3.1607 0.00160
wz[notmissing, selectedcovariates]4 3.4119 0.00065
wz[notmissing, selectedcovariates]5 1.2768 0.20000
wz[notmissing, selectedcovariates]6 2.6686 0.00760
wrua2[notmissing]FRA -1.2559 0.21000
wrua2[notmissing]SWE 1.1569 0.25000
wrua2[notmissing]UNK 0.7882 0.43000
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.482 0.675 0.931 2.358
wz[notmissing, selectedcovariates]1 1.007 0.993 0.803 1.264
wz[notmissing, selectedcovariates]2 0.555 1.802 0.239 1.290
wz[notmissing, selectedcovariates]3 0.883 1.132 0.818 0.954
wz[notmissing, selectedcovariates]4 1.029 0.972 1.012 1.046
wz[notmissing, selectedcovariates]5 1.381 0.724 0.841 2.266
wz[notmissing, selectedcovariates]6 2.914 0.343 1.328 6.391
wrua2[notmissing]FRA 0.392 2.550 0.091 1.690
wrua2[notmissing]SWE 1.576 0.635 0.729 3.406
wrua2[notmissing]UNK 1.608 0.622 0.493 5.244
Rsquare= 0.019 (max possible= 0.46 )
Likelihood ratio test= 40.2 on 10 df, p=1.57e-05
Wald test = 27.5 on 10 df, p=0.00217
Score (logrank) test = 40.9 on 10 df, p=1.19e-05, Robust = 23.2 p=0.00994
(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