response: (death=='1')
covariates: nonhdl hdla bmi bpm dsmoker hisdiab1
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=3441 (65 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.05503 0.946 0.04699 0.06128
wz[notmissing, selectedcovariates]1 -0.04453 0.956 0.02556 0.03372
wz[notmissing, selectedcovariates]2 -0.22388 0.799 0.08973 0.12135
wz[notmissing, selectedcovariates]3 -0.01509 0.985 0.00710 0.00966
wz[notmissing, selectedcovariates]4 0.00946 1.010 0.00202 0.00273
wz[notmissing, selectedcovariates]5 0.59567 1.814 0.06841 0.09421
wz[notmissing, selectedcovariates]6 0.54822 1.730 0.09699 0.14297
wruaco[notmissing]FIN-EAS_03 0.02838 1.029 0.12680 0.16653
wruaco[notmissing]FIN-EAS_24 -0.32906 0.720 0.10732 0.14248
wruaco[notmissing]FIN-WES_03 -0.00890 0.991 0.12086 0.15662
wruaco[notmissing]FIN-WES_24 -0.29475 0.745 0.11175 0.15569
wruaco[notmissing]FRA-LIL_21 -2.81325 0.060 0.71201 0.74661
wruaco[notmissing]FRA-STR_21 -0.27968 0.756 0.19747 0.26027
wruaco[notmissing]FRA-TOU_21 -0.52816 0.590 0.24854 0.29722
wruaco[notmissing]SWE-NSW_02 -0.22436 0.799 0.17266 0.22937
wruaco[notmissing]SWE-NSW_03 -0.64890 0.523 0.18420 0.22918
wruaco[notmissing]UNK-BEL_21 -0.30877 0.734 0.17677 0.21940
z p
allele[notmissing] -0.8979 3.7e-01
wz[notmissing, selectedcovariates]1 -1.3206 1.9e-01
wz[notmissing, selectedcovariates]2 -1.8450 6.5e-02
wz[notmissing, selectedcovariates]3 -1.5625 1.2e-01
wz[notmissing, selectedcovariates]4 3.4639 5.3e-04
wz[notmissing, selectedcovariates]5 6.3230 2.6e-10
wz[notmissing, selectedcovariates]6 3.8345 1.3e-04
wruaco[notmissing]FIN-EAS_03 0.1704 8.6e-01
wruaco[notmissing]FIN-EAS_24 -2.3095 2.1e-02
wruaco[notmissing]FIN-WES_03 -0.0568 9.5e-01
wruaco[notmissing]FIN-WES_24 -1.8932 5.8e-02
wruaco[notmissing]FRA-LIL_21 -3.7680 1.6e-04
wruaco[notmissing]FRA-STR_21 -1.0746 2.8e-01
wruaco[notmissing]FRA-TOU_21 -1.7770 7.6e-02
wruaco[notmissing]SWE-NSW_02 -0.9782 3.3e-01
wruaco[notmissing]SWE-NSW_03 -2.8314 4.6e-03
wruaco[notmissing]UNK-BEL_21 -1.4073 1.6e-01
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.946 1.057 0.8393 1.067
wz[notmissing, selectedcovariates]1 0.956 1.046 0.8953 1.022
wz[notmissing, selectedcovariates]2 0.799 1.251 0.6302 1.014
wz[notmissing, selectedcovariates]3 0.985 1.015 0.9665 1.004
wz[notmissing, selectedcovariates]4 1.010 0.991 1.0041 1.015
wz[notmissing, selectedcovariates]5 1.814 0.551 1.5084 2.182
wz[notmissing, selectedcovariates]6 1.730 0.578 1.3074 2.290
wruaco[notmissing]FIN-EAS_03 1.029 0.972 0.7423 1.426
wruaco[notmissing]FIN-EAS_24 0.720 1.390 0.5443 0.951
wruaco[notmissing]FIN-WES_03 0.991 1.009 0.7292 1.347
wruaco[notmissing]FIN-WES_24 0.745 1.343 0.5489 1.010
wruaco[notmissing]FRA-LIL_21 0.060 16.664 0.0139 0.259
wruaco[notmissing]FRA-STR_21 0.756 1.323 0.4539 1.259
wruaco[notmissing]FRA-TOU_21 0.590 1.696 0.3293 1.056
wruaco[notmissing]SWE-NSW_02 0.799 1.252 0.5097 1.253
wruaco[notmissing]SWE-NSW_03 0.523 1.913 0.3335 0.819
wruaco[notmissing]UNK-BEL_21 0.734 1.362 0.4777 1.129
Rsquare= 0.077 (max possible= 0.998 )
Likelihood ratio test= 276 on 17 df, p=0
Wald test = 142 on 17 df, p=0
Score (logrank) test = 256 on 17 df, p=0, Robust = 149 p=0
(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