response: (stroke2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MIA3_MG_454
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=2109 (32 observations deleted due to missingness)
coef exp(coef) se(coef) robust se z
allele[notmissing] 0.0162 1.016 0.2115 0.22220 0.073
wz[notmissing, selectedcovariates]1 0.0169 1.017 0.1001 0.11194 0.151
wz[notmissing, selectedcovariates]2 -0.6673 0.513 0.3747 0.44493 -1.500
wz[notmissing, selectedcovariates]3 -0.1172 0.889 0.0324 0.03833 -3.058
wz[notmissing, selectedcovariates]4 0.0303 1.031 0.0076 0.00824 3.678
wz[notmissing, selectedcovariates]5 0.4483 1.566 0.2388 0.24883 1.802
wz[notmissing, selectedcovariates]6 0.8828 2.418 0.3537 0.40508 2.179
wrua2[notmissing]FRA -0.9234 0.397 0.7752 0.74328 -1.242
wrua2[notmissing]SWE 0.2911 1.338 0.3896 0.41313 0.705
wrua2[notmissing]UNK 0.3466 1.414 0.6136 0.62448 0.555
p
allele[notmissing] 0.94000
wz[notmissing, selectedcovariates]1 0.88000
wz[notmissing, selectedcovariates]2 0.13000
wz[notmissing, selectedcovariates]3 0.00220
wz[notmissing, selectedcovariates]4 0.00024
wz[notmissing, selectedcovariates]5 0.07200
wz[notmissing, selectedcovariates]6 0.02900
wrua2[notmissing]FRA 0.21000
wrua2[notmissing]SWE 0.48000
wrua2[notmissing]UNK 0.58000
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.016 0.984 0.6575 1.571
wz[notmissing, selectedcovariates]1 1.017 0.983 0.8167 1.267
wz[notmissing, selectedcovariates]2 0.513 1.949 0.2145 1.227
wz[notmissing, selectedcovariates]3 0.889 1.124 0.8250 0.959
wz[notmissing, selectedcovariates]4 1.031 0.970 1.0143 1.048
wz[notmissing, selectedcovariates]5 1.566 0.639 0.9614 2.550
wz[notmissing, selectedcovariates]6 2.418 0.414 1.0929 5.348
wrua2[notmissing]FRA 0.397 2.518 0.0925 1.705
wrua2[notmissing]SWE 1.338 0.747 0.5954 3.007
wrua2[notmissing]UNK 1.414 0.707 0.4159 4.809
Rsquare= 0.017 (max possible= 0.462 )
Likelihood ratio test= 37.0 on 10 df, p=5.76e-05
Wald test = 28.6 on 10 df, p=0.00146
Score (logrank) test = 37.6 on 10 df, p=4.45e-05, Robust = 23.1 p=0.0104
(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