ARTIGO ARTICLE
Infant mortality trends in the State of Rio Grande do Sul, Brazil, 1994-2004: a multilevel analysis of individual and community risk factors
Tendência da mortalidade infantil no Rio Grande do Sul, Brasil, 1994-2004: uma análise multinível de fatores de risco individuais e contextuais
Roselaine Ruviaro ZaniniI, II; Anaelena Bragança de MoraesI, II; Elsa Regina Justo GiuglianiII; João RiboldiII, III
ICentro de Ciências Naturais e Exatas, Universidade Federal de Santa Maria, Santa Maria, Brasil
IIPrograma de Pós-graduação em Epidemiologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brasil
IIIInstituto de Matemática, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brasil
ABSTRACT
The aim of this study was to analyze the trend in infant mortality rates in the State of Rio Grande do Sul, Brazil, from 1994 to 2004, in a longitudinal ecological study, by means of panel data analysis and multilevel linear regression (two levels: microregion and time) to estimate factors associated with infant mortality. The infant mortality rate decreased from 19.2 (1994) to 13.7 (2004) live births, and the principal causes of death in the last five years were perinatal conditions (54.1%). Approximately 47% of the variation in mortality occurred in the microregions, and a 10% increase in coverage by the Family Health Program was associated with a 1 reduction in infant mortality. A 10% increase in the poverty rate was associated with a 2.1 increase in infant deaths. Infant mortality was positively associated with the proportion of low birthweight newborns and the number of hospital beds per thousand inhabitants and negatively associated with the cesarean rate and number of hospitals per 100 thousand inhabitants. The findings suggest that individual and community variables display significant effects on the reduction of infant mortality rates.
Information Systems; Infant Mortality; Cause of Death
RESUMO
O objetivo deste trabalho foi analisar a tendência das taxas de mortalidade infantil no Rio Grande do Sul, Brasil, de 1994 a 2004, em estudo ecológico longitudinal, por meio de análise de dados de painel e regressão linear multinível (dois níveis: microrregião e tempo) para estimar fatores associados à mortalidade infantil. A taxa de mortalidade infantil reduziu de 19,2 (1994) para 13,7 (2004) nascidos vivos, e a principal causa de óbito, nos últimos cinco anos, foi afecções perinatais (54,1%). Aproximadamente 47% da variação nas taxas de mortalidade ocorreram nas microrregiões, e 10% de acréscimo na cobertura do Programa Saúde da Família esteve associado à redução de 1 na mortalidade infantil. O aumento de 10% na taxa de pobreza esteve associado com um aumento de 2,1 nos óbitos infantis. A mortalidade infantil associou-se positivamente com a proporção de recém-nascidos com baixo peso e número de leitos hospitalares por mil habitantes e negativamente com a proporção de cesarianas e número de hospitais por 100 mil habitantes. Os resultados sugerem que variáveis individuais e contextuais apresentam efeitos significativos na redução das taxas de mortalidade infantil.
Sistemas de Informação; Mortalidade Infantil; Causa de Morte
Introduction
The infant mortality rate (IMR), which expresses the probability of a live born infant dying before completing one year, is considered one of the most efficient indicators of a population's social, economic, and ethical development, and its follow-up allows inferences on the population's quality of life and its determinants 1,2,3.
Each year, some 11 million children die before reaching five years of age, the majority from avoidable causes and in poor countries. Of these, 7 million fail to complete their first year of life, with 70% of the deaths associated with nutritional causes or such diseases as respiratory infections, diarrhea, or malaria 4.
In Brazil, the infant mortality rate decreased by approximately 20% in five years, to 21.2 deaths per thousand live births in 2005 5, leaving the country in 84th place in the United Nations (UN) ranking, led by Singapore (2.3), as opposed to Angola (187.5) 6.
The State of Rio Grande do Sul has consistently shown the lowest infant mortality rates in Brazil, decreasing from 19.2 (1994) to 13.7 (2004) deaths per thousand live births 5. However, this rate is still higher than the acceptable value set by the World Health Organization (WHO), namely 10 deaths/thousand live births, a target already reached by such countries as Costa Rica, Chile, and Hungary 6.
Infant mortality is associated with numerous factors, summarized in a conceptual model proposed by Mosley & Chen 7 for developing countries, based on the assumption that social and economic determinants of infant mortality operate through proximate biological determinants. Thus, the decline in infant mortality is attributed, among other factors, to improvement in basic environmental sanitation, expanded access by the population to basic health services, control of diarrheic diseases, immunization, promotion of breastfeeding, greater coverage of prenatal care, and certain social programs 7.
Despite the decreases in infant mortality rates, there are persistent differences in the patterns of decline in mortality between geographic regions and population subgroups, i.e., differential decreases occur between States, as well as between microregions within a single State 3,8,9.
Infant mortality is divided into two components: neonatal (less than 28 days of life) and post-neonatal (28-364 days), since the two periods involve quite specific causes of death and thus require distinct preventive measures. Neonatal deaths can be considered avoidable through adequate prenatal and childbirth care 3,10,11. Likewise, a portion of post-neonatal deaths, which are simpler to control, can be avoided since they are directly associated with quality of life, like basic environmental sanitation and access to health services 12.
Various studies have investigated maternal and infant characteristics as factors associated with infant mortality 12,13,14,15. However, most of these studies overlook the inherent hierarchical structure in the data (multilevel structure), i.e., they fail to take into account the correlation between individuals in the same group, which can lead to incorrect inferences. Among the studies that consider the hierarchical structure, the majority only investigate grouping at the family level 16,17,18, overlooking the effects of the communities or regions in which the children and their mothers live.
In this sense, multilevel longitudinal studies are essential to expand knowledge on the changes that occur in a population 19. In addition, they allow separately estimating the variability at each level, making possible the analysis of the individual or regional characteristics that can explain these differences, thereby contributing to preventive strategies 20,21,22.
Thus, this study's objectives were to analyze the trend in infant mortality rates and their components in the State of Rio Grande do Sul, Brazil, and to identify the relevant individual and community factors by means of panel data and multilevel analytical models.
Method
This was a longitudinal ecological study of data on live births and infant deaths (0-364 days).
The dependent variable was IMR, expressed as the number of deaths in infants (< 1 year) per thousand live births, among children residing in the State's 35 microregions (groupings of neighboring municipalities or counties, with the aim of integrating the organization, planning, and execution of common public services) in the State of Rio Grande do Sul, from 1994 to 2004.
Data were obtained from the records in the Information System on Live Births (SINASC) and the Mortality Information System (SIM), available at the Information Technology Department of the Unified National Health System (DATASUS; http://www.datasus.gov.br).
Underlying cause of death was defined according to the chapters of the 10th revision of the International Classification of Diseases (ICD-10) 23 for the period from 1996 to 2004.
The study considered predictive variables on a continuous scale (centered on the mean) at two levels. The criteria used to include variables in the models were the theoretical framework and the statistical significance of the parameters' estimates.
The determinants that were considered at the individual level, available for the 11 years of the study, were: time, low birthweight rate, cesarean rate, women with seven or more prenatal visits, premature births, and measles and BCG vaccination.
The demographic characteristics of the microregions were: urbanization rate (proportion of the population residing in urban areas), life expectancy at birth (mean number of years that people live, starting from birth), and fertility (mean number of children per woman during her childbearing period). Socioeconomic indicators were: illiteracy (percentage of persons 15 years or older who could not read or write a simple note), mean schooling (years), poverty (percentage of the population with a per capita household income less than half the minimum wage, or approximately U$100/month), hospitalization due to assaults (percentage of hospitalizations due to assaults), population homicide rate, and labor force participation (economically active population as a percentage of the working-age population), in addition to gross domestic product (per capita GDP), socioeconomic development index, and household crowding (percentage of persons living in households with density greater than 2, as expressed by the number of residents per room, not counting bathrooms or one kitchen). Some coverage indicators were also considered, like: Family Health Program (FHP; percentage of persons enrolled), water (proportion of permanent private dwelling connected to the running water supply), sanitation (percentage of households connected to the general sewerage system), and supplementary health (percentage of beneficiaries of private health healthcare plans covering medical consultations, tests, and hospitalization).
Other health indicators were: number of practicing physicians (including medical residents) and hospital beds per thousand inhabitants, hospitals per 100 thousand inhabitants, and healthcare expenditure per inhabitant.
The variables were obtained from the databases of the Brazilian Institute of Geography and Statistics, or National Census Bureau (IBGE); DATASUS (National Immunization Program, under the Primary Health Care Data System - SIAB); National Supplementary Health Care System (ANS); Institute of Applied Economic Research (IPEA); Foundation for Economics and Statistics (FEE); National System of Urban Indicators (SNIU), and the Rio Grande do Sul State Socioeconomic Atlas (Rio Grande do Sul State Secretariat of Planning and Management), from different periods (2000 census, 1994-2004, 2003, or 2004), based on the availability of data.
Data from a same microregion were observed over the course of 11 years. To analyze trends in each time series, percentage rates in annual variation were used, adjusted by exponential regression of the mortality rates.
For analysis of the panel data via mixed models, the use of the MIXED procedure in SAS 9.1 (SAS Inst., Cary, USA) allowed considering different covariance structures, selecting the most adequate for within-individual variation in measures, adjusting a model whose coefficients were significant (p < 0.05). The best covariance structure was that which presented the lowest values for AIC (Akaike's Information Criterion) 24 and BIC (Schwarz's Bayesian Criterion) 25.
For selection of the community variables, in addition to the theoretical framework, an analysis of correlation was performed to evaluate multicollinearity.
In the multilevel linear regression, a hierarchical structure was considered, according to the microregions in Rio Grande do Sul (level 2) and the period from 1994 to 2004 (level 1), for which an unconditional (empty) model was adjusted, in addition to models with fixed and random coefficients. Models with complex variance were also considered, that is, models that allowed modeling the IMR variance as a function of the explanatory variables.
Estimation by the unconditional model produced information on the outcome's variability at each level, and the intra-microregion correlation coefficient was calculated by : 
. The regression coefficients were estimated by IGLS (iterative generalized least squares).
The statistical significance of the estimates was verified by the Wald test, and the best model was chosen by the goodness-of-fit (deviance) test, comparing the main results. As measures of the model's diagnosis, residues were analyzed for the two levels, as well as measures of leverage and influence.
The multilevel and panel data analyses were developed with the aid of MLwiN 2.02 computational applications (Centre for Multilevel Modelling, Bristol, UK) and SAS 9.1, respectively.
Multilevel models for continuous outcomes
Multilevel modeling considers a set of hierarchical data, with the response variable measured at the individual level and with explanatory variables that can be measured at the individual or community level. Thus, conceptually, the model can be viewed as a hierarchical system of regression equations, allowing estimation of the individual (intra-group) and community (inter-group) effects 19.
Level 1 models are developed separately in each unit of level j, considering the possibility of variation of intercepts and slopes 19,26.
For algebraic simplicity, the two-level model is considered here:
with i = 1, 2, ..., nj and j = 1, 2, ..., J; Yij: outcome of the ith unit of level 1, grouped in the jth of level 2; Xij: predictive variable measured in the ith unit of level 1, grouped in the jth of level 2; β0: model's general intercept; β1: slope coefficient associated with predictive variable X; u0j: random effect of level 2; e0ij: random effect of level 1.
Residues u0j and e0ij are assumed to be independent and normally distributed, with mean zero and variances and respectively.
The residual variance, i.e., the variance conditioned on X, is given by:
The covariance between two individuals (i1, i2) in the same group j is:
The intra-group correlation coefficient, which expresses the fraction of total variability that can be attributed to level 2, is given by:
This measure is an indicator of the degree of the population's grouping corresponding to the correlation between the values for two individuals in a group, controlling for variable X 19,26,27.
The model for equation (1) can be extended such that:
β0j: intercept for the jth unit of level 2; β1j: slope coefficient associated with variable X of the ith unit of level 1, grouped in the jth of level 2; β0: expected value of the intercepts in level 2; β1: expected value of the slopes in level 2; u0j: random effect of the jth unit of level 2 on intercept β0j; u1j: random effect of the jth unit of level 2 on slope β1j.
This model is known as a "random coefficients model", which assumes that each group has a different intercept (β0j) and slope (β1j), presenting a complex variance-covariance structure 19,26, given by:
The intercepts and slopes are estimated for units grouped within level 2, producing reduced residues ordered around an overall mean. In this case, the residues can be estimated, multiplying the crude residues 
by a reduction factor:
Ethical issues
The study was approved by the Institutional Review Board of the School of Medicine at the Federal University in Rio Grande do Sul and complies with the requirements in Ruling nº. 196/96 of the Brazilian National Health Council.
Results
In Rio Grande do Sul, from 1994 to 2004, there were 1,897,002 live births and 31,576 deaths in children under one year of age. Figure 1 shows the spatial distribution of mean infant mortality rates by microregion. The lowest means were in the microregions of Guaporé (11.2) and Montenegro (11.9), and the highest were in the microregions of Campanha Meridional (24) and Campanha Central (22.3).
Table 1 shows the trend in infant mortality rates and the component rates, in addition to the percentage rate of annual variation. There was a decrease of 2.5% in this rate for live births (p < 0.001), 2.2% for infant mortality (p = 0.002), 1.2% for neonatal mortality (p = 0.011), and 4% for post-neonatal mortality (p = 0.003). Early neonatal mortality (0-6 days) decreased by 1.9% (p < 0.001), while for late neonatal mortality (7-27 days) there was a non-significant increase of 1% (p = 0.186).
Table 2 shows the distribution of the principal causes of infant deaths. The leading causes were perinatal conditions (8.52), followed by the group of congenital malformations, deformations, and chromosomal abnormalities (3.33).
In the mixed models analysis (model 1), the most adequate covariance structure was autoregressive moving average (ARMA; 1, 1), with residual variance, covariance, and correlation of 12.26, 0.25, and 0.88, respectively, in addition to AIC = 2,045.2, BIC = 2,049.9, and deviance = 31.51 (p < 0.001). Other results are presented in Table 3.
In the multilevel analysis, considering a random intercepts model (model 2), the residual variance of the microregions was estimated at 10.49 (95%CI: 5.06-15.92; p < 0.001), while the residual variance for the years was 12.01 (95%CI: 10.23-13.79; p < 0.001). The estimated intra-microregion correlation coefficient was 0.4662, suggesting that 46.62% of the variation in the infant mortality rate was due to microregional variability in Rio Grande do Sul.
To investigate the possible determinants of this variability, specific characteristics of the microregions were included in the model. Table 3 shows some results from the multilevel models.
Model 3 (random intercepts), with individual-level variables only, showed deviance = 80.14 (p < 0.001), suggesting that the inclusion of these variables significantly improved the model's fit, helping explain the reduction in variance between years (16.78%) as well as between microregions (43.94%). In model 4 (random intercepts), the inclusion of community variables resulted in deviance = 34.22 (p < 0.001), also indicating a better fit, i.e., these variables also helped explain the reduction in variance between microregions (74.83%).
For model 5 (random intercepts and random slopes as a function of time), the goodness-of-fit test (deviance =13.44) indicated that the inclusion of the complex variance component at the level of microregions was significant (p < 0.001), suggesting that the variability in infant mortality rates between microregions is time-dependent, presenting a decrease until 2001 and increasing from that year onward. According to this model, the estimated infant mortality in the microregions was 18.95 (p < 0.001), with variance 5.68, and the estimated time-related slope coefficient (years) was -0.90, with variance 0.11. Estimated covariance was -0.67, producing a correlation of -0.85, indicating that the microregions with the lowest IMR in 1994 showed the largest decreases in infant mortality rates over the years. In relation to level 1, the variation between IMR values was constant over the 11 years.
Considering the effect of community variables, the 10% increase in FHP coverage was associated with a 1 drop in infant mortality, controlling for the other variables, and the 10% increase in the poverty rate was associated with a 2.1 increase in infant deaths. The number of hospital beds per thousand inhabitants was positively associated with the outcome, i.e., a 1% increase in this rate was associated with a 1.91 increase in infant deaths, while a 10% increase in hospitals per 100 thousand inhabitants was associated with a 3.6 drop in infant mortality.
The other variables were tested as random effects. However, these effects were not significant, and only remained in the model as fixed effects. Possible interactions between the variables were also considered, but none was statistically significant.
The analysis of the residues and predicted values indicated that the assumptions of normal distribution and linearity were met, and no outliers or influential values were identified. Some microregions also differed significantly from the overall mean, with Três Passos and Não-me-Toque showing an estimated intercept below the mean value for the microregions, while Campanha Meridional showed an estimated intercept above the mean intercept. Meanwhile, Campanha Meridional and Cachoeira do Sul showed estimated slopes below the mean, while Três Passos and Não-Me-Toque showed slopes above the mean. For all the other microregions, the estimates did not differ significantly from the mean.
Discussion
Infant mortality rates in Brazil decreased by approximately 70% in the last 60 years, principally in the post-neonatal component and in urban and more developed areas 1,2,9,14,28, despite persistent socioeconomic inequalities over time 9. Considered a multifactor process 15, the reduction in infant mortality in Brazil is associated with a series of improvements in living conditions and maternal and infant health care, related to food and nutritional security, basic environmental sanitation, and immunization 18,29.
In the State of Rio Grande do Sul, although the proportion of low birthweight newborns increased during this same period 5, infant mortality decreased, consistent with studies performed in Belo Horizonte (Minas Gerais) 3, Rio de Janeiro 14, and Londrina (Paraná) 30. As in other States of Brazil 3,14,30,31, post-neonatal mortality rates showed the largest decrease, and with the drop in this component a concentration of deaths occurred in the neonatal period, mostly in the first six days of life, thus highlighting the relationship between the deaths and quality of care for pregnant women.
A study 32 on the trend in infant mortality that analyzed the results of 80 surveys from 1980 to 2000 in 31 countries, including Brazil, concluded that the decline has been more pronounced in late neonatal and post-neonatal mortality. In addition, these decreases were not distributed equally between regions, with more rapid increases occurring in countries of Latin America and the Caribbean, North Africa and the Middle East, South and Southeast Asia, and the Pacific, while the decline was slower in Sub-Saharan Africa, especially in the early neonatal period 32.
In recent decades in Brazil, along with the decrease in the share of infectious and parasitic diseases, respiratory diseases, and malnutrition, perinatal conditions have become the main group of causes of infant death, i.e., problems that arise before completing one week of life, like prematurity, intrapartum asphyxia, and neonatal infections, followed by congenital malformations 33. In Rio Grande do Sul, perinatal conditions were the principal causes of infant mortality, followed by congenital malformations, jointly counting for some 75% of all deaths from 2000 to 2004. In the neonatal period, these two causes accounted for 96% of the deaths. Although more difficult to avoid, they can be reduced by measures that guarantee adequate prenatal care and improvements in childbirth care 33, as well as early diagnosis and treatment of maternal and fetal factors that can increase the risk of death.
Although Rio Grande do Sul had one of the lowest infant mortality rates in Brazil in 2004, 20.3% of infant deaths in the State could have been avoided by adequate prenatal and childbirth care in 2004, representing 470 lives that could have been saved. In addition, there are significant regional differences in infant mortality rates, which vary according to the level where they occur and complex social, economic, and demographic factors 34.
Several studies have reported that individuals immersed in health processes undergo influences at different levels. However, this hierarchy is not always taken into consideration, which can produce biased estimates when traditional regression techniques are used 19. Thus, this study used panel data and multilevel models to evaluate the effect of the levels and the determinants of infant mortality in Rio Grande do Sul, testing numerous variables. The results of the two analyses do not display relevant differences, but it was important to consider the different covariance structures in the panel analysis, since it improved the estimates of the respective parameters 19.
Multilevel analysis allows computing the magnitude of each level's contribution to explaining the outcome 19,20. In this study, practically half of the variability in infant mortality rates was due to the community level, i.e., the microregions' effect was considerable. In addition, inclusion of the complex variance term indicated inequality in the rates' variation between microregions over the years, increasing from 2001 onward. Thus, the variability in infant mortality rates between microregions changed as a function of time, but for each microregion the variations remained constant over time. A multilevel study on mortality in children under five years of age in the State of Rio Grande do Norte, based on data from the 2000 population census (http://www.ibge.gov.br), found a 62% effect for microregions and 38% for municipalities, and the main socioeconomic factors were: adequate basic environmental sanitation (access to bathrooms and running water), education (literacy rate), and dependency ratio greater than 75% 35. The larger effect of the microregions observed in the above-mentioned study can be explained by the greater heterogeneity between microregions in the Northeast of Brazil as compared to the South.
The inclusion of individual and community variables according to a multilevel approach allowed elucidating some aspects of infant mortality, and the significant determinants were: time, proportion of low birthweight, cesarean rate, FHP coverage, poverty rate, and number of hospital beds and hospitals per thousand inhabitants.
Infant mortality is closely related to birthweight. This study showed that in Rio Grande do Sul, a 1% decrease in the proportion of low birthweight newborns would decrease the infant mortality rate by 1.2 points. Thus, in order for Rio Grande do Sul to decrease its infant mortality rates, strategies are needed to prevent low birthweight, such as expanding access to qualified prenatal care, better maternal nutrition, and programs to reduce maternal smoking, among others 15,32.
This analysis suggests that an increase in FHP coverage, controlling for the other variables, reduced the infant mortality rate. This result was expected, since by increasing the health services' proximity to the local reality, the program seeks to plan and implement actions based on the problems within a given area, and such measures have contributed decisively to reducing infant mortality. A longitudinal, non-multilevel study (1990-2002) led by the Ministry of Health for all States of Brazil showed that the FHP had a significant impact on the drop in infant mortality: for each 10% increase in FHP coverage, infant mortality decreased by 4.5% 36. This reduction was probably smaller in Rio Grande do Sul since the rates were already lower, and there was less variability between microregions in the State, which is also more developed and provides greater access to basic health services.
The decrease in the poverty rate was associated with a reduction in infant deaths. A study of panel data on infant mortality in Brazil as a whole and in the Southeast region also found a significant effect by this indicator of monetary poverty, but average years of schooling showed a much stronger effect 34.
The number of hospital beds was positively associated with the outcome. Importantly, all hospital beds were considered, not only pediatric beds. We believe that if only intensive or intermediate-care neonatal beds were included, this association might be inverted, i.e., the more such neonatal beds, the lower the infant mortality. Meanwhile, the increase in the hospital rate was negatively associated with infant mortality. We found no multilevel regression study in the literature that considered these community variables, thus making it difficult to compare our findings with others.
Caution is necessary when interpreting the negative association between cesarean rate and infant mortality. The protective effect of cesareans against infant deaths detected in this study is probably due not to the procedure itself, but to the type of population in which cesareans are more common. Cesarean rates are known to be higher in private and outsourced hospitals, which treat a patient population with higher socioeconomic status, better nutrition, and access to qualified prenatal care, in addition to lower infant mortality.
The study's limitations include the possible presence of a measurement bias due to the use of secondary data, the size effect of microregions, which show great internal variability, and the absence of other determinants of infant mortality in the model, like alcohol consumption, smoking, and breastfeeding rates. However, due to the large number of variables that were tested and the correlation between some of them, the explanatory power of some variables that were not included would probably be small. Despite the limitations, multilevel modeling was able to capture significant effects, characterized at the respective levels.
Although Rio Grande do Sul has achieved the lowest infant mortality rate in Brazil, there are still microregions that require greater attention by government, like Cachoeira do Sul and Campanha Meridional, which showed a below-average reduction in infant mortality during the period analyzed.
Some comprehensive, simple, and relatively low-cost interventions should be prioritized, since they are capable of acting in the health promotion of socially vulnerable population groups and in the primary and secondary prevention of more prevalent diseases. These include sex education programs for adolescents, qualified prenatal care, programs for access to (and humanization of) childbirth care, programs to promote exclusive breastfeeding in the first six months and complemented by two years or more, programs for the identification of infant risk and early treatment of illnesses (Program for Integrated Management of Childhood Illness, Ministry of Health), National Immunization Program, and expanded coverage of the FHP 8,18,29,33.
The multilevel analysis of the trend in infant mortality rates provided a pioneering approach in the State of Rio Grande do Sul. In addition to demonstrating the importance of this methodology for identifying each level's contribution to the data's hierarchy, this study showed important results related to the factors associated with infant mortality, suggesting the need for efficient public measures targeting vulnerable groups.
Contributors
R. R. Zanini and A. B. Moraes contributed with the design, data collection, statistical analysis, and interpretation, and writing of the article J. Riboldi and E. R. J. Giugliani collaborated in the supervision, writing, and revision of the article.
References
1. Victora CG, Grassi PR, Schmidt AA. Situação da saúde da criança em área da região sul do Brasil, 1980-1992: tendências temporais e distribuição espacial. Rev Saúde Pública 1994; 28:423-32.
2. Ribeiro VS, Silva AAM. Tendências da mortalidade neonatal em São Luís, Maranhão, Brasil, 1979 a 1996. Cad Saúde Pública 2000; 16:429-38.
3. Caldeira AP, Franca E, Perpetuo IH, Goulart EM. Trends in avoidable causes of infant mortality in Belo Horizonte, Brazil, 1984 to 1998. Rev Saúde Pública 2005; 39:67-74.
4. United Nations Children's Fund. Child mortality statistics. http://www.childinfo.org/cmr/revis/db2.htm (accessed on 22/Nov/2006).
5. Departamento de Informática do SUS. Informações de saúde - estatísticas vitais: mortalidade e nascidos vivos. http://www.tabnet.datasus.gov.br (accessed on 10/Jan/2008).
6. Central Intelligence Agency. The world factbook: guide to rank order pages. Infant mortality rate. https://www.cia.gov/cia/publications/fact book/rankorder/2091rank.html (accessed on 20/Sep/2006).
7. Mosley WH, Chen L. An analytical framework for the study of child survival in developing countries. Popul Dev Rev 1984; 10 Suppl:25-45.
8. Paim JS, Costa MCN. Decréscimo e desigualdade da mortalidade infantil, Salvador, 1980-1988. Bol Oficina Sanit Panam 1993; 114:415-28.
9. Simões CCS, Monteiro CA. Tendência secular e diferenciais regionais da mortalidade infantil no Brasil. In: Monteiro CA, organizador. Velhos e novos males da saúde no Brasil: a evolução do país e de suas doenças. São Paulo: Editora Hucitec; 1995. p. 153-6.
10. Lansky S, França E, Leal MC. Mortalidade perinatal e evitabilidade: revisão da literatura. Rev Saúde Pública 2002; 36:759-72.
11. Ortiz LP. A mortalidade neonatal no Estado de São Paulo: níveis e tendências segundo a duração da gestação. In: XV Encontro Nacional de Estudos Populacionais. http://www.abep.nepo.unicamp.br/encontro2006/docspdf/ABEP2006_621.pdf (accessed on 12/Nov/2006).
12. Lopes AD. Reducing child mortality. Bull World Health Organ 2000; 78:1173.
13. Hill K. Frameworks for studying the determinants of child survival. Bull World Health Organ 2003; 81:138-9.
14. Leal MC, Szwarcwald CL. Evolução da mortalidade neonatal no Estado do Rio de Janeiro, Brasil (1979-1993): Análise por grupo etário segundo região de residência. Rev Saúde Pública 1996; 30:403-12.
15. Rutstein SO. Factors associated with trends in infant and child mortality in developing countries during the 1990s. Bull World Health Organ 2000; 78:1256-70.
16. Das Gupta M. Death clustering, mothers' education and determinants of child mortality in rural Punjab, India. Popul Stud (Camb) 1990; 44: 489-505.
17. Guo G, Rodrigues G. Estimating a multivariate proportional hazards model for clustered data using the EM algorithm with an application to child survival in Guatemala. J Am Stat Assoc 1992; 87:969-76.
18. Sastry N. A nested frailty model for survival data, with an application to the study of child survival in northeast Brazil. J Am Stat Assoc 1997; 92:426-35.
19. Goldstein H. Multilevel statistical models. 3rd Ed. London: Edward Arnold; 2003.
20. Leyland AH, Groenewegen PP. Multilevel modeling and public health policy. Scand J Public Health 2003; 31:267-74.
21. Wright D, Bobashev GV, Novak SP. Decomposing the total variation in a nested random effects model of neighborhood, household, and individual components when the dependent variable is dichotomous: implications for adolescent marijuana use. Drug Alcohol Depend 2005; 78:195-204.
22. Matteson DW, Burr JA, Marshall JR. Infant mortality: a multi-level analysis of individual and community risk factors. Soc Sci Med 1998; 47:1841-54.
23. Organização Mundial da Saúde. Classificação estatística internacional de doenças e problemas relacionados à saúde, 10a revisão. São Paulo: Centro Colaborador da OMS para a Classificação de Doenças em Português; 1995.
24. Akaike H. A new look at the statistical model identification. IEEE Trans Automat Contr 1974; 19: 716-23.
25. Schwarz G. Estimating the dimension of a model. Ann Stat 1978; 6:461-4.
26. Snijders T, Bosker R. Multilevel analysis: an introduction to basic and advanced multilevel modeling. London: Sage Publications; 1999.
27. Barros AJD. Modelos multinível: primeiros passos. http://www.pjeventos.com.br/epi2002/Index2.htm (accessed on 07/May/2002).
28. Bercini LO. Mortalidade neonatal de residentes em localidade urbana da região sul do Brasil. Rev Saúde Pública 1994; 28:38-45.
29. Fundo das Nações Unidas para Infância. Infância ainda vulnerável. http://www.unicef.org/brazil/sib2001/cap2.htm (accessed on 10/Nov/2006).
30. Ferrari LSL, Brito ASJ, Carvalho ABR, Gonzáles MRC. Mortalidade neonatal no Município de Londrina, Paraná, Brasil, nos anos de 1994, 1999 e 2002. Cad Saúde Pública 2006; 22:1063-71.
31. Holcman MM, Latorre MRDO, Santos JLF. Evolução da mortalidade infantil na região metropolitana de São Paulo, 1980-2000. Rev Saúde Pública 2004; 38:180-6.
32. Hill K, Choi Y. Neonatal mortality in the developing world. Demogr Res 2006; 18:429-52.
33. Fundo das Nações Unidas para Infância. Situação da infância brasileira 2006. http://www.unicef.org/brazil/sib06h.htm (accessed on 05/Dec/2006).
34. Gomes FAR, Araújo Júnior AF, Salvato MA. Mortalidade infantil no Brasil e Sudeste: determinantes e perspectivas para o futuro. In: XV Encontro Nacional de Estudos Populacionais. http://www.abep.nepo.unicamp.br/encontro2006/docspdf/ABEP2006_653.pdf (accessed on 18/Nov/2006).
35. Ramos PCF, Formiga MCC, Costa CCL, Ponce de Leon AC. Diferenciais regionais da mortalidade da população infantil no estado do Rio Grande do Norte evidenciados pela utilização de modelos multiníveis. In: XV Encontro Nacional de Estudos Populacionais. http://www.abep.nepo.unicamp.br/encontro2006/docspdf/ABEP2006_660.pdf (accessed on 18/Nov/2006).
36. Macinko J, Guanais FC, Souza MFM. An evaluation of the impact of the Family Health Program on infant mortality in Brazil, 1990-2002. J Epidemiol Community Health 2006; 60:13-9.
Correspondence:
R. R. Zanini
Departamento de Estatística
Centro de Ciências Naturais e Exatas
Universidade Federal de Santa Maria
Av. Roraima 1000, prédio 13
Santa Maria, RS, 97105-900, Brasil
rrzanini@smail.ufsm.br
Submitted on 15/Oct/2007
Final version resubmitted on 12/Mar/2008
Approved on 13/Mar/2008