response: (stroke1 %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=2556 (43 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.28681 0.751 0.22313 0.25812
wz[notmissing, selectedcovariates]1 0.02891 1.029 0.04043 0.04842
wz[notmissing, selectedcovariates]2 -0.78775 0.455 0.15661 0.18466
wz[notmissing, selectedcovariates]3 -0.01725 0.983 0.01135 0.01309
wz[notmissing, selectedcovariates]4 0.02044 1.021 0.00327 0.00391
wz[notmissing, selectedcovariates]5 0.38391 1.468 0.11065 0.12920
wz[notmissing, selectedcovariates]6 0.67242 1.959 0.14670 0.19149
wruaco[notmissing]FIN-EAS_03 0.21360 1.238 0.19402 0.22214
wruaco[notmissing]FIN-EAS_24 -0.23570 0.790 0.17361 0.20437
wruaco[notmissing]FIN-WES_03 0.08627 1.090 0.18895 0.21310
wruaco[notmissing]FIN-WES_24 -0.17245 0.842 0.18108 0.22075
wruaco[notmissing]FRA-LIL_21 -1.27594 0.279 0.51655 0.56905
wruaco[notmissing]FRA-STR_21 -0.89035 0.411 0.38411 0.42310
wruaco[notmissing]FRA-TOU_21 -1.02669 0.358 0.52336 0.56258
wruaco[notmissing]SWE-NSW_02 0.00371 1.004 0.24902 0.32160
wruaco[notmissing]SWE-NSW_03 0.06053 1.062 0.23520 0.26945
wruaco[notmissing]UNK-BEL_21 -0.69287 0.500 0.31718 0.33922
z p
allele[notmissing] -1.1112 2.7e-01
wz[notmissing, selectedcovariates]1 0.5970 5.5e-01
wz[notmissing, selectedcovariates]2 -4.2660 2.0e-05
wz[notmissing, selectedcovariates]3 -1.3171 1.9e-01
wz[notmissing, selectedcovariates]4 5.2214 1.8e-07
wz[notmissing, selectedcovariates]5 2.9714 3.0e-03
wz[notmissing, selectedcovariates]6 3.5115 4.5e-04
wruaco[notmissing]FIN-EAS_03 0.9615 3.4e-01
wruaco[notmissing]FIN-EAS_24 -1.1533 2.5e-01
wruaco[notmissing]FIN-WES_03 0.4048 6.9e-01
wruaco[notmissing]FIN-WES_24 -0.7812 4.3e-01
wruaco[notmissing]FRA-LIL_21 -2.2422 2.5e-02
wruaco[notmissing]FRA-STR_21 -2.1043 3.5e-02
wruaco[notmissing]FRA-TOU_21 -1.8250 6.8e-02
wruaco[notmissing]SWE-NSW_02 0.0115 9.9e-01
wruaco[notmissing]SWE-NSW_03 0.2246 8.2e-01
wruaco[notmissing]UNK-BEL_21 -2.0425 4.1e-02
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.751 1.332 0.4526 1.245
wz[notmissing, selectedcovariates]1 1.029 0.972 0.9361 1.132
wz[notmissing, selectedcovariates]2 0.455 2.198 0.3167 0.653
wz[notmissing, selectedcovariates]3 0.983 1.017 0.9580 1.008
wz[notmissing, selectedcovariates]4 1.021 0.980 1.0128 1.029
wz[notmissing, selectedcovariates]5 1.468 0.681 1.1396 1.891
wz[notmissing, selectedcovariates]6 1.959 0.510 1.3460 2.851
wruaco[notmissing]FIN-EAS_03 1.238 0.808 0.8011 1.914
wruaco[notmissing]FIN-EAS_24 0.790 1.266 0.5293 1.179
wruaco[notmissing]FIN-WES_03 1.090 0.917 0.7179 1.655
wruaco[notmissing]FIN-WES_24 0.842 1.188 0.5460 1.297
wruaco[notmissing]FRA-LIL_21 0.279 3.582 0.0915 0.852
wruaco[notmissing]FRA-STR_21 0.411 2.436 0.1791 0.941
wruaco[notmissing]FRA-TOU_21 0.358 2.792 0.1189 1.079
wruaco[notmissing]SWE-NSW_02 1.004 0.996 0.5344 1.885
wruaco[notmissing]SWE-NSW_03 1.062 0.941 0.6265 1.802
wruaco[notmissing]UNK-BEL_21 0.500 1.999 0.2572 0.972
Rsquare= 0.053 (max possible= 0.955 )
Likelihood ratio test= 139 on 17 df, p=0
Wald test = 97.4 on 17 df, p=2.69e-13
Score (logrank) test = 137 on 17 df, p=0, Robust = 91.3 p=3.47e-12
(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