response: (stroke2=='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=2121 (32 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] 0.00126 1.001 0.1857 0.20038
wz[notmissing, selectedcovariates]1 0.01579 1.016 0.1005 0.11356
wz[notmissing, selectedcovariates]2 -0.64970 0.522 0.3718 0.43476
wz[notmissing, selectedcovariates]3 -0.11651 0.890 0.0325 0.03823
wz[notmissing, selectedcovariates]4 0.03016 1.031 0.0076 0.00823
wz[notmissing, selectedcovariates]5 0.45107 1.570 0.2386 0.24728
wz[notmissing, selectedcovariates]6 0.87372 2.396 0.3524 0.40270
wrua2[notmissing]FRA -0.92354 0.397 0.7754 0.74397
wrua2[notmissing]SWE 0.31752 1.374 0.3878 0.41089
wrua2[notmissing]UNK 0.34761 1.416 0.6142 0.62228
z p
allele[notmissing] 0.0063 0.99000
wz[notmissing, selectedcovariates]1 0.1390 0.89000
wz[notmissing, selectedcovariates]2 -1.4944 0.14000
wz[notmissing, selectedcovariates]3 -3.0479 0.00230
wz[notmissing, selectedcovariates]4 3.6631 0.00025
wz[notmissing, selectedcovariates]5 1.8241 0.06800
wz[notmissing, selectedcovariates]6 2.1697 0.03000
wrua2[notmissing]FRA -1.2414 0.21000
wrua2[notmissing]SWE 0.7728 0.44000
wrua2[notmissing]UNK 0.5586 0.58000
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.001 0.999 0.6761 1.48
wz[notmissing, selectedcovariates]1 1.016 0.984 0.8132 1.27
wz[notmissing, selectedcovariates]2 0.522 1.915 0.2227 1.22
wz[notmissing, selectedcovariates]3 0.890 1.124 0.8258 0.96
wz[notmissing, selectedcovariates]4 1.031 0.970 1.0141 1.05
wz[notmissing, selectedcovariates]5 1.570 0.637 0.9670 2.55
wz[notmissing, selectedcovariates]6 2.396 0.417 1.0881 5.28
wrua2[notmissing]FRA 0.397 2.518 0.0924 1.71
wrua2[notmissing]SWE 1.374 0.728 0.6140 3.07
wrua2[notmissing]UNK 1.416 0.706 0.4181 4.79
Rsquare= 0.017 (max possible= 0.46 )
Likelihood ratio test= 36.8 on 10 df, p=6.13e-05
Wald test = 28.9 on 10 df, p=0.00132
Score (logrank) test = 37.3 on 10 df, p=5.02e-05, Robust = 23.3 p=0.00968
(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