response: (stroke2=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
SNP: MTHFD1L_MG_445
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 z
allele[notmissing] 0.2154 1.240 0.17320 0.18750 1.149
wz[notmissing, selectedcovariates]1 0.0278 1.028 0.09818 0.11348 0.245
wz[notmissing, selectedcovariates]2 -0.5927 0.553 0.36904 0.43239 -1.371
wz[notmissing, selectedcovariates]3 -0.1208 0.886 0.03238 0.03851 -3.136
wz[notmissing, selectedcovariates]4 0.0280 1.028 0.00751 0.00815 3.436
wz[notmissing, selectedcovariates]5 0.3574 1.430 0.23547 0.24856 1.438
wz[notmissing, selectedcovariates]6 0.9572 2.605 0.34010 0.39920 2.398
wrua2[notmissing]FRA -1.0177 0.361 0.77583 0.74653 -1.363
wrua2[notmissing]SWE 0.3720 1.451 0.37209 0.39102 0.951
wrua2[notmissing]UNK 0.3027 1.353 0.61425 0.61983 0.488
p
allele[notmissing] 0.25000
wz[notmissing, selectedcovariates]1 0.81000
wz[notmissing, selectedcovariates]2 0.17000
wz[notmissing, selectedcovariates]3 0.00170
wz[notmissing, selectedcovariates]4 0.00059
wz[notmissing, selectedcovariates]5 0.15000
wz[notmissing, selectedcovariates]6 0.01600
wrua2[notmissing]FRA 0.17000
wrua2[notmissing]SWE 0.34000
wrua2[notmissing]UNK 0.63000
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.240 0.806 0.8589 1.791
wz[notmissing, selectedcovariates]1 1.028 0.973 0.8231 1.284
wz[notmissing, selectedcovariates]2 0.553 1.809 0.2369 1.290
wz[notmissing, selectedcovariates]3 0.886 1.128 0.8218 0.956
wz[notmissing, selectedcovariates]4 1.028 0.972 1.0121 1.045
wz[notmissing, selectedcovariates]5 1.430 0.699 0.8783 2.327
wz[notmissing, selectedcovariates]6 2.605 0.384 1.1910 5.695
wrua2[notmissing]FRA 0.361 2.767 0.0837 1.561
wrua2[notmissing]SWE 1.451 0.689 0.6741 3.122
wrua2[notmissing]UNK 1.353 0.739 0.4016 4.561
Rsquare= 0.018 (max possible= 0.468 )
Likelihood ratio test= 39.6 on 10 df, p=1.95e-05
Wald test = 32.6 on 10 df, p=0.000312
Score (logrank) test = 41.4 on 10 df, p=9.73e-06, Robust = 25.7 p=0.00418
(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