A CLOSER LOOK AT INDIRECT CAUSES OF DEATH AFTER HURRICANE MARIA USING A SEMIPARAMETRIC MODEL

Oscar Lugo, Roberto Rivera
doi: https://doi.org/10.1101/2023.01.03.23284158

Abstract

The Covid-19 pandemic as well as other recent natural emergencies have put the spotlight on emergency planning. One important aspect is that natural disasters or emergencies often lead to indirect deaths and studying the behavior of indirect deaths during emergencies can guide emergency planning. While many studies have implied a large number of indirect deaths in Puerto Rico due to Hurricane Maria; the specific causes of these deaths have not being carefully studied. In this paper, we use a semiparametric model and mortality data to evaluate cause of death trends. Our model adjusts for cause of death effect potentially varying over time while also inferring on how long excess deaths occurred. From September 2017 to March 2018, after adjusting for intra-annual variability and population displacement, we find evidence of significant excess deaths due to Alzheimer’s/Parkinson, heart disease, sepsis, diabetes, renal failure, and pneumonia & influenza. In contrast, for the same time period we find no evidence of significant excess deaths due to cancer, hypertension, respiratory diseases, cerebrovascular disease, suicide, homicide, falling accidents and traffic accidents.

1 INTRODUCTION

On September 20 2017, Hurricane Maria struck Puerto Rico with 155 miles per hour sustained winds and wind gusts of up to 190 mph (Pasch et al., 2018). Torrential rains of up to 40 inches of water accompanied the winds (Hu and Smith, 2018), which resulted in large floods on the entire island. The hurricane destroyed the island’s electricity grid, leaving 3.4 million inhabitants without electricity (Lu and Alcantara, 2020). In addition, infrastructure was affected: aqueducts and sewers, roads, bridges, ports, airports, cell phone towers, and gas stations were left without fuel for vehicles and generators (Ficek, 2018). It is estimated that the effects of the hurricane caused monetary damages in Puerto Rico for more than the US $82 billion (Álvarez Febles and Félix, 2020). Also, Maria destroyed around 80% of the value of the island’s crops, representing a loss of $780 million in agricultural yield (Rodríguez Cruz and Niles, 2018).

Indirect causes of death have been recognized as an important consequence of national emergencies. For example, damage to the transportation network during an earthquake can affect access to medicine and health services generating serious or even fatal repercussions for affected people (Geronimus et al., 2006; Haentjens et al., 2010; WHO, 2018a). However, the different interpretations of guidelines and varied training of personnel determining the cause of death hinder mortality counts (German et al., 2001). Several studies have looked into excess deaths due to Hurricane Maria (Santos-Lozada and Howard, 2017; Rivera and Rolke, 2018; Kishore et al., 2018; Santos-Burgoa et al., 2018; Rivera and Rolke, 2019; Cruz-Cano and Mead, 2019). While Cruz-Cano and Mead (2019) assessed specific indirect causes of death after Hurricane Maria, it only did so until the end of October 2017. Moreover, the authors limited their analysis to a few causes (heart disease, other causes, diabetes, Alzheimer’s disease, and septicemia), and did not account for population displacement, possibly leading to an underestimation of excess deaths. Importantly, the duration of Hurricane Maria’s excess mortality effect may have differed by cause of death, and modeling this can guide the development of emergency management planning.

For future emergency planning and preparedness purposes, it is important to predict fatalities related to the state of emergency, whether the deaths are directly or indirectly related to the event (Brunkard et al., 2008; Ragan et al., 2008; Rivera et al., 2020; Rosenbaum et al., 2021) and excess mortality is a useful method of assessing deaths due to an emergency (Hiam et al., 2017). This research focuses on building a semi-parametric Poisson model that allows estimating the deaths caused by a state of emergency incorporating in the model the intra-annual variability, the population displacement, and causes of death due to non-communicable diseases, and other causes indirect. At the same time, the model captures the duration of the effect on mortality of an emergency according to each death cause.

2 A Semiparametric Model

We propose a semiparametric Poisson regression model Wood (2017) to estimate whether some deaths due to some causes increased after the passage of Hurricane Maria in Puerto Rico while controlling for the temporal variation in the mortality rate of the causes of death and the population size. Population displacement plays an important role, which can be directly influenced by economic, social, labor, academic, climatic variations, and other factors (Neyra, 2017; Acosta et al., 2020).

Additionally, index functions will be deployed in the model to evaluate if there are effects of changes in the behavior of certain leading causes of death including some Non-Communicable Diseases (NCDs) in different post-emergency periods.

Let Dt = number of deaths at time index t, Nt = population size at time t. We assume that Dt follows a Poisson distribution. For i = 0, …, I, we assign pi,t as an indicator of the period i for time t. Specifically, i = 0 represents the pre-emergency period, i = 1 is a given period after the emergency, and so on. For j = 1, … J, we consider cj,t as an indicator of cause of death j at time t. On the other hand, we define dayt = day of the year, and yeart = a categorical variable that captures the mortality trend through the years. We define the following generalized additive semiparametric Poisson regression model with interaction (Ruppert et al., 2003), Embedded Image where μi,j,t = E(Dt t, Embedded Image, Nt, dayt, yeart). We define Embedded Image and α = (α0, …, αI)′ denotes the coefficients of the model for the I pre and post-emergency periods. Similarly, Embedded Image and γ = (γ1, γ2, …, γJ)′ denotes the coefficients of the model for the J causes of death. While (pt⊗ ct)′ = (p0,tc1,t, …, pI,tcJ,t) is the Kronecker product between vectors Embedded Image and Embedded Image represents the interaction between the I pre or post-emergency periods and the J causes of death at time t and ω = (ω11, …, ωMQ)′ represents the coefficients of the interaction. On the other hand, the smooth function f (dayt) captures within-year variation in deaths while β0 is the contrast of the smooth function and the explanatory variables, and ϵi,j,t contains the variability not explained by the other terms. For post-emergency periods, pi,t = 1 and the interaction coefficient make it possible to adapt the effect of the post-emergency period for cause j. The coefficients of the model (1) are estimated by penalized likelihood maximization, taking into account the population displacement caused by a state of emergency while the smoothing penalty parameters are determined by the restricted maximum likelihood (REML). Model (1) can be used to estimate excess death of cause of death j during period i post-emergency of index time t through the difference between the estimation of the model with pi,t = 1, cj,t = 1 and pi,tcj,t = 1 versus the estimated model with pi,t = 0, cj,t = 1 and pi,tcj,t = 0. For i ≥ 1 and j ≥ 1 we obtain: Embedded Image Embedded Image and êi,j estimate β0, αi, γj and ωij respectively. Also, notice that Embedded Image uses Embedded Image, a population size unaltered by post-emergency displacement because âi and êij are already considering the impact of population displacement. Using it again for Embedded Image would imply accounting for the effect twice. This way, Embedded Image The expression Embedded Image in (1) represents the number of deaths typically seen per and Embedded Image represents percentage changes during the post-emergency periods adjusting for the interaction between the J causes of death and the post-emergency periods. Equation (1) is the maximum likelihood estimator for the excess deaths expected from the J causes of death in post-emergency period i. Using equation (1) the accumulated excess of the J causes of death is, Embedded Image for any period that begins at the index q and ends at r.

In general terms, model (1) allows statistical inference of the duration of the effect in the I post-emergency periods, the effect on the J causes of death during a state of emergency, and the effect caused by the interaction between the I post-emergency periods and the J causes of death in the population displacement.

3 Results

3.1 Data

The data used during this study correspond to the deaths reported daily by the Puerto Rico Demographic and Vital Statistics Registry between January 1, 2014 to December 31, 2018 and includes a total of 148,490 records. This work is exempt from requiring human subjects board review. Death records without autopsy are processed in the software “Automatic Classification of Medical Entry”(ACME), this program selects the underlying cause of death (Variables Coding and ACME Description) using the information based on the codes of the International Classification of Diseases (ICD-10 (2016)) on all medical conditions reported for death along with information on the reported relationship between the various conditions (Lu, 2003).

3.2 Population displacement

Demographic data indicates that Puerto Rico registers a high migratory movement, due to the economic crisis it has faced since 2006 (DeWaard et al., 2020; Castro-Prieto et al., 2017; Abel and Deitz, 2014; Rivera, 2016). After the passage of Hurricane Maria in Puerto Rico, lack of water or water service (DeWaard et al., 2020), and a health service crisis (Panditharatne, 2018; Watkins et al., 2020; Glassman, 2019) lead to an increase in population displacement, which can permanently alter the demographic composition of the affected regions (Acosta et al., 2020).

Like Rivera and Rolke (2019) we estimate population displacement through annual estimates of the resident population of Puerto Rico according to the data US Census Bureau adjusted for passenger movement. Usually, the movement passenger is a biased indirect variable of resident migration since it includes the visitor seasonal movements (Velazquez-Estrada et al., 2018), and tourism has been one of the few growing sectors in the Puerto Rican economy (Rivera, 2016). However, a dramatic drop in visitors (Rivera and Rolke, 2019) occurred the months after the arrival of Hurricane Maria. Other methods like the use of disaster maps (Echenique and Melgar, 2018) did not significantly affect our results.

Data from daily death certificates (Dt) were combined with population estimates (Nt) to produce daily mortality rates (Rt); Embedded Image Figure 1 illustrates the daily mortality rate during the years 2014-2018 in Puerto Rico. This graph indicates that the first and last months of each year generally have the highest death rates a pattern commonly seen in other locations (Phillips et al., 2010). On the other hand, Figure 1 displays a drastic spike in mortality rate in the last days of September 2017 and remained high until the beginning of 2018. It is not clear from the figure whether mortality rates increased over the years, yet our model of excess deaths takes the variable yeart can explain this type of temporal fluctuation.

Figure 1:
Figure 1:

Average daily mortality rate from 2014-2018.

In the preliminary analysis, the models did not indicate over-dispersion, so a Poisson model was used.

3.3 Cumulative Excess Deaths by Cause

According to the World Health Organization (WHO, 2018b), the Institute of Statistics of Puerto Rico (Health Department, 2013), and Cruz-Cano and Mead (2019), the main causes of death in Puerto Rico are cancer, hypertension, cardiovascular, kidney, respiratory diseases, and conditions that affect the nervous and endocrine systems. Other causes of death of interest are traffic accidents, homicides, and suicides (Petrosky et al., 2020). In this study, we consider these causes of death along with sepsis, accidents due to any type of fall, and ‘undefined deaths’ to determine if said causes presented an excess of mortality after the natural disaster.

No interannual pattern is observed in the selected causes of death (Figure 2). On the other hand, deaths due to “Other Causes of Death”, “Heart”, “Cancer” and “Diabetes” had the highest average monthly mortality rate, followed by deaths due to “Alzheimer’s or Parkinson’s”, “Respiratory Diseases”, “Hypertension”, “Cerebrovascular”, “Sepsis”, “Pneumonia and Influenza”, “Renal Insufficiency”, “Homicide”, and “Undefined Causes”. Lastly, “Suicide”, “Fall Accidents” and “Traffic” display low mortality rates. We use the class “Other Causes of Death” and the year “2014” as the reference categories for these variables.

Figure 2:
Figure 2:

Monthly mortality rates for several causes of death during the years 2014-2018 in Puerto Rico.

Table 1 presents the estimates of excess and cumulative excess deaths for the 14 causes of death according to 7 post-Maria periods and their respective 95% confidence intervals based on model (1).

View this table:
Table 1:

Estimates of excess death and cumulative excess death for the different causes of death according to model 1 and their respective 95% confidence intervals.

21 excess deaths during September 20th to 30th, 2017 were estimated for “Cancer”, this period being the only one that presented excess deaths. Consequently, the cumulative total excess death intervals with 95% confidence from October 2017 to March 2018 included zero. For October 2017, 148 excess deaths are estimated for “Heart Diseases”. For the periods from September 20 to 30, December, February, and March, the excess ranged from 80 to 99 deaths. While, from September to March we are 95% confident that (284, 916) excess heart disease deaths occurred. For “Diabetes”, the first two post-Maria periods showed an excess of 93 and 56 deaths, respectively. For the entire post-Maria period, considerable uncertainty is found on cumulative excess deaths: (28, 501). It was found that for “Alzheimer’s or Parkinson’s”, “Pneumonia & Influenza”, “Sepsis” “Renal Insufficiency” and “Undefined” excess deaths during the entire post-Maria period were significant. According to the 95% confidence intervals for causes “Suicide”, “Falling Accidents” and “Traffic”, there is no statistically significant evidence of excess deaths from September 20-30, 2017 to March 2018.

Overall, the causes of death with the largest excess deaths during the seven post-emergency periods were: “Heart” and “Undefined”. It should be emphasized that the excess from November 2017 to January 2018 for “Pneumonia and Influenza” were not necessarily due to Maria as seasonal cases are expected (Paz-Bailey et al., 2020).

Figures 3 and 4 display cumulative excess deaths for all causes; in red is the estimated value while black functions are upper and lower 95%. confidence bounds.

Figure 3:
Figure 3:

Cumulative excess death estimates from September 20 to March 31, 2018 according to: other causes, heart disease, cancer, diabetes, Alzheimer or Parkinson’s, respiratory disease, cardiovascular disease, sepsis, pneumonia, and homicide. Dashed line represents the 64 deaths authorities initially attributed to Hurricane Maria.

Figure 4:
Figure 4:

Cumulative excess death estimates from September 20 to March 31, 2018 according to: suicide, fall, traffic accident, hypertension, undefined, and renal insufficiency. Dashed line represents the 64 deaths authorities initially attributed to Hurricane Maria.

4 Discussion

This study explores the excess mortality due to non-communicable diseases and other indirect causes during the state of emergency generated after Hurricane Maria’s passage over Puerto Rico. The model makes adjustments for intra-annual variability and population displacement. This allows greater precision when estimating excess mortality. Moreover, our method can infer on how long an emergency leads to an increase in deaths. To achieve our objective, the statistical model from Rivera and Rolke (2019) was extended to account for different causes of deaths and included cause of death and post-emergency interaction terms so that the model can adjust the effect during the post-emergency period i for each cause of death.

It is determined that the causes of deaths due to sepsis, and diabetes, presented the highest excess deaths for residents of Puerto Rico during September and October, 2017 compared to the pre-Maria period and relative to “Other Causes of death”. Excess death due to heart disease, Alzheimer’s/Parkinson, pneumonia & influenza and renal failure persisted for months after landfall.

The excess deaths observed from renal insufficiency may be attributed to the dialysis centers being left without electric power and electric generators not being designed to be used for long periods (Michaud and Kates, 2017). On the other hand, an increase in sepsis deaths may be due to health problems caused by a lack of supply of potable water since 20 out of 51 wastewater treatment plants were out of service one month after the arrival of Hurricane Maria (Michaud and Kates, 2017). For its part, the excess mortality from this disease was caused by multiorgan dysfunction syndrome, due to poor access to care in intensive care units and the shortage of drug supply (Driessen et al., 2021). Furthermore, Alzheimer’s excess deaths could be due to people who suffer from this disease tending to live in care centers, and the state of emergency led to a reduction in their medical care, routine treatments (Ruhm, 2021). Difficulty in continuing with adequate nutrition adequate medical care, lack and loss of medicines, and difficulties in monitoring their levels of glucose resulted in an increase in deaths due to diabetes (Quast et al., 2019).

Cruz-Cano and Mead (2019) evaluated some causes of death that were traced and reported by the Department of Health of Puerto Rico but only until October 2017. Our estimates are comparable to theirs, though we study more causes of deaths, find evidence Hurricane Maria lead to excess deaths until March 2018, we find significant renal failure excess deaths and infer there was an increase in deaths where the cause was left undefined.

Estimating excess deaths in Puerto Rico by cause after Hurricane Maria helps us understand the effect of the natural disaster and allows us to assess the weaknesses in the health system to better face these types of events in the future (Perreira et al., 2016). Authorities need tools to reduce the number of deaths due to state emergencies. Some strategies include the use of technology in telemedicine, applications for mobile devices, virtual programs (Hacker et al., 2021) and remote monitoring Bowman et al. (2021); LRM (2020). Taking into account the above, it is expected that these measures will reduce the risk of death in the population and better prepare health professionals for future states of emergency generated by different disasters (e.g. hurricanes, tsunamis, and pandemic).

Competing interests: The author(s) declare none

Data Availability

Available upon request

References

  1. ↵
    Abel, J. R. and Deitz, R. (2014). The causes and consequences of puerto rico’s declining population. Current issues in Economics and Finance, 20(4).
  2. ↵
    Acosta, R. J., Kishore, N., Irizarry, R. A., and Buckee, C. O. (2020). Quantifying the dynamics of migration after hurricane maria in puerto rico. Proceedings of the National Academy of Sciences, 117(51):32772–32778.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    Álvarez Febles, N. and Félix, G. (2020). Hurricane maría, agroecology and climate change resiliency. Climate Justice and Community Renewal: Resistance and Grassroots Solutions, pages 1–14.
  4. ↵
    Bowman, C. A., Nuh, M., and Rahim, A. (2021). Covid-19 telehealth expansion can help solve the health care underutilization challenge. The American Journal of Managed Care, 27(1):9–11.
    OpenUrl
  5. ↵
    Brunkard, J., Namulanda, G., and Ratard, R. (2008). Hurricane Katrina deaths, Louisiana, 2005. Disaster Med Public Health Prep, 2(4):215–223.
    OpenUrlCrossRefPubMedWeb of Science
  6. ↵
    Castro-Prieto, J., Martinuzzi, S., Radeloff, V., Helmers, D., Quiñones, M., and Gould, W. (2017). Declining human population but increasing residential development around protected areas in puerto rico. Biological Conservation, 209:473–481.
    OpenUrl
  7. ↵
    Cruz-Cano, R. and Mead, E. L. (2019). Causes of excess deaths in puerto rico after hurricane maria: a time-series estimation. American Journal of Public Health, 109(7):1050–1052.
    OpenUrlPubMed
  8. ↵
    DeWaard, J., Johnson, J. E., and Whitaker, S. D. (2020). Out-migration from and return migration to puerto rico after hurricane maria: evidence from the consumer credit panel. Population and Environment, pages 1–15.
  9. ↵
    Driessen, R. G., Heijnen, N. F., Hulsewe, R. P., Holtkamp, J. W., Winkens, B., van de Poll, M. C., van der Horst, I. C., Bergmans, D. C., and Schnabel, R. M. (2021). Early icu-mortality in sepsis–causes, influencing factors and variability in clinical judgement: a retrospective cohort study. Infectious Diseases, 53(1):61–68.
    OpenUrl
  10. ↵
    Echenique, M. and Melgar, L. (2018). Mapping puerto rico’s hurricane migration with mobile phone data. Retrieved March, 17:2019.
    OpenUrl
  11. ↵
    Ficek, R. E. (2018). Infrastructure and colonial difference in puerto rico after hurricane maría. Transforming Anthropology, 26(2):102–117.
    OpenUrl
  12. ↵
    German, R., Lee, L., Horan, J., Milstein, R., Pertowski, C., Waller, M., for Disease Control,. G. W. G. C., and (CDC), P. (2001). Updated guidelines for evaluating public health surveillance systems: Recommendations from the guidelines working group. MMWR. Recommendations and Reports: Morbidity and Mortality Weekly Report, 50(13):1–35.
    OpenUrl
  13. ↵
    Geronimus, A., Colen, C., Shochet, T., Ingber, L., and James, S. (2006). Urban-rural differences in ex-cess mortality among high-poverty populations: Evidence from the harlem household survey and the pitt county, north carolina study of african american health. Journal of Health Care for the Poor and Underserved. Bulletin of the American Meteorological Society, 17(3):532–558.
    OpenUrl
  14. ↵
    Glassman, B. (2019). A third of movers from puerto rico to the mainland united states relocated to florida in 2018. census. gov, 26.
  15. ↵
    Hacker, K. A., Briss, P. A., Richardson, L., Wright, J., and Petersen, R. (2021). Peer reviewed: Covid-19 and chronic disease: The impact now and in the future. Preventing Chronic Disease, 18.
  16. ↵
    Haentjens, P., Magaziner, J., Colon-Emeric, C. S., Vanderschueren, D., Milisen, K., Velkeniers, B., and Boonen, S. (2010). Meta-analysis: Excess mortality after hip fracture among older women and men. Annals of Internal Medicine. Bulletin of the American Meteorological Society, 152(6):380–390.
    OpenUrl
  17. ↵
    Health Department, Secretaría Auxiliar de Planificación y Desarrollo, D. d. A. E. (2013). Mortalidad en puerto rico. análisis de situación de la salud en puerto rico. Estado libre asociado de Puerto Rico. Depar-tamento de Salud.
  18. ↵
    Hiam, L., Dorling, D., Harrison, D., and McKee, M. (2017). Why has mortality in england and wales been increasing? an iterative demographic analysis. Journal of the Royal Society of Medicine, 110(4):153–162.
    OpenUrlCrossRefPubMed
  19. ↵
    Hu, T. and Smith, R. B. (2018). The impact of hurricane maria on the vegetation of dominica and puerto rico using multispectral remote sensing. Remote Sensing, 10(6):827.
    OpenUrl
  20. ↵
    ICD-10 (2016). International Statistical Classification of Diseases and Related Health Problems 10th, ICD-10, volume 3.
  21. ↵
    Kishore, N., Marqués-Reyes, D., Mahmud, A., Kiang, M., Rodriguez Rivera, I., Fuller, A., Ebner, P., Sorensen, C., Racy, F., Lemery, J., Maas, L., Leaning, J., Irizarry, R., Balsari, S., and Buckee, C. (2018). Mortality in puerto rico after hurricane maria. New England Journal of Medicine, 379.
  22. ↵
    Lrm, T. L. R. M. (2020). Covid-19 heralds a new era for chronic diseases in primary care. The Lancet. Respiratory Medicine, 8(7):647.
    OpenUrl
  23. ↵
    Lu, D. and Alcantara, C. (2020). Después del paso del huracán maría, puerto rico estuvo a oscuras durante 181 días, 6 horas y 45 minutos. The Washington Post. Noviembre 29.
  24. ↵
    Lu, T.-H. (2003). Using acme (automatic classification of medical entry) software to monitor and improve the quality of cause of death statistics. Journal of Epidemiology & Community Health, 57(6):470–471.
    OpenUrlFREE Full Text
  25. ↵
    Michaud, J. and Kates, J. (2017). Public health in puerto rico after hurricane maria. Kaiser family Foundation. Available: http://files.kff.org/attachment/Issue-Brief-Public-Health-in-Puerto-Rico-after-Hurricane-Maria [Accessed 6 May 2018].
  26. ↵
    Neyra, P. C. D. (2017). Los desplazados ambientales por efectos del cambio climático y los intentos de determinación de responsabilidades por daños climáticos. Revista IUS, 1(13).
  27. ↵
    Panditharatne, M. (2018). Six months after maria: Puerto rico’s growing health crisis. Natural Resources Defense Council, March, 21.
  28. ↵
    Pasch, R. J., Penny, A. B., and Berg, R. (2018). National hurricane center tropical cyclone report: Hurricane maria. TROPICAL CYCLONE REPORT AL152017, National Oceanic And Atmospheric Administration and the National Weather Service, pages 1–48.
  29. ↵
    Paz-Bailey, G., Quandelacy, T. M., Adams, L. E., Olsen, S. J., Blanton, L., Munoz-Jordan, J. L., Lozier, M., Alvarado, L. I., and Johansson, M. A. (2020). Recent influenza activity in tropical puerto rico has become synchronized with mainland us. Influenza and other respiratory viruses, 14(5).
  30. ↵
    Perreira, K., Lallemand, N., Napoles, A., and Zuckerman, S. (2016). Environmental scan of puerto rico’s health care infrastructure. Washington (DC): Urban Institute.
  31. ↵
    Petrosky, E., Ertl, A., Sheats, K. J., Wilson, R., Betz, C. J., and Blair, J. M. (2020). Surveillance for violent deaths—national violent death reporting system, 34 states, four california counties, the district of columbia, and puerto rico, 2017. MMWR Surveillance Summaries, 69(8):1.
    OpenUrlCrossRef
  32. ↵
    Phillips, D., Barker, G. E., and Brewer, K. M. (2010). Christmas and new year as risk factors for death. Social science & medicine, 71(8):1463–1471.
    OpenUrl
  33. ↵
    Quast, T., Andel, R., and Sadhu, A. R. (2019). Long-term effects of disasters on seniors with diabetes: evidence from hurricanes katrina and rita. Diabetes care, 42(11):2090–2097.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    Ragan, P., Schulte, J., Nelson, J., and KT. (2008). Mortality surveillance: 2004 to 2005 Florida hurricane-related deaths. Am J Forensic Med Pathol, 29(2):148–153.
    OpenUrlCrossRefPubMed
  35. ↵
    Rivera, R. (2016). A dynamic linear model to forecast hotel registrations in puerto rico using google trends data. Tourism Management, 57:12–20.
    OpenUrl
  36. ↵
    Rivera, R. and Rolke, W. (2018). Estimating the death toll of hurricane maria. Significance, 15(1):8–9.
    OpenUrl
  37. ↵
    Rivera, R. and Rolke, W. (2019). Modeling excess deaths after a natural disaster with application to hurricane maria. Statistics in Medicine, 38(23):4545–4554.
    OpenUrlCrossRefPubMed
  38. ↵
    Rivera, R., Rosenbaum, J., and Quispe, W. (2020). Excess mortality in the united states during the first three months of the covid-19 pandemic. Epidemiology & Infection, pages 1–20.
  39. ↵
    Rodríguez Cruz, L. and Niles, M. (2018). Hurricane maria’s impacts on puerto rican farmers: Experience, challenges, and perceptions.
  40. ↵
    Rosenbaum, J. E., Stillo, M., Graves, N., and Rivera, R. (2021). Timeliness of us mortality data releases during the covid-19 pandemic: delays are associated with electronic death registration system and elevated weekly mortality. Journal of Public Health Policy, 42:536–549.
    OpenUrl
  41. ↵
    Ruhm, C. J. (2021). Excess deaths in the united states during the first year of covid-19. Technical report, National Bureau of Economic Research.
  42. ↵
    Ruppert, D., Wand, M. P., and Carroll, R. J. (2003). Semiparametric regression. Number 12. Cambridge university press.
  43. ↵
    Santos-Burgoa, C., Sandberg, J., Suárez, E., Goldman-Hawes, A., Zeger, S., Garcia-Meza, A., Pérez, C. M., Estrada-Merly, N., Colón-Ramos, U., Nazario, C. M., et al. (2018). Differential and persistent risk of excess mortality from hurricane maria in puerto rico: a time-series analysis. The Lancet Planetary Health, 2(11):e478–e488.
    OpenUrl
  44. ↵
    Santos-Lozada, A. R. and Howard, J. T. (2017). Estimates of excess deaths in puerto rico following hurricane maria. SocArxiv, pages 1–8. doi:10.31235/osf.io/s7dmu.
    OpenUrlCrossRef
  45. ↵
    Velazquez-Estrada, A., Marazzi, M., Orenstein-Cardona, J., and Romero, N. (2018). Perfil del migrante: 2016. tech. rep. Puerto Rico Institute of Statistics.
  46. ↵
    Watkins, D. J., Zayas, H. R. T., Vega, C. M. V., Rosario, Z., Welton, M., Arroyo, L. D. A., Cardona, N., Reguero, Z. J. D., Rivera, A. S., Huerta-Montanez, G., et al. (2020). Investigating the impact of hurricane maria on an ongoing birth cohort in puerto rico. Population and Environment, pages 1–17.
  47. ↵
    WHO (2018a). Definitions: emergencies. Retrieved from:. World Health Organization.
  48. ↵
    WHO (2018b). The top 10 causes of death. 2018.[internet].[consultado febrero 2020]. Disponible en: https://www.who.int/news-room/fact-sheets/detail/the-top-10-causes-of-death.
  49. ↵
    Wood, S. N. (2017). Generalized additive models: an introduction with R. Chapman and Hall/CRC.
Back to top
Posted January 04, 2023.
Download PDF
Data/Code
Email
A CLOSER LOOK AT INDIRECT CAUSES OF DEATH AFTER HURRICANE MARIA USING A SEMIPARAMETRIC MODEL
Oscar Lugo, Roberto Rivera
medRxiv 2023.01.03.23284158; doi: https://doi.org/10.1101/2023.01.03.23284158
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area