response: (stroke3 %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= 2485
coef exp(coef) se(coef) robust se z
allele[notmissing] 0.0241 1.024 0.0839 0.0948 0.254
wz[notmissing, selectedcovariates] 0.3108 1.365 0.1223 0.1326 2.345
wruaco[notmissing]FIN-EAS_03 0.2262 1.254 0.2144 0.2287 0.989
wruaco[notmissing]FIN-EAS_24 -0.1274 0.880 0.1828 0.2066 -0.617
wruaco[notmissing]FIN-WES_03 -0.2039 0.816 0.2315 0.2459 -0.829
wruaco[notmissing]FIN-WES_24 -0.2268 0.797 0.1974 0.2258 -1.005
wruaco[notmissing]FRA-LIL_21 -0.8762 0.416 0.4710 0.5026 -1.744
wruaco[notmissing]FRA-STR_21 -0.9123 0.402 0.4402 0.4646 -1.963
wruaco[notmissing]FRA-TOU_21 -1.1642 0.312 0.5254 0.5565 -2.092
wruaco[notmissing]SWE-NSW_02 -0.0321 0.968 0.2770 0.3100 -0.103
wruaco[notmissing]SWE-NSW_03 0.0789 1.082 0.2568 0.2786 0.283
wruaco[notmissing]UNK-BEL_21 -0.7862 0.456 0.3587 0.3760 -2.091
p
allele[notmissing] 0.800
wz[notmissing, selectedcovariates] 0.019
wruaco[notmissing]FIN-EAS_03 0.320
wruaco[notmissing]FIN-EAS_24 0.540
wruaco[notmissing]FIN-WES_03 0.410
wruaco[notmissing]FIN-WES_24 0.320
wruaco[notmissing]FRA-LIL_21 0.081
wruaco[notmissing]FRA-STR_21 0.050
wruaco[notmissing]FRA-TOU_21 0.036
wruaco[notmissing]SWE-NSW_02 0.920
wruaco[notmissing]SWE-NSW_03 0.780
wruaco[notmissing]UNK-BEL_21 0.037
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 1.024 0.976 0.851 1.234
wz[notmissing, selectedcovariates] 1.365 0.733 1.052 1.769
wruaco[notmissing]FIN-EAS_03 1.254 0.798 0.801 1.963
wruaco[notmissing]FIN-EAS_24 0.880 1.136 0.587 1.320
wruaco[notmissing]FIN-WES_03 0.816 1.226 0.504 1.320
wruaco[notmissing]FIN-WES_24 0.797 1.255 0.512 1.241
wruaco[notmissing]FRA-LIL_21 0.416 2.402 0.155 1.115
wruaco[notmissing]FRA-STR_21 0.402 2.490 0.162 0.998
wruaco[notmissing]FRA-TOU_21 0.312 3.203 0.105 0.929
wruaco[notmissing]SWE-NSW_02 0.968 1.033 0.527 1.778
wruaco[notmissing]SWE-NSW_03 1.082 0.924 0.627 1.868
wruaco[notmissing]UNK-BEL_21 0.456 2.195 0.218 0.952
Rsquare= 0.015 (max possible= 0.924 )
Likelihood ratio test= 36.7 on 12 df, p=0.000252
Wald test = 31.4 on 12 df, p=0.00168
Score (logrank) test = 35.1 on 12 df, p=0.000451, Robust = 32.9 p=0.00101
(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