response: (chd2=='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=2363 (45 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.07215 1.075 0.10085 0.11809
wz[notmissing, selectedcovariates]1 0.10357 1.109 0.04763 0.05315
wz[notmissing, selectedcovariates]2 -0.73193 0.481 0.19917 0.23525
wz[notmissing, selectedcovariates]3 0.00615 1.006 0.01399 0.01660
wz[notmissing, selectedcovariates]4 0.01570 1.016 0.00406 0.00447
wz[notmissing, selectedcovariates]5 0.73095 2.077 0.12662 0.14927
wz[notmissing, selectedcovariates]6 0.88455 2.422 0.17452 0.21059
wrua2[notmissing]FRA -1.32906 0.265 0.40427 0.43310
wrua2[notmissing]SWE -0.54661 0.579 0.27820 0.29702
wrua2[notmissing]UNK -0.10261 0.902 0.31436 0.34797
z p
allele[notmissing] 0.611 5.4e-01
wz[notmissing, selectedcovariates]1 1.949 5.1e-02
wz[notmissing, selectedcovariates]2 -3.111 1.9e-03
wz[notmissing, selectedcovariates]3 0.371 7.1e-01
wz[notmissing, selectedcovariates]4 3.513 4.4e-04
wz[notmissing, selectedcovariates]5 4.897 9.7e-07
wz[notmissing, selectedcovariates]6 4.200 2.7e-05
wrua2[notmissing]FRA -3.069 2.1e-03
wrua2[notmissing]SWE -1.840 6.6e-02
wrua2[notmissing]UNK -0.295 7.7e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.075 0.930 0.853 1.355
wz[notmissing, selectedcovariates]1 1.109 0.902 0.999 1.231
wz[notmissing, selectedcovariates]2 0.481 2.079 0.303 0.763
wz[notmissing, selectedcovariates]3 1.006 0.994 0.974 1.039
wz[notmissing, selectedcovariates]4 1.016 0.984 1.007 1.025
wz[notmissing, selectedcovariates]5 2.077 0.481 1.550 2.783
wz[notmissing, selectedcovariates]6 2.422 0.413 1.603 3.659
wrua2[notmissing]FRA 0.265 3.777 0.113 0.619
wrua2[notmissing]SWE 0.579 1.727 0.323 1.036
wrua2[notmissing]UNK 0.902 1.108 0.456 1.785
Rsquare= 0.052 (max possible= 0.877 )
Likelihood ratio test= 125 on 10 df, p=0
Wald test = 91.3 on 10 df, p=3e-15
Score (logrank) test = 125 on 10 df, p=0, Robust = 80 p=4.84e-13
(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