response: (stroke3 %in% c('1','2'))
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
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=2443 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.215406 0.806 0.24762 0.28046
wz[notmissing, selectedcovariates]1 0.063435 1.065 0.04470 0.05293
wz[notmissing, selectedcovariates]2 -0.750162 0.472 0.17739 0.20161
wz[notmissing, selectedcovariates]3 -0.000882 0.999 0.01246 0.01382
wz[notmissing, selectedcovariates]4 0.016279 1.016 0.00372 0.00429
wz[notmissing, selectedcovariates]5 0.371067 1.449 0.12572 0.13874
wz[notmissing, selectedcovariates]6 0.738812 2.093 0.16131 0.19933
wruaco[notmissing]FIN-EAS_03 0.251508 1.286 0.21832 0.23733
wruaco[notmissing]FIN-EAS_24 -0.204327 0.815 0.19341 0.22229
wruaco[notmissing]FIN-WES_03 -0.196992 0.821 0.23271 0.24922
wruaco[notmissing]FIN-WES_24 -0.230898 0.794 0.20671 0.24322
wruaco[notmissing]FRA-LIL_21 -1.008506 0.365 0.52058 0.56968
wruaco[notmissing]FRA-STR_21 -1.128979 0.323 0.47823 0.51909
wruaco[notmissing]FRA-TOU_21 -0.776865 0.460 0.52982 0.57127
wruaco[notmissing]SWE-NSW_02 0.012094 1.012 0.28414 0.33841
wruaco[notmissing]SWE-NSW_03 0.046054 1.047 0.26234 0.29240
wruaco[notmissing]UNK-BEL_21 -0.692441 0.500 0.36109 0.38292
z p
allele[notmissing] -0.7681 0.44000
wz[notmissing, selectedcovariates]1 1.1984 0.23000
wz[notmissing, selectedcovariates]2 -3.7208 0.00020
wz[notmissing, selectedcovariates]3 -0.0638 0.95000
wz[notmissing, selectedcovariates]4 3.7929 0.00015
wz[notmissing, selectedcovariates]5 2.6745 0.00750
wz[notmissing, selectedcovariates]6 3.7065 0.00021
wruaco[notmissing]FIN-EAS_03 1.0598 0.29000
wruaco[notmissing]FIN-EAS_24 -0.9192 0.36000
wruaco[notmissing]FIN-WES_03 -0.7904 0.43000
wruaco[notmissing]FIN-WES_24 -0.9494 0.34000
wruaco[notmissing]FRA-LIL_21 -1.7703 0.07700
wruaco[notmissing]FRA-STR_21 -2.1749 0.03000
wruaco[notmissing]FRA-TOU_21 -1.3599 0.17000
wruaco[notmissing]SWE-NSW_02 0.0357 0.97000
wruaco[notmissing]SWE-NSW_03 0.1575 0.87000
wruaco[notmissing]UNK-BEL_21 -1.8083 0.07100
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.806 1.240 0.465 1.397
wz[notmissing, selectedcovariates]1 1.065 0.939 0.960 1.182
wz[notmissing, selectedcovariates]2 0.472 2.117 0.318 0.701
wz[notmissing, selectedcovariates]3 0.999 1.001 0.972 1.027
wz[notmissing, selectedcovariates]4 1.016 0.984 1.008 1.025
wz[notmissing, selectedcovariates]5 1.449 0.690 1.104 1.902
wz[notmissing, selectedcovariates]6 2.093 0.478 1.416 3.094
wruaco[notmissing]FIN-EAS_03 1.286 0.778 0.808 2.048
wruaco[notmissing]FIN-EAS_24 0.815 1.227 0.527 1.260
wruaco[notmissing]FIN-WES_03 0.821 1.218 0.504 1.338
wruaco[notmissing]FIN-WES_24 0.794 1.260 0.493 1.279
wruaco[notmissing]FRA-LIL_21 0.365 2.742 0.119 1.114
wruaco[notmissing]FRA-STR_21 0.323 3.092 0.117 0.894
wruaco[notmissing]FRA-TOU_21 0.460 2.175 0.150 1.409
wruaco[notmissing]SWE-NSW_02 1.012 0.988 0.521 1.965
wruaco[notmissing]SWE-NSW_03 1.047 0.955 0.590 1.857
wruaco[notmissing]UNK-BEL_21 0.500 1.999 0.236 1.060
Rsquare= 0.043 (max possible= 0.922 )
Likelihood ratio test= 107 on 17 df, p=5.22e-15
Wald test = 82.8 on 17 df, p=1.19e-10
Score (logrank) test = 107 on 17 df, p=3.77e-15, Robust = 73.8 p=4.66e-09
(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