response: (death=='1')
covariates: dsmoker
SNP: MTHFD1L_MG_449
Call:
coxph(formula = y[notmissing] ~ allele[notmissing] + wz[notmissing,
selectedcovariates] + wruaco[notmissing] + strata(wsex[notmissing]) +
cluster(wdata$id[notmissing]), weights = weight[notmissing],
robust = TRUE)
n= 3507
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.1652 0.8477 0.1260 0.165 -1.000
wz[notmissing, selectedcovariates] 0.5881 1.8005 0.0667 0.091 6.460
wruaco[notmissing]FIN-EAS_03 0.0576 1.0593 0.1242 0.162 0.356
wruaco[notmissing]FIN-EAS_24 -0.1927 0.8247 0.1015 0.135 -1.427
wruaco[notmissing]FIN-WES_03 0.0430 1.0439 0.1198 0.155 0.277
wruaco[notmissing]FIN-WES_24 -0.2213 0.8015 0.1067 0.148 -1.498
wruaco[notmissing]FRA-LIL_21 -2.4551 0.0859 0.5828 0.614 -3.996
wruaco[notmissing]FRA-STR_21 -0.2554 0.7746 0.1968 0.254 -1.006
wruaco[notmissing]FRA-TOU_21 -0.5220 0.5933 0.2313 0.279 -1.874
wruaco[notmissing]SWE-NSW_02 -0.2433 0.7841 0.1672 0.217 -1.121
wruaco[notmissing]SWE-NSW_03 -0.6169 0.5396 0.1804 0.218 -2.825
wruaco[notmissing]UNK-BEL_21 -0.3061 0.7363 0.1747 0.213 -1.435
p
allele[notmissing] 3.2e-01
wz[notmissing, selectedcovariates] 1.0e-10
wruaco[notmissing]FIN-EAS_03 7.2e-01
wruaco[notmissing]FIN-EAS_24 1.5e-01
wruaco[notmissing]FIN-WES_03 7.8e-01
wruaco[notmissing]FIN-WES_24 1.3e-01
wruaco[notmissing]FRA-LIL_21 6.4e-05
wruaco[notmissing]FRA-STR_21 3.1e-01
wruaco[notmissing]FRA-TOU_21 6.1e-02
wruaco[notmissing]SWE-NSW_02 2.6e-01
wruaco[notmissing]SWE-NSW_03 4.7e-03
wruaco[notmissing]UNK-BEL_21 1.5e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.8477 1.180 0.6132 1.172
wz[notmissing, selectedcovariates] 1.8005 0.555 1.5063 2.152
wruaco[notmissing]FIN-EAS_03 1.0593 0.944 0.7714 1.454
wruaco[notmissing]FIN-EAS_24 0.8247 1.213 0.6329 1.075
wruaco[notmissing]FIN-WES_03 1.0439 0.958 0.7704 1.415
wruaco[notmissing]FIN-WES_24 0.8015 1.248 0.6000 1.071
wruaco[notmissing]FRA-LIL_21 0.0859 11.647 0.0258 0.286
wruaco[notmissing]FRA-STR_21 0.7746 1.291 0.4709 1.274
wruaco[notmissing]FRA-TOU_21 0.5933 1.685 0.3437 1.024
wruaco[notmissing]SWE-NSW_02 0.7841 1.275 0.5124 1.200
wruaco[notmissing]SWE-NSW_03 0.5396 1.853 0.3517 0.828
wruaco[notmissing]UNK-BEL_21 0.7363 1.358 0.4847 1.118
Rsquare= 0.06 (max possible= 0.998 )
Likelihood ratio test= 216 on 12 df, p=0
Wald test = 116 on 12 df, p=0
Score (logrank) test = 196 on 12 df, p=0, Robust = 121 p=0
(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