response: (chd1 %in% c('1','2'))
covariates: dsmoker
SNP: SMAD3_MG_452
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= 3076
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.0207 0.979 0.0533 0.067 -0.3093
wz[notmissing, selectedcovariates] 0.5005 1.650 0.0754 0.101 4.9582
wruaco[notmissing]FIN-EAS_03 0.7552 2.128 0.1288 0.162 4.6721
wruaco[notmissing]FIN-EAS_24 0.3301 1.391 0.1139 0.144 2.2861
wruaco[notmissing]FIN-WES_03 0.2554 1.291 0.1419 0.176 1.4482
wruaco[notmissing]FIN-WES_24 0.1657 1.180 0.1269 0.164 1.0117
wruaco[notmissing]FRA-LIL_21 -0.1245 0.883 0.2207 0.295 -0.4214
wruaco[notmissing]FRA-STR_21 0.0242 1.024 0.1953 0.251 0.0964
wruaco[notmissing]FRA-TOU_21 0.0532 1.055 0.1985 0.260 0.2044
wruaco[notmissing]SWE-NSW_02 0.4736 1.606 0.2052 0.255 1.8601
wruaco[notmissing]SWE-NSW_03 -0.1169 0.890 0.2942 0.341 -0.3429
wruaco[notmissing]UNK-BEL_21 0.5439 1.723 0.1440 0.186 2.9231
p
allele[notmissing] 7.6e-01
wz[notmissing, selectedcovariates] 7.1e-07
wruaco[notmissing]FIN-EAS_03 3.0e-06
wruaco[notmissing]FIN-EAS_24 2.2e-02
wruaco[notmissing]FIN-WES_03 1.5e-01
wruaco[notmissing]FIN-WES_24 3.1e-01
wruaco[notmissing]FRA-LIL_21 6.7e-01
wruaco[notmissing]FRA-STR_21 9.2e-01
wruaco[notmissing]FRA-TOU_21 8.4e-01
wruaco[notmissing]SWE-NSW_02 6.3e-02
wruaco[notmissing]SWE-NSW_03 7.3e-01
wruaco[notmissing]UNK-BEL_21 3.5e-03
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.979 1.021 0.859 1.12
wz[notmissing, selectedcovariates] 1.650 0.606 1.353 2.01
wruaco[notmissing]FIN-EAS_03 2.128 0.470 1.550 2.92
wruaco[notmissing]FIN-EAS_24 1.391 0.719 1.048 1.85
wruaco[notmissing]FIN-WES_03 1.291 0.775 0.914 1.82
wruaco[notmissing]FIN-WES_24 1.180 0.847 0.856 1.63
wruaco[notmissing]FRA-LIL_21 0.883 1.133 0.495 1.58
wruaco[notmissing]FRA-STR_21 1.024 0.976 0.627 1.68
wruaco[notmissing]FRA-TOU_21 1.055 0.948 0.633 1.76
wruaco[notmissing]SWE-NSW_02 1.606 0.623 0.975 2.64
wruaco[notmissing]SWE-NSW_03 0.890 1.124 0.456 1.74
wruaco[notmissing]UNK-BEL_21 1.723 0.580 1.196 2.48
Rsquare= 0.027 (max possible= 0.995 )
Likelihood ratio test= 84 on 12 df, p=7.04e-13
Wald test = 49.7 on 12 df, p=1.60e-06
Score (logrank) test = 86.2 on 12 df, p=2.68e-13, Robust = 46.4 p=5.87e-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