response: (stroke3 %in% c('1','2'))
covariates: dsmoker
SNP: CELSR2_MG_438
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= 2494
coef exp(coef) se(coef) robust se z
allele[notmissing] 0.0108 1.011 0.0865 0.0969 0.1115
wz[notmissing, selectedcovariates] 0.3098 1.363 0.1218 0.1327 2.3349
wruaco[notmissing]FIN-EAS_03 0.2130 1.237 0.2138 0.2276 0.9360
wruaco[notmissing]FIN-EAS_24 -0.0945 0.910 0.1814 0.2053 -0.4604
wruaco[notmissing]FIN-WES_03 -0.1758 0.839 0.2283 0.2435 -0.7221
wruaco[notmissing]FIN-WES_24 -0.1877 0.829 0.1937 0.2216 -0.8469
wruaco[notmissing]FRA-LIL_21 -0.9507 0.386 0.4707 0.5005 -1.8995
wruaco[notmissing]FRA-STR_21 -0.9980 0.369 0.4391 0.4627 -2.1570
wruaco[notmissing]FRA-TOU_21 -1.2195 0.295 0.5260 0.5577 -2.1866
wruaco[notmissing]SWE-NSW_02 -0.0085 0.992 0.2760 0.3089 -0.0275
wruaco[notmissing]SWE-NSW_03 0.0817 1.085 0.2563 0.2792 0.2926
wruaco[notmissing]UNK-BEL_21 -0.7847 0.456 0.3573 0.3751 -2.0918
p
allele[notmissing] 0.910
wz[notmissing, selectedcovariates] 0.020
wruaco[notmissing]FIN-EAS_03 0.350
wruaco[notmissing]FIN-EAS_24 0.650
wruaco[notmissing]FIN-WES_03 0.470
wruaco[notmissing]FIN-WES_24 0.400
wruaco[notmissing]FRA-LIL_21 0.057
wruaco[notmissing]FRA-STR_21 0.031
wruaco[notmissing]FRA-TOU_21 0.029
wruaco[notmissing]SWE-NSW_02 0.980
wruaco[notmissing]SWE-NSW_03 0.770
wruaco[notmissing]UNK-BEL_21 0.036
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.011 0.989 0.836 1.222
wz[notmissing, selectedcovariates] 1.363 0.734 1.051 1.768
wruaco[notmissing]FIN-EAS_03 1.237 0.808 0.792 1.933
wruaco[notmissing]FIN-EAS_24 0.910 1.099 0.608 1.360
wruaco[notmissing]FIN-WES_03 0.839 1.192 0.520 1.352
wruaco[notmissing]FIN-WES_24 0.829 1.206 0.537 1.280
wruaco[notmissing]FRA-LIL_21 0.386 2.588 0.145 1.031
wruaco[notmissing]FRA-STR_21 0.369 2.713 0.149 0.913
wruaco[notmissing]FRA-TOU_21 0.295 3.386 0.099 0.881
wruaco[notmissing]SWE-NSW_02 0.992 1.009 0.541 1.816
wruaco[notmissing]SWE-NSW_03 1.085 0.922 0.628 1.875
wruaco[notmissing]UNK-BEL_21 0.456 2.192 0.219 0.952
Rsquare= 0.015 (max possible= 0.925 )
Likelihood ratio test= 37.9 on 12 df, p=0.000157
Wald test = 32.4 on 12 df, p=0.00120
Score (logrank) test = 36 on 12 df, p=0.000324, Robust = 34.4 p=0.000592
(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