response: (chd1=='2')
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=2742 (41 observations deleted due to missingness)
coef exp(coef) se(coef) robust se
allele[notmissing] -0.01806 0.982 0.06359 0.08139
wz[notmissing, selectedcovariates]1 0.23235 1.262 0.02891 0.05198
wz[notmissing, selectedcovariates]2 -0.96396 0.381 0.13524 0.16988
wz[notmissing, selectedcovariates]3 0.00972 1.010 0.00946 0.01133
wz[notmissing, selectedcovariates]4 0.01441 1.015 0.00269 0.00346
wz[notmissing, selectedcovariates]5 0.44942 1.567 0.08978 0.12345
wz[notmissing, selectedcovariates]6 0.79413 2.213 0.11917 0.17442
wruaco[notmissing]FIN-EAS_03 1.10909 3.032 0.15013 0.18925
wruaco[notmissing]FIN-EAS_24 0.72266 2.060 0.13893 0.17351
wruaco[notmissing]FIN-WES_03 0.54626 1.727 0.16758 0.20457
wruaco[notmissing]FIN-WES_24 0.47357 1.606 0.15632 0.21932
wruaco[notmissing]FRA-LIL_21 0.49218 1.636 0.23454 0.34950
wruaco[notmissing]FRA-STR_21 0.42137 1.524 0.21728 0.29178
wruaco[notmissing]FRA-TOU_21 0.92412 2.520 0.22537 0.29917
wruaco[notmissing]SWE-NSW_02 0.73476 2.085 0.23874 0.31558
wruaco[notmissing]SWE-NSW_03 0.22839 1.257 0.32773 0.39309
wruaco[notmissing]UNK-BEL_21 0.94975 2.585 0.16723 0.21788
z p
allele[notmissing] -0.222 8.2e-01
wz[notmissing, selectedcovariates]1 4.470 7.8e-06
wz[notmissing, selectedcovariates]2 -5.675 1.4e-08
wz[notmissing, selectedcovariates]3 0.858 3.9e-01
wz[notmissing, selectedcovariates]4 4.160 3.2e-05
wz[notmissing, selectedcovariates]5 3.640 2.7e-04
wz[notmissing, selectedcovariates]6 4.553 5.3e-06
wruaco[notmissing]FIN-EAS_03 5.861 4.6e-09
wruaco[notmissing]FIN-EAS_24 4.165 3.1e-05
wruaco[notmissing]FIN-WES_03 2.670 7.6e-03
wruaco[notmissing]FIN-WES_24 2.159 3.1e-02
wruaco[notmissing]FRA-LIL_21 1.408 1.6e-01
wruaco[notmissing]FRA-STR_21 1.444 1.5e-01
wruaco[notmissing]FRA-TOU_21 3.089 2.0e-03
wruaco[notmissing]SWE-NSW_02 2.328 2.0e-02
wruaco[notmissing]SWE-NSW_03 0.581 5.6e-01
wruaco[notmissing]UNK-BEL_21 4.359 1.3e-05
exp(coef) exp(-coef) lower .95 upper .95
allele[notmissing] 0.982 1.018 0.837 1.152
wz[notmissing, selectedcovariates]1 1.262 0.793 1.139 1.397
wz[notmissing, selectedcovariates]2 0.381 2.622 0.273 0.532
wz[notmissing, selectedcovariates]3 1.010 0.990 0.988 1.032
wz[notmissing, selectedcovariates]4 1.015 0.986 1.008 1.021
wz[notmissing, selectedcovariates]5 1.567 0.638 1.231 1.996
wz[notmissing, selectedcovariates]6 2.213 0.452 1.572 3.114
wruaco[notmissing]FIN-EAS_03 3.032 0.330 2.092 4.393
wruaco[notmissing]FIN-EAS_24 2.060 0.485 1.466 2.894
wruaco[notmissing]FIN-WES_03 1.727 0.579 1.156 2.579
wruaco[notmissing]FIN-WES_24 1.606 0.623 1.045 2.468
wruaco[notmissing]FRA-LIL_21 1.636 0.611 0.825 3.245
wruaco[notmissing]FRA-STR_21 1.524 0.656 0.860 2.700
wruaco[notmissing]FRA-TOU_21 2.520 0.397 1.402 4.529
wruaco[notmissing]SWE-NSW_02 2.085 0.480 1.123 3.870
wruaco[notmissing]SWE-NSW_03 1.257 0.796 0.582 2.715
wruaco[notmissing]UNK-BEL_21 2.585 0.387 1.687 3.962
Rsquare= 0.107 (max possible= 0.986 )
Likelihood ratio test= 309 on 17 df, p=0
Wald test = 151 on 17 df, p=0
Score (logrank) test = 338 on 17 df, p=0, Robust = 171 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