Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 26, Number 8—August 2020
Research Letter

Severe Acute Respiratory Syndrome Coronavirus 2 Transmission Potential, Iran, 2020

Kamalich Muniz-Rodriguez1, Isaac Chun-Hai Fung1Comments to Author , Shayesteh R. Ferdosi, Sylvia K. Ofori, Yiseul Lee, Amna Tariq, and Gerardo ChowellComments to Author 
Author affiliations: Georgia Southern University, Statesboro, Georgia, USA (K. Muniz-Rodriguez, I.C.-H. Fung, S.K. Ofori); The Translational Genomics Research Institute, Phoenix, Arizona, USA (S.R. Ferdosi); Georgia State University School of Public Health, Atlanta, Georgia, USA (Y. Lee, A. Tariq, G. Chowell)

Cite This Article


To determine the transmission potential of severe acute respiratory syndrome coronavirus 2 in Iran in 2020, we estimated the reproduction number as 4.4 (95% CI 3.9–4.9) by using a generalized growth model and 3.5 (95% CI 1.3–8.1) by using epidemic doubling time. The reproduction number decreased to 1.55 after social distancing interventions were implemented.

Since early 2020, Iran has been experiencing a devastating epidemic of coronavirus disease (COVID-19) (1). To determine the transmission potential of severe acute respiratory syndrome coronavirus 2 and thereby guide outbreak response efforts, we calculated basic reproduction numbers (R0). During the early transmission phase, R0 quantifies the average number of secondary cases generated by a primary case in a completely susceptible population, absent interventions or behavioral changes. R0>1 indicates the possibility of sustained transmission; R0<1 implies that transmission chains cannot sustain epidemic growth. As the epidemic continues, the effective reproduction number (Re) offers a time-dependent record of the average number of secondary cases per case as the number of susceptible persons becomes depleted and control interventions take effect. We used 2 methods to quantify the reproduction number by using the curve of reported COVID-19 cases in Iran and its 5 regions (Appendix Table 1). The Georgia Southern University Institutional Review Board made a non–human subjects determination for this project (H20364), under the G8 exemption category.

For method 1, we used a generalized growth model (2) with the growth rate and its scaling factor to characterize the daily reported incidence. Next, we simulated the calibrated generalized growth model by using a discretized probability distribution of the serial interval and assuming a Poisson error structure (Appendix).

We based method 2 on calculation of the epidemic’s doubling times, which correspond to the times when the cumulative incidence doubles and are estimated by using the curve of cumulative daily reported cases. To quantify parameter uncertainty, we used parametric bootstrapping with a Poisson error structure around the number of new reported cases to derive 95% CIs (35). Assuming exponential growth, the epidemic growth rate is equal to ln(2)/doubling time. Assuming that the preinfectious and infectious periods follow an exponential distribution, R0 ≈ (1 + growth rate × serial interval) (Appendix) (6).

For both methods, the serial interval was assumed to follow a gamma distribution; mean (± SD) = 4.41 (± 3.17) days (7; C. You et al., unpub. data, We used MATLAB version R2019b ( and R version 3.6.2 ( for data analyses and creating figures. We determined a priori that α = 0.05.

Using Wikipedia as a starting point, we double-checked the daily reported new cases during February 19–March 19, 2020 (the day before the Iranian New Year) against official Iran press releases and other credible news sources and corrected the data according to official data (Appendix Tables 2, 3, Figure 1). Incident cases for the 5 regions were missing for 2 days (March 2–3), which we excluded from our analysis. Because the reported national number of new cases did not match the sum of new cases reported in Iran’s 5 regions on March 5, we treated each time series as independent and used the data as reported. Using method 1, we estimated R0 data for February 19–March 1, 2020. Using method 2, we estimated R0 from the early transmission phase (February 19–March 1, 2020) and Re based on the growth rate estimated during March 6–19, 2020, when the epidemic slowed, probably reflecting the effect of social distancing.


Thumbnail of Estimates of transmission potential for severe acute respiratory syndrome coronavirus 2 in Iran, 2020. A) Growth rate, r; B) scaling of the growth rate parameter, p; C) mean basic reproduction number, R0; and D) fit of the generalized growth model (method 1) to the Iran data, assuming Poisson error structure as of March 1, 2020. Dashed lines indicate 95% CIs.

Figure. Estimates of transmission potential for severe acute respiratory syndrome coronavirus 2 in Iran, 2020. A) Growth rate, r; B) scaling of the growth rate parameter, p; C) mean basic reproduction number,...

Using method 1, we estimated an R0 of 4.4 (95% CI 3.9–4.9) for COVID-19 in Iran. We estimated a growth rate of 0.65 (95% CI 0.56–0.75) and a scaling parameter of 0.96 (95% CI 0.93–1.00) (Appendix Table 4). The scaling parameter indicated near-exponential epidemic growth (Figure). Using method 2, we found that during February 19–March 1, the cumulative incidence of confirmed cases in Iran had doubled 8 times. The estimated epidemic doubling time was 1.20 (95% CI 1.05–1.45) days, and the corresponding R0 estimate was 3.50 (95% CI 1.28–8.14). During March 6–19, the cumulative incidence of confirmed cases doubled 1 time; doubling time was 5.46 (95% CI 5.29–5.65) days. The corresponding Re estimate was 1.55 (95% CI 1.06–2.57) (Appendix Table 5, Figures 7, 8). Our results are robust and consistent with Iran’s COVID-19 R0 estimates of 4.7 (A. Ahmadi et al., unpub. data, and 4.86 (E. Sahafizadeh, unpub. data, but higher than the R0 of 2.72 estimated by N. Ghaffarzadegan and H. Rahmandad (unpub. data,

Our study has limitations. Our analysis is based on the number of daily reported cases, whereas it would be ideal to analyze case counts by dates of symptoms onset, which were not available. Case counts could be underreported because of underdiagnosis, given subclinical or asymptomatic cases or limited testing capacity to test persons with mild illness. The rapid increase in case counts might represent a belated realization of epidemic severity and rapid catching up and testing many persons with suspected cases. If the reporting ratio remains constant over the study period, and given the near-exponential growth of the epidemic’s trajectory, our estimates would remain reliable. Although data are not stratified according to imported and local cases, we assumed that persons were infected locally because transmission has probably been ongoing in Iran for some time (8).

Although the COVID-19 epidemic in Iran has slowed substantially, the situation remains dire. Tighter social distancing interventions are needed to bring this epidemic under control.

Ms. Muniz-Rodriguez is a doctoral student in epidemiology and Dr. Fung is an associate professor of epidemiology at Jiann-Ping Hsu College of Public Health, Georgia Southern University. Their research interests include infectious disease epidemiology, digital health, and disaster emergency responses.



G.C. received support from NSF grant 1414374 as part of the joint NSF-NIH-USDA Ecology and Evolution of Infectious Diseases program. I.C.-H.F. received salary support from the Centers for Disease Control and Prevention (CDC) (19IPA1908208). This article is not part of I.C.-H.F’s CDC-sponsored projects.



  1. Wood  G. Coronavirus could break Iranian society [cited 2020 Feb 29].
  2. Viboud  C, Simonsen  L, Chowell  G. A generalized-growth model to characterize the early ascending phase of infectious disease outbreaks. Epidemics. 2016;15:2737. DOIPubMed
  3. Banks  HT, Hu  S, Thompson  WC. Modeling and inverse problems in the presence of uncertainty: CRC Press; 2014.
  4. Chowell  G, Ammon  CE, Hengartner  NW, Hyman  JM. Transmission dynamics of the great influenza pandemic of 1918 in Geneva, Switzerland: Assessing the effects of hypothetical interventions. J Theor Biol. 2006;241:193204. DOIPubMed
  5. Chowell  G, Shim  E, Brauer  F, Diaz-Dueñas  P, Hyman  JM, Castillo-Chavez  C. Modelling the transmission dynamics of acute haemorrhagic conjunctivitis: application to the 2003 outbreak in Mexico. Stat Med. 2006;25:184057. DOIPubMed
  6. Vynnycky  E, White  RG. An introduction to infectious disease modelling. Oxford (UK): Oxford University Press; 2010.
  7. Nishiura  H, Linton  NM, Akhmetzhanov  AR. Serial interval of novel coronavirus (COVID-19) infections. Int J Infect Dis. 2020;93:2846. DOIPubMed
  8. Tuite  AR, Bogoch  II, Sherbo  R, Watts  A, Fisman  D, Khan  K. Estimation of coronavirus disease 2019 (COVID-19) burden and potential for international dissemination of infection from Iran. Ann Intern Med. 2020.




Cite This Article

DOI: 10.3201/eid2608.200536

Original Publication Date: April 22, 2020

1These first authors contributed equally to this article.

Table of Contents – Volume 26, Number 8—August 2020


Please use the form below to submit correspondence to the authors or contact them at the following address:

Isaac Chun-Hai Fung, Department of Biostatistics, Epidemiology and Environmental Health Sciences, Jiann-Ping Hsu College of Public Health, Georgia Southern University, 
PO Box 7989, Statesboro, GA 30460-7989 USA; ; or Gerardo Chowell, Department of Population Health Sciences, School of Public Health, Georgia State University, Suite 662, Office 640B, Atlanta, GA 30303, USA

Send To

character(s) remaining.

Comment submitted successfully, thank you for your feedback.


Page created: April 22, 2020
Page updated: July 19, 2020
Page reviewed: July 19, 2020
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.