Abstract
The COVID-19 pandemic has left other pathologies commonly present in the population in a secondary context. Therefore, it is necessary to study the evolution of these diseases in the presence or absence of COVID-19.
Objectives The present study had the following objectives: 1. To evaluate the relationship between the COVID-19 epidemic and the possible decrease in deaths from respiratory diseases in Chile. 2 Study the relationship between meteorological variables and severity of COVID-19 with respect to deaths from respiratory diseases from January 2018 to February 2021.
Methods The variable number of deaths due to respiratory diseases in Chile was analyzed considering the monthly records of meteorological variables (temperature, precipitation and humidity) in each Region of Chile and severity of COVID-19, to evaluate the mortality trend before and after the pandemic. For this, different Non-Observable Heterogeneity Models for Panel Data were used
Results Our findings show that among the variables that affect the mortality rate from respiratory diseases, there are the number of deaths from COVID-19 that has a negative effect, the number of patients with COVID-19 in intensive care unit (ICU) that has a positive effect and the minimum temperature with a negative effect. These results are supported by the application of the panel regression with one-way random-effects models.
Conclusion This study revealed that there is an unexpected reduction in deaths from respiratory diseases in Chile in the post-pandemic period. Therefore, it can be concluded that this variable decreased with the appearance of the COVID-19 pandemic.
Introduction
In December 2019, the People’s Republic of China reported a cluster of pneumonia cases of unknown etiology, subsequently identified as a new coronavirus by the Chinese Center for Disease Control and Prevention. Later time, in February 2020, the World Health Organization (WHO) named the disease COVID-19, short for “coronavirus disease 2019” and, in March of the same year, declared an international public health emergency for its rapid community, regional and global spread. From the confirmation of the first COVID-19 cases to March 2021, 116.736.437 confirmed cases were reported, including 2.593.285 deaths worldwide, of which 45 % of confirmed cases and 48 % of deaths were contributed by the Region of the Americas. Specifically, South America contributed with a higher proportion of deaths (85 %) of the total, surpassing North America (14.5 %) [1].
In March 2021, the Ministry of Health of Chile (MINSAL for the Spanish acronym for Ministry of Health) published in its epidemiological report [2], that up to that date a total of 1.018.677 cases (confirmed and probable) and a total of 28.756 deaths from COVID-19 had been registered. Another report [3] also prepared by the MINSAL, describes the analysis of the causes of deaths from COVID-19, where Chile applied rules for coding mortality and simulating the profile of deaths from COVID-19, not considering the existence of the virus as a cause death in medical certification. This evidenced that respiratory causes would have the first place in deaths in the country, observing that deaths from Influenza and Pneumonia represent 47 %, and deaths from Other Diseases of the Respiratory System represent 5 %, displacing even circulatory diseases and tumors, which in previous years occupied the first places in the ranking of deaths in the country.
In Chile, in 2018 (pre-pandemic) respiratory diseases ranked third in the national ranking of cause of death, only surpassed by diseases of the circulatory system and cancer; observing a total of 12.228 deaths from this cause [4]. However, as indicated in the previous paragraph, this situation changed with the arrival of COVID-19 in the country.
A recent study published by [5] reveals that 545 million people currently live with respiratory diseases (Rd hereafter), which represents 7.4 % of the world population and that 3.9 million people die each year from Rd, which constitutes a great health burden, due to its severity and mortality figures worldwide. The associated healthcare costs are an increasing burden on the economies of all countries, and if the loss of family or caregiver productivity of people with respiratory diseases is considered, the cost to society is much higher [6]. Likewise, at the national level in people over 65 years of age, Rd are the second cause of hospital discharge [7]. The previous discussion reflects the need to carry out a study on Rd in the pre and post pandemic period so that health entities are aware of a quantitative measure of how COVID-19 has affected the statistics on deaths from respiratory diseases, in such a way that this information can contribute to the effective distribution of health resources.
There are several studies that have addressed the health impact relating mortality from respiratory diseases and mortality from COVID-19. In this context, a Brazilian study made comparisons between deaths from respiratory diseases and deaths from COVID-19 with a time series statistical methodology for monthly registry data, in order to describe the peak of mortality during the summer seasons. This study revealed a latitudinal gradient in the country, with 27 peak deaths occurring in April in most of the northern and northeastern states, and gradually later in the southernmost states, with peak mortality in June and July (season of cold) in the Southeast and South regions, respectively [8].
On the other hand, a Spanish study used four regression models to explore the associations of the average daily temperature and air quality (PM 2.5) with the new daily cases of COVID-19 in the four main regions of Spain (Castilla y León, Castilla-La Mancha, Catalonia and Madrid) [9]. Among these models are, Panel Regression Models, Quantile Regression Model, Clustered OLS and Fixed Effects Regression Model known as Individual Heterogeneity Models. Recently, a study used Individual Heterogeneity Models to examine the effect between COVID-19 using the variables confirmed cases and deaths, which reflect the severity of the pandemic, and the meteorological factors and the air pollutant (PM 2.5), in six countries in South Asia. The authors have concluded that high temperature and humidity increase the transmission of COVID-19 infections and that it can also be applied to regions with higher transmission rates, where the minimum temperature is mostly above 21◦ C and the Humidity hovers around 80 % for months. Additionally, the air pollutant (PM 2.5) exhibits significant negative and positive effects on the number of confirmed COVID-19 cases [10]. Other investigations have also used Individual Heterogeneity Models as statistical methodology to examine the relationship between meteorological variables and environmental pollution with respect to confirmed cases of COVID-19 and mortality, [11],[12] and their references.
Another study evaluated the co-infection of COVID-19 in the presence of other respiratory viruses, analyzing the statistically significant correlations of the variables. In particular, [13] evaluated the presence of influenza A/B virus, human metapneumovirus, bocavirus, adenovirus, respiratory syncytial virus, and parainfluenza virus in 105 patients who died from COVID-19; finding co-infection with influenza virus in 22.3 %, respiratory syncytial virus and bocavirus in 9.7 %, parainfluenza virus in 3.9 %, human metapneumovirus in 2.9 % and finally adenovirus in 1.9 % of COVID-19 positive dead cases. In the same line [14] suggested higher rates of co-infection between SARS-CoV-2 and other respiratory pathogens. There are also studies that have explored the possibility that COVID-19 has affected the seasonal trend of other viruses, for example, the circulation of influenza has been decreased see [15]. However, it circulates together with COVID-19, being able to generate a greater impact on the evolution and management of respiratory diseases [16]. An important aspect that has aroused great interest among different scientific groups is the relationship between influenza and the COVID-19 outbreak. In this context, a study [17] analyzed whether influenza vaccination minimizes the spread of COVID-19 in the Italian population aged 65 to 48 years or older; concluding that the influenza vaccination coverage rate in people 65 years of age and older was associated with reduced spread and less severe clinical expression of COVID-19. The reader can deepen this topic by reviewing the studies carried out in Europe and the USA. [18],[19] and [20].
Among the various topics about the impact of COVID-19 and its relationship with respiratory diseases, there are studies that compared the rates of transmissibility, hospitalization and mortality of COVID-19 with respect to other pathogens that cause respiratory diseases [21]. The review of the research mentioned in the previous paragraphs highlights the effects that the COVID-19 pandemic has had on the health of the population. So much so, that other pathologies that cause great damage or deterioration of life and even that can cause death have been neglected. It is in this context, it seems necessary to study the effect of the COVID-19 pandemic outbreak that has had on deaths from respiratory diseases in Chile, considering the regional effect and climatic variables that predominate throughout the country.
In this study we focus on respiratory diseases, classified in J00-J99 according to the International Classification of Diseases (ICD 10); and that constitute those pathologies of the respiratory tract and other structures of the lung, associated with the main causes of morbidity and mortality worldwide [22]. In particular, this research proposes to quantify the relationship between deaths from respiratory diseases with the pandemic generated by COVID-19 and with meteorological variables in Chile. To do this, we generate a data panel with monthly records from different regions of Chile from January 2018 to February 2021 and Individual Heterogeneity Models are proposed, such as; Fixed and Random Effects Regression Models. To the best of our knowledge, this of the few studies to use panel data regression analysis to determine the effect of deaths from respiratory diseases and COVID-19 outbreak. The objective of this work is to provide an empirical study to determine the variables that have affected mortality from respiratory diseases in Chile in the period of the COVID-19 pandemic. In this way, it is intended to contribute useful statistical information for government authorities on the evolution of deaths from respiratory diseases in Chile with respect to the COVID-19 pandemic and meteorological variables. Another aspect that we must highlight from this study is that once COVID-19 is controlled in the population, chronic diseases that have been neglected in this period will reappear with force, and therefore, have information on the effect of COVID-19 on these chronic diseases is crucial to effectively allocate both available human and financial resources. In this context, this study has focused on deaths from respiratory diseases, without prejudice that it can be extended to other chronic diseases which must be continuously monitored by the Health authorities.
The rest of the article is organized as follows: Section 1 introduces the data, defines the variables, along with the empirical models; Section 2 reports the empirical results. Section 3 presents the strengths and weaknesses of this study through a brief discussion, and the last section concludes the study.
1. Materials and methods
Data Sources and the Definition of Variables
To carry out this study, the monthly records of the variables (Table 1) were considered for 14 regions of Chile, from January 2018 to February 2021. These variables include the data on the number of deaths from respiratory diseases obtained from the Ministry of Health https://deis.minsal.cl/, whereas the data on the variables regarding the severity of COVId-19 were obtained from the Ministry of Science and Technology https://www.minciencia.gob.cl/covid19. Finally, the data of the temperature and precipitation variables were obtained from the updated report of the Chilean Meteorological Directorate DGAC, which are available on the website https://climatologia.meteochile.gob.cl. With this information, a balanced data panel was generated with a total of 14 regions and 38 temporal observations. Figure (1) shows a map of Chile with the regions where data recording was possible (gray color). In the regions highlighted in white, it was not possible to build the observation panel since there are no meteorological stations in these geographical areas that allow obtaining information on temperature and precipitation. Table (1) describes the variables used in this study. In particular, the dependent variable is the number of deaths from respiratory diseases (DRd). In addition, to investigate the impact of COVID-19 on the number of DRd, we consider three variables, namely: the monthly number of confirmed deaths from COVID-19 (CD), the incidence rate (IC), and the monthly number of patients in the intensive care unit diagnosed with COVID-19 (ICU).
Chile’ regiones. Dataset
Varibales under study
Empirical Models
To quantify the impact that the COVID-19 pandemic has had on the number of DRd in Chile, we used a panel data regression. In this regard, we can point out that the panel data models take into account time and regions simultaneously, while other models have the limitation of only expressing these heterogeneities across units or over time. Furthermore, panel data models are better at capturing the heterogeneity involved in both cross-sectional units and time dimensions, to reduce estimation bias and multicollinearity, and are better suited for studying the dynamics of change and complex behavioral models [23] and [24]. In this study, we consider two types of panel data models. The first, denoted by Model 1, is a two-way fixed or random effects model, that is, this model captures the individual and temporal effect of the data. While the second model, called Model 2, is a one-way random or fixed effects model, which only considers the temporal effect of the data.
▪ Two-way fixed-effects model
First, we consider a two-way fixed-effects model with all the explanatory variables as follows:
Model 1
where DRd denotes the number of deaths from respiratory diseases, CD is the number of monthly deaths from COVID-19, the log (CD) is with respect to the observations obtained at the beginning of the pandemic, that is, observed from March 2020, CI is the incidence rate, which is the number of monthly confirmed cases divided by the gross permanent resident population (number confirmed/population (10.000 per unit)), maxT minT denotes temperature (measured in Celsius), P denotes precipitation (measured in hPa), ICU is the number of COVID-19 patients in intensive care, i denotes the region, t denotes the period under study, the δi and ηt are the region and time fixed effects, respectively, and εit denotes the error term. In the case of considering a two-way random-effects model, the individual and temporal effects satisfy the following,
and
For the random effects models, there are three estimation methodologies for the parameters, which are available in the free R statistical software [25], through the plmtest command of the plm package. This command provides the following methods, “swar” [26] (default), “amemiya” [27], and “walhus” [28].
▪ One-way fixed-effects model
The data come from 14 different regions, each with peculiar characteristics, which could suggest that the observations from the same region share some common characteristics, for example, climatic conditions, economic development, among other characteristics. Therefore, such characteristics could be related to the regressors, so it is advisable to incorporate the 14 regions as factors in the model and thus control by region. For these reasons, a model is proposed where the errors uit consider only a time effect to control the contemporary correlation in the cross section as follows
where sit = δiRiIi(t) with Ri represent the i th region and Ii(t) 0, 1 an indicator variable, for i = 1, …, 14, t = 1, …, T, defined by
2. Results
In order to know the behavior of the variable of DRd in all the regions registered to priori and posterior of the pandemic, we consider the prior and posterior mean-variance tests, where the data from January 2020 are posterior to the pandemic and other months as prior to the pandemic. Table 2 reports the results obtained from applying the mean-variance tests. We have differentiated the cases of non-homogeneity of variance, in which case we use the Welch’s test. From the Table 2 we can say that deaths from respiratory diseases in all the study regions are significantly lower than before the outbreak, except in the Aysén region, which has a non-significant difference.
Mean-variance tests. This table reports the results of the mean-variance test of number of deaths due to respiratory diseases prior and posterior to the pandemic. The symbols ***, **, and * represent significance levels of 1 %, 5 %, and 10 %, respectively.
To complement the previous analysis, we consider one-way model with the post and pre-pandemic dummy variable to study its effect on Chile’s deaths from respiratory diseases, that is:
where Pdummyt is a dummy variable, assuming a value of 1 for records after the COVID-19 pandemic and 0 for the months before the COVID-19 pandemic. This dummy variable can be used to directly compare the differences between the DRd before and after COVID-19. Among the models used to describe the relationship of equation (3) are the Pooling Model; which groups all the cross sections, without considering any form of individual effect and among the non-observable heterogeneity models are the Within Model, which is also called fixed effects models, and the random effects models defined above. Table (3) shows the estimated results of the aforementioned models. Note that in all models the Pdummyt covariate is statistically significant. In particular, Column 4 shows that the COVID-19 pandemic has significantly reduced the DRd in a 34.12 %. It should be taken into consideration that, in the previous analysis, and as indicated in Eq (3), we only the regional fixed effect was controlled in the regression. Now let’s go a step further to consider panel data regression as Eqs (1)-(2). Next, the models proposed in this study (1 and 2) will be studied to analyze which are the variables that affect the DRd in Chile.
The individual effect model.
2.1. Using two-way fixed-effects model: Model 1
In this case, the performance of Model 1 was evaluated to determine the degree of impact of the pandemic on deaths from respiratory diseases in the regions of Chile. For this, it was controlled for both the regional fixed effect and the temporal effect in the panel data. Table (4) reports the benchmark estimation results of Model 1. From column 5, it shows that the number of deaths from Covid-19 has a significant negative effect on the response variable. On the other hand, the number of beds in ICU used by Covid-19 patients has a significant positive effect. Finally, among the climatic variables, the minimum temperature has a significant negative effect on the number of deaths from respiratory diseases in Chile, while precipitation has a positive effect. For example, for each percentage unit of increase in deaths from Covid-19, the number of deaths from respiratory diseases decreases by 7.6 %, while for each unit of increase in the number of ICU beds due to Covid-19 the deaths from respiratory diseases increased by 0.055 %. Another important aspect is that if the minimum temperature decreases by one degree the number of deaths from respiratory diseases rises 4.942 % and, a unit increase in the precipitation would result in a go up of 0.089 % in the number deaths from respiratory diseases. It is important to note that the performance of Model 1 with respect to adjusted R2 is low, which suggests that this model would not be the most appropriate to adjust the data under study.
The twoways effect model.
Figure 2 display the fitted values for each regression model of the data panel, where the red line represents the actual data, while the colored lines represent the fitted values for Model 1. From these figures, we can conclude that the best model that captures the trend structure is a Swar Type Random Effect Model.
Model 1: Number of deaths from respiratory diseases (red lines) versus fitted values (colored lines).
2.2. Using one-way fixed-effects model: Model 2
Chile is one of the most centralized countries in Latin America and of the OCDE [30], that is, there is a high concentration of resources and development in the Metropolitan Region. Specifically, the capital of the country concentrates the best schools, hospitals, clinics, universities, cultural events and social development, among others. To eliminate the “centralization” effect, the regional effect was used as a covariate to study the impact of COVID-19 on the DRd based on the one-way fixed-effects model. Table (5) shows the estimated results with respect to Model 2. From Column 5, it is observed that COVID-19 has significantly reduced the number of deaths from respiratory diseases. For example, a unit percentage increase in deaths-Covid-19 would result in a 5.987 % drop in the number of deaths from respiratory diseases. On the other hand, a unit increase in the number of COVID-19 patients in intensive care leads to a 0.039 % increase in the number of deaths from respiratory diseases. Another important aspect is that if the minimum temperature increases by one degree, the number of deaths from respiratory diseases increases by 2,681 %. Unlike Model 1, the performance of Model 2 with respect to the adjusted R2 is high, this is supported by Fig 3 which shows the adjusted values of the panel regression model defined in Table (5). This suggests that Model 2 would be the most suitable for fitting the data in the study. A Hausman Type Test will determine whether individual effects should be treated as fixed or can be assumed incorrect with regressors, using a more efficient random effects specification. Under the hypothesis of no correlation between the regressors and the individual effects, this hypothesis, with a p-value of 18 %, is not rejected at the 5 % confidence level. Therefore, a random effect model with the Amemiya Method is more appropriate for this model, with an adjusted R2 = 0.95521. Finally, we can observe that the Regions of the extreme south (Region 11 and 12) and that of the extreme north (Region 15) show a decrease in the number of deaths due to respiratory diseases compared to the Tarapacá region (Region 1). While the rest of the regions present a significant increase in deaths from respiratory diseases when compared to the same region (Region 1). The foregoing is due to the fact that Regions 11, 12 and 15 have a lower number of population than the Tarapacá Region.
Time effect model.
Model 2: Number of deaths from respiratory diseases (red lines) versus fitted values (colored lines).
3. Discussion
An important strength of this study is the general design used to assess the possible relationship between the COVID-19 pandemic with respect to climatic conditions and mortality from respiratory diseases. In all statistical analyzes performed in this study, the magnitude and direction of the observed phenomena were found to be consistent. Importantly, it was also revealed that by incorporating geographic location in Model 2 as an explanatory variable, the magnitude and significance of the observed results increased considerably.
Regarding the behavior of the variable number of deaths from respiratory causes (DRd) in all the regions registered after the pandemic, it was significantly lower than the records prior to the COVID-19 pandemic. These findings confirm the hypothesis of a study conducted in Brazil where the authors estimate an average underreporting of 40.68 % (range 25.9 % -62.7 %) for deaths related to COVID-19 in this country [8].
In this study it was also observed (Fig 3) that in the extreme North and South regions of Chile, Model 2 has a low performance to explain the number of deaths from respiratory diseases. Among the possible causes of the above, it can be mentioned that; the population density of these regions is lower compared to other areas. Another possible cause is the large difference in PIB per capita between the extreme regions of the country and the Metropolitan region with the rest of the regions of the national territory [30], which is associated with compliance with the different confinement measures, and the preventive health care that the inhabitants of these regions have had.
The findings of the present study have limitations. First, this study had a longitudinal type design. Therefore, in these types of studies, variables recorded both in time and by locality are needed, which are difficult to obtain. In this context, there are sociodemographic variables whose monthly records were not available by region, and therefore, they were not present in this study to analyze the association between COVID-19 and cultural differences to avoid deaths from respiratory diseases. It is proposed to carry out a subsequent study, which not only considers the geographical differences that affect self-protection habits against respiratory diseases and COVID-19, but also cultural differences.
Second, our findings could be affected by the number of individuals vaccinated against influenza during the COVID-19 pandemic. Another possible limitation of this study could be the existence of errors in the identification/recoding of the true cause of death in the database and errors in the recording/calculation of the climatic conditions analyzed in the database. It is presumed that, if such inaccuracies occurred in the data, they would have occurred with a similar frequency among the different databases examined. Therefore, if such phenomena were present in the data examined, the statistical power of this study would have been reduced. Still, the present study’s results provide some insight into the association between the COVID-19 pandemic and the number of deaths from respiratory diseases.
4. Conclusion
Our findings show that among the variables that affect the mortality rate from respiratory diseases are; the number of deaths from COVID-19 which has a negative effect, the number of patients with COVID-19 in the ICU, which has a positive effect, and the minimum temperature with a negative effect. These results are supported by the application of the panel regression for one-way random-effect model. We can conclude that the number of deaths from respiratory diseases in Chile decreased with the appearance of the COVID-19 pandemic.
Data Availability
These variables 127 include the data on the number of deaths from respiratory diseases obtained from the 128 Ministry of Health https://deis.minsal.cl/, whereas the data on the variables 129 regarding the severity of COVId-19 were obtained from the Ministry of Science and 130 Technology https://www.minciencia.gob.cl/covid19. Finally, the data of the 131 temperature and precipitation variables were obtained from the updated report of the 132 Chilean Meteorological Directorate DGAC, which are available on the website 133 https://climatologia.meteochile.gob.cl.