response: (chd2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: SMAD3_MG_452
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=2357 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.02228 1.023 0.0960 0.10863
wz[notmissing, selectedcovariates]1 0.09929 1.104 0.0479 0.05269
wz[notmissing, selectedcovariates]2 -0.73981 0.477 0.2000 0.23534
wz[notmissing, selectedcovariates]3 0.00511 1.005 0.0140 0.01665
wz[notmissing, selectedcovariates]4 0.01629 1.016 0.0041 0.00453
wz[notmissing, selectedcovariates]5 0.76992 2.160 0.1278 0.15035
wz[notmissing, selectedcovariates]6 0.87881 2.408 0.1743 0.21003
wrua2[notmissing]FRA -1.29038 0.275 0.4038 0.43371
wrua2[notmissing]SWE -0.46935 0.625 0.2707 0.29192
wrua2[notmissing]UNK -0.17176 0.842 0.3232 0.35870
z p
allele[notmissing] 0.205 8.4e-01
wz[notmissing, selectedcovariates]1 1.885 5.9e-02
wz[notmissing, selectedcovariates]2 -3.144 1.7e-03
wz[notmissing, selectedcovariates]3 0.307 7.6e-01
wz[notmissing, selectedcovariates]4 3.592 3.3e-04
wz[notmissing, selectedcovariates]5 5.121 3.0e-07
wz[notmissing, selectedcovariates]6 4.184 2.9e-05
wrua2[notmissing]FRA -2.975 2.9e-03
wrua2[notmissing]SWE -1.608 1.1e-01
wrua2[notmissing]UNK -0.479 6.3e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.023 0.978 0.826 1.265
wz[notmissing, selectedcovariates]1 1.104 0.905 0.996 1.225
wz[notmissing, selectedcovariates]2 0.477 2.096 0.301 0.757
wz[notmissing, selectedcovariates]3 1.005 0.995 0.973 1.038
wz[notmissing, selectedcovariates]4 1.016 0.984 1.007 1.025
wz[notmissing, selectedcovariates]5 2.160 0.463 1.608 2.900
wz[notmissing, selectedcovariates]6 2.408 0.415 1.595 3.634
wrua2[notmissing]FRA 0.275 3.634 0.118 0.644
wrua2[notmissing]SWE 0.625 1.599 0.353 1.108
wrua2[notmissing]UNK 0.842 1.187 0.417 1.701
Rsquare= 0.052 (max possible= 0.875 )
Likelihood ratio test= 126 on 10 df, p=0
Wald test = 92 on 10 df, p=2.11e-15
Score (logrank) test = 126 on 10 df, p=0, Robust = 80.3 p=4.42e-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