Estimates were made using the survey commands functions of the STATA program (STATA Corp, College Station, Texas, USA), which enabled us to incorporate the sampling design and weights into all statistical calculations (descriptive, χ2, logistic regression). Statistical significance was set at a 2-tailed α < 0.05.