response: (stroke3 %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= 2486
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.2029 0.816 0.241 0.275 -0.738
wz[notmissing, selectedcovariates] 0.3138 1.369 0.122 0.133 2.365
wruaco[notmissing]FIN-EAS_03 0.2175 1.243 0.214 0.228 0.952
wruaco[notmissing]FIN-EAS_24 -0.1307 0.877 0.183 0.207 -0.633
wruaco[notmissing]FIN-WES_03 -0.2043 0.815 0.232 0.246 -0.830
wruaco[notmissing]FIN-WES_24 -0.2213 0.801 0.197 0.227 -0.977
wruaco[notmissing]FRA-LIL_21 -0.8594 0.423 0.470 0.499 -1.721
wruaco[notmissing]FRA-STR_21 -0.9144 0.401 0.440 0.464 -1.973
wruaco[notmissing]FRA-TOU_21 -1.1735 0.309 0.525 0.557 -2.109
wruaco[notmissing]SWE-NSW_02 -0.0290 0.971 0.276 0.308 -0.094
wruaco[notmissing]SWE-NSW_03 0.0790 1.082 0.257 0.279 0.283
wruaco[notmissing]UNK-BEL_21 -0.7772 0.460 0.358 0.376 -2.070
p
allele[notmissing] 0.460
wz[notmissing, selectedcovariates] 0.018
wruaco[notmissing]FIN-EAS_03 0.340
wruaco[notmissing]FIN-EAS_24 0.530
wruaco[notmissing]FIN-WES_03 0.410
wruaco[notmissing]FIN-WES_24 0.330
wruaco[notmissing]FRA-LIL_21 0.085
wruaco[notmissing]FRA-STR_21 0.049
wruaco[notmissing]FRA-TOU_21 0.035
wruaco[notmissing]SWE-NSW_02 0.930
wruaco[notmissing]SWE-NSW_03 0.780
wruaco[notmissing]UNK-BEL_21 0.038
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.816 1.225 0.476 1.399
wz[notmissing, selectedcovariates] 1.369 0.731 1.055 1.775
wruaco[notmissing]FIN-EAS_03 1.243 0.805 0.794 1.945
wruaco[notmissing]FIN-EAS_24 0.877 1.140 0.585 1.315
wruaco[notmissing]FIN-WES_03 0.815 1.227 0.503 1.321
wruaco[notmissing]FIN-WES_24 0.801 1.248 0.514 1.250
wruaco[notmissing]FRA-LIL_21 0.423 2.362 0.159 1.127
wruaco[notmissing]FRA-STR_21 0.401 2.495 0.162 0.994
wruaco[notmissing]FRA-TOU_21 0.309 3.233 0.104 0.921
wruaco[notmissing]SWE-NSW_02 0.971 1.029 0.531 1.776
wruaco[notmissing]SWE-NSW_03 1.082 0.924 0.626 1.871
wruaco[notmissing]UNK-BEL_21 0.460 2.175 0.220 0.960
Rsquare= 0.015 (max possible= 0.924 )
Likelihood ratio test= 37.3 on 12 df, p=0.000198
Wald test = 31.8 on 12 df, p=0.00146
Score (logrank) test = 35.7 on 12 df, p=0.000360, Robust = 33.2 p=0.000887
(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