Mortality due to garbage codes in Brazilian municipalities: differences in rate estimates by the direct and Bayesian methods from 2015 to 2017

Renato Azeredo Teixeira Lenice Harumi Ishitani Elisabeth França Pedro Cisalpino Pinheiro Marina Martins Lobato Deborah Carvalho Malta About the authors

ABSTRACT:

Objective:

To generate estimates of mortality rates due to garbage codes (GC) for Brazilian municipalities by comparing the direct and the Bayesian methods, based on deaths registered in the Mortality Information System (SIM) between 2015 and 2017.

Methods:

Data from the SIM were used. The analysis was performed in groups of GC levels 1 and 2, levels 3 and 4, and total GC. Mortality rates were estimated directly and also according to the Bayesian method by applying the Empirical Bayesian Estimator.

Results:

About 38% of GC were estimated and regional differences in mortality rates were observed, higher in the Northeast and Southeast and lower in the South and Midwest regions. The Southeast presented similar rates for the two analyzed groups of GC. The smallest differences between direct and Bayesian method estimates were observed in large cities with a population over 500 thousand inhabitants. Municipalities in the north of the state of Minas Gerais and those in the states of Rio de Janeiro, São Paulo, and Bahia presented high rates at levels 1 and 2.

Conclusion:

There are differences in the quality of the definition of the underlying causes of death, even with the use of Bayesian methodology, which assists in smoothing the rates. The quality of the definition of causes of death is important, as they are associated with the access to and quality of healthcare services and support health planning.

Keywords:
Data accuracy; Mortality registries; Cause of death; Small-area analysis

INTRODUCTION

In health planning, the availability of good-quality data is essential. Mortality data are one of the best known and most used, as they enable to know the health condition by age, sex, place of residence, and causes of death11. Mathers CD, Fat DM, Inoue M, Rao C, Lopez AD. Counting the dead and what they died from: An assessment of the global status of cause of death data. Bull World Health Organ 2005; 83(3): 171-7.. Based on these data, it is possible to identify populations at risk, more frequent causes of death, premature or preventable deaths, temporal trends and, thus, define priorities and interventions. However, for the generated information to adequately support the planning and direction of healthcare actions, with a consequent impact on the population's health profile, it is essential for the databases to have good coverage and quality11. Mathers CD, Fat DM, Inoue M, Rao C, Lopez AD. Counting the dead and what they died from: An assessment of the global status of cause of death data. Bull World Health Organ 2005; 83(3): 171-7..

In this sense, the Brazilian Ministry of Health carried out actions to improve the information on national vital statistics. They are practical for:

  • improving the capture of deaths by the Mortality Information System (Sistema de Informações sobre Mortalidade – SIM), such as investments in the training of healthcare teams, expansion of codification related to causes of death, search and legalization of clandestine cemeteries, awareness of managers regarding the underreport of deaths, expansion of family health teams, and hiring doctors for the inland of the country via Programa Mais Médicos para o Brasil (More Doctors Program), among others22. Malta DC, Teixeira R, Oliveira GMM de, Ribeiro AL. Mortalidade por Doenças Cardiovasculares Segundo o Sistema de Informação sobre Mortalidade e as Estimativas do Estudo Carga Global de Doenças no Brasil, 2000-2017. Arq Bras Cardiol 2020; 115(2): 152-60. https://doi.org/10.36660/abc.20190867
    https://doi.org/10.36660/abc.20190867...
    ,33. Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância de Doenças e Agravos Não Transmissíveis e Promoção da Saúde. Saúde Brasil 2015/2016: uma análise da situação de saúde e da epidemia pelo vírus Zika e por outras doenças transmitidas pelo Aedes aegypti. Brasil: Ministério da Saúde; 2017. 388 p.;

  • reducing deaths from ill-defined causes (chapter 18 of the 10th edition of the International Statistical Classification of Diseases and Related Health Problems – ICD-10) and from other garbage codes (GC)44. França EB, Ishitani LH, Teixeira RA, da Cunha CC, Marinho MF. Improving the usefulness of mortality data: Reclassification of ill-defined causes based on medical records and home interviews in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): e1900010.supl.3. https://doi.org/10.1590/1980-549720190010.supl.3
    https://doi.org/10.1590/1980-54972019001...
    and groups due to underlying cause of death (UCD) deemed incorrect or nonspecific such as UCD declared as sepsis and cardiac arrest.

Thus, the group of causes called “GC” is considered an indicator of the quality of health information. The lower the incidence, the better the quality of these data44. França EB, Ishitani LH, Teixeira RA, da Cunha CC, Marinho MF. Improving the usefulness of mortality data: Reclassification of ill-defined causes based on medical records and home interviews in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): e1900010.supl.3. https://doi.org/10.1590/1980-549720190010.supl.3
https://doi.org/10.1590/1980-54972019001...
66. Marinho MF, França EB, Teixeira RA, Ishitani LH, da Cunha CC, dos Santos MR, et al. Data for health: Impact on improving the quality of cause-of-death information in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): 1-13. https://doi.org/10.1590/1980-549720190005.supl.3
https://doi.org/10.1590/1980-54972019000...
.

The magnitude of GC has been analyzed according to regions, states, capitals, and groups of municipalities44. França EB, Ishitani LH, Teixeira RA, da Cunha CC, Marinho MF. Improving the usefulness of mortality data: Reclassification of ill-defined causes based on medical records and home interviews in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): e1900010.supl.3. https://doi.org/10.1590/1980-549720190010.supl.3
https://doi.org/10.1590/1980-54972019001...
,77. Teixeira RA, Naghavi M, Guimarães MDC, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: Garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol 2019; 22(Supl. 3): e190002.supl.3. https://doi.org/10.1590/1980-549720190002.supl.3
https://doi.org/10.1590/1980-54972019000...
. Knowing its distribution at the municipal level is important for planning local actions and reducing the occurrence of deaths from these causes. Studies from the Global Burden of Disease (GBD) have pointed out that poorer countries and locations have worse health indicators and lower quality of databases as well as a higher proportion of GC88. Naghavi M, Makela S, Foreman K, O’Brien J, Pourmalek F, Lozano R. Algorithms for enhancing public health utility of national causes-of-death data. Popul Health Metr 2010; 8: 9. https://doi.org/10.1186/1478-7954-8-9
https://doi.org/10.1186/1478-7954-8-9...
. Also in Brazil, the Busca Ativa de Óbitos study (Proactive Search of Deaths) has already identified worse SIM quality in small municipalities of the North and Northeast regions of the country99. Almeida W da S, Szwarcwald CL. Adequação das informações de mortalidade e correção dos óbitos informados a partir da Pesquisa de Busca Ativa. Ciênc Saúde Coletiva 2017; 22(10): 3193-203. https://doi.org/10.1590/1413-812320172210.12002016
https://doi.org/10.1590/1413-81232017221...
.

It is assumed that in smaller municipalities and in poorer regions the proportion of GC is higher. Considering that, of the 5,570 municipalities in the country, 88% had a population of less than 50 thousand inhabitants in 20101010. Instituto Brasileiro de Geografia e Estatística. Perfil dos municípios [Internet]. Instituto Brasileiro de Geografia e Estatística; 2020 [access on Mar. 15, 2020]. Available from: https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
https://www.ibge.gov.br/estatisticas/dow...
, it is a great challenge to directly calculate mortality rates in these locations due to instability and great variability in the estimates. Thus, as an alternative, some authors have been using Bayesian methods, such as the Empirical Bayesian Estimator, in order to estimate mortality rates in municipalities1111. Pinheiro PC, Queiroz BL, Teixeira RA, Ribeiro ALP, Malta DC. A mortalidade feminina por acidentes de motocicleta nos municípios brasileiros, 2005, 2010 e 2015. Rev Bras Epidemiol 2020; 23(Supl. 1): e200010.SUPL.1. https://doi.org/10.1590/1980-549720200010.supl.1
https://doi.org/10.1590/1980-54972020001...
,1212. Assunção RM, Barreto SM, Guerra HL, Sakurai E. Mapas de taxas epidemiológicas: uma abordagem Bayesiana. Cad Saúde Pública 1998; 14(4): 713-23. https://doi.org/10.1590/S0102-311X1998000400013
https://doi.org/10.1590/S0102-311X199800...
.

Therefore, this study aimed at generating estimates of mortality rates due to GC for Brazilian municipalities by comparing the direct and the Bayesian methods, based on deaths registered in SIM between 2015 and 2017.

METHODS

Descriptive study using public data from SIM from 2015 to 2017. The analysis considered the municipality of residence and, in order to minimize fluctuations in the number of deaths at the municipal level, where small numbers generate high variability in the rates, it was decided to work with the three-year period.

The selection of GC was based on the GBD 2017 study1313. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
https://doi.org/10.1016/S0140-6736(18)32...
, which classified the defined causes into three major groups: communicable, maternal, neonatal, and nutritional; non-communicable diseases (NCD); and external causes. In addition to these groups, the GBD study defines four levels of GC, with level 1 being the worst scenario in relation to the quality of the definition of causes in this order. They are as follows1414. Naghavi M, Richards N, Chowdhury H, Eynstone-Hinkins J, Franca E, Hegnauer M, et al. Improving the quality of cause of death data for public health policy: Are all “garbage” codes equally problematic? BMC Med 2020; 18(1): 55. https://doi.org/10.1186/s12916-020-01525-w
https://doi.org/10.1186/s12916-020-01525...
:

  • Level 1: GCs can be redistributed to any of the large groups of defined causes in the GBD study, as aforementioned. For instance, a GC can be coded as sepsis by any group, which can result from a death from transport accident, from an infectious disease, such as pneumonia, or from a chronic disease such as cancer;

  • Level 2: GCs are redistributed to a large group, or at most to a second group (for instance, UCD defined as gastrointestinal bleeding, unspecified, should be redistributed only to the group of non-communicable diseases);

  • Level 3: They refer to causes that are likely to be in the same chapter as ICD-10. For example: unspecified cancer, although requiring greater specificity of type or organ, this GC will be attributed to the disease and redistributed to the same group of specific causes of cancer;

  • Level 4: The UCD probably refers to a single disease, such as unspecified stroke, which may be ischemic or hemorrhagic; diabetes, which can be redistributed as type I or type II.

Thus, it is considered that GC levels 1 and 2 are the most problematic and may have greater impact on the quality of statistics on causes of death, as they contain little information on the actual UCD for being highly nonspecific1313. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
https://doi.org/10.1016/S0140-6736(18)32...
,1414. Naghavi M, Richards N, Chowdhury H, Eynstone-Hinkins J, Franca E, Hegnauer M, et al. Improving the quality of cause of death data for public health policy: Are all “garbage” codes equally problematic? BMC Med 2020; 18(1): 55. https://doi.org/10.1186/s12916-020-01525-w
https://doi.org/10.1186/s12916-020-01525...
. For the present study, the analyses were made for the total GC and groups of levels 1 and 2 and 3 and 4.

For the purposes of a more adequate comparison, the limitations of the local quality of UCD were considered when using the SIM at the municipal level as a treatment for improving the quality of data to level the limitations. The treatment of missing data consisted in step 1. A proportional redistribution was applied in the case of missing data on age, sex, or municipality of residence, and the last variable comprised the Federative Unit (FU)22. Malta DC, Teixeira R, Oliveira GMM de, Ribeiro AL. Mortalidade por Doenças Cardiovasculares Segundo o Sistema de Informação sobre Mortalidade e as Estimativas do Estudo Carga Global de Doenças no Brasil, 2000-2017. Arq Bras Cardiol 2020; 115(2): 152-60. https://doi.org/10.36660/abc.20190867
https://doi.org/10.36660/abc.20190867...
.

Correction was applied for unregistered deaths, i.e., correction of underreporting, taking into account the heterogeneity of SIM coverage in the country1515. Szwarcwald CL, de Frias PG, de Souza Júnior PRB, da Silva de Almeida W, de Morais Neto OL. Correction of vital statistics based on a proactive search of deaths and live births: Evidence from a study of the North and Northeast regions of Brazil. Popul Health Metrics 2014; 12: 16. https://doi.org/10.1186/1478-7954-12-16
https://doi.org/10.1186/1478-7954-12-16...
. GBD 2017 correction was used according to sex, age, and FU. The correction coefficient was generated by the ratio between deaths estimated by the GBD and those observed in the SIM according to states, but it was applied in the municipalities with a general mortality rate of less than five deaths per 100 thousand inhabitants, avoiding overestimations in municipalities whose death rate was classified as of good quality1616. Abouzahr C, Mikkelsen L, Rampatige R, Lopez A. Mortality statistics: a tool to improve understanding and quality. Health Information Systems; 2010.. Municipalities whose value was less than 1 were not corrected.

The estimates of mortality rates were prepared with the Empirical Bayesian Estimator (EBE)1212. Assunção RM, Barreto SM, Guerra HL, Sakurai E. Mapas de taxas epidemiológicas: uma abordagem Bayesiana. Cad Saúde Pública 1998; 14(4): 713-23. https://doi.org/10.1590/S0102-311X1998000400013
https://doi.org/10.1590/S0102-311X199800...
,1717. de Carvalho AXY, da Silva GDM, de Almeida Júnior GR, de Albuquerque PHM. Taxas bayesianas para o mapeamento de homicídios nos municípios brasileiros. Cad Saúde Pública 2012; 28(7): 1249-62. https://doi.org/10.1590/S0102-311X2012000700004
https://doi.org/10.1590/S0102-311X201200...
1919. Marshall RJ. Mapping Disease and Mortality Rates Using Empirical Bayes Estimators. J R Stat Soc Ser C Appl Stat 1991; 40(2): 283-94.. This methodology considers the neighborhood distribution, allowing to minimize the effects caused by the small numbers in the denominator analyzed in small populations. In these municipalities, one death can considerably impact the mortality rates estimated in the localities. In addition, EBE allows estimating rates in places without death incidence, which enables to calculate risks in which the observed event is zero, through information from neighboring regions. In this study, the distribution of the eight closest neighbors of the evaluated municipalities was considered.

Similar to the direct method, for the calculation of mortality rates, Bayesian estimators have as parameters the number of deaths and the population; however, it is known that advanced ages have high mortality rates due to GC. For minimizing the effect of the age distribution of the municipal population, the standardized rates for age were calculated, considering the standard population of the 2010 Census2020. Instituto Brasileiro de Geografia e Estatística. Censo 2010. Atlas censo demográfico. Brasil: Instituto Brasileiro de Geografia e Estatística; 2013. and using the absolute values of expected deaths. Inhabitants of the municipalities, according to sex and age, were estimated according to the demographic cohort-component method for population projections, with an empirical Bayesian contraction estimator, to minimize the instability in the estimates of differential factors of the method growth in smaller areas2121. Freire FHMA, Gonzaga MR, Queiroz BL. Projeção populacional municipal com estimadores bayesianos, Brasil 2010-2030. In: Sawyer DO, editor. Seguridade Social Municipais. Projeto Brasil 3 Tempos. Brasília: Secretaria Especial de Assuntos Estratégicos da Presidência da República (SAE/SG/PR), United Nations Development Programme, Brazil (UNDP) and International Policy Centre for Inclusive Growth; 2019.. Thus, the rates estimated by the direct method (crude rate) and the Bayesian method (Bayesian rate) used the expected number of deaths considering the municipal standardized rates and the respective population, during the three-year period from 2015 to 2017.

Municipal descriptive analyses for GC, total, and groups were generated according to regions of Brazil. Histograms were used to show the absolute differences between the rates estimated by the Bayesian method, with the use of EBE, and the direct method, division of the expected number of deaths and population, in the municipalities and according to population size. Finally, maps were prepared containing the spatial distributions of Brazilian municipalities with the estimated Bayesian mortality rates.

The current study was approved by the Ethics Committee for Research involving Human Beings of Universidade Federal de Minas Gerais, part of the project Desigualdades em Pequenas Áreas Geográficas dos Indicadores de Doenças Crônicas Não Transmissíveis, Violências e seus Fatores de Risco [Inequalities in Small Geographic Areas of the Indicators of Chronic Non-communicable Diseases, Violence, and their Risk Factors].

RESULTS

1,508,021 deaths from GC were estimated in Brazil in the three-year period from 2015 to 2017, representing 38.6% of the 3,909,339 deaths registered in the SIM. Among the GC, 45.1% were of levels 1 and 2. The descriptive Bayesian mortality rates shown in Table 1 demonstrate significant regional differences. The Northeast and Southeast regions had the highest Bayesian rates by total GC: 262.2 (260.1; 264.3) and 253.8 (251.8; 255.9) deaths per 100 thousand inhabitants, respectively. On the other hand, the lowest total GC rates were observed in the South and Midwest regions, with 192.3 (189.9; 194.6) and 197.8 (194.2; 201.4), respectively.

Table 1
Descriptive statistics, mean and respective 95% confidence interval, standard deviation, median, minimum and maximum of Bayesian mortality rates due to garbage codes, total, and groups. Brazilian municipalities and regions, 2015 to 2017.

Among the GC subgroups, GC levels 1 and 2 in the Southeast stand out, with 125.7 (123.8; 127.5), for having presented a rate value similar to that of levels 3 and 4, with 128.6 (127.7; 129.5). It was also verified that the means of rates estimated by the Bayesian method reached higher values than those estimated by the direct method, but with less variability, as the standard deviation and amplitude were lower (descriptive not shown in Table 1). The states with the lowest mortality rates for GC per region were: Amapá (177.1), Rio Grande do Norte (224.3), Espírito Santo (169.3), Rio Grande do Sul (176.8), and the Federal District (140.5), followed by North, Northeast, Southeast, Midwest, and South regions, which are not in the table.

Figure 1 shows the means and 95% confidence intervals of the differences in Bayesian and direct rates in Brazil and its regions. The Midwest, the North, and the South regions had the highest means. GC levels 3 and 4 showed the greatest mean differences.

Figure 1
Mean of the absolute differences between the Bayesian mortality rates and the direct method and respective 95% confidence intervals. Brazil and regions, 2015 to 2017.

Figure 2 shows the histograms of the direct and Bayesian mortality rates by GC and groups. The decrease in the variability of mortality rates after the correction is noteworthy, considering that the frequency is more concentrated in the center of the histogram for this indicator.

Figure 2
Histogram of mortality rates, by estimation methods, for (A) GC levels 1 and 2; (B) GC levels 3 and 4; and (C) Total GC. Brazilian municipalities, 2015 to 2017.

Figure 3 indicates the histograms of the absolute differences between the Bayesian rates and the crude rates for the total GC and groups according to Brazilian regions, in addition to differentiating the size of the municipalities by color. Large cities have differences close to zero, especially those with a population over 500 thousand inhabitants. Conversely, smaller municipalities have distributions throughout the x-axis, i.e., a more heterogeneous distribution. There is a large volume of municipalities with a population of less than 10 thousand inhabitants in the South of the country, which stand out with values that are more distant from zero, especially for this region. In the Northeast and Southeast regions, high frequencies above zero are observed in cities with 10 to 50 thousand inhabitants for the differences in the rates of total GC.

Figure 3
Histogram of the absolute differences in mortality rates according to the Bayesian and direct methods, per size of municipalities, by estimation methods, for (A) GC levels 1 and 2; (B) GC levels 3 and 4; (C) Total GC. Brazilian municipalities, 2015 to 2017.

The municipal analysis of the geographic distributions of mortality rates by total GC and investigated subgroups is shown in Figure 1. In the caption it is possible to observe that the darker the color, the higher the mortality rate in the municipality. When analyzing the total GC (Figure 1C), the highest concentration of dark colors in the Southeast region is highlighted, mainly in the north of Minas Gerais and in the states of Rio de Janeiro and São Paulo, and in the Northeast region, mainly in the state of Bahia. On the other hand, the state of Espírito Santo and the South and Midwest regions stand out for the presence of lighter colors, i.e., lower rates. When observing the map of GC rates concerning levels 1 and 2, a geographic distribution similar to that of total GC is verified (Figure 1A).

Finally, when analyzing the distribution of mortality rates of GC levels 3 and 4 (Figure 1B), it is possible to observe a random distribution throughout the country, with no visual patterns being identified, not even in Espírito Santo, which stood out for having presented a standard of lowest rates for total GC and GC levels 1 and 2.

DISCUSSION

The results of the present study highlight the high proportion of GC in the country, with more than a third of deaths thus classified, being distributed with 12.9, 4.5, 4.1, and 17.1% among the respective levels of 1 to 4 in the three-year period from 2015 to 2017. The Northeast and Southeast regions had the highest Bayesian death rates from total GC, and the lowest rates of total GC were observed in the South and Midwest regions. Smaller municipalities concentrate higher rates of GC. The heterogeneity of the quality of mortality data in the period from 2015 to 2017, considering GC as quality indicators, is evident.

The importance of correcting and using the EBE method to calculate mortality rates for GC in Brazilian municipalities is emphasized. These results may be related to difficulties of access to health and the scarcity of resources in health care, including the quality of provided services and diagnoses, as these are factors that negatively interfere in the accuracy of the definition of UCD2222. Mahapatra P, Shibuya K, Lopez AD, Coullare F, Notzon FC, Rao C, et al. Civil registration systems and vital statistics: successes and missed opportunities. Lancet 2007; 370(9599): 1653-63. https://doi.org/10.1016/s0140-6736(07)61308-7
https://doi.org/10.1016/s0140-6736(07)61...
. Furthermore, the use of EBE for small areas minimizes fluctuations, considering the observed regional realities, as the neighboring municipalities are taken into account to calculate the estimates of final rates.

Although SIM has been considered as a source of good-quality data in recent years, this characteristic is regionally differentiated. In addition, the analysis of mortality rate in small populations can generate high variability, as small numbers can considerably change the mortality rate. Thus, some methodological aspects of the present study should be highlighted.

The treatment of the raw SIM data and the methodology used to work with small areas allowed the analysis of mortality rates, minimizing the random fluctuations in the spatial distribution of rates between municipalities. First, the redistribution of missing data and the treatment of SIM underreporting of deaths brought the level of quality of raw data closer to the municipal level. As the use of correction coefficients developed for states in GBD studies1313. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
https://doi.org/10.1016/S0140-6736(18)32...
can generate values that do not correspond to the municipal reality, corrections were applied in this study only to cities where the overall mortality rates were considered lower than expected1616. Abouzahr C, Mikkelsen L, Rampatige R, Lopez A. Mortality statistics: a tool to improve understanding and quality. Health Information Systems; 2010.. In addition, the effect of differences in age distribution was removed by using age-standardized rates. Finally, the use of EBE, when considering information from neighbors in the rate estimates, enabled estimates without random spatial changes.

Smoothing in the estimates of mortality rates by the Bayesian method, which takes into account the neighboring municipalities to generate the estimates, applied in this study, proved to be adequate, considering the heterogeneity in the quality of the mortality data and the large number of small cities in Brazil1010. Instituto Brasileiro de Geografia e Estatística. Perfil dos municípios [Internet]. Instituto Brasileiro de Geografia e Estatística; 2020 [access on Mar. 15, 2020]. Available from: https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
https://www.ibge.gov.br/estatisticas/dow...
.

As expected, Bayesian estimates showed less variability than the direct ones in the methodology for calculating rates, as in smaller municipalities significant correction is expected due to the weight of larger neighboring cities1919. Marshall RJ. Mapping Disease and Mortality Rates Using Empirical Bayes Estimators. J R Stat Soc Ser C Appl Stat 1991; 40(2): 283-94.. Furthermore, when dealing with small numbers, in the direct method it was observed the nonoccurrence of deaths from GC, which means that the risk of mortality equals zero. This may not be in line with the local reality, considering that fluctuations caused by small numbers can interfere with rate estimates with the use of direct methods1212. Assunção RM, Barreto SM, Guerra HL, Sakurai E. Mapas de taxas epidemiológicas: uma abordagem Bayesiana. Cad Saúde Pública 1998; 14(4): 713-23. https://doi.org/10.1590/S0102-311X1998000400013
https://doi.org/10.1590/S0102-311X199800...
.

Based on estimates that use the values of the neighboring municipalities for the size of the population, greater differences between the rates estimated by the Bayesian and direct methods were observed, especially in smaller municipalities. This can be verified in Figure 3, in which there is considerable frequency of the difference between the rates of the direct and Bayesian methods far from zero in the municipalities with smaller population sizes.

The analysis per regions also enabled to verify that the quality of the mortality data, considering total GC, overall, presents consubstantial regional differences. There were lower rates in the South and Midwest compared with other regions. These results corroborate the hypothesis that more developed areas have better quality in the most appropriate definition of causes of death. Access to health and the quality of these services are factors that can contribute to the better definition of UCD2323. Abreu DMX de, Sakurai E, Campos LN. A evolução da mortalidade por causas mal definidas na população idosa em quatro capitais brasileiras, 1996-2007. Rev Bras Estud Popul 2010; 27(1): 75-88. https://doi.org/10.1590/S0102-30982010000100006
https://doi.org/10.1590/S0102-3098201000...
,2424. Kanso S, Romero DE, da Costa Leite I, de Moraes EN. Diferenciais geográficos, socioeconômicos e demográficos da qualidade da informação da CB de morte dos idosos no Brasil. Cad Saúde Pública 2011; 27(7): 1323-39. https://doi.org/10.1590/S0102-311X2011000700008
https://doi.org/10.1590/S0102-311X201100...
.

An unexpected finding was that the Southeast region had the second highest Bayesian mortality rate due to total GC, only lower than that of the Northeast. Moreover, it was noted that the rates of GC groups levels 1 and 2 and 3 and 4 showed very close values. The GC levels 1 and 2 group represents deaths with little information for an adequate definition of UCD. In this group, according to GBD 20171313. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
https://doi.org/10.1016/S0140-6736(18)32...
, there are deaths to be redistributed among all defined causes such as, for example, R98 – unattended death – and R99 – other ill-defined and unspecified causes of mortality. Causes like these do not provide information to support health management, as they do not enable to target preventive actions for them. Considering that the Southeast region is one of the most economically developed in the country2525. Instituto de Pesquisa Econômica Aplicada. Atlas do Desenvolvimento Humano. PNUD; Fundação João Pinheiro; Instituto de Pesquisa Econômica Aplicada; 2013., with higher number of larger municipalities and, therefore, with better access to healthcare services on the part of the population2626. Stopa SR, Malta DC, Monteiro CN, Szwarcwald CL, Goldbaum M, Cesar CLG. Use of and access to health services in Brazil, 2013 National Health Survey. Rev Saúde Pública 2017; 51(Supl. 1): 1S-10S. https://doi.org/10.1590/s1518-8787.2017051000074
https://doi.org/10.1590/s1518-8787.20170...
, together with the South region, a hypothesis to be raised is that the high population density is preventing a more accurate diagnosis for the definition of UCD.

Regional differences can affect the comparability of mortality indicators for specific groups of causes of death that are incorrectly classified as GC. Such differences may still be due to the different types of GC, because, depending on the location, there may be a predominance of GC more related to the group of communicable diseases, chronic non-communicable diseases, or external causes2727. GBD 2016 Causes of Death Collaborators. Global, regional, and national age-sex specifc mortality for 264 causes of death, 1980-2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet 2017; 390(10100): 1151-210. https://doi.org/10.1016/S0140-6736(17)32152-9
https://doi.org/10.1016/S0140-6736(17)32...
. In addition, it is very likely that, depending on the FU or municipality, there will be variations in the certification and coding of the causes of death2828. Danilova I, Shkolnikov VM, Jdanov DA, Meslé F, Vallin J. Identifying potential differences in cause-of-death coding practices across Russian regions. Popul Health Metrics 2016; 14: 8. https://doi.org/10.1186/s12963-016-0078-0
https://doi.org/10.1186/s12963-016-0078-...
.

According to the GBD 2017 study1313. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
https://doi.org/10.1016/S0140-6736(18)32...
, Brazil accounts for an estimate of 17.2% of major GC, whereas Latin American countries, such as Argentina, Peru and Uruguay, have rates of 29.2, 20.4, and 22.2%, respectively, and even in Europe, such as France (22.3%), there are worse scenarios. Conversely, countries with better proportions are also observed such as Chile (10.2%) and Mexico (11.9%).

Despite the improvement already observed in the quality of mortality data66. Marinho MF, França EB, Teixeira RA, Ishitani LH, da Cunha CC, dos Santos MR, et al. Data for health: Impact on improving the quality of cause-of-death information in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): 1-13. https://doi.org/10.1590/1980-549720190005.supl.3
https://doi.org/10.1590/1980-54972019000...
,2929. Balieiro PC da S, Silva LCF da, Sampaio V de S, Monte EX do, Pereira EM dos S, Queiroz LAF de, et al. Fatores associados à mortalidade por causas inespecíficas e mal definidas no estado do Amazonas, Brasil, de 2006 a 2012. Ciênc Saúde Coletiva 2020; 25(1): 339-52. https://doi.org/10.1590/1413-81232020251.27182017
https://doi.org/10.1590/1413-81232020251...
, the analysis per municipalities shows that mortality rates by GC are heterogeneous in the country. This analysis (Figure 4) shows spatial distribution containing groups of regionalized cities and points out areas with high rates, demonstrating intraregional inequalities. Municipalities with the highest rates were concentrated in northern Minas Gerais and southern Bahia, areas that presented groups of municipalities with values classified in the same category as the range of causes by total GC, represented by a darker color. These areas are characterized as of poor socioeconomic development3030. Albuquerque MV de, Viana AL d’Ávila, Lima LD de, Ferreira MP, Fusaro ER, Iozzi FL. Desigualdades regionais na saúde: mudanças observadas no Brasil de 2000 a 2016. Ciênc Saúde Coletiva 2017; 22(4): 1055-64. https://doi.org/10.1590/1413-81232017224.26862016
https://doi.org/10.1590/1413-81232017224...
. Possibly, these regions, in addition to being more distant from the capitals, which are references, are located in areas distant from regional healthcare centers, which increases the chance of nonspecific diagnoses2929. Balieiro PC da S, Silva LCF da, Sampaio V de S, Monte EX do, Pereira EM dos S, Queiroz LAF de, et al. Fatores associados à mortalidade por causas inespecíficas e mal definidas no estado do Amazonas, Brasil, de 2006 a 2012. Ciênc Saúde Coletiva 2020; 25(1): 339-52. https://doi.org/10.1590/1413-81232020251.27182017
https://doi.org/10.1590/1413-81232020251...
. On the contrary, municipalities in the South region, where places with higher socioeconomic development and high supply/complexity of healthcare services are concentrated, showed lower rates, represented by a lighter color3030. Albuquerque MV de, Viana AL d’Ávila, Lima LD de, Ferreira MP, Fusaro ER, Iozzi FL. Desigualdades regionais na saúde: mudanças observadas no Brasil de 2000 a 2016. Ciênc Saúde Coletiva 2017; 22(4): 1055-64. https://doi.org/10.1590/1413-81232017224.26862016
https://doi.org/10.1590/1413-81232017224...
.

Figure 4
Geographic distribution of the Bayesian mortality rates by total GC, GC levels 1 and 2, and GC levels 3 and 4 in Brazilian municipalities, 2015 to 2017.

In the analysis by types of GC (Figures 4A and 4B), it is verified that GC levels 3 and 4 have a heterogeneous distribution throughout the country, whereas GC levels 1 and 2 (the most serious ones) have few, but important points with high rates in the South and Midwest regions. These findings enable to consider the need to locally prioritize a plan aiming at reducing deaths certified as GC such as better access to healthcare services and diagnosis and the improvement of death surveillance.

Although there was no correlation analysis, visually, it is clear that other factors may be related to the generated clusters of municipalities such as the coverage of Programa Saúde da Família (Family Health Program), socioeconomic level, among others. This suggests an in-depth study, aiming to evaluate the relationship between the findings and the aforementioned characteristics of these places. To better understand this situation, it is suggested to use more analytical methodologies, using different data sources to identify significantly associated factors.

The results of this study show differences observed in the quality of the definition of UCD with the use of an adequate methodology for the analysis of small areas. When investigating differences between the crude and estimated rates, with the aid of the Bayesian methodology, it was possible to verify that the situation of GC in Brazil, even with methodology for smoothing it, is still regionally differentiated. The quality of the definition of UCD is extremely important for public health, considering that it is associated with the access to and quality of healthcare services and supports health planning. Therefore, an analysis of small areas is very important for the actors responsible for health management in the country.

  • Financial support: Small Areas Project, Department of Health Surveillance, Ministry of Health (TED 148-2018)

REFERENCES

  • 1
    Mathers CD, Fat DM, Inoue M, Rao C, Lopez AD. Counting the dead and what they died from: An assessment of the global status of cause of death data. Bull World Health Organ 2005; 83(3): 171-7.
  • 2
    Malta DC, Teixeira R, Oliveira GMM de, Ribeiro AL. Mortalidade por Doenças Cardiovasculares Segundo o Sistema de Informação sobre Mortalidade e as Estimativas do Estudo Carga Global de Doenças no Brasil, 2000-2017. Arq Bras Cardiol 2020; 115(2): 152-60. https://doi.org/10.36660/abc.20190867
    » https://doi.org/10.36660/abc.20190867
  • 3
    Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância de Doenças e Agravos Não Transmissíveis e Promoção da Saúde. Saúde Brasil 2015/2016: uma análise da situação de saúde e da epidemia pelo vírus Zika e por outras doenças transmitidas pelo Aedes aegypti. Brasil: Ministério da Saúde; 2017. 388 p.
  • 4
    França EB, Ishitani LH, Teixeira RA, da Cunha CC, Marinho MF. Improving the usefulness of mortality data: Reclassification of ill-defined causes based on medical records and home interviews in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): e1900010.supl.3. https://doi.org/10.1590/1980-549720190010.supl.3
    » https://doi.org/10.1590/1980-549720190010.supl.3
  • 5
    Murray CJL, Lopez, Alan D. The Global burden of disease: a comprehensive assessment of mortality and disability from diseases, injuries, and risk factors in 1990 and projected to 2020. World Health Organization; 1996.
  • 6
    Marinho MF, França EB, Teixeira RA, Ishitani LH, da Cunha CC, dos Santos MR, et al. Data for health: Impact on improving the quality of cause-of-death information in Brazil. Rev Bras Epidemiol 2019; 22(Supl. 3): 1-13. https://doi.org/10.1590/1980-549720190005.supl.3
    » https://doi.org/10.1590/1980-549720190005.supl.3
  • 7
    Teixeira RA, Naghavi M, Guimarães MDC, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: Garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol 2019; 22(Supl. 3): e190002.supl.3. https://doi.org/10.1590/1980-549720190002.supl.3
    » https://doi.org/10.1590/1980-549720190002.supl.3
  • 8
    Naghavi M, Makela S, Foreman K, O’Brien J, Pourmalek F, Lozano R. Algorithms for enhancing public health utility of national causes-of-death data. Popul Health Metr 2010; 8: 9. https://doi.org/10.1186/1478-7954-8-9
    » https://doi.org/10.1186/1478-7954-8-9
  • 9
    Almeida W da S, Szwarcwald CL. Adequação das informações de mortalidade e correção dos óbitos informados a partir da Pesquisa de Busca Ativa. Ciênc Saúde Coletiva 2017; 22(10): 3193-203. https://doi.org/10.1590/1413-812320172210.12002016
    » https://doi.org/10.1590/1413-812320172210.12002016
  • 10
    Instituto Brasileiro de Geografia e Estatística. Perfil dos municípios [Internet]. Instituto Brasileiro de Geografia e Estatística; 2020 [access on Mar. 15, 2020]. Available from: https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
    » https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
  • 11
    Pinheiro PC, Queiroz BL, Teixeira RA, Ribeiro ALP, Malta DC. A mortalidade feminina por acidentes de motocicleta nos municípios brasileiros, 2005, 2010 e 2015. Rev Bras Epidemiol 2020; 23(Supl. 1): e200010.SUPL.1. https://doi.org/10.1590/1980-549720200010.supl.1
    » https://doi.org/10.1590/1980-549720200010.supl.1
  • 12
    Assunção RM, Barreto SM, Guerra HL, Sakurai E. Mapas de taxas epidemiológicas: uma abordagem Bayesiana. Cad Saúde Pública 1998; 14(4): 713-23. https://doi.org/10.1590/S0102-311X1998000400013
    » https://doi.org/10.1590/S0102-311X1998000400013
  • 13
    GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018; 392(10159): 1736-88. https://doi.org/10.1016/S0140-6736(18)32203-7
    » https://doi.org/10.1016/S0140-6736(18)32203-7
  • 14
    Naghavi M, Richards N, Chowdhury H, Eynstone-Hinkins J, Franca E, Hegnauer M, et al. Improving the quality of cause of death data for public health policy: Are all “garbage” codes equally problematic? BMC Med 2020; 18(1): 55. https://doi.org/10.1186/s12916-020-01525-w
    » https://doi.org/10.1186/s12916-020-01525-w
  • 15
    Szwarcwald CL, de Frias PG, de Souza Júnior PRB, da Silva de Almeida W, de Morais Neto OL. Correction of vital statistics based on a proactive search of deaths and live births: Evidence from a study of the North and Northeast regions of Brazil. Popul Health Metrics 2014; 12: 16. https://doi.org/10.1186/1478-7954-12-16
    » https://doi.org/10.1186/1478-7954-12-16
  • 16
    Abouzahr C, Mikkelsen L, Rampatige R, Lopez A. Mortality statistics: a tool to improve understanding and quality. Health Information Systems; 2010.
  • 17
    de Carvalho AXY, da Silva GDM, de Almeida Júnior GR, de Albuquerque PHM. Taxas bayesianas para o mapeamento de homicídios nos municípios brasileiros. Cad Saúde Pública 2012; 28(7): 1249-62. https://doi.org/10.1590/S0102-311X2012000700004
    » https://doi.org/10.1590/S0102-311X2012000700004
  • 18
    Justino JR, Freire FHM de A, Lucio PS. Estimação de sub-registros de óbitos em pequenas áreas com os métodos bayesiano empírico e algoritmo EM. Rev Bras Estud Popul 2012; 29(1): 87-100. https://doi.org/10.1590/S0102-30982012000100006
    » https://doi.org/10.1590/S0102-30982012000100006
  • 19
    Marshall RJ. Mapping Disease and Mortality Rates Using Empirical Bayes Estimators. J R Stat Soc Ser C Appl Stat 1991; 40(2): 283-94.
  • 20
    Instituto Brasileiro de Geografia e Estatística. Censo 2010. Atlas censo demográfico. Brasil: Instituto Brasileiro de Geografia e Estatística; 2013.
  • 21
    Freire FHMA, Gonzaga MR, Queiroz BL. Projeção populacional municipal com estimadores bayesianos, Brasil 2010-2030. In: Sawyer DO, editor. Seguridade Social Municipais. Projeto Brasil 3 Tempos. Brasília: Secretaria Especial de Assuntos Estratégicos da Presidência da República (SAE/SG/PR), United Nations Development Programme, Brazil (UNDP) and International Policy Centre for Inclusive Growth; 2019.
  • 22
    Mahapatra P, Shibuya K, Lopez AD, Coullare F, Notzon FC, Rao C, et al. Civil registration systems and vital statistics: successes and missed opportunities. Lancet 2007; 370(9599): 1653-63. https://doi.org/10.1016/s0140-6736(07)61308-7
    » https://doi.org/10.1016/s0140-6736(07)61308-7
  • 23
    Abreu DMX de, Sakurai E, Campos LN. A evolução da mortalidade por causas mal definidas na população idosa em quatro capitais brasileiras, 1996-2007. Rev Bras Estud Popul 2010; 27(1): 75-88. https://doi.org/10.1590/S0102-30982010000100006
    » https://doi.org/10.1590/S0102-30982010000100006
  • 24
    Kanso S, Romero DE, da Costa Leite I, de Moraes EN. Diferenciais geográficos, socioeconômicos e demográficos da qualidade da informação da CB de morte dos idosos no Brasil. Cad Saúde Pública 2011; 27(7): 1323-39. https://doi.org/10.1590/S0102-311X2011000700008
    » https://doi.org/10.1590/S0102-311X2011000700008
  • 25
    Instituto de Pesquisa Econômica Aplicada. Atlas do Desenvolvimento Humano. PNUD; Fundação João Pinheiro; Instituto de Pesquisa Econômica Aplicada; 2013.
  • 26
    Stopa SR, Malta DC, Monteiro CN, Szwarcwald CL, Goldbaum M, Cesar CLG. Use of and access to health services in Brazil, 2013 National Health Survey. Rev Saúde Pública 2017; 51(Supl. 1): 1S-10S. https://doi.org/10.1590/s1518-8787.2017051000074
    » https://doi.org/10.1590/s1518-8787.2017051000074
  • 27
    GBD 2016 Causes of Death Collaborators. Global, regional, and national age-sex specifc mortality for 264 causes of death, 1980-2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet 2017; 390(10100): 1151-210. https://doi.org/10.1016/S0140-6736(17)32152-9
    » https://doi.org/10.1016/S0140-6736(17)32152-9
  • 28
    Danilova I, Shkolnikov VM, Jdanov DA, Meslé F, Vallin J. Identifying potential differences in cause-of-death coding practices across Russian regions. Popul Health Metrics 2016; 14: 8. https://doi.org/10.1186/s12963-016-0078-0
    » https://doi.org/10.1186/s12963-016-0078-0
  • 29
    Balieiro PC da S, Silva LCF da, Sampaio V de S, Monte EX do, Pereira EM dos S, Queiroz LAF de, et al. Fatores associados à mortalidade por causas inespecíficas e mal definidas no estado do Amazonas, Brasil, de 2006 a 2012. Ciênc Saúde Coletiva 2020; 25(1): 339-52. https://doi.org/10.1590/1413-81232020251.27182017
    » https://doi.org/10.1590/1413-81232020251.27182017
  • 30
    Albuquerque MV de, Viana AL d’Ávila, Lima LD de, Ferreira MP, Fusaro ER, Iozzi FL. Desigualdades regionais na saúde: mudanças observadas no Brasil de 2000 a 2016. Ciênc Saúde Coletiva 2017; 22(4): 1055-64. https://doi.org/10.1590/1413-81232017224.26862016
    » https://doi.org/10.1590/1413-81232017224.26862016

Publication Dates

  • Publication in this collection
    16 Apr 2021
  • Date of issue
    2021

History

  • Received
    22 Sept 2020
  • Reviewed
    24 Nov 2020
  • Accepted
    30 Nov 2020
  • Preprint
    14 Dec 2020
Associação Brasileira de Pós -Graduação em Saúde Coletiva São Paulo - SP - Brazil
E-mail: revbrepi@usp.br