Next Article in Journal
Opioid Misuse: A Review of the Main Issues, Challenges, and Strategies
Next Article in Special Issue
The Association between COVID-19 Pandemic and the Quality of Life of Medical Students in Silesian Voivodeship, Poland
Previous Article in Journal
Alcohol Recognition and Desire to Drink of Extended Alcohol Brand Logos
Previous Article in Special Issue
Impact of COVID-19 on Health-Related Quality of Life: A Longitudinal Study in a Spanish Clinical Sample
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

COVID-19, Non-Communicable Diseases, and Behavioral Factors in the Peruvian Population ≥ 15 Years: An Ecological Study during the First and Second Year of the COVID-19 Pandemic

Facultad de Ciencias de la Salud, Universidad Tecnológica del Peru, Lima 15046, Peru
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(18), 11757; https://doi.org/10.3390/ijerph191811757
Submission received: 16 August 2022 / Revised: 10 September 2022 / Accepted: 13 September 2022 / Published: 17 September 2022

Abstract

:
A range of health-related and behavioral risk factors are associated with COVID-19 incidence and mortality. In the present study, we assess the association between incidence, mortality, and case fatality rate due to COVID-19 and the prevalence of hypertension, obesity, overweight, tobacco and alcohol use in the Peruvian population aged ≥15 years during the first and second year of the COVID-19 pandemic. In this ecological study, we used the prevalence rates of hypertension, overweight, obesity, tobacco, and alcohol use obtained from the Encuesta Demográfica y de Salud Familiar (ENDES) 2020 and 2021. We estimated the crude incidence and mortality rates (per 100,000 habitants) and case fatality rate (%) of COVID-19 in 25 Peruvian regions using data from the Peruvian Ministry of Health that were accurate as of 31 December 2021. Spearman correlation and lineal regression analysis was applied to assess the correlations between the study variables as well as multivariable regression analysis adjusted by confounding factors affecting the incidence and mortality rate and case fatality rate of COVID-19. In 2020, adjusted by confounding factors, the prevalence rate of obesity (β = 0.582; p = 0.037) was found to be associated with the COVID-19 mortality rate (per 100,000 habitants). There was also an association between obesity and the COVID-19 case fatality rate (β = 0.993; p = 0.014). In 2021, the prevalence of obesity was also found to be associated with the COVID-19 mortality rate (β = 0.713; p = 0.028); however, adjusted by confounding factors, including COVID-19 vaccination coverage rates, no association was found between the obesity prevalence and the COVID-19 mortality rate (β = 0.031; p = 0.895). In summary, Peruvian regions with higher obesity prevalence rates had higher COVID-19 mortality and case fatality rates during the first year of the COVID-19 pandemic. However, adjusted by the COVID-19 vaccination coverage, no association between the obesity prevalence rate and the COVID-19 mortality rate was found during the second year of the COVID-19 pandemic.

1. Introduction

The incidence and mortality due to COVID-19 is associated with multiple risk factors, such as sociodemographic and health indexes, and genetic and immunology factors [1]. Among these, hypertension and obesity are significant risk factors for severe illness and death [2,3,4,5]. A recent study suggested that countries with a higher body mass index (BMI) or cholesterol level have a higher rate of COVID-19 incidence and mortality [6]. In addition, obesity is associated with the risk of hospitalization and mortality due to COVID-19 [3,5]. Behavioral risk factors have also been found to be associated with COVID-19 outcomes, above all in some regions where the prevalence of tobacco and alcohol use is high [7]. A systematic review suggested that smokers appear to be at a low risk of SARS-CoV-2 infection, while former smokers appear to be at increased risk of hospitalization, increased disease severity, and mortality from COVID-19 [8]. Another study examined whether social and behavioral risk factors were associated with mortality among U.S. veterans with COVID-19, but no association was found [9]. Despite these findings, we remain unclear about whether behavioral risk factors are associated with incidence and mortality due to COVID-19.
Peru has been deeply affected the COVID-19 pandemic. By the end of December 2021, more than 4 million confirmed COVID-19 cases and more than 215,000 deaths had been reported [10]. In addition, recent studies revealed that more than 100,000 excess all-cause deaths occurred in 2020, and there was a strong excess in geographical and temporal mortality patterns according to region [11]. This excess mortality can be explained by social factors, individual factors, the health care system, and other direct and indirect factors related to COVID-19 infection. Several Peruvian studies have shown an association between obesity and COVID-19 mortality, both in the overall population and in patients infected with SARS-CoV-2 [12,13]. To control the COVID-19 pandemic and mitigate mortality associated with its complications, Peru launched the COVID-19 vaccination on 9 February 2021, which initially included immunization of health workers and subsequently encompassed the overall population [14]. During the post-vaccination months, there has been a reduction in the COVID-19 mortality rate in the Peruvian population [15,16]. In this setting, we assessed the association between incidence, mortality, and case fatality rate due to COVID-19 and the prevalence rates of hypertension, obesity, overweight, tobacco, and alcohol use in the Peruvian population aged >15 years during the first and second year (since the launch of the COVID-19 vaccination) of the COVID-19 pandemic. These findings will help to determine the indirect impact of the SARS-CoV-2 vaccination on the COVID-19 mortality risk in regions where the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use are high.

2. Materials and Methods

2.1. Design and Area Study

We conducted an ecological study to analyze the association between the prevalence rates of hypertension, overweight, obesity, tobacco use, and alcohol use and the incidence rate, mortality rate, and case fatality rate due to COVID-19 in the Peruvian population between 1 March 2020 and 31 December 2021. Peru has an estimated 33 million inhabitants, according to the National Institute of Statistics and Informatics (INEI), divided into 24 geographic regions and one constitutional province, Callao [17]. The study was performed following the Strengthening the Reporting of Observational Studies in Epidemiology (Appendix A Table A1. Checklist STROBE Statement) reporting guidelines [18].

2.2. Risk Factors in Peru

For this study, we retrieved secondary data of each geographic region on the prevalence rates of hypertension, overweight, obesity, tobacco use, and alcohol use in Peruvian individuals aged ≥ 15 years from the Encuesta Demografica y de Salud Familiar (ENDES) 2020 and 2021 of the INEI [19,20]. ENDES 2020 included to 32,197 men and women aged ≥ 15 years from 25 Peruvian regions, from January to December 2020, and ENDES 2021 includes to 32,124 men and women aged ≥15 years [19,20]. Hypertension was defined as an elevation above the normal values of systolic blood pressure over diastolic blood pressure of 140/90 mm Hg [19]. Overweight and obesity were defined as the disproportionate deposition of fat in the body that can damage health. Tobacco and alcohol use were defined as risk factors for developing various chronic diseases that alter cardiovascular function, liver function, and renal function, among others [19].

2.3. COVID-19 Data Collection

COVID-19 case data were obtained from the freely accessible National Open Data Platform website, Peru (COVID-19 Open Data, https://www.datosabiertos.gob.pe/group/datos-abiertos-de-covid-19 (accessed on 20 January 2022)), accurate as of 31 December 2021. A COVID-19 confirmed case was defined as a case with a positive result for SARS-CoV-2 from an RT-PCR assay or antigen test or immunochromatographic test. In the National Open Data Platform, COVID-19 Open Data (“COVID-19 positive cases” and “COVID-19 deaths”), and data on cases of and deaths from COVID-19 include information on sex, age, region, and date of notification (cases or deaths). Other information on cases and deaths from COVID-19 is not available in the National Open Data Platform, COVID-19 Open Data. COVID-19 deaths were obtained from the Sistema Informatico Nacional de Defunciones (SINADEF) from the Peruvian Ministry of Health (MINSA), from March 2020 to December 2021 [21]. The SINADEF database records all deaths that occur in Peru. We included the deaths by COVID-19 as the underlying cause of death.

2.4. Confounding Factors

Confounding factors included the mean age in the region (cases or deaths), mean income and gender balance (cases or deaths), and the number of intensive care unit (ICU) beds. These confounding factors were used in a previous observational study [12]. We also used COVID-19 vaccination coverage rates in individuals ≥ 18 years in 2021 as a confounding factor. Data on the mean age in the region, mean income, and gender balance were retrieved obtained from the INEI [17]. The number of ICU beds was obtained by way of App. F500.2 from the Superintendencia Nacional de Salud, Perú (SUSALUD) [22]. COVID-19 vaccination coverage rates in 2021 were obtained from the Single National Health Information Repository (REUNIS), MINSA, Peru [23].

2.5. Statistical Analysis

Crude incidence and mortality rates (per 100,000 habitants) among individuals ≥ 15 years old were calculated by dividing the number of COVID-19 cases and deaths per region in 2020 and 2021 by the population of each region. We also calculated the COVID-19 case fatality rate among those aged ≥15 years in each region. Estimated populations for calculating incidence and mortality rates were obtained from the INEI, Peru. Spearman’s test correlation and linear regression models were used to estimate the associations between the incidence, mortality, and case fatality rates due to COVID-19 and the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use. Non-parametric approaches were applied for non-normally distributed data. Additionally, multivariable regression models were applied to assess the association between COVID-19 incidence, mortality, and case fatality rates and the prevalence rates of hypertension, overweight, obesity, tobacco and alcohol use. Confounding factors to multivariable regression models in 2020 and 2021 included the mean age in the region (years), mean monthly income (PEN), gender balance in deaths, and number of ICU beds in 2020 and 2021. We also used COVID-19 vaccination coverage rates in 2021 as a confounding factor. p-values < 0.05 were considered significant. Scatter plots were provided to indicate the correlation between the study variables. Statistical analyses were done using StataSE 17.0 Software.

2.6. Ethics

As this study used publicly available cases and death data, and there was no direct patient involvement.

3. Results

A total of 2,206,897 COVID-19 cases from 6 March 2020 to 31 December 2021 in individuals aged ≥15 years were included, of whom 970,590 were in the incidence analysis for 2020, and 1,236,307 for 2021 (Figure 1A). During this period, 198,447 deaths from COVID-19 in individuals aged ≥15 years were also included, of whom 94,225 were included in the mortality analysis for 2020, and 108,351 were included in the mortality analysis for 2021 (Figure 1B).
Table 1 describes the prevalence, incidence, mortality, and fatality rates in the Peruvian population stratified by region in 2020. The hypertension prevalence rate was higher in Callao, and the prevalence rates of overweight and tobacco use were higher in Madre de Dios. The prevalence rates of obesity and alcohol use were higher in Moquegua and Piura, respectively. The COVID-19 incidence rate was higher in Moquegua, while the mortality rates were highest in Callao, Moquegua, Ica, and Lambayeque. The fatality case rates were highest in Piura, Lambayeque, La Libertad, Ica, and Callao regions (Table 1).
In 2021, the prevalence rates of hypertension and obesity were higher in Tacna region, and the prevalence rate of overweight was higher in Pasco region. The prevalence rates of tobacco and alcohol use were higher in Madre de Dios and Moquegua regions, respectively (Table 2). The COVID-19 incidence rates were highest in Moquegua, Callao and Lima regions. The mortality rates were highest in Ica, Callao, and Lima regions; and the fatality case rates were highest in Ica, Ucayali, and Lambayeque regions (Table 2).

3.1. Correlation between the Study Variables and COVID-19 Measures in 2020

There was an association between the prevalence rates of obesity (r = 0.60; p = 0.001), tobacco use (r = 0.52; p = 0.006), and alcohol use (r = 0.53; p = 0.005) and COVID-19 incidence rates (per 100,000 inhabitants) (Figure 2A). There was also an association between obesity (r = 0.80; p = 0.0001) and alcohol use prevalence (r = 0.72; p = 0.0001) and COVID-19 mortality rates (per 100,000 inhabitants) (Figure 2B). Finally, there was an association between the prevalence of obesity and the COVID-19 case fatality rate (r = 0.45; p = 0.023) (Figure 2C).

3.2. Regression Analysis in 2020

Table 3 shows the results of a multiple regression analysis adjusted on the association between obesity, smoking, and alcohol use prevalence and the incidence, mortality, and case fatality rate of COVID-19. Results revealed that, adjusted by confounding factors, the prevalence rate of obesity (β = 0.582; p = 0.037) was associated with the COVID-19 mortality rate (per 100,000 inhabitants). There was also an association between obesity and the COVID-19 case fatality rate (β = 0.993; p = 0.014) (Table 3).

3.3. Correlation between the Study Variables and COVID-19 Measures in 2021

There was no association between hypertension, obesity, smoking and alcohol use prevalence and the COVID-19 incidence rate (per 100,000 inhabitants), except for the overweight (r = 0.451; p = 0.023) (Figure 3A). The overweight (r = 0.454; p = 0.022) and obesity prevalence (r = 0.619; p = 0.001) were associated with the COVID-19 mortality rate (per 100,000 inhabitants) (Figure 3B). There was no association between hypertension, obesity, overweight, smoking, and alcohol use prevalence and the COVID-19 case fatality rate (Figure 3C).

3.4. Regression Analysis in 2021

Table 4 shows the results of multiple regression analysis adjusted on the association between obesity and the COVID-19 mortality rate (per 100,000 inhabitants). When adjusted by confounding factors, the prevalence rate of obesity was associated with COVID-19 mortality rate (per 100,000 inhabitants) (β = 0.713; p = 0.028). However, when adjusted by confounding factors including COVID-19 vaccination coverage rates, there was no association between the obesity prevalence and the COVID-19 mortality rate (per 100,000 inhabitants) (β = 0.031; p = 0.895) (Table 4).

4. Discussion

In this nationally based study, we examined the association between health status (prevalence of hypertension, overweight, and obesity) and behavioral factors (tobacco and alcohol use prevalence) and the COVID-19 incidence and mortality rates using an ecological design. Despite relatively high levels of hypertension, overweight and obesity, and behavioral risk factors (tobacco and alcohol use prevalence), no association with the incidence of COVID-19 was found among individuals in the Peruvian population aged ≥15 years during the first and second year of the COVID-19 pandemic. Instead, we identified obesity as being associated with the COVID-19 mortality or fatality rate, which was consistent with findings reported in other studies.

4.1. Potential Explanations and Implications

Initially, we found a positive correlation between the prevalence rates of obesity, tobacco and alcohol use, and the incidence rate due to COVID-19; however, after adjusting for potential confounding factors, there was no association. With regard to obesity, our findings were different from a previous study indicating an association between BMI and cholesterol and the COVID-19 incidence rate in 193 countries [6].
Similar to the previous population-level studies that suggested an association between obesity and COVID-19 mortality in several countries, including a study in Peru [6,12], we found a strong correlation between obesity and mortality and fatality due to COVID-19 in 2020 even after adjusting for confounding factors. The variability in mortality rates observed among different Peruvian regions can be partly explained by the higher obesity prevalence rate. Therefore, we observed that as the obesity prevalence rate increased, the COVID-19 mortality rate increased in the Peruvian population. Obesity is a factor for hospitalization and adverse outcomes from COVID-19 [2,3,5]. Several studies explain the association of obesity and the adverse outcomes and mortality from COVID-19. First, obesity rate is associated with a higher prevalence of hypertension and diabetes, diseases that are associated with worse outcomes from COVID-19 [5,24]. Second, the SARS-CoV-2 virus has a high affinity for angiotensin-converting enzyme 2, which is essential for host cell infection [1,24]. This enzyme is expressed in adipose tissue; therefore, its presence in adipose tissue may exacerbate the severity of SARS-CoV-2 infection and worse outcomes from COVID-19 disease [25]. Third, another potential explanation is that the higher obesity grade may be associated to a more severe COVID-19 among obese individuals [2,3,4,5]. However, these data were not available in our study. Fourth, the association between the obesity rate and COVID-19 mortality rate also could be explained by external factors such as the capacity of the Peruvian health system (lack of ICU beds, oxygen, and medicines), social factors (sex, age or population density) [26,27], and other individual factors such as multimorbidity and frailty.
During the COVID-19 pandemic, concerns were widespread that behavioral risk factors (prevalence of tobacco and alcohol use) may be associated with a severity greater and death due to COVID-19, particularly during the first year of the COVID-19 pandemic. In that setting, we found a positive correlation between alcohol use and COVID-19 mortality in 2020, however, after adjusting for confounding factors in the study population, there was correlation between alcohol use and COVID-19 mortality. Although we found no ecological studies to explain this association, a cohort study among a population of U.S. veterans examined whether behavioral risk factors (tobacco and alcohol use) were associated with COVID-19 mortality, but researchers found no association [9]. The association found in our study may be difficult to explain since we did not investigate the frequency of alcohol consumption and complications among the alcohol-consuming population. Rather, the correlation between alcohol use and COVID-19 mortality may be attributed to external factors such as the overload of the health system (infrastructure, lack of ICU beds, oxygen, and medicines) and other individual factors such as multimorbidity and frailty in the Peruvian population.
Among all the comparisons made in 2021, only a correlation between obesity and COVID-19 mortality was found even after adjusting for confounding factors. Curiously, when adjusted by confounding factors, including COVID-19 vaccination coverage rates in 2021, there was no association between the obesity prevalence rate and the COVID-19 mortality rate. Vaccination against SARS-CoV-2 is a more effective intervention in addressing the COVID-19 pandemic [1,4,28]. Despite patients with obesity showing a reduced immune response to the influenza A (H1N1) vaccine [29], the phase 3 clinical trials of Moderna COVID-19 and Pfizer-BioNTech COVID-19 vaccines included subgroups of individuals with obesity and severe obesity and showed effective results [28,30,31]. The efficacy of the Pfizer-BioNTech vaccine in individuals with obesity from 7 days after the second dose was 95.4% (95% CI: 86.0–99.1%) in preventing COVID-19 [30]. The efficacy of the Moderna vaccine in obese individuals from 14 days after the second dose was 94.1% (95% CI: 89.3–96.8%) [31]. In addition, a study investigated the vaccine’s effectiveness against SARS-CoV-2 infection and severe outcomes among patients with immune-mediated inflammatory diseases and found two vaccine doses to be effective against SARS-CoV-2 infection and severe COVID-19 outcomes in patients with rheumatoid arthritis, ankylosing spondylitis, inflammatory diseases, and psoriasis [32]. Peru implemented the SARS-CoV-2 vaccination in February 2021, initially including immunization of health workers, the population with comorbidities, and subsequently, the overall population. The vaccination led to a reduction in the COVID-19 mortality rate in the Peruvian population [15,16]. In this setting, our findings can help to explain the indirect impact of the SARS-CoV-2 vaccination on the COVID-19 mortality risk reduction in regions where the prevalence of obesity is high.

4.2. Limitations

Our study has several limitations. First is the ecological design; thus, our study shows findings at the population level and should not be generalizable to the individual level. Second, we used several information sources (ENDES 2020 and 2021, INEI, MINSA, SINADEF, and SUSALUD); therefore, our findings may have in a possible bias. Third, our estimates could have been overestimated owing to unmeasured variates, such as obesity grade, multimorbidity, and frailty in the Peruvian population. Fourth, in the multivariable regression analysis adjusted in 2021, we used the COVID-19 vaccination coverage rates in individuals ≥18 years, when the study population were individuals aged ≥15 years. Despite these limitations, our study is robust and large since it included confounding factors in the analysis.

4.3. Implications for Public Health

Our findings have important implications for public health. First, current data from the ENDES in Peru highlight that obesity prevalence rates in individuals aged ≥15 years vary by region, from 9.6 to 35.8% in 2020, and from 10.4 to 37.4% in 2021 [19]. Therefore, these individuals are more likely to have severe disease and higher risk to death, compared with individuals who have a normal weight. Second, our results can help us to explain the high mortality rates from COVID-19 in the Peruvian population. Therefore, although we cannot guarantee that weight loss will reduce COVID-19 mortality in the Peruvian population, interventions for preventing obesity and other chronic diseases are necessary to reduce the impact of future pandemics.

5. Conclusions

Peruvian regions with higher prevalence rates of obesity had higher COVID-19 mortality and case fatality rates in 2020. However, adjusted by COVID-19 vaccination coverage, no association was found between the obesity prevalence and the COVID-19 mortality rate in 2021. These findings support the need for interventions and public health messaging regarding COVID-19 prevention in obese individuals aged ≥15 years. In addition, these findings can help to explain the indirect impact of SARS-CoV-2 vaccination on COVID-19, informing future vaccine strategies.

Author Contributions

Conceptualization, J.C., F.S. and M.C.R.-S.; methodology, J.C. and F.S.; software, M.C.R.-S.; validation M.C.R.-S.; formal analysis, J.C., F.S. and M.C.R.-S.; investigation, J.C., F.S. and M.C.R.-S.; data curation, J.C. and F.S.; writing—original draft preparation, M.C.R.-S.; visualization, M.C.R.-S.; supervision, M.C.R.-S. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by the Universidad Tecnologica del Peru, Lima, Peru.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data underlying the results presented in the study are available from: National Institute of Statistics and Informatics (INEI). Peruvian Population. Available online: https://www.inei.gob.pe/estadisticas/indice-tematico/population-estimates-andprojections/ (accessed on 3 June 2021). National Institute of Statistics and Informatics (INEI); Peru—National Demographic and Health Survey (ENDES) 2020. Noncommunicable and Communicable Diseases. 2020. Available online: https://proyectos.inei.gob.pe/endes/2020/SALUD/ENFERMEDADES_ENDES_2020.pdf (accessed on 20 January 2021). Noncommunicable and Communicable Diseases. 2020. Available online: Available online: https://proyectos.inei.gob.pe/endes/2021/SALUD/ENFERMEDADES_ENDES_2021.pdf (accessed on 22 July 2021). Peruvian Ministry of Health (MINSA). COVID-19 deaths. National System of Deaths (SINADEF). Available online: https://www.datosabiertos.gob.pe/dataset/fallecidos-porcovid-19-ministerio-de-salud-minsa (accessed on 3 June 2021). Superintendencia Nacional de Salud, Peru (SUSALUD). Daily Report on Form F500.2, App. for Centralized Management of the Availability of Hospitalization and ICU Beds at the National Level and of All Subsystems (Application F500.2); SUSALUD: Lima, Peru, 2020; Available online: http://portal.susalud.gob.pe/wp-content/uploads/archivo/registro-camas/2020/diciembre/31/20201231_turnonoche-7.-%20Reporte%20Ejecutivo%20-%20Comando%20COVID-19%20-%20F500.2%20(Camas)-%20ZC.pdf (accessed on 12 August 2022). Daily Report on Form F500.2, App. for Centralized Management of the Availability of Hospitalization and ICU Beds at the National Level and of All Subsystems (Application F500.2); SUSALUD: Lima, Peru, 2021; Available online: http://portal.susalud.gob.pe/wp-content/uploads/archivo/registro-camas/2021/diciembre/31/20211231_turnodia-1.-%20Reporte%20Ejecutivo%20-%20Comando%20COVID-19%20-%20F500.2%20(Camas)-%20ZC_v2.pdf (accessed on 12 August 2022). Peruvian Ministry of Health (MINSA). Single National Health Information Repository (REUNIS). COVID-19 vaccine in Peru. Available online: https://www.minsa.gob.pe/reunis/data/vacunas-covid19.asp (accessed on 31 December 2021).

Acknowledgments

The authors would like to thank the Universidad Tecnológica del Peru for the APC for this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. STROBE Statement—checklist of items that should be included in reports of observational studies.
Table A1. STROBE Statement—checklist of items that should be included in reports of observational studies.
Item NoRecommendation
Title and abstract1(a) Indicate the study’s design with a commonly used term in the title or the abstract
(b) Provide in the abstract an informative and balanced summary of what was done and what was found
Introduction
Background/rationale1Explain the scientific background and rationale for the investigation being reported
Objectives2State specific objectives, including any prespecified hypotheses
Methods
Study design2Present key elements of study design early in the paper
Setting2Describe the setting, locations, and relevant dates, including periods of recruitment, exposure, follow-up, and data collection
Participants2–3(a) Ecological study
Case-control study—Give the eligibility criteria and the sources and methods of case ascertainment and control selection. Give the rationale for the choice of cases and controls
Cross-sectional study—Give the eligibility criteria and the sources and methods of selection of participants
(b) Cohort study—For matched studies, give matching criteria and number of exposed and unexposed
Case-control study—For matched studies, give matching criteria and the number of controls per case
Variables3Clearly define all outcomes, exposures, predictors, potential confounders, and effect modifiers. Give diagnostic criteria, if applicable
Data sources/measurement3 *For each variable of interest, give sources of data and details of methods of assessment (measurement). Describe comparability of assessment methods if there is more than one group
BiasNADescribe any efforts to address potential sources of bias
Study sizeNAExplain how the study size was arrived at
Quantitative variables3Explain how quantitative variables were handled in the analyses. If applicable, describe which groupings were chosen and why
Statistical methods3(a) Describe all statistical methods, including those used to control for confounding
(b) Describe any methods used to examine subgroups and interactions
(c) Explain how missing data were addressed
(d) Cohort study—If applicable, explain how loss to follow-up was addressed
Case-control study—If applicable, explain how matching of cases and controls was addressed
Cross-sectional study—If applicable, describe analytical methods taking account of sampling strategy
(e) Describe any sensitivity analyses
Results
Participants3 *(a) Report numbers of individuals at each stage of study—e.g., numbers potentially eligible, examined for eligibility, confirmed eligible, included in the study, completing follow-up, and analyzed
(b) Give reasons for non-participation at each stage
(c) Consider use of a flow diagram
Descriptive data4 *(a) Give characteristics of study participants (e.g., demographic, clinical, social) and information on exposures and potential confounders
(b) Indicate number of participants with missing data for each variable of interest
(c) Cohort study—Summarize follow-up time (e.g., average and total amount)
Outcome data4–9 *Cohort study—Report numbers of outcome events or summary measures over time
Case-control study—Report numbers in each exposure category or summary measures of exposure
Cross-sectional study—Report numbers of outcome events or summary measures
Main results4–9(a) Give unadjusted estimates and, if applicable, confounder-adjusted estimates and their precision (e.g., 95% confidence interval). Make clear which confounders were adjusted for and why they were included
(b) Report category boundaries when continuous variables were categorized
(c) If relevant, consider translating estimates of relative risk into absolute risk for a meaningful time period
Other analyses4–9Report other analyses done—e.g., analyses of subgroups and interactions and sensitivity analyses
Discussion
Key results9–10Summarize key results with reference to study objectives
Limitations11Discuss limitations of the study, taking into account sources of potential bias or imprecision. Discuss both direction and magnitude of any potential bias
Interpretation9–11Give a cautious overall interpretation of results considering objectives, limitations, multiplicity of analyses, results from similar studies, and other relevant evidence
Generalizability9–11Discuss the generalizability (external validity) of the study results
Other information
Funding11Give the source of funding and the role of the funders for the present study and, if applicable, for the original study on which the present article is based
* Give information separately for cases and controls in case-control studies and, if applicable, for exposed and unexposed groups in cohort and cross-sectional studies. Note: An Explanation and Elaboration article discusses each checklist item and gives methodological background and published examples of transparent reporting. The STROBE checklist is best used in conjunction with this article (freely available on the Web sites of PLoS Medicine at http://www.plosmedicine.org/, Annals of Internal Medicine at http://www.annals.org/, and Epidemiology at http://www.epidem.com/). Information on the STROBE Initiative is available at www.strobe-statement.org.

References

  1. Dhama, K.; Khan, S.; Tiwari, R.; Sircar, S.; Bhat, S.; Malik, Y.S.; Singh, K.P.; Chaicumpa, W.; Bonilla-Aldana, D.K.; Rodriguez-Morales, A.J. Coronavirus Disease 2019–COVID-19. Clin. Microbiol. Rev. 2020, 33, 1–48. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, J.; Zhu, L.; Liu, L.; Zhao, X.; Zhang, Z.; Xue, L.; Yan, X.; Huang, S.; Li, Y.; Cheng, J.; et al. Overweight and Obesity are Risk Factors of Severe Illness in Patients with COVID-19. Obesity 2020, 28, 2049. [Google Scholar] [CrossRef] [PubMed]
  3. Gao, M.; Piernas, C.; Astbury, N.M.; Hippisley-Cox, J.; O’Rahilly, S.; Aveyard, P.; Jebb, S.A. Associations between body-mass index and COVID-19 severity in 6·9 million people in England: A prospective, community-based, cohort study. Lancet Diabetes Endocrinol. 2021, 9, 350–359. [Google Scholar] [CrossRef]
  4. Westheim, A.J.F.; Bitorina, A.V.; Theys, J.; Shiri-Sverdlov, R. COVID-19 infection, progression, and vaccination: Focus on obesity and related metabolic disturbances. Obes Rev. 2021, 22, e13313. [Google Scholar] [CrossRef]
  5. Hendren, N.S.; de Lemos, J.A.; Ayers, C.; Das, S.R.; Rao, A.; Carter, S.; Rosenblatt, A.; Walchok, J.; Omar, W.; Khera, R.; et al. Association of Body Mass Index and Age with Morbidity and Mortality in Patients Hospitalized With COVID-19. Circulation 2021, 143, 135–144. [Google Scholar] [CrossRef]
  6. Sarmadi, M.; Ahmadi-Soleimani, S.M.; Fararouei, M.; Dianatinasab, M. COVID-19, body mass index and cholesterol: An ecological study using global data. BMC Public Health. 2021, 21, 1–14. [Google Scholar] [CrossRef]
  7. Rosoff, D.B.; Yoo, J.; Lohoff, F.W. A genetically-informed study disentangling the relationships between tobacco smoking, cannabis use, alcohol consumption, substance use disorders and respiratory infections, including COVID-19. medRxiv 2021. [CrossRef]
  8. Simons, D.; Shahab, L.; Brown, J.; Perski, O. The association of smoking status with SARS-CoV-2 infection, hospitalization and mortality from COVID-19: A living rapid evidence review with Bayesian meta-analyses (version 7). Addiction 2021, 116, 1319. [Google Scholar] [CrossRef]
  9. Kelly, J.D.; Bravata, D.M.; Bent, S.; Wray, C.M.; Leonard, S.J.; Boscardin, W.J.; Myers, L.J.; Keyhani, S. Association of Social and Behavioral Risk Factors with Mortality Among US Veterans With COVID-19. JAMA Netw. Open 2021, 4, e2113031. [Google Scholar] [CrossRef]
  10. Peruvian Ministry of Health (MINSA). Available online: https://covid19.minsa.gob.pe/ (accessed on 8 September 2022).
  11. Ramírez-Soto, M.C.; Ortega-Cáceres, G.; Arroyo-Hernández, H. Original research: Excess all-cause deaths stratified by sex and age in Peru: A time series analysis during the COVID-19 pandemic. BMJ Open 2022, 12, e057056. [Google Scholar] [CrossRef]
  12. Ramírez-Soto, M.C.; Alarcón-Arroyo, M.; Chilcon-Vitor, Y.; Chirinos-Pérez, Y.; Quispe-Vargas, G.; Solsol-Jacome, K.; Zavaleta, E.Q. Association between Obesity and COVID-19 Mortality in Peru: An Ecological Study. Trop Med Infect Dis. 2021, 6, 182. [Google Scholar] [CrossRef] [PubMed]
  13. Acosta, G.; Escobar, G.; Bernaola, G.; Alfaro, J.; Taype, W.; Marcos, C.; Amado, J. Description of patients with severe COVID-19 treated in a national referral hospital in Peru. Rev. Peru. Med. Exp. Salud Publica 2020, 37, 253–258. [Google Scholar] [CrossRef] [PubMed]
  14. The Peruvian. Peru Starts Vaccination Plan against COVID-19; The Peruvian: Lima, Peru, 2021. [Google Scholar]
  15. Escobar-Agreda, S.; Vargas Herrera, J.; Rojas-Mezarina, L. Preliminary evidence on the effect of vaccination against COVID-19 in Peru. Bol. Inst. Nac. Salud 2021, 27, 35–44. [Google Scholar]
  16. Escobar-Agreda, S.; Silva-Valencia, J.; Rojas-Mezarina, L.; Vargas-Herrera, J. Survival of health workers infected by SARS-CoV-2 in the context of vaccination against COVID-19 in Peru. An. Fac. Med. 2021, 82, 106–112. [Google Scholar]
  17. National Institute of Statistics and Informatics (INEI). Peruvian Population. Available online: https://www.inei.gob.pe/estadisticas/indice-tematico/population-estimates-and-projections/ (accessed on 3 June 2021).
  18. Vandenbroucke, J.P.; von Elm, E.; Altman, D.G.; Gøtzsche, P.C.; Mulrow, C.D.; Pocock, S.J.; Poole, C.; Schlesselman, J.J.; Egger, M. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): Explanation and elaboration. Int. J. Surg. 2014, 12, 1500–1524. [Google Scholar] [CrossRef]
  19. National Institute of Statistics and Informatics (INEI). Peru—National Demographic and Health Survey (ENDES) 2020. Noncommunicable and Communicable Diseases. 2020. Available online: https://proyectos.inei.gob.pe/endes/2020/SALUD/ENFERMEDADES_ENDES_2020.pdf (accessed on 8 September 2022).
  20. National Institute of Statistics and Informatics (INEI). Peru—National Demographic and Health Survey (ENDES) 2021. Noncommunicable and Communicable Diseases. 2021. Available online: https://proyectos.inei.gob.pe/endes/2021/SALUD/ENFERMEDADES_ENDES_2021.pdf (accessed on 8 September 2022).
  21. Peruvian Ministry of Health (MINSA). COVID-19 deaths. National System of Deaths (SINADEF). Available online: https://www.datosabiertos.gob.pe/dataset/fallecidos-por-covid-19-ministerio-de-salud-minsa (accessed on 3 June 2021).
  22. Superintendencia Nacional de Salud, Peru (SUSALUD). Daily Report on Form F500.2, App. for Centralized Management of the Availability of Hospitalization and ICU Beds at the National Level and of All Subsystems (Application F500.2); SUSALUD: Lima, Peru, 2020; Available online: http://portal.susalud.gob.pe/seguimiento-del-registro-de-camas-f500-2/ (accessed on 26 December 2021).
  23. Peruvian Ministry of Health (MINSA). Single National Health Information Repository (REUNIS). COVID-19 vaccine in Peru. Available online: https://www.minsa.gob.pe/reunis/data/vacunas-covid19.asp (accessed on 8 September 2022).
  24. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: A retrospective cohort study. Lancet 2020, 395, 1054–1062. [Google Scholar] [CrossRef]
  25. Al-Benna, S. Association of high level gene expression of ACE2 in adipose tissue with mortality of COVID-19 infection in obese patients. Obes. Med. 2020, 19, 100283. [Google Scholar] [CrossRef]
  26. Schwalb, A.; Seas, C. The COVID-19 Pandemic in Peru: What Went Wrong? Am. J. Trop. Med. Hyg. 2021, 104, 1176–1178. [Google Scholar] [CrossRef]
  27. Taylor, L. COVID-19: Why Peru suffers from one of the highest excess death rates in the world. BMJ 2021, 372, n611. [Google Scholar] [CrossRef]
  28. Townsend, M.J.; Kyle, T.K.; Stanford, F.C. COVID-19 Vaccination and Obesity: Optimism and Challenges. Obesity 2021, 29, 634. [Google Scholar] [CrossRef]
  29. Egawa, Y.; Ohfuji, S.; Fukushima, W.; Yamazaki, Y.; Morioka, T.; Emoto, M.; Maeda, K.; Inaba, M.; Hirota, Y. Immunogenicity of influenza A(H1N1)pdm09 vaccine in patients with diabetes mellitus: With special reference to age, body mass index, and HbA1c. Hum. Vaccin. Immunother. 2014, 10, 1187. [Google Scholar] [CrossRef] [PubMed]
  30. Polack, F.P.; Thomas, S.J.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Perez, J.L.; Pérez Marc, G.; Moreira, E.D.; Zerbini, C.; et al. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. N. Engl. J. Med. 2020, 383, 2603–2615. [Google Scholar] [CrossRef] [PubMed]
  31. Baden, L.R.; El Sahly, H.M.; Essink, B.; Kotloff, K.; Frey, S.; Novak, R.; Diemert, D.; Spector, S.A.; Rouphael, N.; Creech, C.B.; et al. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. N. Engl. J. Med. 2021, 384, 403–416. [Google Scholar] [PubMed]
  32. Widdifield, J.; Kwong, J.C.; Chen, S.; Eder, L.; Benchimol, E.I.; Kaplan, G.G.; Hitchon, C.; Aviña-Zubieta, J.A.; Lacaille, D.; Chung, H. Vaccine effectiveness against SARS-CoV-2 infection and severe outcomes among individuals with immune-mediated inflammatory diseases tested between March 1 and Nov 22, 2021, in Ontario, Canada: A population-based analysis. Lancet Rheumatol. 2022, 4, e430–e440. [Google Scholar] [CrossRef]
Figure 1. Study flow chart using the STROBE reporting guidelines. COVID-19 cases (A) and deaths (B) included in the analysis for 2020 and 2021.
Figure 1. Study flow chart using the STROBE reporting guidelines. COVID-19 cases (A) and deaths (B) included in the analysis for 2020 and 2021.
Ijerph 19 11757 g001
Figure 2. Linear regression models and Spearman correlation between the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use and the incidence rate (A), mortality rate (B), and case fatality rate due to COVID-19 (C) in 2020.
Figure 2. Linear regression models and Spearman correlation between the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use and the incidence rate (A), mortality rate (B), and case fatality rate due to COVID-19 (C) in 2020.
Ijerph 19 11757 g002
Figure 3. Linear regression models and Spearman correlation between the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use and the incidence rate (A), mortality rate (B), and case fatality rate due to COVID-19 (C) in 2021.
Figure 3. Linear regression models and Spearman correlation between the prevalence rates of hypertension, overweight, obesity, smoking, and alcohol use and the incidence rate (A), mortality rate (B), and case fatality rate due to COVID-19 (C) in 2021.
Ijerph 19 11757 g003
Table 1. Prevalence rates of hypertension, overweight and obesity, smoking, and alcohol use and the incidence, mortality, and fatality rates due to COVID-19 in ≥15 year in 25 Peruvian regions in 2020.
Table 1. Prevalence rates of hypertension, overweight and obesity, smoking, and alcohol use and the incidence, mortality, and fatality rates due to COVID-19 in ≥15 year in 25 Peruvian regions in 2020.
RegionHypertension Prevalence (%)Overweight Prevalence (%)Obesity Prevalence (%)Prevalence of Tobacco Use (%)Prevalence of Alcohol Use (%)Population ≥ 15 Years in 2020Incidence Rate (per 100,000 Habitants)Mortality Rate (per 100,000 Habitants)Case Fatality Rate (%)Gender
Balance in Deaths (Men/Women) *
Mean Age in Deaths (Years) *Mean
Monthly
Income
(PEN) *
No. of ICU Beds in 2020 *†
Amazonas7.936.115.814.316.8289,8025852.62201.523.442.064.8992.914
Ancash7.135.621.812.214.1876,7033351.65328.59.82.168.21057.244
Apurímac10.633.214.711.723.5300,3952384.53119.845.031.568.11004.526
Arequipa10.636.828.814.824.31,187,9313962.01358.19.041.268.91530.370
Ayacucho1033.915.513.419.7464,1363162.87182.065.762.166.31095.420
Cajamarca8.638.613.99.813.61,016,7922348.66159.826.82.168.1850.226
Callao1335.331.821.730.2902,6094655.17584.5312.562.167.41355.684
Cusco9.836.616.811.522.8988,8972504.61160.896.421.366.0963.117
Huancavelica9.530.89.610.811.8236,9553083.29161.215.232.066.6669.010
Huánuco836.915.913.616.6524,3713701.39205.015.541.667.5892.449
Ica9.536.533.514.528.8725,6104025.72520.6712.931.767.11478.248
Junín5.839.31717.621.7982,1992728.06245.168.991.865.71082.743
La Libertad10.238.527.812.121.71,531,6682366.96319.4613.51.967.11167.263
Lambayeque10.139.5257.923.8991,1213280.23484.714.781.767.11159.639
Lima8.93728.917.3526.28,750,4174967.75482.949.721.866.91653.5739
Loreto10.429.622.118.625.7680,9273396.69415.0212.222.065.61180.48
Madre de Dios8.643.932.422.829.4135,4286460.26324.95.032.863.01399.98
Moquegua7.937.735.812.625.3155,5459861.45568.325.762.668.81693.718
Pasco6.137.717.216.812.5195,1143198.64167.595.241.763.7834.812
Piura10.134.2257.630.81,535,4332694.29420.415.61.966.8992.681
Puno10.837.920.49.314.1904,2672088.54170.418.161.763.4809.820
San Martín11.336.519.914.925639,5333682.53245.966.682.165.8983.317
Tacna10.538.734.4926.7303,7014613.42262.15.682.365.51259.926
Tumbes11.44027.614.525.8191,8504501.43401.368.921.965.91142.68
Ucayali7.437.72217.530.7416,9324454.68383.528.611.964.31203.118
Abbreviation: COVID-19, coronavirus disease 2019. * Confounding factors to multivariable regression models in 2020 included the mean age in deaths (years), mean monthly income (PEN), gender balance, and number of ICU beds. † SICOVID App. F500.2, SUSALUD, 2020 (accessed on 12 August 2022).
Table 2. Prevalence rates of hypertension, overweight and obesity, smoking, and alcohol use and the incidence, mortality, and fatality rates due to COVID-19 in ≥15 year in 25 Peruvian regions in 2021.
Table 2. Prevalence rates of hypertension, overweight and obesity, smoking, and alcohol use and the incidence, mortality, and fatality rates due to COVID-19 in ≥15 year in 25 Peruvian regions in 2021.
RegionHypertension Prevalence (%)Overweight Prevalence (%)Obesity Prevalence (%)Prevalence of Tobacco Use (%)Prevalence of Alcohol Use (%)Population ≥ 15 Years in 2021Incidence Rate (per 100,000 Habitants)Mortality Rate (per 100,000 Habitants)Case Fatality Rate (%)Gender
Balance in Deaths (Men/Women) *
Mean Age in Deaths (Years) *No. of ICU Beds in 2021 *†Vaccination Coverage in Individuals ≥ 18 Years (%) *
Amazonas8.834.815.613.019.1340,7174084.33203.394.981.766.42068.0
Ancash6.840.221.412.819.5953,8165092.07405.217.961.767.34988.4
Apurímac6.035.113.712.329.6351,0745305.15330.416.231.669.53682.1
Arequipa11.338.428.517.330.51,215,1795560.17451.298.121.665.95384.5
Ayacucho8.133.617.715.523.4519,6563685.52254.596.911.668.22070.6
Cajamarca8.634.615.611.319.31,201,6973387.38212.626.281.768.05475.4
Callao11.638.730.821.425.9877,1616741.64553.388.211.765.49090.0
Cusco10.537.618.113.126.71,111,8684482.73290.056.471.767.53479.0
Huancavelica8.831.110.413.113.3330,5662563.78235.359.181.767.82173.3
Huánuco7.334.218.113.216.8645,2002445.44255.2710.441.967.13169.7
Ica8.235.135.014.429.4700,3944047.44686.0416.951.565.39592.8
Junín5.437.717.322.823.91,063,8495538.47437.667.91.866.26481.2
La Libertad8.337.128.210.823.11,544,9773533.45353.610.011.666.99983.0
Lambayeque9.939.223.88.231.81,050,9822935.35382.513.031.866.45480.8
Lima11.637.630.615.232.08,796,3476206.02536.178.641.765.566288.8
Loreto10.637.020.818.230.8785,3012301.03181.717.91.564.43664.5
Madre de Dios10.639.631.927.926.2129,4653332.17251.037.532.262.02557.8
Moquegua9.738.934.810.837.0157,7048552.1410.264.81.866.92483.7
Pasco6.940.916.317.816.0222,4114259.23323.727.61.464.02981.4
Piura10.035.827.010.035.91,521,2513209.07387.1212.061.566.412283.4
Puno9.936.520.111.119.1989,1252248.05273.4712.162.164.54057.5
San Martín8.536.421.614.930.1701,5173178.97206.556.51.666.22973.6
Tacna12.338.237.48.229.6301,1485434.54393.167.232.063.83575.2
Tumbes10.238.931.013.424.8181,7695291.88451.128.521.965.92186.8
Ucayali7.535.724.916.431.5436,0452579.09342.6313.281.864.52864.3
Abbreviation: COVID-19, coronavirus disease 2019. * Confounding factors to multivariable regression models in 2021 included the mean age in the region (years), mean monthly income (PEN), gender balance, number of ICU beds, and COVID-19 vaccination coverage rates in 2021. † SICOVID App. F500.2, SUSALUD, 2021 (accessed on 12 August 2022).
Table 3. Multiple regression analysis of the prevalence rates of obesity, smoking, and alcohol use with the incidence, mortality, and case fatality rates due to COVID-19 in 2020 *.
Table 3. Multiple regression analysis of the prevalence rates of obesity, smoking, and alcohol use with the incidence, mortality, and case fatality rates due to COVID-19 in 2020 *.
No Adjusted Analysis Full Adjusted Analysis
ModelCoef.SEBeta (β)tpCoef.SEBeta (β)tp
Crude incidence rate per 100,000 habitants
Obesity prevalence130.8138.360.5793.410.00228.3172.890.1250.390.7
Smoking prevalence144.9682.190.3451.760.091−54.8680.7−0.13−0.680.51
Prevalence of alcohol use96.1354.240.3461.770.09−50.7355.3−0.182−0.920.37
Crude mortality rate per 100,000 habitants
Obesity prevalence 15.142.460.7876.130.000111.184.980.5822.250.037
Prevalence of alcohol use15.953.510.6874.540.00018.094.470.3481.810.087
Fatality case rate (%)
Obesity prevalence 0.1850.0890.3972.080.0490.4630.1710.9932.700.014
Abbreviation: COVID-19, coronavirus disease 2019; SE, standard error. * Model is adjusted for the following confounders: mean age in the region (cases or deaths), mean income and gender balance (cases or deaths), and the number of intensive care unit beds.
Table 4. Multiple regression analysis of the prevalence rates of overweight and obesity with the COVID-19 mortality rate in 2021.
Table 4. Multiple regression analysis of the prevalence rates of overweight and obesity with the COVID-19 mortality rate in 2021.
Crude Mortality Rate per 100,000 Habitants (No Adjusted)Crude Mortality Rate per 100,000 Habitants (Full Adjusted)
Coef.SEBeta (β)tpCoef.SEBeta (β)tp
Overweight prevalence17.9210.570.3331.690.1048.8411.070.1640.800.434 *
Obesity prevalence10.672.630.6444.050.00111.804.90.7132.370.028 *
Overweight prevalence **NANANANANA−18.76.35−0.349−2.960.008
Obesity prevalence **NANANANANA0.523.940.0310.130.895
Abbreviation: COVID-19, coronavirus disease 2019; SE, standard error; NA, not applicable. * Model is adjusted for the following confounders: mean age in the region (cases or deaths), mean income and gender balance (cases or deaths), and the number of intensive care unit beds. ** Model is adjusted for the confounders, including COVID-19 vaccination coverage rates in 2021.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Canorio, J.; Sánchez, F.; Ramírez-Soto, M.C. COVID-19, Non-Communicable Diseases, and Behavioral Factors in the Peruvian Population ≥ 15 Years: An Ecological Study during the First and Second Year of the COVID-19 Pandemic. Int. J. Environ. Res. Public Health 2022, 19, 11757. https://doi.org/10.3390/ijerph191811757

AMA Style

Canorio J, Sánchez F, Ramírez-Soto MC. COVID-19, Non-Communicable Diseases, and Behavioral Factors in the Peruvian Population ≥ 15 Years: An Ecological Study during the First and Second Year of the COVID-19 Pandemic. International Journal of Environmental Research and Public Health. 2022; 19(18):11757. https://doi.org/10.3390/ijerph191811757

Chicago/Turabian Style

Canorio, Jordan, Flor Sánchez, and Max Carlos Ramírez-Soto. 2022. "COVID-19, Non-Communicable Diseases, and Behavioral Factors in the Peruvian Population ≥ 15 Years: An Ecological Study during the First and Second Year of the COVID-19 Pandemic" International Journal of Environmental Research and Public Health 19, no. 18: 11757. https://doi.org/10.3390/ijerph191811757

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop