Sugar-sweetened beverages consumption and BMI in Mexican adolescents. Mexican National Health and Nutrition Survey 2006

Consumo de bebidas azucaradas y su relación con el IMC en adolescentes mexicanos. Encuesta Nacional de Salud y Nutrición 2006

Alejandra Jiménez-Aguilar, MSc; Mario Flores, MD, MSc; Teresa Shamah-Levy, MSc

Centro de Investigación en Nutrición y Salud, Instituto Nacional de Salud Pública. Cuernavaca, Morelos, México

ABSTRACT

OBJECTIVE: To evaluate the association between the consumption of sugar-sweetened beverages (SSBs) and body mass index (BMI) in Mexican adolescents.
MATERIAL AND METHODS: We analyzed the data of 10 689 adolescents (ages 10 to 19 years old) who participated in the Mexican National Health and Nutrition Survey 2006 (ENSANUT 2006). Consumption of SSBs (i.e. sodas, fruit beverages and sugar beverages) was evaluated by means of a semi-quantitative food frequency questionnaire. BMI was calculated (kg/m²).
RESULTS: Mean age was 13.8 ± 2.7 years. Fifty percent were females. Mean BMI was 21.7 ± 4.5. Thirty percent of adolescents were overweight or obese. Ninety percent of adolescents consumed at least one SSB during the 7 days before the interview. The median consumption of SSBs was 0.89 portion per day. Multiple-linear regression analysis showed that for each portion of sodas consumed, a 0.17-point increase in BMI was observed in boys after adjusting for confounders (95% CI; 0.02-0.32, p 0.03). Positive interactions of SSB consumption with age and time watching TV were observed in boys.
CONCLUSIONS: Consumption of sodas was positively associated with BMI in Mexican boys.

Key words: beverages; body mass index; overweight; obesity; adolescents; Mexico

RESUMEN

OBJETIVO: Examinar la asociación entre el consumo de bebidas refrescantes azucaradas (BRA) y el índice de masa corporal (IMC) en adolescentes mexicanos.
MATERIAL Y MÉTODOS: Se analizaron datos de 10 689 adolescentes (10 a 19 años de edad) de la Encuesta Nacional de Salud y Nutrición 2006 (ENSANUT 2006). El consumo de bebidas refrescantes azucaradas (BRA: refrescos, bebidas de fruta y bebidas endulzadas) se evaluó con un cuestionario semicuantitativo de frecuencia de consumo de alimentos. Se calculó el índice de masa corporal (kg/m²). ]]> RESULTADOS: La media de edad fue de 13.8 ± 2.7 años. El 50.4% fueron mujeres. La media de IMC fue de 21.7 ± 4.5. Un 30% de los adolescentes presentó sobrepeso u obesidad. El 90% de los adolescentes consumieron al menos una BRA en los 7 días previos a la encuesta, con una mediana de 0.89 porciones/día. El análisis de regresión lineal mostró que por cada porción consumida de refrescos, el IMC en los adolescentes varones aumentó 0.17 unidades, después de ajustar por variables confusoras (IC 95%: 0.02, 0.32, p= 0.03). Se observaron interacciones entre el consumo de BRA con la edad y el tiempo viendo televisión en los varones adolescentes.
CONCLUSIONES: El consumo de refrescos se asoció positivamente con el IMC en varones adolescentes mexicanos.

Palabras clave: bebidas; índice de masa corporal; sobrepeso; obesidad; adolescentes; México

The prevalence of overweight and obesity have dramatically increased in Mexico.¹ The most recent Mexican National Health and Nutrition Survey (ENSANUT 2006) found that one out of three adolescents is overweight or obese.²At the same time an increase in consumption of sugar-sweetened beverages (SSBs) –such as sodas, fruit beverages and sweet-flavored drinks– has been observed. Mexico is the second largest consumer of soft drinks, after the US.³ According to data from National Income-Expenditure Surveys (NIES) household expenditure on soft drinks has risen by 37.2% between 1986 and 1998.¹ At present, beverages contribute a fifth of all calories consumed by Mexicans. Because of this increase, recent Mexican guidelines suggest that energy from beverages must not exceed 10% of the recommendations of total energy.⁴

Recent evidence has shown that a high consumption of SSBs was associated with body weight gain and obesity in different populations.⁵ In Mexico, there is scarce information about this problem.

The objective of the present study is to examine the association between consumption of SSBs and body mass index (BMI) in Mexican adolescents using data from the ENSANUT 2006.

Material and Methods

Briefly, the ENSANUT 2006 was a probabilistic survey with a complex, stratified sampling design, carried out between October 2005 and May 2006. The ENSANUT 2006 is representative of four regions (north, center, south and Mexico City)^* and of rural (< 2 500 population) and urban (> 2 500 population) areas in the Mexican Republic. It obtained information from 48 304 households.^‡

For determining sample size it was considered that the minimum interest prevalence should have an accuracy of 8.1%. It was also considered that state estimators obtained through the survey should have a 25% maximum relative error, a 95% confidence interval, a 20% non response rate, and a design effect of 1.7, determining a sample size of at least 1 476 households by state. A subsample corresponding to a third of the whole sample was selected for dietary information.²

The study sample was limited to adolescents (boys and girls 10 to 19 years old) who had dietary data. Prior to the interview, informed consent was obtained from the family head and from the adolescent. Protocol was reviewed and approved by the Ethics, Bio-safety and Research Committees at the National Institute of Public Health (INSP).

Dietary information. Trained personnel conducted an adapted version of the food-frequency semi-quantitative questionnaire (FFSQ), for its application in digital format; found in the Procedure Handbook for Nutrition Projects, published by the INSP.⁶ The questionnaire included 101 food and beverage items classified into 14 groups. Data were entered into HP Compaq nx 6120 and Dell Latitude D510 laptop computers. Personnel from the INSP developed and validated the electronic questionnaires for input of data into the Fox Pro program version 7.

Non beverages energy. Dietitians converted the reported consumptions into grams or ml of food and beverages items. Energy from sources other than SSBs was then calculated using a comprehensive database compiled from diverse sources.^** For data cleanup, observations with an energy intake less than 500 kcal or more than 7000 kcal were eliminated. ⁷

Consumption of sugar-sweetened beverages. The standard portion of beverages was a 240 ml glass. The FFSQ collected information about the consumption of sodas, fruit drinks, sugar beverages and diet beverages. Participants were asked for the numbers of days (never, one, two to four, five to six, and seven days), times a day (one, two to three, four to five, and six times per day), and total of glasses consumed for each beverage during a seven-day period before the interview. We obtained the total consumption of each beverage from the number of servings per day. The consumption of diet beverages was categorized as low or high intake with respect to the median of its distribution.

Body mass index. Standardized interviewers^8,9 took measurements of weight with a portable electronic scale (Tanita) with a precision of 100 g. Height was obtained with a Dynatop stadiometer, with a capacity of 2 meters and an accuracy of 1 mm.

Body mass index was calculated as the weight in kilograms divided by the height in meters squared. BMI data were considered valid when they were between 10 and 58 units.² The classification of overweight and obesity was done in accordance with age and sex-specific BMI cut points recommended by the International Obesity Task Force (IOTF).^2,10 The 15^th percentile of NCHS' weight-for-age z-scores was used to identify individuals with low weight.¹¹

Socieconomic status. An indicator of socioeconomic status (SES) was constructed using a principal components analysis.^12,13 It included variables related to housing conditions (such as flooring and roofing materials), home appliances (including refrigerators, stoves, washing machines, TV sets, radios, videoplayers, telephones, and computers) and number of rooms (other than bathrooms, kitchens and corridors). The first component accounted for 46% of the total variance. The standardized factor was divided into tertiles.

Sex and age. Sex and age was obtained by the survey team.

Presence of menarche. Presence of menarche was obtained by interviews and used as an indicator of sexual maturation in girls.¹⁶

Data Analysis. Descriptive statistics (medians, 25th and 75th percentiles) of servings per day of SSBs were obtained. Median intakes of SSBs across biological and socioeconomic categories were compared. Differences between categories were obtained by median regressions.

To assess the relationship between consumption of SSBs and BMI, linear regression models were developed. After examination of multi-collinearity, beverages as servings consumed per day were included simultaneously as independent variables (sodas, fruit drinks and sugar beverages). The variable SSBs as a group was tested in another model.

We sequentially adjusted for sets of variables that might confound the association between SSBs and BMI: non-beverage energy, diet drinks, region (urban or rural), age, SES, and time watching TV. Models were stratified by sex, and for female adolescents presence of menarche was included.

Two-way interaction terms were tested in the linear regression models for consumption of SSBs and age, and for TV-viewing and menarche.

Sampling weights (expansion factors) were calculated as the inverse of the selection probability for each individual, accordingly to the sampling scheme used.² All calculations were weighted by expansion factors and adjusted for sampling (clustering) effects using the SVY STATA 9.0 module (Stata Corporation, TX, USA)^*** for complex surveys. A statistically significant level of 0.05 was used.

Results

The comparison between subjects included and excluded from analysis showed no significant differences with respect to age, sex or BMI (p > 0.05).

Fifty percent of the subjects were female adolescents. The mean age was 13.8 ± 2.7 years. More than 30% of the male adolescents were in the group of 13 to 15 years, and nearly 30% of the females were in the oldest age group (16 to 19 years). Mean BMI was 21.4 ± 4.4. On average, female adolescents had a higher BMI and a higher prevalence of overweight in relation to males (p< 0.05). The greater proportion of adolescents had a low SES (49.9%). Adolescents were concentrated in the center (39.8%) and southern (40.6%) regions of the country. The greater proportion of adolescents lived in urban areas (57.5%). The average time watching TV was 7.5 ± 5.5 h/w. Fifty-two percent of adolescents watched TV less than seven hours a week. A higher proportion of female adolescents watched TV for more than 21 hours a week as compared to males (p< 0.05).The average age at menarche was 12.1 ± 1.2 years and 70.8% of the females already had it (Table I).

During the seven days prior to the survey, 77.2% of the adolescents consumed sodas, 39.7% consumed fruit beverages and 27.7% consumed sugar beverages. Ninety percent consumed one or more SSBs. On average, a somewhat higher percentage of male adolescents consumed sodas (79%) compared to females (75.5%) (p< 0.01). Moreover, approximately 91% of males and 89% of females (p< 0.05) consumed one or more SSBs (data not shown).

At the national level, the daily median consumption of SSBs was 0.89 (percentile 25: 0.43, p75: 1.96). Median SSB consumption was higher for male adolescents than for females (p< 0.001). Median SSB consumption was higher among older adolescents (13 to 19 years) compared to the youngest adolescents (10 to 12 years). In general, overweight adolescents consumed more SSBs compared to adolescents who had a normal weight. The consumption of SSBs was positively associated with SES. A higher consumption of SSBs was observed in the northern region and in urban areas (p< 0.05). Adolescents who already had menarche consumed, on average, more SSBs than those who had not yet had it (p< 0.05). Adolescents who watched less than 7 hrs of TV per week consumed less SSBs (p< 0.05) (Table II).

For male adolescents, the unadjusted analysis showed that consumption of sodas was positively associated with BMI. For each serving of soda consumed, BMI increased 0.54 units (95% CI: 0.39-0.70, p< 0.001), whereas consumption of fruit drinks and sugar beverages were not significantly associated with BMI. It was noted that before adjusting for non-beverage energy the coefficient related to the consumption of soft drinks decreased (β= 0.16, 95% CI: 0.01-0.31, p= 0.04), while after adjustment for non-beverage energy the coefficient increased slightly (β= 0.17, 95% CI: 0.02-0.32, p= 0.03). None showed significant associations between consumption of fruit drinks and sugar beverages with BMI before or after adjusting for non-beverage energy. After adjusting for non-beverage energy, the consumption of all SSBs were not significantly associated with BMI in male adolescents (β= 0.06, 95% CI: -0.04-0.17, p= 0.25) (Table III).

In female adolescents, the unadjusted model showed that for every serving of soda consumed, BMI increased 0.19 units (95% CI: 0.03-0.35, p= 0.03). For fruit drinks and sugar beverages, no significant associations were observed. The adjusted analysis showed that before and after adjusting for non beverages energy, no drink studied, including the set of SSBs was significantly associated with BMI in females (Table III).

Figure 1 presents the adjusted mean BMI for male adolescents derived from the regression models that included interaction terms between the consumption of SSBs and age, and between the consumption of SSBs and TV-viewing time. In panel A, it is observed that in each age group, males with a higher consumption of SSBs had a higher BMI compared to those with a lower consumption of SSBs (p< 0.05).

In panel B, it is observed that in the category of 7 hrs of TV per week, there were no significant differences between male adolescents with a lower consumption of SSBs compared to males with a higher consumption of SSBs (p > 0.05), while in the categories of more than 7 hrs of TV per week it was observed that males with a greater consumption of SSBs had a higher BMI compared to those with a lower consumption of SSBs (p< 0.05) (Figure 1).

Figure 2 presents the adjusted mean BMI for girls resulting from the regression models that included interaction terms between the consumption of SSBs and age, and between the consumption of SSBs and TV-viewing time. In panel A, it is observed for the group of female adolescents 10 to 12 years old that the BMI was higher in those who had a higher consumption of SSBs compared to those who had a lower consumption of SSBs (p< 0.05). The oldest females (16 to 19 years) who had a higher intake had a BMI lower than those who had a lower consumption of SSBs (p< 0.05). In panel B, it is observed in the categories of more than 7 hrs of TV per week that the females who had a greater consumption of SSBs had a lower BMI with respect to those who had a lower consumption of SSBs, but the differences were not significant (p > 0.05) (Figure 2).

Discussion

We analyzed the relationship between consumption of sugar-sweetened beverages (SSBs) and BMI in a sample of adolescents who participated in a nationally representative Mexican survey (ENSANUT 2006).

The results showed that consumption of sodas were positively associated with BMI for male adolescents, after adjusting for potential confounders. Moreover, positive interactions of SSB consumption with age and TV-viewing time were observed in males. Giammattei et al. showed a similar interaction between SSB consumption with TV-viewing time, finding a significant positive association in adolescents between BMI with hours of television watched per evening and daily soft drink consumption (p< 0.05).¹⁷

Among female adolescents, no significant associations were found. This finding is consistent with previous studies, which have shown more consistent associations in males compared to females. However, the results are quite diverse.^18,19

It is important to mention that the results in female adolescents with respect to age group show a different direction than expected; we observed that in the oldest group of females (16 to 19 years), those who had a higher intake had a lower BMI compared with those females who had a lower consumption of SSBs (p< 0.05). This situation could be due to reverse causation, since it is possible that the oldest group of females has a greater concern about body weight than the youngest group. Thus, female adolescents 16 to 19 years with a greater BMI tended to consume less SSBs than those 16 to 19 years who had a lower BMI.

In general terms, the prevalence of overweight and obesity, consumption of different types of beverages, as well as the effect of the consumption of SSBs on BMI found in our study coincide with other studies. For example, Berkey et al. found that for male adolescents who increased their soda intake by > 2 servings per day from the previous year gained weight (0.14 increase in BMI).¹⁹ Ludwig et al. found at baseline consumption of SSBs that in the fully adjusted model, BMI increased by 0.18 for each serving consumed per day (95% CI: 0.09-0.27; p= 0.02).⁸ Phillips et al. found that subjects in the third and fourth quartiles of percentage of calories from soda had BMI _Z scores of 0.17 units higher, on average, than those of subjects in the fist quartile (p< 0.001).²⁰ In another study conducted among Mexican adolescents (10 to 19 years), the association between the consumption of SSBs and BMI was greater than that observed in this study. Denova et al. found that for each additional SSBs serving consumed daily, adolescents BMI increased by 0.33 (p< 0.001).²¹

The review by Malik⁵ has confirmed the notion about a positive association between consumption of SSBs and weight gain in different populations. The biological plausibility of how the consumption of SSBs can influence body weight is still unknown. It has been suggested that the consumption of energy through liquid food tends to generate a lower sense of satiety as compared to the intake of solid foods, resulting in a higher consumption of food in subsequent meals and, therefore, a larger total energy consumption.²² However, in our study we controlled for non-beverage energy intake and this did not affect the results. In one study among obese adolescents that compared satiety after consuming the same amount of energy through liquid foods versus solid foods, no difference in satiety was observed.²³ Among other assumptions, it has been proposed that the sugars contained in the SSBs could prove addictive. A study in rats showed that after withdrawal of a sugar-rich diet, the animals showed qualitatively clear signs of dependence on sugars, similar to those presented in morphine and nicotine addiction, including teeth chattering and anxiety.²⁴ On the other hand, it is possible that consumption of SSBs is only one component of a feeding pattern that favors body weight gain, or perhaps a reflection of an "obesogenic" lifestyle. A relationship was tested in an intervention study aimed to reduce the consumption of SSBs by American adolescents, and a significant reduction of body weight was observed among overweight and obese adolescents regardless of lifestyle variables.²⁵

The interpretation of our results should consider certain limitations. The most important limitation deals with the cross-sectional nature of our study, which does not allow for making causal inferences between consumption of SSBs and increased BMI. However, our results coincide with previous longitudinal studies that suggest a cause-effect association (that is, consumption of SSBs precedes and has a positive association with weight gain). Thus, one could rule out the idea of a reverse causality for male adolescents (being overweight or obese promotes the intake of SSBs).

Some of the strengths of the present study consist of the use of a large, nationally representative sample of Mexican adolescents, as well as the use of standardized personnel and methods. Additionally, the subjects and interviewers were blinded to the hypothesis of the study. Finally, we controlled in the multivariate analyses for potential confounders.

In conclusion, there was a positive association between the consumption of sodas and BMI in a nationally representative sample of Mexican adolescents. However this association was observed in males only. While it is necessary to conduct further studies on the subject with longitudinal design, it is also important to take immediate actions in this regard, such as encouraging the decrease in consumption of SSBs, as well as promoting the consumption of plain water for preventing and controlling the epidemic of overweight and obesity in Mexico. The Secretary of Health has recently undertaken an effort in this regard when designing the "recommendations for beverages' consumption for a healthy lifestyle for the mexican population".

Acknowledgments

We would like to acknowledge Dr. Bernardo Hernández for his advice in writing the manuscript and Alfonso Mendoza for database management.

References

1. Rivera J, Barquera S, Campirano F, Campos I, Safdie M, Tovar V. Epidemiological and nutritional transition in Mexico: rapid increase of non-communicable chronic diseases and obesity. Public Health Nutrition 2002;5:113-122.         [ Links ]

2. Fernández O, Rivera J, Shamah T, Villalpando S, Hernández M, Sepúlveda J. Encuesta Nacional de Salud y Nutrición 2006. Cuernavaca, Mexico: Instituto Nacional de Salud Pública, 2006.         [ Links ]

3. Gómez A. Cada año, 152 litros de refresco por mexicano. La prensa. Organización Editorial Mexicana. 2008 June. [Consulted 2008 November 21st.]. Available at: http://www.oem.com.mx/laprensa/notas/n752335.htm.         [ Links ]

4. Rivera JA, Muñoz O, Rosas M, Aguilar CA, Popkin BM, Willett WC. Consumo de bebidas para una vida saludable: recomendaciones para la población mexicana. Salud Publica Mex 2008; 50:173-195.         [ Links ]

5. Malik VS, Schulze MB, Hu FB. Intake of sugar-sweetened beverages and weight gain: a systematic review. Am J Clin Nutr 2006;84:274-288.         [ Links ]

6. Ramírez I, Mundo V, Rodríguez S, Vizuet I, Hernández N, Jiménez A. Encuestas dietéticas. In: Manual de procedimientos para proyectos de nutrición. Cuernavaca, Mexico: Instituto Nacional de Salud Pública, 2006; 27-108.         [ Links ]

7. Ludwig DS, Peterson KE, Gortmaaker S. Relation between consumption of sugar-sweetened drinks and childhood obesity: a prospective, observational analysis. Lancet 2001;357:505-508.         [ Links ]

8. Lohman TG, Roche AF, Martorell R. Anthropometric standardization reference manual. Champaign (IL): Human Kinetics Books, 1998.         [ Links ]

9. Habicht JP. Standarization of anthopometric methods in the field. PAHO Bull 1974;76:375-384.         [ Links ]

10. Cole T J, Bellizzi MC, Flegal KM, Dietz WH. Establishing a Standard definition for child overweight and obesity worldwide: international survey. BMJ 2000;320:1240.         [ Links ]

11. Must A, Dallal G, Dietz W. Reference data for obesity: 85th and 95th percentiles of body mass index (wt/ht2) and triceps skin fold thickness. Am J Clin Nutr 1991;53:839-846.         [ Links ]

12. Reyment R, Jöreskog K. Applied factor analysis in the natural sciences. Cambridge, UK: Cambridge University Press, 1996: 371.         [ Links ]

13. Hair JF, Anderson R, Tatham R, Black W. Multivariate data analysis with readings. 3rd ed. New York (NY): Macmillan Publishing Company, 1992: 223-253.         [ Links ]

14. Hernández B, Gortmaker S, Laird N, Colditz G, Parra-Cabrera S, Peterson K. Validez y reproducibilidad de un cuestionario de actividad e inactividad física para escolares de la ciudad de México. Salud Publica Mex 2000;42:315-323.         [ Links ]

15. American Academy of Pediatrics. Committee on Public Education Pediatrics. Pediatrics 2001;108:1222-1225.         [ Links ]

16. World Health Organization. El estado físico: uso e interpretación de la antropometría. Informe de un Comité de Expertos de la OMS. Geneva: WHO, 1995.         [ Links ]

17 Giammattei J, Blix G, Marshak HH, Wollitzer AO, Pettitt DJ. Television watching and soft drink consumption: associations with obesity in 11- to 13- year-old schoolchildren. Arch Pediatr Adolesc Med 2003;157:882-886.         [ Links ]

18. Gillis LJ, Bar-Or O. Food away from home, sugar-sweetened drink consumption and juvenile obesity. J Am Coll Nutr 2003;22:539-545.         [ Links ]

19. Berkey CS, Rockett HR, Field AE, Gillman MW, Coditz GA. Sugar added beverage and adolescent weight change. Obes Res 2004;12:778-788.         [ Links ]

20. Phillips SM, Bandini LG, Naumova EN, Cyr H, Colclough S, Dietz WH, et al. Energy-dense snack food intake in adolescents: longitudinal relationship to weight and fatness. Obes Res 2004; 12:461-472.         [ Links ]

21. Denova E, Jiménez A, Halley E, Huitrón G, Talavera JO, Pineda D, et al. Association between Sweetened Beverage Consumption and Body Mass Index, Proportion of Body Fat and Body Fat Distribution in Mexican Adolescents. Ann Nutr Metab 2008;53:245-251.         [ Links ]

22. Mattes RD. Dietary compensation by humans for supplemental energy provided as ethanol or carbohydrate in fluids. Physiol Behav 1996;59:179-187.         [ Links ]

23. Almiron E, Flores SY, Drewnowski A. No difference in satiety or in subsequent energy intakes between a beverage and a solid food. Physiol Behav 2004;82:671-677.         [ Links ]

24. Colantuoni C, Rada P, McCarthy J, Petten C, Avena NM, Chadeayne A, et al. Evidence that intermittent, excessive sugar intake cause endogenous opiate dependence. Obes Res 2002; 10:478-488.         [ Links ]

25. Ebbeling CB, Feldman HA, Osganian SK, Chomitz VR, Ellenbogen SJ, Ludwig DS. Effects of decreasing sugar-sweetened beverage consumption on body weight in adolescents: a randomized, controlled pilot study. Pediatrics 2006;117:673-680.         [ Links ]

Received on: April 11, 2008 ]]> Accepted on: February 25, 2009

Address reprint requests to: Dr. Mario Flores. Instituto Nacional de Salud Pública. Av. Universidad 655, col. Santa María Ahuacatitlán. 62100 Cuernavaca, Morelos, México. E-mail: mflores@correo.insp.mx
* The regions were established according to the definition by INEGI (National Institute of Statistics, Geography and Informatics): north (Baja California, Baja California Sur, Coahuila, Chihuahua, Durango, Nuevo Leon, Sonora and Tamaulipas), center (Aguascalientes, Colima, Guanajuato, Jalisco, Mexico (excluding urbanized municipalities and localities adjacent to Mexico City), Michoacan, Morelos, Nayarit, Queretaro, San Luis Potosi, Sinaloa, and Zacatecas), Mexico City (including the Federal District and urbanized municipalities in the State of Mexico), and south (Campeche, Chiapas, Guerrero, Hidalgo, Oaxaca, Puebla, Quintana Roo, Tabasco, Tlaxcala, Veracruz and Yucatan).
‡ A household unit was defined as all persons who usually sleep in the same house and benefit from a common income, whether or not they are connected by kinship.
** Safdie M, Barquera S, Porcayo M, Rodríguez S, Ramírez C, Rivera J, et al. Bases de datos del valor nutritivo de los alimentos. Compilación del Instituto Nacional de Salud Pública, 2004. (Unpublished document).
*** Stata Corp. Intercooled Stata 9.2 College Station Texas, USA, 2006.