Volume 18, Number 1—January 2012
The 1918–19 Influenza Pandemic in Boyacá, Colombia
To quantify age-specific excess-mortality rates and transmissibility patterns for the 1918–20 influenza pandemic in Boyacá, Colombia, we reviewed archival mortality records. We identified a severe pandemic wave during October 1918–January1919 associated with 40 excess deaths per 10,000 population. The age profile for excess deaths was W shaped; highest mortality rates were among infants (<5 y of age), followed by elderly persons (>60 y) and young adults (25–29 y). Mean reproduction number was estimated at 1.4–1.7, assuming 3- or 4-day generation intervals. Boyacá, unlike cities in Europe, the United States, or Mexico, experienced neither a herald pandemic wave of deaths early in 1918 nor a recrudescent wave in 1920. In agreement with reports from Mexico, our study found no death-sparing effect for elderly persons in Colombia. We found regional disparities in prior immunity and timing of introduction of the 1918 pandemic virus across populations.
Quantitative analyses of age-specific death rates, transmissibility, and dissemination patterns of the 1918 influenza pandemic in the United States (1,2), Mexico (3), Peru (4), Japan (5), Europe (6,7), Taiwan (8), and Singapore (9) have shed light on the epidemiology of the most devastating pandemic in recent history (10). These studies revealed the pandemic’s unusual severity in young adults, occurrence in multiple waves, and higher transmission potential than that of seasonal epidemics (11). However, quantitative historical studies remain scarce for Latin America, Africa, and Asia, where our understanding of influenza disease patterns remains particularly weak.
The emergence of the pandemic influenza A (H1N1) 2009 virus in Mexico (12,13) reinforced the need to understand the epidemiology of past pandemics in the Americas to inform preparedness plans. We therefore analyzed death patterns for the 1918 influenza pandemic in Boyacá, a rural area in central Colombia, where influenza seasonality is less defined than in temperate regions (14). By using archival records, we quantified the age-specific excess-death rates and transmission potential of the 1918–19 pandemic in Boyacá and compared these findings with those reported for other locations, especially Mexico City, Mexico.
Materials and Methods
Boyacá is located in the central part of Colombia within the Andes Mountains at latitude ≈5.5°N (Figure 1). In 1918, the population of Boyacá was 659,947 and <50% of the area was occupied. Hygienic conditions were poor. A centralized disease notification system was lacking; however, death records were maintained by parishes.
The climate in Boyacá varies from high humidity and high mean temperature (≈40°C) in low areas near the Magdalena River (altitude 600 m) to cold mean temperature (<6°C) and permanent snow in the Cocuy Mountains (altitude 5,500 m). The 2 rainy seasons, April–May and October–November, produce ≈1,000 mm3/rainfall/year.
Historical Death Records
A total of 32,843 death records, written mostly by Catholic priests and corresponding to January 1917–December 1920, were manually retrieved from the parish archives of 78 municipalities in the department of Boyacá. From these archival records, we extracted age, cause, and exact date of death. To estimate mortality rates, we compiled weekly numbers of deaths from all causes and from respiratory illness, stratified into 5-year age groups (Figure 2, Figure 3). To obtain precise estimates of the transmission potential, we compiled daily death time series, combining all age groups.
We obtained age-specific estimates of population size for the department of Boyacá from a 1918 census report (15). In 1918, ≈70% of Boyacá’s population was located in rural areas. During 1912–1918, the average annual population growth rate in Boyacá was 1.7%; during 1918–1920, it was 3.8%.
Estimation of Excess Deaths
For characterization of mortality rates for the Boyacá pandemic, influenza-associated mortality rates must be separated from background mortality rates (deaths from respiratory illness other than influenza) and considered separately for each age group and cause of death (respiratory or all causes). To estimate pandemic mortality rate, we can define a discrete period of pandemic influenza activity and estimate the number of deaths in excess of background deaths that occurred during the pandemic period. Because mortality rates tend to vary seasonally throughout the year, our background estimate must also vary seasonally. To find the best estimate for baseline mortality rate in the absence of influenza activity, we applied regression methods, using harmonic terms and time trends, to mortality rate data (6,16,17) (Technical Appendix [PDF - 186 KB - 4 pages]).
The regression model determines the extent to which observed weekly mortality rate fit the expectation of background mortality rate. Periods of poor fit indicate that observed mortality rate exceeds typical baseline levels, presumably because of increased influenza activity.
We defined pandemic periods as the weeks when deaths from respiratory illness exceeded the upper limit of the 95% CI of the background model. To estimate the mortality rate during the pandemic, for each age group we summed the weekly number of deaths from respiratory illness and from all causes that exceeded model baseline rates during each pandemic period during 1918–20.
To ensure that our estimates were not sensitive to modeling assumptions, we also estimated excess deaths by using an alternative approach to calculate background deaths. In this approach, background mortality rates for a given week are obtained by averaging mortality rates during the same week in previous years (Technical Appendix [PDF - 186 KB - 4 pages]).
Finally, we estimated a relative measure of the effects of pandemic-associated deaths for each age group, which considers the typical mortality rate experienced by that age group. We calculated relative risk for pandemic-associated death, defined as the ratio of excess deaths during pandemic periods to expected baseline deaths. Relative risk has been shown to facilitate comparison between age groups or countries, which have different background risks for death (17,18).
Comparing Patterns of Age-specific Deaths
We compared patterns of age-specific excess deaths from the 1918–19 Boyacá pandemic with those recently published for Mexico City (3). The estimates for Mexico City were based on excess-death rates obtained from monthly pneumonia and influenza records (1916–1920), stratified by 6 age groups (<5, 5–19, 20–29, 30–49, 50–69, and >70 years). Excess-death rates for Mexico City were calculated with a method similar to that used in this study.
We also reviewed key epidemiologic features of the pandemic in various locations as recently reported (1,3,4,6–9,19,20), focusing on comparisons of overall excess-death rates associated with the pandemic. We also reviewed epidemiologic evidence for early (herald) waves occurring before September 1918 and for death sparing among elderly persons. We limited the review to studies that provided monthly or weekly historical death data because such data enable identification of herald waves and precise estimation of excess-death rates.
Estimation of Transmission Potential
Transmissibility of an infectious pathogen is measured by the basic reproduction number (R0), which is the average number of secondary infections generated by an infectious person in an entirely susceptible population (21). A related quantity is the reproduction number, R, which can be used for partially immune populations who have been vaccinated or previously exposed to similar pathogens (21).
We estimated R for the 1918 pandemic virus in Boyacá by using a simple method that relies on the epidemic growth rate, a measure of how fast the number of cases increases over time (Technical Appendix [PDF - 186 KB - 4 pages]). Briefly, in the early ascending phase of an epidemic, the daily number of cases (or deaths) should follow an exponential function. By taking the log of daily deaths in the ascending phase, a straight line can be fit to the data. R can be derived from the growth rate estimate r by a simple equation involving the duration of the latency and infectious periods (22) (Technical Appendix [PDF - 186 KB - 4 pages]).
Because of the uncertainty associated with duration of the latency and infectious periods for influenza, we considered periods of 1.5 and 2 days each (23,24). Latency and infectious periods can be summed into a single statistic called the generation interval, which measures the interval between disease onset in 2 successive cases. The generation intervals considered in this study were 3 and 4 days (23,24).
We defined the ascending phase as the period between the day of pandemic onset (defined as the first day of the period of steadily increasing deaths) and the day immediately before the epidemic peak. We tested the robustness of R estimates to the choice of death indicator (deaths from respiratory illness or from all causes). We also compared estimates derived from crude numbers of deaths and excess deaths from respiratory illness that were above background rates.
The same approach and assumptions have been used to quantify Rs associated with the 1918 pandemic in Copenhagen, Denmark, and Mexico City, and hence the Boyacá estimates are directly comparable to estimates from these studies (3,6). For comparison with Boyacá, we also reviewed the literature for published estimates of R associated with the 1918 pandemic in the Americas (2–4).
Timing of Pandemic Waves and Age-specific Patterns of Death
The age-stratified time series of deaths from respiratory illness or all causes in Boyacá indicated that a severe pandemic wave occurred during a 15-week period, October 20, 1918–January 26, 1919 (Figure 2, Figure 3). The profile of age-specific excess deaths from respiratory illness associated with the pandemic period formed a W-shaped pattern; peak mortality rates among infants (<5 years of age) were followed by peak rates among elderly persons (>60 years) and young adults (25–29 years) (Table 1). Excess deaths were lowest among children 5–14 years of age and adults 50–59 years of age. Similar age patterns were found for all-cause deaths (Figure 4); the correlation coefficient between respiratory and all-cause excess-death rates was >0.99 (p<0.01). Excess deaths from respiratory illness captured most influenza-related all-cause excess deaths across all age groups (95% on average, range 81%–100%). Confidence intervals were larger for the most extreme age groups.
To facilitate the comparison between population age groups with different background risks for death, we calculated the risk for excess-death rates relative to baseline rates (Table 1). Although absolute excess-death rates were highest for young children (0–4 years of age) and elderly persons (>60 years), during the pandemic the relative risks were lowest for these age groups. Relative risk was highest for young adults 25–29 years of age; excess-death rates increased 51-fold above background death rates for respiratory causes and 6-fold for all causes.
Comparison of Boyacá and Mexico City shows that age-specific excess-death rates produced a W-shaped pattern for both locations (Figure 5). However, excess-death rates among young adults (20–29 years) were substantially higher for Mexico City than for Boyacá. By contrast, excess-death rates among infants were 2-fold lower for Mexico City than for Boyacá. Excess-death rates for elderly persons were similar for both cities. Overall, we estimate that the October 1918–January 1919 pandemic period was associated with 47 and 40 excess respiratory deaths per 10,000 population in Mexico City and Boyacá, respectively.
A broader comparison of epidemiologic patterns associated with the pandemic at 12 locations on different continents highlights substantial variations in the timing, number of pandemic waves, and age-specific death rates (Table 2). Europe and the United States generally experienced herald waves in during March–August 1918 (except for Paris) and low excess-death rates among elderly populations. In contrast, there was no evidence of death sparing among elderly populations in Latin America or Asia, and herald waves occurred at 4 of the 7 locations studied in these regions. Excess-death rates from respiratory illness were high for Iquitos, Peru; Toluca, Mexico; and Basque Provinces, Spain (121–288 deaths/10,000 population); intermediate in Taiwan (78–180 deaths/10,000 population); and lower elsewhere, including in Boyacá (29–67 deaths/10,000 population).
Reproduction Number Estimates
Table 3 provides summary estimates for the R for the 1918 influenza pandemic in Boyacá, based on growth in daily rates for death from respiratory illness. R was estimated to be 1.4, assuming a short generation interval of 3 days, and 1.5–1.6, assuming a longer interval of 4 days. A sensitivity analysis, based on excess deaths from respiratory illness occurring above a background of expected deaths, generated slightly higher R estimates (1.4–1.5 for a generation interval of 3 days and 1.6–1.7 for a generation interval of 4 days). Different approaches for estimating background deaths resulted in R estimate differences of <0.06 (4%).
Comparison of estimates derived from different locations in the Americas revealed some geographic variations in the transmission potential of the 1918–19 pandemic wave (Table 4). Although R estimates were 1.3–1.8 in most locations in the Americas, assuming a 3-day generation interval, the transmissibility of influenza during the autumn wave might have been particularly high in Toluca, Mexico (estimated R = 2.0–2.5).
Our study makes use of extensive archival death records covering before and during the 1918–19 influenza pandemic in Boyacá, Colombia, and confirms the substantial number of deaths caused by the pandemic in this region. The main epidemiologic features of the pandemic in Boyacá include a single wave of excess deaths during October 1918–January 1919; high excess-death rates among infants and elderly persons; and a moderate R (estimated at 1.4–1.5, assuming a 3-day generation interval).
We did not identify a herald wave of deaths from pandemic influenza in the early part of 1918 in Boyacá. According to epidemiologic data, herald waves of mild pandemic activity have been reported for the spring and summer of 1918 in other regions of the world, including New York City (1), Mexico (3), Lima (4) Geneva (25,26), Copenhagen (6), military camps in the United States (6), the United Kingdom (27), and Singapore (9). The absence of a herald wave in Boyacá could be explained by late introduction of the pandemic influenza virus; alternatively, a mild first wave may have occurred without causing many deaths. Thus, we cannot rule out early pandemic activity, which might have been associated with mild illnesses, before October 1918 in Boyacá. For instance, the summer pandemic wave of 1918 in Denmark was clearly evident only from time-series case data (6). These epidemiologic findings suggesting early pandemic virus activity have recently been confirmed by sequencing of pandemic influenza virus specimens isolated from Army camp populations in the United States as early as May 1918 (28).
Although substantial postpandemic waves have been reported for 1919–20 in New York City (1), Mexico City (3), Lima (4), Japan (5), and Taiwan (8), we could not identify a clear recrudescent pandemic wave in 1920 in Boyacá. A 3-week period in January 1920 and a 4-week period in April–May 1920 were associated with a small increase in deaths from respiratory illness, mostly affecting elderly persons, but we cannot with certainty attribute these deaths to pandemic influenza. Early public health warnings and effective implementation of control interventions in large cities such as New York City, Mexico City, Lima, and Taiwan, could have contributed to maintaining a large pool of susceptible persons, which could fuel subsequent pandemic waves (29). In Japan, postpandemic waves were somewhat limited to regions that escaped earlier waves (5). Given that Boyacá was a relatively small rural area, pandemic activity in 1918 might have proceeded unabated, with no particular interventions, medical or nonmedical. Alternatively, Boyacá could have escaped the recrudescent pandemic wave in 1920 because of its remote location. Overall, the main wave of deaths from pandemic influenza that occurred during October 1918–January 1919 in Boyacá is reminiscent of the single wave of pandemic influenza A (H1N1) 2009 wave that occurred in the Southern Hemisphere during the winter of 2009 (e.g., Chile , Australia , and New Zealand ). Additional data from the 1918 pandemic in other Southern Hemisphere locations are warranted before these findings can be generalized.
The W-shaped age-specific pattern of deaths during the 1918–19 pandemic wave in Boyacá is in agreement with recent reports from the Mexico City area (3) and Peru (4). These reports suggest a lack of death sparing among elderly populations of urban and rural areas of Latin America, although data from additional locations would be useful for generalizing these conclusions. This pattern is also in agreement with anecdotal evidence from aboriginal populations in Alaska in 1918 (32). In contrast to reports for Latin America and Alaska, reports for the United States and Europe suggest that elderly populations were substantially protected from influenza-associated death in 1918 (1,5,6). Previous studies have hypothesized that childhood exposure to influenza A (H1N1) viruses before 1870 might account for prior immunity among elderly persons during the 1918 pandemic. A similar phenomenon has been noted for pandemic (H1N1) 2009, during which risk for clinical infection and death was lower during the pandemic than during seasonal epidemics for persons >60 years of age (13,33).
Regional differences in prior immunity to influenza might result from heterogeneous circulation of influenza viruses during the 19th century, when long-distance travel was much less common than it is today (3). In 1918, Colombia’s population of 5.8 million was heterogeneously distributed and relatively isolated from the rest of the world (34); this isolation could explain the lack of exposure to influenza viruses during the middle of the 19th century. Also in 1918, transportation was underdeveloped in Colombia, consisting mostly of horse- or mule-drawn street cars, waterways, and sparse railroads that did not connect with Boyacá (34). Remoteness could have affected the probability of introduction and of local dissemination of influenza viruses in the Boyacá region. A similar phenomenon could also explain the apparent lack of a herald pandemic wave in the spring of 1918, when pandemic virus activity was not yet globally widespread. Of note, the capital city, Bogota, was the first area in Colombia to report increased influenza activity in October 1918; the virus quickly spread to other Colombia locations (34).
Excess-death rates among young adults were lower in Boyacá than in Mexico City (3). The reasons for this difference are unclear but could be associated with a more sporadic distribution of the population in Boyacá, resulting in lower overall influenza attack rates; however, we do not have epidemiologic evidence to support this assumption. Alternatively, the unidentified factors that made young adults particularly susceptible to influenza-related death in Europe, the United States, and Mexico in 1918 (1,5,6) might have been less common among young adults in Colombia. Despite these geographic differences in absolute risk for death from pandemic influenza, in all locations with sufficient data the relative risk for death consistently peaked among adults 20–29 years of age when compared with baseline death rates during nonpandemic years. Hence, our study confirms the universal atypical severity of this virus in young adults, as previously reported for the United States (1), Mexico (3), Europe (6,7), and Taiwan (8). We also note that data from Boyacá and Mexico City do not support the pessimistic hypothesis that populations lacking prior immunity to the 1918 virus would experience a V-shaped age-associated risk for death, in which risk would rapidly and continuously rise past teenage years (35).
Our excess-deaths approach warrants some caveats. The regression model used to estimate background deaths poorly fit the Boyacá data during the nonpandemic period, probably because of weak seasonality. However, our estimates of excess deaths from pandemic influenza based on deaths from respiratory illness and all causes were highly correlated, similar to those from other temperate countries, where baseline death rates are more seasonal (1,3,6). The sensitivity analysis that we conducted by using an alternative approach to estimate background deaths did not make assumptions about seasonality (20). This analysis produced excess-death estimates highly correlated with those derived from the regression approach (correlation = 0.97; p<0.01; mean difference 4%–7%).
Transmissibility estimates derived from 1918–20 pandemic illness and death data are 1.5–5.4 for community-based settings in several regions of the world (2,6,36,37) (Table 4). Our transmissibility estimates for Boyacá, Colombia, assuming a generation interval of 3 days, are in close agreement with those reported for the wave in autumn in Mexico City (3), Lima (4), England and Wales (27), and Copenhagen (6) and slightly lower than estimates reported for the city of Toluca, Mexico (3), and US cities (2,38). Boyacá’s sparsely distributed population could explain why the estimated disease transmissibility is relatively low. It remains unclear whether differences in reproduction number estimates across locations and pandemic waves reflect true differences attributable to variation in attack rates or local factors affecting transmission or merely illustrate difficulties in measuring this parameter with precision (38). In previous studies focused on reproduction number estimates in which we used similar data and approaches, we have shown that inclusion of a delay between disease onset and death has little effect on the estimates (39).
In conclusion, historical studies from understudied areas are especially helpful for documenting the global death rates and transmission patterns of the 1918 pandemic and for revealing substantial variations among locations. In particular, the lack of death sparing for elderly persons in Colombia and Mexico differs markedly from contemporaneous observations in the United States and Europe. During the 19th century, the Latin American region was relatively isolated (and still is today) (40), which would affect the circulation of historical influenza viruses and baseline population immunity to influenza. We believe that this finding suggests recycling of influenza viruses as the best explanation for death sparing among elderly persons in the United States and Europe in 1918. Preservation and interpretation of archival epidemiologic data are crucial for a better understanding of past pandemics and for better preparedness against future pandemics.
Dr Chowell is an associate professor in the School of Human Evolution and Social Change at Arizona State University and a research fellow at the Fogarty International Center, National Institutes of Health. His research interests include mathematical and statistical modeling of infectious disease transmission and control interventions, with a focus on seasonal and pandemic influenza and quantitative characterization of past influenza pandemics.
We thank Tanya Wilcox for editorial assistance.
This research was conducted in the context of the Multinational Influenza Seasonal Mortality Study, which is an ongoing international collaborative effort for understanding influenza epidemiological and evolutionary patterns, and which is led by the Fogarty International Center, National Institutes of Health (www.origem.info/misms/index.php). Funding for this project comes in part (to L.S.) from the Research and Policy for Infectious Disease Dynamics program of the Science and Technology Directorate, Department of Homeland Security and Fogarty International Center, and from the Office of Global Affairs, International Influenza Unit, in the Office of the Secretary of the Department of Health and Human Services.
- Olson DR, Simonsen L, Edelson PJ, Morse SS. Epidemiological evidence of an early wave of the 1918 influenza pandemic in New York City. Proc Natl Acad Sci U S A. 2005;102:11059–63.
- Mills CE, Robins JM, Lipsitch M. Transmissibility of 1918 pandemic influenza. Nature. 2004;432:904–6.
- Chowell G, Viboud C, Simonsen L, Miller MA, Acuna-Soto R. Mortality patterns associated with the 1918 influenza pandemic in Mexico: evidence for a spring herald wave and lack of preexisting immunity in older populations. J Infect Dis. 2010;202:567–75.
- Chowell G, Viboud C, Simonsen L, Miller MA, Hurtado J, Soto G, The 1918–1920 influenza pandemic in Peru. Vaccine. 2011;29(Suppl 2):B21–6.
- Richard SA, Sugaya N, Simonsen L, Miller MA, Viboud C. A comparative study of the 1918–1920 influenza pandemic in Japan, USA and UK: mortality impact and implications for pandemic planning. Epidemiol Infect. 2009;137:1062–72.
- Andreasen V, Viboud C, Simonsen L. Epidemiologic characterization of the 1918 influenza pandemic summer wave in Copenhagen: implications for pandemic control strategies. J Infect Dis. 2008;197:270–8.
- Erkoreka A. The Spanish influenza pandemic in occidental Europe (1918–1920) and victim age. Influenza Other Respi Viruses. 2010;4:81–9.
- Hsieh YH. Excess deaths and immunoprotection during 1918–1920 influenza pandemic, Taiwan. Emerg Infect Dis. 2009;15:1617–9.
- Lee VJ, Chen MI, Chan SP, Wong CS, Cutter J, Goh KT, Influenza pandemics in Singapore, a tropical, globally connected city. Emerg Infect Dis. 2007;13:1052–7.
- Morens DM, Taubenberger JK, Harvey HA, Memoli MJ. The 1918 influenza pandemic: lessons for 2009 and the future. Crit Care Med. 2010;38(Suppl):e10–20.
- Miller MA, Viboud C, Balinska M, Simonsen L. The signature features of influenza pandemics–implications for policy. N Engl J Med. 2009;360:2595–8.
- Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Pandemic potential of a strain of influenza A (H1N1): early findings. Science. 2009;324:1557–61.
- Chowell G, Bertozzi SM, Colchero MA, Lopez-Gatell H, Alpuche-Aranda C, Hernandez M, Severe respiratory disease concurrent with the circulation of H1N1 influenza. N Engl J Med. 2009;361:674–9.
- Viboud C, Alonso WJ, Simonsen L. Influenza in tropical regions. PLoS Med. 2006;3:e89.
- de Boyacá D. Colombia, Departamento Administrativo Nacional de Estadistica. Censo de 1918. Boyacá (Colombia): The Department; 1918.
- Serfling RE. Methods for current statistical analysis of excess pneumonia–influenza deaths. Public Health Rep. 1963;78:494–506.
- Viboud C, Grais RF, Lafont BA, Miller MA, Simonsen L. Multinational impact of the 1968 Hong Kong influenza pandemic: evidence for a smoldering pandemic. J Infect Dis. 2005;192:233–48.
- Serfling RE, Sherman IL, Houseworth WJ. Excess pneumonia–influenza mortality by age and sex in three major influenza A2 epidemics, United States, 1957–58, 1960 and 1963. Am J Epidemiol. 1967;86:433–41.
- Erkoreka A. Spanish influenza in the heart of Europe. A study of a significant sample of the Basque population. Gesnerus. 2008;65:30–41.
- Murray CJ, Lopez AD, Chin B, Feehan D, Hill KH. Estimation of potential global pandemic influenza mortality on the basis of vital registry data from the 1918–20 pandemic: a quantitative analysis. Lancet. 2006;368:2211–8.
- Anderson RM, May RM. Infectious diseases of humans. Oxford (UK): Oxford University Press; 1991.
- Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci. 2007;274:599–604.
- Ferguson NM, Cummings DA, Fraser C, Cajka JC, Cooley PC, Burke DS. Strategies for mitigating an influenza pandemic. Nature. 2006;442:448–52.
- Cauchemez S, Carrat F, Viboud C, Valleron AJ, Boelle PYA. Bayesian MCMC approach to study transmission of influenza: application to household longitudinal data. Stat Med. 2004;23:3469–87.
- Ammon CE. Spanish flu epidemic in 1918 in Geneva, Switzerland. Euro Surveill. 2002;7:190–2.
- Chowell G, Ammon CE, Hengartner NW, Hyman JM. Estimation of the reproductive number of the Spanish flu epidemic in Geneva, Switzerland. Vaccine. 2006;24:6747–50.
- Chowell G, Bettencourt LM, Johnson N, Alonso WJ, Viboud C. The 1918–1919 influenza pandemic in England and Wales: spatial patterns in transmissibility and mortality impact. Proc Biol Sci. 2008;275:501–9.
- Sheng Z-M, Chertow DS, Ambroggio X, McCall S, Przygodzki RM, Cunningham RE, Autopsy series of 68 cases dying before and during the 1918 influenza pandemic peak. Proc Natl Acad Sci U S A. 2011;108:16416–21.
- Hsieh YH, Chan CH. Excess healthcare burden during 1918–1920 influenza pandemic in Taiwan: implications for post-pandemic preparedness. BMC Public Health. 2011;11:41.
- Ministerio de Salud de Chile. Situacion de influenza A(H1N1). Reporte 01/26/2010 [cited 2010 Feb 15]. http://www.pandemia.cl/index.php?id=61
- ANZIC Investigators, Webb SA, Pettilä V, Seppelt I, Bellomo R, Bailey M, et al. Critical care services and 2009 H1N1 influenza in Australia and New Zealand. N Engl J Med. 2009;361:1925–34.
- Crosby AW. Epidemic and peace, 1918. Westport (CT): Greenwood Press; 1976.
- Vaillant L, La Ruche G, Tarantola A, Barboza P. Epidemiology of fatal cases associated with pandemic H1N1 influenza 2009. Euro Surveill. 2009;14:pii:19309.
- Ospina Díaz JM, Martínez Martín AF, Herrán Falla OF. Impact of the flu pandemic of 1918–1919 on the general mortality profile in Boyacá, Colombia [in Spanish]. Hist Cienc Saude Manguinhos. 2009;16:53–81.
- Palese P. Influenza: old and new threats. Nat Med. 2004;10(Suppl):S82–7.
- Vynnycky E, Trindall A, Mangtani P. Estimates of the reproduction numbers of Spanish influenza using morbidity data. Int J Epidemiol. 2007;36:881–9.
- White LF, Pagano M. Transmissibility of the influenza virus in the 1918 pandemic. PLoS ONE. 2008;3:e1498.
- Chowell G, Nishiura H, Bettencourt LM. Comparative estimation of the reproduction number for pandemic influenza from daily case notification data. J R Soc Interface. 2007;4:155–66.
- Chowell G, Viboud C, Simonsen L, Miller M, Alonso WJ. The reproduction number of seasonal influenza epidemics in Brazil, 1996–2006. Proc Biol Sci. 2010;277:1857–66.
- Russell CA, Jones TC, Barr IG, Cox NJ, Garten RJ, Gregory V, The global circulation of seasonal influenza A (H3N2) viruses. Science. 2008;320:340–6.
Suggested citation for this article: Chowell G, Viboud C, Simonsen L, Miller MA, Acuna-Soto R, Ospina Díaz JM, et al. The 1918–19 influenza pandemic in Boyacá, Colombia. Emerg Infect Dis [serial on the Internet]. 2012 Jan [date cited]. http://dx.doi.org/10.3201/eid1801.101969