response: (chd1 %in% c('1','2'))
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= 3078
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.0820 0.921 0.1418 0.181 -0.4537
wz[notmissing, selectedcovariates] 0.5053 1.658 0.0753 0.101 5.0022
wruaco[notmissing]FIN-EAS_03 0.7490 2.115 0.1288 0.162 4.6361
wruaco[notmissing]FIN-EAS_24 0.3293 1.390 0.1139 0.145 2.2769
wruaco[notmissing]FIN-WES_03 0.2574 1.294 0.1420 0.176 1.4597
wruaco[notmissing]FIN-WES_24 0.1675 1.182 0.1269 0.164 1.0223
wruaco[notmissing]FRA-LIL_21 -0.1265 0.881 0.2200 0.294 -0.4307
wruaco[notmissing]FRA-STR_21 0.0201 1.020 0.1952 0.251 0.0803
wruaco[notmissing]FRA-TOU_21 0.0517 1.053 0.1985 0.260 0.1986
wruaco[notmissing]SWE-NSW_02 0.4681 1.597 0.2050 0.254 1.8421
wruaco[notmissing]SWE-NSW_03 -0.1190 0.888 0.2940 0.340 -0.3498
wruaco[notmissing]UNK-BEL_21 0.5490 1.732 0.1430 0.185 2.9596
p
allele[notmissing] 6.5e-01
wz[notmissing, selectedcovariates] 5.7e-07
wruaco[notmissing]FIN-EAS_03 3.6e-06
wruaco[notmissing]FIN-EAS_24 2.3e-02
wruaco[notmissing]FIN-WES_03 1.4e-01
wruaco[notmissing]FIN-WES_24 3.1e-01
wruaco[notmissing]FRA-LIL_21 6.7e-01
wruaco[notmissing]FRA-STR_21 9.4e-01
wruaco[notmissing]FRA-TOU_21 8.4e-01
wruaco[notmissing]SWE-NSW_02 6.5e-02
wruaco[notmissing]SWE-NSW_03 7.3e-01
wruaco[notmissing]UNK-BEL_21 3.1e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.921 1.085 0.646 1.31
wz[notmissing, selectedcovariates] 1.658 0.603 1.360 2.02
wruaco[notmissing]FIN-EAS_03 2.115 0.473 1.541 2.90
wruaco[notmissing]FIN-EAS_24 1.390 0.719 1.047 1.85
wruaco[notmissing]FIN-WES_03 1.294 0.773 0.916 1.83
wruaco[notmissing]FIN-WES_24 1.182 0.846 0.858 1.63
wruaco[notmissing]FRA-LIL_21 0.881 1.135 0.495 1.57
wruaco[notmissing]FRA-STR_21 1.020 0.980 0.624 1.67
wruaco[notmissing]FRA-TOU_21 1.053 0.950 0.632 1.75
wruaco[notmissing]SWE-NSW_02 1.597 0.626 0.971 2.63
wruaco[notmissing]SWE-NSW_03 0.888 1.126 0.456 1.73
wruaco[notmissing]UNK-BEL_21 1.732 0.578 1.204 2.49
Rsquare= 0.027 (max possible= 0.995 )
Likelihood ratio test= 85 on 12 df, p=4.42e-13
Wald test = 49.9 on 12 df, p=1.47e-06
Score (logrank) test = 87.1 on 12 df, p=1.81e-13, Robust = 46.6 p=5.39e-06
(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