Next Article in Journal
Hemophagocytic Lymphohistiocytosis (HLH) in Patients with Tick-Borne Illness: A Scoping Review of 98 Cases
Previous Article in Journal
Respiratory Syncytial Virus (RSV) and Intention to Recommend RSV Vaccination: A Cross-Sectional Survey of Cardiologists and Cardiac Nurses in Southern Italy
Previous Article in Special Issue
Effects of COVID-19-Associated Infection Control on the Pattern of Infections Imported by German Soldiers and Police Officers Returning from Predominantly Tropical Deployment Sites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Excess Mortality Stratified by Age and Sex for Croatia and Croatian Counties during the 2020–2021 COVID-19 Pandemic

1
Department of Clinical Epidemiology, University Hospital of Split, 21000 Split, Croatia
2
Department of Medical Physics and Biophysics, University of Split School of Medicine, 21000 Split, Croatia
3
Division for Epidemiology and Prevention of Noncommunicable Chronic Diseases, Croatian Institute of Public Health, 10000 Zagreb, Croatia
4
Department of Public Health, University of Split School of Medicine, 21000 Split, Croatia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Infect. Dis. Rep. 2024, 16(2), 142-153; https://0-doi-org.brum.beds.ac.uk/10.3390/idr16020011
Submission received: 29 December 2023 / Revised: 7 February 2024 / Accepted: 16 February 2024 / Published: 20 February 2024
(This article belongs to the Special Issue Emerging Infections: Epidemiology, Diagnostics, Clinics and Evolution)

Abstract

:
Excess mortality is often used to estimate the effect of a certain crisis on the population. It is defined as the number of deaths during a crisis exceeding the expected number based on historical trends. Here, we calculated excess mortality due to the COVID-19 pandemic for Croatia in the 2020–2021 period. The excess was calculated on the national and county level for different age and sex categories. In addition to the absolute number, the excess mortality was also expressed as a ratio of excess deaths to the predicted baseline and excess mortality rate. We showed that using both measures is necessary to avoid incorrect conclusions. The estimated excess mortality on the national level was 14,963, corresponding to an excess percentage of 14.3%. With respect to sex, there was a higher excess mortality rate for men compared to women. An exponential relationship was observed between age and the excess mortality rate.These trends wee representative of most counties as well, with large variations in the magnitude of the effect. However, there were also exceptions to the general rule. The reasons for these deviations were discussed in terms of between-county differences in demographic structure, population density and special events that took place during the pandemic.

1. Introduction

The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, was first identified in Wuhan, China, towards the end of 2019. It quickly escalated into a global health emergency, with the World Health Organization declaring it a pandemic in March 2020. This crisis tested public health systems worldwide, disrupted international travel and led to significant economic and social challenges.
The impact of an infectious disease is traditionally estimated based on the number of cases and the number of deaths. However, both of these metrics can be greatly influenced by the limited availability of testing and different definitions of COVID-19 deaths used by different countries [1]. For example, some countries only count PCR-confirmed COVID-19 deaths, while others, such as Croatia, include suspected COVID-19 deaths as well. Additionally, even if counted perfectly, the official number of deaths from COVID-19 only provides a measure of the direct impact of the pandemic and does not take into account the indirect effects on total country mortality. These effects might include deaths caused by the collapse of the medical system, by the psychological effects of isolation, etc.
This is why excess mortality emerged as a crucial indicator during the COVID-19 pandemic [1]. This metric shows the number of deaths exceeding reference levels calculated based on historical mortality data. Assuming the incidence of other diseases remained stable over the period of interest, excess mortality can give us an insight into the number of deaths directly or indirectly related to COVID-19.
Excess mortality can be displayed just as an absolute number of excess deaths, but more commonly, it is presented as relative excess mortality, or specific excess mortality rate. Relative excess mortality is defined as the change in mortality relative to an expected mortality value calculated based on data from some reference time frame. It is often used to express excess mortality since it can be calculated solely from the mortality data and provides information about the relative change in mortality due to an event. However, it may not be suitable in certain situations, such as those where the baseline mortality rates are extremely low, so small changes can result in very high relative excess mortality values that may be misleading. Additionally, it does not explicitly account for the population size, making it harder to compare the absolute impact of an event in different locations. Excess mortality rate addresses this issue by standardizing the absolute excess mortality using population size, but does not include information on the relative impact of those additional deaths compared to the baseline. Consequently, this study will use both metrics in order to obtain a more complete picture, capturing both the relative change in mortality provided by the excess percentage and the absolute impact of the event provided by the excess rate.
To provide a global perspective on excess mortality during the COVID-19 pandemic, excess mortality calculations for the total country populations (including Croatia) have already been performed by other studies such as the one by Karlinsky and Kobak [2] and the COVID-19 Excess Mortality Collaborators [3]. However, both of those studies were based on the World Mortality Dataset, which possibly included preliminary data subject to future revisions, potentially impacting the overall accuracy. Moreover, since exact population data are not available for each year for every country, the Karlinsky and Kobak study relied on the population size estimates for 2020 from the United Nations World Population Prospect (WPP) dataset [4]. Finally, because of their global scope, the above studies do not include calculations over different age–sex subgroups (on national and subnational levels), and these differences must be considered in order to fully understand the differences in the effects of the COVID-19 pandemic on a certain population. We also found a single study on excess mortality focusing on Croatia [5], but the excess mortality was calculated there only for the year 2020, and only at the national level, without accounting for variations across different counties.
The COVID-19 pandemic had varying impact on different segments of the population, particularly when considering factors like age and sex. Older adults were at much higher risk of severe outcomes compared to younger people [6,7], with an exponential relationship being observed between age and infection fatality rate [7]. With respect to age, there seems to be a significant male bias in mortality from COVID-19 [8,9,10], with an estimated average female/male mortality incidence ratio of 0.7 across European countries.
The goal of this retrospective observational population-based study is to calculate excess mortality in Croatia for the 2020–2021 period based on historical data from the 2016–2019 period. The excess will be calculated separately for different counties and age and sex groups. The variations in excess mortality will be discussed in terms of between-county differences in demographic structure, population density and special events that took place during the pandemic.

2. Materials and Methods

The data used in this study were obtained from the Croatian Bureau of Statistics and Croatian Institute of Public Health [11]. The epidemiological indicator used was excess mortality, defined as the difference between the actual and expected number of deaths in a certain time period. In addition to the cumulative absolute excess mortality, we also expressed the excess as relative excess mortality and excess mortality rate per 1000 people. The specific mortality rates per 1000 people were obtained by dividing the absolute excess mortality by the number of people belonging to a certain age and sex subcategory and multiplying by thousand. Expressed mathematically, R e l a t i v e   e x c e s s   m o r a t l i t y = E x c e s s   m o r t a l i t y E x p e c t e d   m o r t a l i t y , and E x c e s s   m o r t a l i t y   r a t e   p e r   1000   p e o p l e = E x c e s s   m o r t a l i t y C a t e g o r y   p o p u l a t i o n   s i z e × 1000 . The number of people in each subcategory was obtained from the 2021 Census of Population performed by the Croatian Bureau of Statistics [12]. The population was stratified into four age (0–64, 65–74, 75–84 and 85+) and two sex (male, female) categories.
The excess mortality was calculated for the 2020–2021 period based on the baseline predictions obtained using the mortality data from the 2016–2019 period. For estimation of the baseline, we adopted the methodology from Karlinsky and Kobak [2]. Briefly, the baseline for different counties was modeled using a linear regression model with a common slope coefficient, but a different intercept for each month, or expressed mathematically— D m o n t h ,   y e a r = β y e a r + α m o n t h + ϵ . This is the simplest model that takes into account both the between-year (through the slope coefficient) and seasonal trends (through different monthly intercepts). Since the mortality data for age and sex subgroups was available only on a yearly basis, the baseline was calculated as D c a t e g o r y ,   y e a r = β y e a r + α c a t e g o r y + ϵ in those cases.
The variance of excess mortality was calculated by summing the variance due to the baseline uncertainty and variance due to Gaussian noise present even without the pandemic. See the Materials and Methods section in the article by Karlinsky and Kobak for more details [2]. Excess mortality prediction was deemed statistically significant if the z-score was larger than 2, where the z-score was calculated as z = e x c e s s   m o r t a l i t y   V a r ( e x c e s s   m o r a t l i t y ) . All data modeling, data analysis and data visualization was performed using the R programming language (R language for statistical computing version 4.2.1, R Foundation for Statistical Computing, Vienna, Austria) [13].

3. Results

3.1. Excess Mortality on the National Level

3.1.1. Excess Mortality for the Whole Population

On the country level, the total mortality in the 2020–2021 period was 119,706. Subtracting the modeled expected mortality from the recorded deaths, an excess mortality prediction of 14,963 was obtained. Translated into relative excess mortality and excess mortality death rate per 1000 people, the obtained values were 14.3% and 3.86, respectively.

3.1.2. Excess Mortality with Respect to Age and Sex Differences

In addition to calculating excess mortality for the entire country population, we also estimated the excess across different age–sex subgroups. We observed an increased excess mortality rate in men (Figure 1, Table 1), with similar female/male excess mortality rates over different age groups.
In terms of age, the greatest relative increase in mortality appeared in the 65–74 years age group for both sexes. The excess mortality rates for both sexes indicate a continuous increase, with age resembling an exponential function (Figure 1, Table 1).

3.1.3. Excess Mortality vs. Official COVID-19 Mortality

The ratio of expected mortality and official COVID-19 mortality can inform us on the quality of COVID-19 mortality tracking and the magnitude of indirect COVID-19 effects. The Croatian Institute of Public Health reports a total of 13035 deaths due to COVID-19 [14], resulting in an excess mortality/COVID-19 mortality ratio of 1.15. This points to an undercount of COVID-19 induced deaths.

3.2. Excess Mortality over Different Counties

3.2.1. Excess Mortalities over Entire County Populations

This section analyzes the total excess mortality across Croatian counties. To reduce clutter on the charts, the county names are abbreviated, but the full names are available in Table A1. The results indicate large variability in the effect of the pandemic across counties (Figure 2). Both metrics unanimously show that the pandemic hit the KK and SK counties the most, and the effect was lowest for the PS and SM counties.
Figure 2 also shows that using only a single metric is often not enough to reach the correct conclusion on the impact of the pandemic. Using only the relative excess mortality or excess mortality rate would be enough if county rankings were always equal for both metrics (yielding horizontal connection lines in the middle). However, we can see that this is often not the case, and there are many counties for which the rankings differ strongly across the two metrics.
Counties for which the relative excess mortality is low and excess mortality rate is high (large positive slopes of gray lines in the middle, for example, the LS and KA counties) could indicate a high historical mortality baseline. Therefore, even a significant increase in excess mortality might not markedly change the relative excess mortality. A high historical mortality rate could be attributed to an older population demographic for such counties (Table 2). Conversely, counties with high relative excess mortality but low excess mortality rate suggest a low mortality baseline. Such examples are the GZ, SD and DN counties.

3.2.2. Excess Mortality across Counties with Respect to Age and Sex Differences

Just like with national-level data, we calculated excess mortality stratified by age and sex subgroups across different counties. In terms of excess mortality differences across sexes, the results again indicate a male bias in most counties, in agreement with the national trend (Figure 3). Stratifying according to age and sex also reveals that the main reason for the high excess mortality rate in the KK county is the large excess mortality rate in the 75–84 years age group, where it is almost twice as high as that of the county with the second highest value for that age group.
With respect to age, we again see an approximate exponential dependance of excess mortality rate on age for both sexes in most counties (Figure 4 and Figure 5). A notable exception to the general rule is the LS county, which displays negative excess percentages and mortality rates in the oldest age category.
Both Figure 4 and Figure 5 again point to a large difference in the behavior of the two excess mortality metrics, highlighting the necessity of simultaneously using both to avoid misconclusions.

4. Discussion

This study presents the excess mortality numbers for Croatia during the 2020–2021 COVID-19 pandemic. The excess was calculated both at the national and county levels, taking into account the age and sex structure of the population.
An excess mortality prediction of 14,963 was obtained for the whole country. Translated into cumulative relative excess mortality and excess mortality rate per 1000 people, that corresponds to values of 14.3% and 3.86, respectively. Other studies calculating excess mortality for Croatia over the same time span offer a wide range of estimates—from 12,205 to 22,900 [2,3,15,16,17,18]. These differences can be attributed to the usage of different statistical models and/or different time spans for the estimation of excess mortality. Also, all of those studies based their estimates on the Human Mortality Database [19,20] which (depending on the time of access and the time span in question) might include preliminary data subject to further revision. We obtained our mortality data from the Croatian Health Statistics Yearbook published by the Croatian Institute of Public Health [11]. Furthermore, their population estimates for Croatia were obtained from the United Nations World Population Prospects [4], and we used the 2021 Population Census data [12] for better population estimate accuracy.
In terms of comparison to other European countries, estimates of relative excess mortality based on the methodology by Karlinsky and Kobak show that Albania had the highest score of 34%, and Denmark and Iceland the lowest score of only 1% [17]. Croatia was positioned somewhere in the middle with a value of 16%, outperforming most Southeastern and Eastern European countries and scoring worse than most Western European countries. Since Western European countries have, on average, a higher per capita GDP than Eastern European countries, this indicates a possible role of a country’s economic status on the impact of COVID-19 [21,22].
A greater excess mortality rate was seen in men, in line with results from other studies for other countries [9,10] and Croatia [5]. The female/male excess mortality rates over different age groups on the national level were in agreement with another study on excess mortality across different countries [8]. An exponential relationship of excess mortality rate with age was observed, in accordance with the results from recent meta-analyses [6,7]. The relative increase in mortality was greatest for the 65–74 years age group, similar to the findings from another study focusing on excess mortality in Croatia, but only for the year 2020 [5].
Calculating the ratio of expected mortality and official COVID-19 mortality can inform us on the quality of COVID-19 mortality tracking and the magnitude of indirect COVID-19 effects. The obtained ratio of 1.15 is similar to estimates obtained by Karlinsky and Kobak using the same methodology, but accounting only for the year 2020 [2]. Such a value points to an undercount of COVID-19 induced deaths—either due to the imperfect counting of direct COVID-19 deaths, or due to a large increase in mortality from other causes induced by the COVID-19 pandemic [23,24].
Since the baseline for calculation of excess mortality includes the effect of mortality increase in winter due to influenza, the obtained estimate of excess mortality should be interpreted as excess mortality above the one observed in the presence of other seasonal infectious disease (such as influenza). Consequently, the excess mortality estimates presented here should be taken as the lower limit of excess mortality solely due to COVID-19. This is corroborated by recent studies indicating a substantial decrease of influenza levels during the beginning of the pandemic in the European region [25]. The recently published report by the Croatian Institute of Public Health on mortality causes in 2021 confirms the above statement, as there was only a single death due to influenza recorded during that entire year [14]. Nevertheless, we think that it is sensible to keep the effect of seasonal influenza in the baseline, as this approach gives us a number representing an additional mortality burden countries had to deal with on top of seasonal fluctuations.
The results on the county level underscore the necessity of simultaneously using both the relative excess mortality and excess mortality rate metrics when interpreting the impact of the pandemic. Using just one of them might be enough if county rankings were always equal for both metrics. However, this was often not the case, and there were many counties where the relative rankings differed strongly (Figure 2).
Counties for which the relative excess mortality was low and excess mortality rate was high could indicate a high historical mortality baseline in those counties. Conversely, a high excess percentage coupled with a low excess mortality rate could suggest a low county mortality baseline. All such counties have a younger demographic structure, probably because they contain large urban centers (cities of Zagreb, Split and Dubrovnik) which attract younger people due to a higher chance of finding work and better infrastructure for young families. Furthermore, such large cities have easier access to hospital care and better quality of hospitals, and both of these factors have been shown to decrease excess mortality rates [26,27,28].
Very low excess mortality values for the SM county deserve a special mention, since in December 2020, the county was hit by an earthquake of magnitude 6.3 on the Richter scale. The destruction of many homes, and/or fear of another disaster and aftershocks triggered migrations to other parts of the country [29,30]. The official number of people who left their homes is 2861, although the actual number is probably higher [31]. This could have contributed to low excess mortality in the region, since it decreased the population size, making the projected baseline overly high.
Just like with national-level data, we calculated excess mortality stratified by age and sex across different counties. In terms of excess mortality differences across sexes, the results indicate a male bias in most counties, in agreement with the national trend (Figure 3).
With respect to age, we again see an approximate exponential dependance of excess mortality on age for both sexes in most counties (Figure 4 and Figure 5). A notable exception to the general rule is the LS county, which displays negative excess percentages and mortality rates in the oldest age category.
A possible factor driving this deviation from expected behavior could be the effect of population density, as higher density means a greater probability of disease transmission. However, this parameter is not easily quantifiable, as it cannot be assessed by simply dividing the population count by the county area, since some counties may spread over large surfaces, in spite of most people living in fairly dense areas [32]. Consequently, “microscopic” effective densities have to be assessed, which, in addition to the local population density, also take into account additional factors such as the rate of daily migrations and percentage of people who use the public transportation system [32,33,34,35].
The LS County spans nearly 10% of Croatia’s surface, but contains only 1% of the total population, yielding a population density of just 7.99 people per square kilometer. Additionally, most people there live in houses, as opposed to large residential buildings containing many apartments. Consequently, in addition to low general density, the LS county also has a low effective (local) density, making it easier to impose social distancing. This could explain the unexpected low excess mortality rate in the oldest age group, in spite of this county having the highest proportion of people in that age group (largest number of people at risk) across all counties (Table 2).
Although excess mortality is a great metric for assessment of the overall magnitude of a crisis event, studies calculating excess mortality are faced with several methodological limitations. Estimating excess mortality requires defining a baseline level of expected mortality. The baseline value can be highly influenced by the choice of statistical model and/or the duration of the historical reference period, impacting the final excess mortality value [15]. The estimates are also influenced by how detailed the mortality data are (weekly, monthly or yearly mortality records). Although the total mortality counts across counties in our study were available on a monthly basis, the mortality data for different age and sex subgroups across counties were available only on a yearly basis, affecting the excess estimates. Finally, without detailed information on specific causes of death, excess mortality studies cannot capture the full complexity of how different factors contribute to mortality during a crisis. Further studies involving such information could elucidate the underlying causes for the difference between the recorded and estimated mortality during the COVID-19 pandemic.

5. Conclusions

Excess mortality emerged as a crucial indicator during the COVID-19 pandemic. Assuming the incidence of other diseases remained stable over the period of interest, it can give us an insight into the number of deaths directly or indirectly related to COVID-19. We calculated excess mortality on the national and subnational (county) levels across different age–sex subgroups. An excess mortality prediction of 14,963 was obtained for the whole country, in agreement with estimates from other studies ranging from 12,205 to 22,900.
The relationships of excess mortality with age and sex on the national level were in line with previous research, showing a greater mortality rate for men and an exponential increase in excess mortality rate with increasing age. The largest increase in mortality (largest relative excess mortality value) was observed in the 65–74 years age group.
On the county level, we showed that both relative excess mortality and excess mortality rate have to be used simultaneously in order to fully grasp the local effect of the pandemic. Different county rankings for each of these metrics were explained through differences in baseline mortality rates due to the distinct demographic structure of each county.
Although the age–sex trends across counties mostly followed the national trends, there were also exceptions to the general rule. These exceptions were discussed in terms of between-county differences in demographic structure, population density and special events that took place during the pandemic.
In summary, this study offers a comprehensive analysis of the effect of the COVID-19 pandemic on excess mortality in Croatia during the 2020–2021 period. Unlike previous studies involving Croatian data, it takes into account the age and sex differences on both the national and county level. Although the impact of the pandemic can be assessed using only absolute or relative measures of excess mortality, our results underline the necessity of using both metrics in order to be able to make correct and informed decisions.

Author Contributions

Conceptualization: M.Š., Z.B., M.E. and N.B.; methodology, Z.B., M.E. and N.B.; software, Z.B.; validation: M.Š., Z.B., M.E. and N.B.; formal analysis: Z.B.; investigation: M.Š., Z.B. and N.B.; resources: M.E.; data curation: Z.B. and M.E.; writing—original draft preparation: M.Š., Z.B. and N.B.; writing—review and editing: M.Š., Z.B., M.E. and N.B.; visualization, Z.B.; supervision, M.E. and N.B. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded privately by the authors.

Institutional Review Board Statement

Ethical review and approval were waived for this study due to using publicly available mortality statistics data.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data and code are available upon reasonable request from the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Full names and abbreviations of Croatian counties.
Table A1. Full names and abbreviations of Croatian counties.
CountyAbbreviation
Bjelovar-Bilogora CountyBB
Brod-Posavina CountyBP
Dubrovnik-Neretva CountyDN
Istria CountyIS
Karlovac CountyKA
Koprivnica-Križevci CountyKK
Krapina-Zagorje CountyKZ
Lika-Senj CountyLS
Međimurje CountyME
Osijek-Baranja CountyOB
Požega-Slavonia CountyPS
Primorje-Gorski Kotar CountyPG
Sisak-Moslavina CountySM
Split-Dalmatia CountySD
Varaždin CountyVA
Virovitica-Podravina CountyVP
Vukovar-Srijem CountyVS
Zadar CountyZD
ZagrebGZ
Zagreb CountyZG
Šibenik-Knin CountySK

References

  1. Beaney, T.; Clarke, J.M.; Jain, V.; Golestaneh, A.K.; Lyons, G.; Salman, D.; Majeed, A. Excess Mortality: The Gold Standard in Measuring the Impact of COVID-19 Worldwide? J. R. Soc. Med. 2020, 113, 329–334. [Google Scholar] [CrossRef]
  2. Karlinsky, A.; Kobak, D. Tracking Excess Mortality across Countries during the COVID-19 Pandemic with the World Mortality Dataset. eLife 2021, 10, e69336. [Google Scholar] [CrossRef] [PubMed]
  3. COVID-19 Excess Mortality Collaborators. Estimating Excess Mortality due to the COVID-19 Pandemic: A Systematic Analysis of COVID-19-Related Mortality, 2020–2021. Lancet 2022, 399, 1513–1536. [Google Scholar] [CrossRef] [PubMed]
  4. United Nations. World Population Prospects 2019: Department of Economic and Social Affairs. World Population Prospects 2019. Available online: https://population.un.org/wpp/ (accessed on 15 February 2024).
  5. Cipin, I.; Mustac, D.; Medjimurec, P. Impact of COVID-19 on Mortality in Croatia. Stanovnistvo 2021, 59, 1–16. [Google Scholar] [CrossRef]
  6. Starke, K.R.; Reissig, D.; Petereit-Haack, G.; Schmauder, S.; Nienhaus, A.; Seidler, A. The Isolated Effect of Age on the Risk of COVID-19 Severe Outcomes: A Systematic Review with Meta-Analysis. BMJ Glob. Health 2021, 6, e006434. [Google Scholar] [CrossRef]
  7. Levin, A.T.; Hanage, W.P.; Owusu-Boaitey, N.; Cochran, K.B.; Walsh, S.P.; Meyerowitz-Katz, G. Assessing the Age Specificity of Infection Fatality Rates for COVID-19: Systematic Review, Meta-Analysis, and Public Policy Implications. Eur. J. Epidemiol. 2020, 35, 1123–1138. [Google Scholar] [CrossRef]
  8. Nielsen, J.; Nørgaard, S.K.; Lanzieri, G.; Vestergaard, L.S.; Moelbak, K. Sex-Differences in COVID-19 Associated Excess Mortality Is Not Exceptional for the COVID-19 Pandemic. Sci. Rep. 2021, 11, 20815. [Google Scholar] [CrossRef]
  9. Scully, E.P.; Haverfield, J.; Ursin, R.L.; Tannenbaum, C.; Klein, S.L. Considering How Biological Sex Impacts Immune Responses and COVID-19 Outcomes. Nat. Rev. Immunol. 2020, 20, 442–447. [Google Scholar] [CrossRef]
  10. Peckham, H.; de Gruijter, N.M.; Raine, C.; Radziszewska, A.; Ciurtin, C.; Wedderburn, L.R.; Rosser, E.C.; Webb, K.; Deakin, C.T. Male Sex Identified by Global COVID-19 Meta-Analysis as a Risk Factor for Death and ITU Admission. Nat. Commun. 2020, 11, 6317. [Google Scholar] [CrossRef]
  11. Croatian Institute of Public Health. Croatian Health Statistics Yearbooks 2017–2021. Available online: https://www.hzjz.hr/cat/hrvatski-zdravstveno-statisticki-ljetopis/ (accessed on 15 February 2024).
  12. Croatian Bureau of Statistics, Population Census 2021. Available online: https://dzs.gov.hr/vijesti/objavljeni-konacni-rezultati-popisa-2021/1270 (accessed on 15 February 2024).
  13. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computin: Vienna, Austria, 2022; Available online: https://www.R-Project.org/ (accessed on 15 February 2024).
  14. Croatian Institute of Public Health. Mortality Statistics. Available online: https://www.hzjz.hr/tag/mortalitet-statistika/ (accessed on 15 February 2024).
  15. Levitt, M.; Zonta, F.; Ioannidis, J.P.A. Comparison of Pandemic Excess Mortality in 2020–2021 across Different Empirical Calculations. Environ. Res. 2022, 213, 113754. [Google Scholar] [CrossRef] [PubMed]
  16. The Economist. Tracking COVID-19 Excess Deaths. Available online: https://www.economist.com/graphic-detail/coronavirus-excess-deaths-tracker (accessed on 15 February 2024).
  17. Estimated 2020–2021 Excess Death Count Based on Karlinsky and Kobak Methods. Available online: https://ourworldindata.org/excess-mortality-covid (accessed on 15 February 2024).
  18. Estimated 2020–2021 Excess Death Count from the Economist. Available online: https://ourworldindata.org/excess-mortality-covid (accessed on 15 February 2024).
  19. Barbieri, M.; Wilmoth, J.R.; Shkolnikov, V.M.; Glei, D.; Jasilionis, D.; Jdanov, D.; Boe, C.; Riffe, T.; Grigoriev, P.; Winant, C. Data Resource Profile: The Human Mortality Database (HMD). Int. J. Epidemiol. 2015, 44, 1549–1556. [Google Scholar] [CrossRef]
  20. Human Mortality Database: University of California Berkeley (USA) and Max Planck Institute for Demographic Research (Germany). Available online: https://www.mortality.org (accessed on 15 February 2024).
  21. López-Gigosos, R.; Vegas-Romero, M.; Mariscal, A.; Mariscal-López, E.; Fang, J.; Gutiérrez-Bedmar, M. Excess Mortality in 22 European Countries in 2020 and 2021: Relationship with Socioeconomic Indicators. Perspect. Public Health 2023, 17579139231180800. [Google Scholar] [CrossRef] [PubMed]
  22. Ioannidis, J.P.A.; Zonta, F.; Levitt, M. Variability in Excess Deaths across Countries with Different Vulnerability during 2020–2023. Proc. Natl. Acad. Sci. USA 2023, 120, e2309557120. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, W.-E.; Woo Park, S.; Weinberger, D.M.; Olson, D.; Simonsen, L.; Grenfell, B.T.; Viboud, C. Direct and Indirect Mortality Impacts of the COVID-19 Pandemic in the United States, March 1, 2020 to January 1, 2022. eLife 2023, 12, e77562. [Google Scholar] [CrossRef]
  24. Riou, J.; Hauser, A.; Fesser, A.; Althaus, C.L.; Egger, M.; Konstantinoudis, G. Direct and Indirect Effects of the COVID-19 Pandemic on Mortality in Switzerland. Nat. Commun. 2023, 14, 90. [Google Scholar] [CrossRef] [PubMed]
  25. Adlhoch, C.; Mook, P.; Lamb, F.; Ferland, L.; Melidou, A.; Amato-Gauci, A.J.; Pebody, R. Very Little Influenza in the WHO European Region during the 2020/21 Season, Weeks 40 2020 to 8 2021. Eurosurveillance 2021, 26, 2100221. [Google Scholar] [CrossRef] [PubMed]
  26. Janke, A.T.; Mei, H.; Rothenberg, C.; Becher, R.D.; Lin, Z.; Venkatesh, A.K. Analysis of Hospital Resource Availability and COVID-19 Mortality Across the United States. J. Hosp. Med. 2021, 16, 211–214. [Google Scholar] [CrossRef]
  27. Tchicaya, A.; Lorentz, N.; Leduc, K.; de Lanchy, G. COVID-19 Mortality with Regard to Healthcare Services Availability, Health Risks, and Socio-Spatial Factors at Department Level in France: A Spatial Cross-Sectional Analysis. PLoS ONE 2021, 16, e0256857. [Google Scholar] [CrossRef]
  28. Kunz, J.S.; Propper, C. JUE Insight: Is Hospital Quality Predictive of Pandemic Deaths? Evidence from US Counties. J. Urban Econ. 2023, 133, 103472. [Google Scholar] [CrossRef]
  29. Čivljak, R.; Markotić, A.; Capak, K. Earthquake in the Time of COVID-19: The Story from Croatia (CroVID-20). J. Glob. Health 2020, 10, 3–6. [Google Scholar] [CrossRef]
  30. Maršanić, V.B.; Dobrović, N.; Tadić, M.F.; Flander, G.B. 2020 Double Crisis in Croatia: Earthquakes in the Time of COVID-19. Eur. Child Adolesc. Psychiatry 2021, 30, 1309–1313. [Google Scholar] [CrossRef] [PubMed]
  31. Civil Protection Headquarter, Earthquake Information. Available online: https://potresinfo.gov.hr/sisacko-moslavacka-zupan (accessed on 15 February 2024).
  32. Desmet, K.; Wacziarg, R. JUE Insight: Understanding Spatial Variation in COVID-19 across the United States. J. Urban Econ. 2022, 127, 103332. [Google Scholar] [CrossRef] [PubMed]
  33. Pascoal, R.; Rocha, H. Population Density Impact on COVID-19 Mortality Rate: A Multifractal Analysis Using French Data. Phys. A Stat. Mech. Appl. 2022, 593, 126979. [Google Scholar] [CrossRef] [PubMed]
  34. Parkes, B.; Stafoggia, M.; Fecht, D.; Davies, B.; Bonander, C.; De’ Donato, F.; Michelozzi, P.; Piel, F.B.; Strömberg, U.; Blangiardo, M. Community Factors and Excess Mortality in the COVID-19 Pandemic in England, Italy and Sweden. Eur. J. Public Health 2023, 33, 695–703. [Google Scholar] [CrossRef]
  35. Harris, J.E.; Boarnet, M.; Boynton, K.; Brodsky, G.; Cardone, P.; Cohen-gould, L.; Cooley, P.; Cragg, M.; Diamond, P.; Everett, D.; et al. The Subways Seeded the Massive Coronavirus Epidemic in New York City; National Bureau of Economic Research: Cambridge, MA, USA, 2020. [Google Scholar] [CrossRef]
Figure 1. Country-level relative excess mortalities and excess mortality rates per 1000 people for different sex–age group combinations.
Figure 1. Country-level relative excess mortalities and excess mortality rates per 1000 people for different sex–age group combinations.
Idr 16 00011 g001
Figure 2. (A) County-level relative excess mortalities (left) and excess mortality rates per 1000 people (right) with lines in between connecting the positions representing the same counties for easier comparison of different metrics’ values. (B) Maps of relative excess mortalities and excess mortality rates across different counties.
Figure 2. (A) County-level relative excess mortalities (left) and excess mortality rates per 1000 people (right) with lines in between connecting the positions representing the same counties for easier comparison of different metrics’ values. (B) Maps of relative excess mortalities and excess mortality rates across different counties.
Idr 16 00011 g002
Figure 3. Comparison of specific excess mortality rates per 1000 people across counties and sexes for each of the age categories.
Figure 3. Comparison of specific excess mortality rates per 1000 people across counties and sexes for each of the age categories.
Idr 16 00011 g003
Figure 4. Comparison of two excess mortality metrics across counties for men in different age categories. Relative excess mortality is given in violet, and excess mortality rate per 1000 people in blue. Correspondingly colored dashed lines denote the zero value for each metric to facilitate comparison.
Figure 4. Comparison of two excess mortality metrics across counties for men in different age categories. Relative excess mortality is given in violet, and excess mortality rate per 1000 people in blue. Correspondingly colored dashed lines denote the zero value for each metric to facilitate comparison.
Idr 16 00011 g004
Figure 5. Comparison of two excess mortality metrics across counties for women in different age categories. Relative excess mortality is given in violet, and excess mortality rate per 1000 people in blue. Correspondingly colored dashed lines denote the zero value for each metric to facilitate comparison.
Figure 5. Comparison of two excess mortality metrics across counties for women in different age categories. Relative excess mortality is given in violet, and excess mortality rate per 1000 people in blue. Correspondingly colored dashed lines denote the zero value for each metric to facilitate comparison.
Idr 16 00011 g005
Table 1. Country-level excess mortalities for the entire population and different age–sex subgroups.
Table 1. Country-level excess mortalities for the entire population and different age–sex subgroups.
AgeSexPopulationExcess Mortalityz-Score 1
0–64M1,503,6118781.98
0–64F1,498,9835252.15
65–74M230,234290023.7
65–74F277,48419497.27
75–84M103,73726474.82
75–84F165,83722432.23
85+M27,54714775.36
85+F64,40023522.28
All agesBoth sexes3,871,83314,9635.09
1 Excess mortality prediction is deemed statistically significant if the z-score is larger than 2.
Table 2. Total county populations with percentages of population in each age category.
Table 2. Total county populations with percentages of population in each age category.
CountyPopulation0–6465–7475–8485+
BB101,87976.813.96.92.4
BP130,26777.413.17.22.3
DN115,56477.513.06.92.6
GZ767,13179.311.66.82.3
IS195,23775.814.17.22.8
KA112,19575.513.97.63.1
KK101,22178.013.06.92.1
KZ120,70279.212.16.42.2
LS42,74873.713.98.73.6
ME105,25079.412.26.32.1
OB258,02678.013.36.72.0
PG265,41974.215.47.62.8
PS64,08477.113.37.32.2
SD423,40778.212.86.52.4
SK96,38172.615.58.53.3
SM139,60375.314.77.62.5
VA159,48779.212.16.62.1
VP70,36878.013.26.62.1
VS143,11377.213.67.12.2
ZD159,76676.014.07.52.5
ZG299,98578.712.96.42.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Šošić, M.; Boban, Z.; Erceg, M.; Boban, N. Excess Mortality Stratified by Age and Sex for Croatia and Croatian Counties during the 2020–2021 COVID-19 Pandemic. Infect. Dis. Rep. 2024, 16, 142-153. https://0-doi-org.brum.beds.ac.uk/10.3390/idr16020011

AMA Style

Šošić M, Boban Z, Erceg M, Boban N. Excess Mortality Stratified by Age and Sex for Croatia and Croatian Counties during the 2020–2021 COVID-19 Pandemic. Infectious Disease Reports. 2024; 16(2):142-153. https://0-doi-org.brum.beds.ac.uk/10.3390/idr16020011

Chicago/Turabian Style

Šošić, Mara, Zvonimir Boban, Marijan Erceg, and Nataša Boban. 2024. "Excess Mortality Stratified by Age and Sex for Croatia and Croatian Counties during the 2020–2021 COVID-19 Pandemic" Infectious Disease Reports 16, no. 2: 142-153. https://0-doi-org.brum.beds.ac.uk/10.3390/idr16020011

Article Metrics

Back to TopTop