response: (chd1 %in% c('1','2'))
covariates: dsmoker
SNP: MIA3_MG_454
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= 3061
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.1630 0.850 0.0659 0.0807 -2.0193
wz[notmissing, selectedcovariates] 0.5053 1.658 0.0757 0.1016 4.9712
wruaco[notmissing]FIN-EAS_03 0.7901 2.204 0.1284 0.1617 4.8855
wruaco[notmissing]FIN-EAS_24 0.3260 1.385 0.1139 0.1446 2.2545
wruaco[notmissing]FIN-WES_03 0.2272 1.255 0.1429 0.1774 1.2807
wruaco[notmissing]FIN-WES_24 0.1572 1.170 0.1269 0.1639 0.9590
wruaco[notmissing]FRA-LIL_21 -0.1229 0.884 0.2199 0.2921 -0.4207
wruaco[notmissing]FRA-STR_21 0.0122 1.012 0.1975 0.2546 0.0479
wruaco[notmissing]FRA-TOU_21 0.0370 1.038 0.1986 0.2604 0.1422
wruaco[notmissing]SWE-NSW_02 0.4526 1.572 0.2053 0.2554 1.7721
wruaco[notmissing]SWE-NSW_03 -0.1313 0.877 0.2941 0.3399 -0.3862
wruaco[notmissing]UNK-BEL_21 0.5226 1.686 0.1436 0.1868 2.7976
p
allele[notmissing] 4.3e-02
wz[notmissing, selectedcovariates] 6.7e-07
wruaco[notmissing]FIN-EAS_03 1.0e-06
wruaco[notmissing]FIN-EAS_24 2.4e-02
wruaco[notmissing]FIN-WES_03 2.0e-01
wruaco[notmissing]FIN-WES_24 3.4e-01
wruaco[notmissing]FRA-LIL_21 6.7e-01
wruaco[notmissing]FRA-STR_21 9.6e-01
wruaco[notmissing]FRA-TOU_21 8.9e-01
wruaco[notmissing]SWE-NSW_02 7.6e-02
wruaco[notmissing]SWE-NSW_03 7.0e-01
wruaco[notmissing]UNK-BEL_21 5.1e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.850 1.177 0.725 0.995
wz[notmissing, selectedcovariates] 1.658 0.603 1.358 2.023
wruaco[notmissing]FIN-EAS_03 2.204 0.454 1.605 3.026
wruaco[notmissing]FIN-EAS_24 1.385 0.722 1.044 1.839
wruaco[notmissing]FIN-WES_03 1.255 0.797 0.886 1.777
wruaco[notmissing]FIN-WES_24 1.170 0.855 0.849 1.614
wruaco[notmissing]FRA-LIL_21 0.884 1.131 0.499 1.568
wruaco[notmissing]FRA-STR_21 1.012 0.988 0.615 1.667
wruaco[notmissing]FRA-TOU_21 1.038 0.964 0.623 1.729
wruaco[notmissing]SWE-NSW_02 1.572 0.636 0.953 2.594
wruaco[notmissing]SWE-NSW_03 0.877 1.140 0.450 1.707
wruaco[notmissing]UNK-BEL_21 1.686 0.593 1.169 2.432
Rsquare= 0.03 (max possible= 0.995 )
Likelihood ratio test= 94.1 on 12 df, p=7.88e-15
Wald test = 55.9 on 12 df, p=1.27e-07
Score (logrank) test = 96.7 on 12 df, p=2.44e-15, Robust = 52.3 p=5.48e-07
(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