response: (stroke4 %in% c('1','2','3','4','5','6'))
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= 2482
coef exp(coef) se(coef) robust se z
allele[notmissing] -0.0568 0.945 0.0957 0.107 -0.530
wz[notmissing, selectedcovariates] 0.3061 1.358 0.1213 0.134 2.287
wruaco[notmissing]FIN-EAS_03 0.2166 1.242 0.2132 0.229 0.948
wruaco[notmissing]FIN-EAS_24 -0.1701 0.844 0.1822 0.207 -0.820
wruaco[notmissing]FIN-WES_03 -0.0758 0.927 0.2200 0.238 -0.319
wruaco[notmissing]FIN-WES_24 -0.2867 0.751 0.1983 0.227 -1.262
wruaco[notmissing]FRA-LIL_21 -0.9059 0.404 0.4694 0.499 -1.816
wruaco[notmissing]FRA-STR_21 -0.9110 0.402 0.4399 0.465 -1.959
wruaco[notmissing]FRA-TOU_21 -1.2103 0.298 0.5250 0.557 -2.174
wruaco[notmissing]SWE-NSW_02 -0.0691 0.933 0.2763 0.313 -0.221
wruaco[notmissing]SWE-NSW_03 0.0473 1.048 0.2564 0.281 0.168
wruaco[notmissing]UNK-BEL_21 -0.8152 0.443 0.3572 0.377 -2.163
p
allele[notmissing] 0.600
wz[notmissing, selectedcovariates] 0.022
wruaco[notmissing]FIN-EAS_03 0.340
wruaco[notmissing]FIN-EAS_24 0.410
wruaco[notmissing]FIN-WES_03 0.750
wruaco[notmissing]FIN-WES_24 0.210
wruaco[notmissing]FRA-LIL_21 0.069
wruaco[notmissing]FRA-STR_21 0.050
wruaco[notmissing]FRA-TOU_21 0.030
wruaco[notmissing]SWE-NSW_02 0.830
wruaco[notmissing]SWE-NSW_03 0.870
wruaco[notmissing]UNK-BEL_21 0.031
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.945 1.058 0.766 1.166
wz[notmissing, selectedcovariates] 1.358 0.736 1.045 1.766
wruaco[notmissing]FIN-EAS_03 1.242 0.805 0.793 1.943
wruaco[notmissing]FIN-EAS_24 0.844 1.185 0.562 1.267
wruaco[notmissing]FIN-WES_03 0.927 1.079 0.581 1.478
wruaco[notmissing]FIN-WES_24 0.751 1.332 0.481 1.172
wruaco[notmissing]FRA-LIL_21 0.404 2.474 0.152 1.074
wruaco[notmissing]FRA-STR_21 0.402 2.487 0.162 1.001
wruaco[notmissing]FRA-TOU_21 0.298 3.354 0.100 0.888
wruaco[notmissing]SWE-NSW_02 0.933 1.072 0.506 1.722
wruaco[notmissing]SWE-NSW_03 1.048 0.954 0.604 1.820
wruaco[notmissing]UNK-BEL_21 0.443 2.260 0.211 0.926
Rsquare= 0.015 (max possible= 0.928 )
Likelihood ratio test= 38.6 on 12 df, p=0.000124
Wald test = 31.9 on 12 df, p=0.00141
Score (logrank) test = 36.7 on 12 df, p=0.000246, Robust = 33.7 p=0.00076
(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