Haplotype frequencies were estimated using a Bayesian approach implemented with PHASE (Stephens and Donnelly, 2003). Haploview version 2.04 Software (Whitehead Institute for Biomedical Research, USA) was used to produce linkage disequilibrium (LD) matrices. Since rare and uncommon haplotypes are subject to estimation errors because of increased sampling variance, all analyses were conducted with haplotypes ≥5% frequency. These haplotypes were included in models for CRHR1 and for CRH together with the predictor variables listed above and backward stepwise regression was performed.
