response: (stroke3 %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= 2471
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.0623 0.940 0.0971 0.109 -0.5720
wz[notmissing, selectedcovariates] 0.3199 1.377 0.1228 0.134 2.3891
wruaco[notmissing]FIN-EAS_03 0.2200 1.246 0.2162 0.229 0.9587
wruaco[notmissing]FIN-EAS_24 -0.1355 0.873 0.1833 0.208 -0.6527
wruaco[notmissing]FIN-WES_03 -0.1547 0.857 0.2283 0.244 -0.6350
wruaco[notmissing]FIN-WES_24 -0.2562 0.774 0.1992 0.227 -1.1266
wruaco[notmissing]FRA-LIL_21 -0.8570 0.424 0.4702 0.500 -1.7150
wruaco[notmissing]FRA-STR_21 -0.8680 0.420 0.4409 0.466 -1.8629
wruaco[notmissing]FRA-TOU_21 -1.1684 0.311 0.5258 0.558 -2.0957
wruaco[notmissing]SWE-NSW_02 -0.0268 0.974 0.2766 0.311 -0.0861
wruaco[notmissing]SWE-NSW_03 0.0771 1.080 0.2572 0.281 0.2742
wruaco[notmissing]UNK-BEL_21 -0.7771 0.460 0.3585 0.378 -2.0548
p
allele[notmissing] 0.570
wz[notmissing, selectedcovariates] 0.017
wruaco[notmissing]FIN-EAS_03 0.340
wruaco[notmissing]FIN-EAS_24 0.510
wruaco[notmissing]FIN-WES_03 0.530
wruaco[notmissing]FIN-WES_24 0.260
wruaco[notmissing]FRA-LIL_21 0.086
wruaco[notmissing]FRA-STR_21 0.062
wruaco[notmissing]FRA-TOU_21 0.036
wruaco[notmissing]SWE-NSW_02 0.930
wruaco[notmissing]SWE-NSW_03 0.780
wruaco[notmissing]UNK-BEL_21 0.040
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.940 1.064 0.759 1.163
wz[notmissing, selectedcovariates] 1.377 0.726 1.059 1.790
wruaco[notmissing]FIN-EAS_03 1.246 0.803 0.795 1.954
wruaco[notmissing]FIN-EAS_24 0.873 1.145 0.581 1.312
wruaco[notmissing]FIN-WES_03 0.857 1.167 0.531 1.381
wruaco[notmissing]FIN-WES_24 0.774 1.292 0.496 1.209
wruaco[notmissing]FRA-LIL_21 0.424 2.356 0.159 1.130
wruaco[notmissing]FRA-STR_21 0.420 2.382 0.168 1.046
wruaco[notmissing]FRA-TOU_21 0.311 3.217 0.104 0.927
wruaco[notmissing]SWE-NSW_02 0.974 1.027 0.529 1.790
wruaco[notmissing]SWE-NSW_03 1.080 0.926 0.623 1.874
wruaco[notmissing]UNK-BEL_21 0.460 2.175 0.219 0.965
Rsquare= 0.015 (max possible= 0.924 )
Likelihood ratio test= 36.6 on 12 df, p=0.000259
Wald test = 30.9 on 12 df, p=0.00201
Score (logrank) test = 35.0 on 12 df, p=0.00046, Robust = 32.4 p=0.00119
(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