Transmission potential of mpox in Mainland China, June-July 2023: estimating reproduction number during the initial phase of the epidemic ========================================================================================================================================= * Andrei R. Akhmetzhanov * Pei-Hsuan Wu ## Abstract Despite reporting very few mpox cases in early 2023, Mainland China observed a surge of over 500 cases during the summer. Amid ambiguous prevention strategies and stigma surrounding mpox transmission, the epidemic silently escalated. This study aims to quantify the scale of the mpox epidemic and assess the transmission dynamics of the virus by estimating the effective reproduction number (*R**e*) during its early phase. Publicly available data were aggregated to obtain daily mpox incidences in Mainland China, and the *R**e* was estimated using an exponential growth model. The mean *R**e* value was found to be 1.76 (95% credible interval: 1.51–2.06), suggesting a case doubling time of approximately 2 weeks. This estimate was compared with *R**e* values from 16 other countries’ national outbreaks in 2022 that had cumulative incidences exceeding 700 symptomatic cases by the end of that year. The *R**e* estimates for these outbreaks ranged from 1.17 for Portugal to 2.88 for Colombia. The pooled mean *R**e* was 1.66 (95% credible interval: 1.40–1.92), which aligns closely with the *R**e* for Mainland China. These findings underscore the need for immediate and effective control measures including targeted vaccination campaigns to mitigate the further spread and impact of the epidemic. Keywords * mpox * global outbreak * reproduction number ## 1. Introduction On 9 August 2023, the China Centers for Disease Control and Prevention (CDC) announced 491 cases in Mainland China in July, representing a marked escalation from 106 cases in June and a five-fold increase (1). Although concerning, this surge was largely anticipated by health experts (2, 3). Amid ambiguous prevention strategies and the stigmatization of specific patient groups—especially men who have sex with men (MSM)—the mpox epidemic intensified unnoticed. Mpox is a viral infection mainly transmitted through direct skin-to-skin contact (4). The global mpox outbreak of 2022–2023 was primarily driven by sexual transmission within MSM networks (5). With limited spillover to the broader general population, the 2022 outbreaks predominantly affected high-risk MSM communities (6, 7), making these communities primary targets for subsequent vaccination campaigns (8, 9). The global transmission of mpox was evident in 2022, but it was not until 2023 that Western Pacific nations, including Japan, South Korea, and Taiwan, saw significant case counts (10-13). By the end of 2022, many of these countries had minimal or no local mpox transmissions, possibly due to ongoing coronavirus disease 2019 (COVID-19)-related travel and social restrictions (14). Given the considerable MSM population in Mainland China, the uncontrollable spread of mpox presents a significant public health concern. An effective mpox vaccine exists, however, challenges such as limited accessibility and disease-related stigma have hampered its broad dissemination (2), underscoring the need to understand mpox transmission dynamics in Mainland China and predict the evolution of future outbreaks. Epidemics typically follow a consistent transmission cycle: introduction, escalation, peak, and eventual decline (15). The initial stages, often marked by exponential growth, offer important insights into understanding of disease spread dynamics. It is crucial to estimate the effective reproduction number, *R**e*, which indicates the average number of secondary transmissions per primary case, during the period in which interventions, behavioral shifts, or acquired immunity have not yet made an impact. Prior *R**e* evaluations for the 2022 national mpox outbreaks have revealed variability (16-19). Using a sub-exponential growth model, Du and colleagues (16) estimated that the *R**e* ranged from 1.02 in Portugal to 1.95 in the United States. Conversely, other reports have noted higher *R**e* values. For example, for the United Kingdom, *R**e* values of 2.32 and 3.14 were reported (17, 18), compared with a value of 1.18 estimated by Du et al. (16). Whereas these studies (17, 18) applied variations of a simple susceptible–exposed–infected–recovered model, they also relied on case counts that were classified by the date of reporting. Such an approach can result in an upward bias in estimates because it may reflect not only actual transmission dynamics but also the recognition of an outbreak and the discovery of infection clusters. Administrative batching in case reporting can further bias the estimation of *R**e*. Nevertheless, a recent UK mpox study (19) employed a detailed, more complex susceptible–exposed– infected–recovered model that factored in potential biases, generating *R**e* values that ranged from 5.16 to 3.02 during the initial stage of the outbreak. Given diverse factors, including cultural differences, differences in outbreak response, and challenges with data collection, identifying how these estimates align with *R**e* in Mainland China is crucial. In this study, we analyzed data from the World Health Organization (WHO) and other governmental sources to enhance our understanding of mpox transmission in Mainland China (20-22). By focusing on case counts solely attributed to Mainland China and estimating the *R**e* during the initial phase of the outbreak, we aimed to compare counts with estimates from 2022 and project potential case counts for the near future. Our ultimate goal was to inform ongoing public health interventions, and potentially optimize efficiency of these interventions, while also contributing to global epidemic preparedness. ## 2. Methods ### Data collection This study used three primary data sources. First, daily incidences of confirmed mpox cases in Mainland China, Taiwan, and the Hong Kong Special Administrative Region were extracted from the WHO global mpox dashboard (2022-2023) as of 15 August 2023 (20). These incidences were categorized by the date of symptom onset, followed by the date of diagnosis and the date of reporting. Second, a dataset of the daily incidences that was confirmed in Taiwan, categorized either by the date of symptom onset or the date of reporting, was retrieved from the National Infectious Disease Statistics System of the Taiwan Centers for Disease Control (CDC) on 23 August 2023 (21). Third, a line list of cases that were confirmed in Hong Kong was compiled using press releases from the Center for Health Protection of the Department of Health of the Government of Hong Kong Special Administrative Region on 23 August 2023 (22). The WHO data recorded 859 total cases in Mainland China, Taiwan, and Hong Kong. These cases were classified as follows: 615 by the date of symptom onset, 225 by the date of reporting, and 19 by the date of diagnosis. The datasets from Taiwan and Hong Kong recorded 292 and 33 cases, respectively; these cases were all classified by the date of symptom onset. Notably, one patient in Hong Kong lacked a definitive date for symptom onset. The report indicated that this patient began showing symptoms at the end of June. The patient was reported in August 2023 and was unlikely to have been accounted in the WHO dataset; therefore, the patient was omitted from our analysis. To extract cases that were confirmed solely in Mainland China, the combined daily counts from the Taiwan and Hong Kong datasets were subtracted from the incidences recorded in the WHO dataset. In the process, two assumptions regarding data decomposition in the WHO dataset were investigated. First, we assumed that all counts from Taiwan were reported by the date of symptom onset in the WHO dataset. However, upon subtracting these counts from the WHO data, the resulting daily incidence yielded predominantly negative values before 1 June 2023. Subsequently, only one date (21 June 2023) appeared with a negative incidence, when the WHO dataset recorded a daily count of two cases, whereas Taiwan and Hong Kong reported three and no cases, respectively. Therefore, we deduced that our initial assumption regarding the integration of Taiwan incidence data into the WHO dataset was likely different. Thus, we hypothesized that all case counts from Taiwan were incorporated into the WHO dataset by their respective dates of reporting to the WHO. Upon examining the WHO reporting pattern, it became evident that most cases that were defined by their reporting dates—specifically May 2, 9, 16, 23, 30; June 6, 20, 27; and July 4—were recorded at equidistant intervals of 7 days. A single deviation from this pattern was noted: the inclusion of June 8 with the exclusion of June 13 (**Supplementary Figure 1**). Assuming that the reporting date in Taiwan had a one-day lag to its confirmation date, the case counts in Taiwan aligned almost perfectly with the counts of cases by reporting date in the WHO dataset (**Supplementary Figure 2**). Although incidences from the Taiwan dataset predominantly fell below the WHO dataset counts, matches were observed for the last five dates: June 8, 20, and 27, and July 4 and 11. This observation led us to the conclusion that the latter approach of imputing data from Taiwan to the WHO dataset was plausible. Moreover, the upper bound of the exponential growth phase was set to *t**up* = 5 july 2023, indicating the peak of observed daily incidence of symptomatic cases at 21. The exponential growth phase was investigated by considering a 45-day period and placing a lower bound at *t**lo*= 21 May 2023. The robustness of the estimates was explored in the sensitivity analysis by varying the dates for both boundaries, which ranged from 16 to 26 May 2023 for the lower bound and 30 June to 10 July 2023 for the upper bound. The duration of the exponential growth phase that was considered in the sensitivity analysis was thus varied between 15 and 55 days. ### Statistical framework #### Reporting delay distribution Given that all cases with unidentified symptom onset dates were back-projected to their estimated date of symptom onset in the subsequent analysis, estimating the reporting delay distribution was crucial. This distribution measured the period from the date of symptom onset to either the date of reporting or the date of diagnosis. For this estimation, we utilized the Hong Kong line list data, which is composed of 33 records, and a cutoff time of *T* = 23 August 2023. Each of the reporting delays (time to diagnosis or time to reporting) was fit to the generalized gamma distribution (GGD), which is defined by three parameters: shape (*Q*), location (*M*), and scale (*S*) .The choice of GGD was guided by its flexibility and its capability to represent three commonly used distributions: gamma (when *Q* = *S*), Weibull (when *Q* = 1), and log-normal (when *Q* = 0) (23). The mathematical inference for the GGD was based on the following: if the random variable follows a GGD, then the transformation where *Q* exp (*Q* − 2*z*), Where *z* = (log(*x*) − *M*)/*S*, adheres to a gamma distribution with shape *Q* − 2 and scale 1. In particular, *O**i* represented the date of symptom onset, and *D**i* (where *O**i*≤ *D**i*) was the date of reporting or diagnosis, extracted from the Hong Kong line list data, with *I*= 1,…,33. We assumed that the respective times of symptom onset (*o**i*) and reporting or diagnosis (*d**i*) were uniformly distributed within their intervals as follows: ![Formula][1] Where the symbol “∼” implies “is distributed as.” The difference *d**i* − *o**i*, representing the observed delay, conforms to the following truncated gamma distribution: ![Formula][2] where: *y**i* = *Q* exp (*Q* −2 *Z**i*), *Y**i* = *Q* exp (*Q*−2*Z**i*), and *z**i* = (log(*d**i* − *o**i*) − *M*)/*S*.*Z**i*=(log(*T*−o*i*)−*M*)/*S*. The likelihood is represented by the form: ![Formula][3] where *λ* and Г are the probability density function and cumulative density function of the gamma distribution, respectively. The parameters are assumed to have weakly-informative priors: log *Q, M*, log *S* ∼ Normal(mean = 0,SD =1). #### Reconstructed incidence Given that not all of the cases were categorized by their date of symptom onset, it was essential to back-project the cases identified solely by their date of diagnosis or reporting to ascertain their dates of symptom onset. Therefore, *n**t* was the number of cases that were symptomatic on day *t*, counting from day 1 on *t**lo*. Similarly,![Graphic][4] and ![Graphic][5] represented the cases with unknown symptom onset dates but with identified reporting and diagnosis dates, respectively. Any non-zero count ![Graphic][6] (where ∗: = {*r,d*}) was then back-projected using sampling from a multinomial distribution as follows: ![Formula][7] In the above, ![Graphic][8] describes the discretized reporting delay distribution as follows: ![Formula][9] The final equation in set (5) indicates the likelihood of a case being reported on day but manifesting symptoms before day 1 (i.e., before *t**lo*). To compute the overall incidence, the following summation was employed: ![Formula][10] Effective reproduction number (*R**e*) during the initial phase of the epidemic Earlier, we defined the exponential growth phase using the date range *t**lo* ≤ *t* ≤ *t**up*. The *R**e* for the initial phase of the epidemic was estimated by analyzing the case counts (*c**t*) with the presumption of the exponential growth. The growth rate was estimated by incorporating the negative binomial likelihood as follows (24): ![Formula][11] where *t**lo* ≤ *t* ≤ *t**up* and the parameters { *i*,*r*, ϕ } were supported by the following weakly-informative priors: ![Formula][12] *R**e* was computed based on Wallinga and Lipsitch’s formula (25), which links *R**e*,*r*, and gamma-distributed generation time with the mean *μ*and SD *σ*: ![Formula][13] For this study, we adopted the values *μ*=12.5 days and *σ*=5.7 days, as estimated in the study by Guzzetta et al. (26). #### Extrapolation of case counts beyond 5 July 2023 Acknowledging the potential for variations in reported cases due to potential delays or modifications in reporting, we extrapolated the case counts from 5 July to 1 September 2023 under the assumption of sustained exponential growth. The case counts ![Graphic][14], which encompass both previously reported and extrapolated cases based on symptom onset date followed the following equation: ![Formula][15] Here, “∼”denotes sampling from the negative binomial distribution. The cumulative counts of cases for July and August were derived by aggregating ![Graphic][16] values across their respective date ranges. For computing case counts based on the date of reporting (*n**t*), the counts ![Graphic][17] were adjusted by considering the reporting delay. For every time, the following sampling was conducted: ![Formula][18] where ![Graphic][19] refers to (5). Thereafter: ![Formula][20] Finally, the aggregated count of the number of cases (*n**t*) for July and August was determined using summation over the respective date ranges. ### Global trends across the 2022 outbreaks Mpox incidence data for 2022 was collated from the WHO dashboard (20). From this dataset, we selected 16 countries that confirmed over 700 symptomatic cases by the end of the year. Most of these nations provided daily incidence counts categorized by the date of symptom onset, the date of diagnosis, or the date of reporting. However, six counties in the WHO Region of Americas (AMRO)—Argentina, Brazil, Chile, Colombia, Mexico, and Peru—presented only aggregated case counts without any subtyping. For each country, we identified the exponential growth phase as a 45-day window that led up to the epidemic’s peak,which was represented by the highest daily case count for 2022. In situations in which multiple peaks were recorded, we selected the earliest occurrence. In our sensitivity analysis, we considered alternating windows of 30 and 60 days. For simplicity, we applied the reporting delay distribution from the Hong Kong data—a reasonable approach given its close alignment to a globally reported median of 6 days, an inter-quartile range spanning 4–9 days (20), and the absence of more detailed global data. Using the methods described earlier, we then estimated the *R**e* for each of the selected countries. ### 2.1. Technical details We adopted the Bayesian framework to fit each realization of the daily incidence classified by the date of symptom onset (6) to models (7)–(9) using Markov Chain Monte Carlo (MCMC) sampling techniques. The Bayesian estimation process was performed using Stan software (27). The inherent structure of the reconstructed incidence was challenging given its reliance on sampling back-projected counts from a multinomial distribution (4). Because these counts, which were represented as integers, were not immediately compatible with the efficient use of Hamilton Monte Carlo sampling in Stan, we employed an approach similar to that previously described (28). For each realization of (4), we simulated the MCMC chain using 1,000 iterations for the tuning-in process and retained only a single posterior for further analysis. To construct the posterior distribution, we used 4,000 different chains that resulted in 4,000 posterior draws. The convergence of the MCMC simulation was inspected visually and was checked implicitly by fitting the model only for symptomatic cases. In doing so, the back-projection was not needed. One thousand iterations that were used for tuning-in and four chains that were composed of 1,000 posterior draws resulted in an R-hat statistic below 1.1 (29). The code is made available on [http://github.com/aakhmetz/Mpox-in-MainlandChina-2023](http://github.com/aakhmetz/Mpox-in-MainlandChina-2023). ## 3. Results Using the Hong Kong data, we estimated a mean delay of 3.4 days (95% credible interval [CrI]: 2.7–4.4) from symptom onset to case diagnosis and a mean delay of 6.0 days (95% CrI: 5.1–7.0) from symptom onset to case reporting. Upon integrating these delays into the case back-projection (**Figure 1**), the *R**e* was estimated at 1.76 (95% CrI: 1.49–2.06). Correspondingly, the doubling time was 15 days (95% CrI: 11–21). Projecting the incidence of mpox past 5 July 2023 to cover the months of July and August, the mean count of cases that were categorized by the symptom onset date was 1,080 (95% CrI: 650–1,720) for July and 5,200 (95% CrI: 1,820–12,060) for August. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F1) Figure 1: The estimation of the case back-projection for Mainland China, integrating reporting delays. The exponential growth phase (dashed lines) represents the lower bound on 21 May 2023 and the upper bound on 5 July 2023. Adjusting these numbers for the reporting delay, the mean numbers of reported cases were 800 (95% CrI: 540–1,180) and 3,850 (95% CrI: 1,490–8,290) for July and August, respectively. The analysis of national outbreaks from selected countries in 2022 revealed a range in *R**e* values. Portugal had the lowest mean *R**e* value at 1.17 (95% CrI: 1.07–1.30), whereas Colombia had the highest value at 2.88 (95% CrI: 2.46–3.34), as illustrated in **Figure 2**. The pooled mean *R**e* was 1.71 (95% CrI: 1.48–1.95), with a between-country variance (heterogeneity) of τ 2=0.19. However, upon evaluating the data quality provided to the WHO by individual countries, it became evident that the nations within the WHO AMRO—with the exceptions of the United States and Canada—provided only aggregated case counts without subtyping them by date of symptom onset, date of diagnosis, or date of reporting. Such reporting likely introduced higher variability and potential bias in the *R**e* values. By excluding these countries from the analysis, the pooled mean *R**e* dropped slightly to 1.66 (95% CrI: 1.40–1.92), and the between-country variance declined to τ2=0.15.Notably, the *R**e*, value for Mainland China remained similar to the pooled average regardless of whether AMRO countries were included in the analysis. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F2) Figure 2: Estimated effective reproduction number (*R**e*) for Mainland China (red) compared with selected national 2022 outbreaks (black) and their pooled mean (black diamond). The ***R******e*** values are shown in the first column on the right, and a period of 45 days is indicated in the last column. Countries that did not report case counts subtyped by the date of symptom onset, diagnosis date, and reporting date were excluded from the calculation of the pooled mean and are shown in dashed black. The sensitivity analysis confirmed the robustness of our findings but highlighted some variability. When considering only cases with a known symptom onset date, the *R**e* for Mainland China remained consistent at 1.75 (95% CrI: 1.48–2.05), suggesting that most cases classified by date of reporting were predominantly linked to Taiwan. Varying the time window for the identified exponential growth phase—in which the lower bound (*t**up*) varied between 16 May 2023 and 26 June 2023, and the upper bound (*t**lo*) varied between 30 June 2023 and 10 July 2023—the median posterior *R**e* ranged from 1.58 to 1.72 (95% CrI 1.38–2.05; **Supplementary Figure 3**). The posterior predictions for reported cases ranged from 650 to 840 (95% CrI: 480–1,400) for July and 2,150 to 3,540 (95% CrI: 1,110–9,690) for August (**Supplementary Figure 4, Supplementary Figure 5**). For the selected 2022 outbreaks, adjusting the length of the exponential growth phase to either 30 or 60 days did not critically alter the prior estimates (**Supplementary Figure 6**). Importantly, when we focused solely on cases with a known date of symptom onset, the estimated *R**e* remained unchanged, indicating the minimal impact of case records with unknown symptom onset on the *R**e* for the selected countries. ## 4. Discussion In our analysis, we estimated that the *R**e* of the mpox epidemic in Mainland China was within the range of 1.6–1.8. This range suggests a doubling time of approximately 2 weeks. Importantly, our estimate was based on the early phase of the epidemic, in which an evident exponential growth of cases was observed and control measures had likely not yet been implemented. Given the relatively low incidence of mpox and no deaths reported during the epidemic’s early phase, significant behavioral changes in at-risk populations were likely not yet in effect, implying that mpox transmission occurred largely unhindered within affected communities during this early phase. The *R**e* value estimated in our study closely aligns with values derived from selected national outbreaks in 2022, suggesting the potentially universal character of the outbreaks. However, an accurate estimation of infections remains challenging given the likely underreporting of mpox cases. Several factors might contribute to this under-ascertainment. One reason may be the societal stigmatization associated with the disease that leads to individuals’ reluctance to disclose their infection for fear of facing societal prejudice or discrimination. Another possibility is that infections can result from intimate encounters outside of a committed relationship, leading individuals to withhold their health status because of the fear of repercussions in their personal lives. Consequently, most cases in the dataset may be either severe or detected through active surveillance measures such as those implemented at entry points into Mainland China. Our study has several limitations. First, we did not differentiate between indigenous and imported infections in our dataset, a factor that may have influenced the estimation of *R**e*. However, by focusing exclusively on the exponential phase in which case growth is most evident, the influence of imported cases on our estimates was likely minimal. Second, the actual method of imputing cases from Taiwan in the WHO dataset may have differed from ours or could have been modified by the WHO over time. For example, the observed pattern of data imputations by date of reporting that were aggregated every 7 days was no longer feasible at the end of July and August. Third, the reporting delay distribution was estimated using Hong-Kong line list data, which may differ from the data from Mainland China given the differences in healthcare systems, case ascertainment, and case investigation protocols. Fourth, the results of comparisons of *R**e* values across different countries must be interpreted with caution. The proactive control measures adopted by some nations versus the more passive strategies employed by others can introduce variances. For instance, the outbreak in Portugal evolved at a slower pace given the government’s proactive position. The case counts were almost steady over 2–3 months before a decline was observed. Although the *R**e* values in our study and that of Du et al. (16) were close to 1, a granular look at the case counts combined with a time-varying *R**e* revealed a brief *R**e* peak of 2.7 on 10 May 2022 that then gradually declined to 1 and dropped below this threshold by the end of June 2022 (30). Our mathematical framework successfully captured this average trend across Portugal and other countries, and a pooled mean *R**e* of 1.66 summarized our overall findings. Finally, the WHO dataset retrieved on 15 August 2023 contained incomplete case counts given that some cases were yet to be reported. At the time of writing the results of this study, a snapshot for 12 September 2023 became available. The *R**e* estimate only slightly changed to 1.63 (95% CrI: 1.43–1.86); however, the case counts for May and June were updated and revealed possible long reporting delays (**Supplementary Figure 1A** and **Supplementary Figure 7**). Although mpox may be associated with a lower mortality rate compared with diseases such as smallpox, its morbidity is notably high, and mpox is associated with health complications such as pneumonia, confusion, and eye infections (31). Mpox disproportionally affects people living with human immunodeficiency virus (HIV) (32). The transmission of mpox throughout Mainland China in the summer of 2023 was largely unconstrained given that regulatory guidelines were first published only at the end of July. This timeline put vulnerable populations at heightened risk of infections and subsequent complications. Nevertheless, the potential impact of efficient vaccination campaigns and information programs aimed at raising awareness of the disease cannot be overlooked; such interventions could significantly reduce case numbers and prevent a larger future outbreak. In contrast to Mainland China, Taiwan managed to stabilize its mpox situation in 2023, maintaining an *R**e* value of approximately one. Drawing from its previous experience—specifically, its management of the hepatitis A virus (HAV) outbreak among MSM in 2015–2017 (33)—the Taiwan government promptly initiated a vaccination campaign, targeted at high-risk groups that was bolstered by awareness campaign. The strategy not only ensured low vaccine hesitancy but also prevented an escalating outbreak. Mainland China could adopt a similar approach with vaccination campaigns focused on high-risk social groups coupled with the dissemination of accurate health messages. Avoiding stigmatization of the disease, an oversight during the early stages of the oubreak (2, 3), could further contribute to higher vaccination rates. ## 5. Conclusions In this study, the mean *R**e* value estimated for the early phase of the mpox outbreak in Mainland China was between 1.6 and 1.8, closely aligning with the pooled mean *R**e* from selected outbreaks in 2022. This alignment suggests that Mainland China could experience a similar trajectory of case growth as observed in other countries previously affected by mpox. These findings underscore the urgency for implementing immediate and effective control measures, including targeted vaccination coupled with information campaigns to reduce stigmatization of vulnerable groups associated with the disease. ## Data Availability All data produced are available online at the github [https://github.com/aakhmetz/Mpox-in-MainlandChina-2023](https://github.com/aakhmetz/Mpox-in-MainlandChina-2023) ## Financial support This study was supported by the National Science and Technology Council, Taiwan (NSTC #111-2314-B-002-289). ## Conflict of interest The authors declare that there are no known competing financial or personal relationships that could have appeared to influence the work reported in this paper. ## Ethical standards This study was based on publicly available data and did not require ethical approval. ## Data availability statement All data used for this study can be found at: [https://github.com/aakhmetz/Mpox-in-MainlandChina-2023/tree/main/data](https://github.com/aakhmetz/Mpox-in-MainlandChina-2023/tree/main/data) ## Figures ![Supplementary Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F3.medium.gif) [Supplementary Figure 1:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F3) Supplementary Figure 1: **Daily incidences reported by the World Health Organization (A), Taiwan Centers for Disease Control (B), and the Government of the Hong Kong Special Administrative Region (C).** The snapshot dates are indicated in the legends. Each case count is subtyped by the date of symptom onset, date of diagnosis, and date of reporting. ![Supplementary Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F4.medium.gif) [Supplementary Figure 2:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F4) Supplementary Figure 2: **The concordance in case counts between the World Health Organization dataset (blue; as of 15 August 2023) and the Taiwan Centers for Disease Control dataset (pink; as of 23 August 2023) that were aggregated by date of reporting**. The dates of reporting for Taiwan were aggregated to the dates of reporting from the WHO dataset. Taiwan cases beyond the latest date (11 July 2023) were omitted. ![Supplementary Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F5.medium.gif) [Supplementary Figure 3:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F5) Supplementary Figure 3: **Effective reproduction number (*R******e*** **) during the initial phase of the epidemic in Mainland China, depending on the time window that identified the exponential phase**. (A) shows the posterior median ***R******e*** value, and (B) and (C) show the lower (2.5th percentile) and upper (97.5th percentile) bounds of the 95% credible interval, respectively. The horizontal axis indicates the lower bound, and the vertical axis indicates the upper bound that were used for the exponential phase window. The date is indicated as a month/day in 2023. ![Supplementary Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F6.medium.gif) [Supplementary Figure 4:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F6) Supplementary Figure 4: Predicted number of reported cases in July, depending on the time window that identified the exponential phase. (A) shows the posterior median of number of cases, and (B) and (C) show the lower (2.5th percentile) and upper (97.5th percentile) bounds of the 95% credible interval, respectively. The horizontal axis indicates the lower bound date, and the vertical axis indicates the upper bound that were used for the exponential phase window. The date is indicated as a month/day in 2023. ![Supplementary Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F7.medium.gif) [Supplementary Figure 5:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F7) Supplementary Figure 5: Predicted number of reported cases in August, depending on the time window that identified the exponential phase. (A) shows the posterior median of number of cases, and (B) and (C) show the lower (2.5th percentile) and upper (97.5th percentile) bounds of the 95% credible interval, respectively. The horizontal axis indicates the lower bound date, and the vertical axis indicates the upper bound that were used for the exponential phase window. The date is indicated as a month/day in 2023. ![Supplementary Figure 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F8.medium.gif) [Supplementary Figure 6:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F8) Supplementary Figure 6: Effective reproduction number (*R**e*) across selected national 2022 mpox outbreaks. The ***R******e*** values are shown in the first column on the right, whereas a period of 30 days (A) or 60 days (B) is indicated in the last column. Countries that did not report case counts subtyped by the date of symptom onset, diagnosis date, or reporting date, were excluded from the calculation of the pooled mean and are shown in dashed black. ![Supplementary Figure 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/10/05/2023.09.23.23296017/F9.medium.gif) [Supplementary Figure 7:](http://medrxiv.org/content/early/2023/10/05/2023.09.23.23296017/F9) Supplementary Figure 7: Daily incidences reported by the World Health Organization with the snapshot date of 12 September 2023. Each case count is subtyped by the date of symptom onset, date of diagnosis, and date of reporting. ## Acknowledgements We thank Anahid Pinchis from Edanz ([www.edanz.com/ac](http://www.edanz.com/ac)) for editing a draft of this manuscript. ## Footnotes * **E-mails:** s8111127{at}gmail.com (P.W.) * Correcting the typographical mistakes with figures' references (Word issue) + updating the abstract * Received September 23, 2023. * Revision received October 5, 2023. * Accepted October 5, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.China Centers for Disease Control. Mpox epidemic monitoring situation in July 2023 [in Chinese]. 2023 [cited 12 September 2023]; Available from: [https://www.chinacdc.cn/jkzt/crb/qt/szkb\_13037/gwjszl\_13092/202308/t20230809\_268502.html](https://www.chinacdc.cn/jkzt/crb/qt/szkb_13037/gwjszl_13092/202308/t20230809_268502.html) 2. 2.Murphy F. Mpox: China’s health authorities fight surge in cases of unknown origin. BMJ. 2023;382:p1939. doi:10.1136/bmj.p1939 [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzODIvYXVnMjJfNy9wMTkzOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzA1LzIwMjMuMDkuMjMuMjMyOTYwMTcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. 3.Yang Z. China is suddenly dealing with another public health crisis: mpox. MIT Technology Review. 2023 [cited 1 August 2023]; Available from: [https://www.technologyreview.com/2023/08/01/1077047/china-public-health-crisis-mpox-outbreak/](https://www.technologyreview.com/2023/08/01/1077047/china-public-health-crisis-mpox-outbreak/) 4. 4.The United States Centers for Disease Control and Prevention. Mpox: how it spreads. 2023 [cited 12 September 2023]; Available from: [https://www.cdc.gov/poxvirus/mpox/if-sick/transmission.html](https://www.cdc.gov/poxvirus/mpox/if-sick/transmission.html) 5. 5.Thornhill JP, Barkati S, Walmsley S, Rockstroh J, Antinori A, Harrison LB, et al. Monkeypox Virus Infection in Humans across 16 Countries — April–June 2022. New England Journal of Medicine. 2022;387(8):679–91. doi:10.1056/NEJMoa2207323 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2207323&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35866746&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F05%2F2023.09.23.23296017.atom) 6. 6.Endo A, Murayama H, Abbott S, Ratnayake R, Pearson CAB, Edmunds WJ, et al. Heavy-tailed sexual contact networks and monkeypox epidemiology in the global outbreak, 2022. Science. 2022 Oct 7;378(6615):90–4. doi:10.1126/science.add4507 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.add4507&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36137054&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F05%2F2023.09.23.23296017.atom) 7. 7.Murayama H, Pearson CAB, Abbott S, Miura F, Jung S-m, Fearon E, et al. Accumulation of Immunity in Heavy-Tailed Sexual Contact Networks Shapes Mpox Outbreak Sizes. The Journal of Infectious Diseases. 2023. doi:10.1093/infdis/jiad254 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiad254&link_type=DOI) 8. 8.Chitwood MH, Kwon J, Savinkina A, Walker J, Bilinski A, Gonsalves G. Estimated Testing, Tracing, and Vaccination Targets for Containment of the US Mpox Outbreak. JAMA Network Open. 2023;6(1):e2250984–e. doi:10.1001/jamanetworkopen.2022.50984 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2022.50984&link_type=DOI) 9. 9.Vairo F, Leone S, Mazzotta V, Piselli P, De Carli G, Lanini S, et al. The possible effect of sociobehavioral factors and public health actions on the mpox epidemic slowdown. International Journal of Infectious Diseases. 2023 2023/05/01/;130:83–5. doi:10.1016/j.ijid.2023.03.005 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2023.03.005&link_type=DOI) 10. 10.Nishiura H, Hayashi K, Anzai A, Kobayashi T, Miyama T, Ueda M, et al. Upcoming major epidemic of mpox in Japan, 2022-23. ResearchGate Preprint. 2023. doi:10.13140/RG.2.2.34892.08327 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.13140/RG.2.2.34892.08327&link_type=DOI) 11. 11.Kim T, Park E, Eun JS, Lee EY, Mun JW, Choi Y, et al. Detecting mpox infection in the early epidemic: an epidemiologic investigation of the third and fourth cases in Korea. Epidemiol Health. 2023;45:e2023040. doi:10.4178/epih.e2023040 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4178/epih.e2023040&link_type=DOI) 12. 12.Chan PC, Lee CC. Response to Mpox transmission - Taiwan, Jan 1 -May 11, 2023. J Formos Med Assoc. 2023 Sep;122(9):819–21. doi:10.1016/j.jfma.2023.06.019 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jfma.2023.06.019&link_type=DOI) 13. 13.Endo A, Jung S-m, Miura F. Mpox emergence in Japan: ongoing risk of establishment in Asia. The Lancet. 2023;401(10392):1923–4. doi:10.1016/S0140-6736(23)00766-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(23)00766-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37269865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F05%2F2023.09.23.23296017.atom) 14. 14.Akhmetzhanov AR, Ponce L, Thompson RN. Emergence potential of monkeypox in the Western Pacific Region, July 2022. International Journal of Infectious Diseases. 2022 2022/09/01/;122:829–31. doi:10.1016/j.ijid.2022.07.044 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2022.07.044&link_type=DOI) 15. 15.Linton NM, Lovell-Read FA, Southall E, Lee H, Akhmetzhanov AR, Thompson RN, Nishiura H. When Do Epidemics End? Scientific Insights from Mathematical Modelling Studies. Centaurus. 2022;64(1):31–60. doi:10.1484/j.Cnt.5.130125 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1484/j.cnt.5.130125)&link_type=DOI) 16. 16.Du Z, Shao Z, Bai Y, Wang L, Herrera-Diestra JL, Fox SJ, et al. Reproduction number of monkeypox in the early stage of the 2022 multi-country outbreak. Journal of Travel Medicine. 2022;29(8). doi:10.1093/jtm/taac099 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taac099&link_type=DOI) 17. 17.Liu T, Yang S, Luo B, Fan X, Zhuang Y, Gao GF, et al. Anticipating the transmissibility of the 2022 mpox outbreak. Journal of Medical Virology. 2023;95(3):e28683. doi:10.1002/jmv.28683 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.28683&link_type=DOI) 18. 18.Yang S, Guo X, Zhao Z, Abudunaibi B, Zhao Y, Rui J, et al. Possibility of mpox viral transmission and control from high-risk to the general population: a modeling study. BMC Infectious Diseases. 2023 2023/02/24;23(1):119. doi:10.1186/s12879-023-08083-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-023-08083-5&link_type=DOI) 19. 19.Brand SPC, Cavallaro M, Cumming F, Turner C, Florence I, Blomquist P, et al. The role of vaccination and public awareness in forecasts of Mpox incidence in the United Kingdom. Nat Commun. 2023 Jul 11;14(1):4100. doi:10.1038/s41467-023-38816-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-023-38816-8&link_type=DOI) 20. 20.World Health Organization. 2022-23 Mpox (Monkeypox) Outbreak: Global Trends. 2023 [cited 23 August 2023]; Available from: [https://worldhealthorg.shinyapps.io/mpx_global/](https://worldhealthorg.shinyapps.io/mpx_global/) 21. 21.Taiwan Centers for Disease Control. Taiwan National Infectious Disease Statistics System. 2023 [cited 23 August 2023]; Available from: [https://nidss.cdc.gov.tw/en/nndss/disease?id=MPXV](https://nidss.cdc.gov.tw/en/nndss/disease?id=MPXV) 22. 22.The Government of the Hong Kong Special Administrative Region. Center for Health Protection, Department of Health. Mpox (also known as monkeypox). 2023 [cited 23 August 2023]; Available from: [https://www.chp.gov.hk/en/miniweb/press_release/100155.html](https://www.chp.gov.hk/en/miniweb/press_release/100155.html) 23. 23.Hoffmann T, Alsing J. Faecal shedding models for SARS-CoV-2 RNA among hospitalised patients and implications for wastewater-based epidemiology. Journal of the Royal Statistical Society Series C: Applied Statistics. 2023;72(2):330–45. doi:10.1093/jrsssc/qlad011 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jrsssc/qlad011&link_type=DOI) 24. 24.Ma J. Estimating epidemic exponential growth rate and basic reproduction number. Infect Dis Model. 2020;5:129–41. doi:10.1016/j.idm.2019.12.009 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.idm.2019.12.009&link_type=DOI) 25. 25.Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci. 2007 Feb 22;274(1609):599–604. doi:10.1098/rspb.2006.3754 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2006.3754&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17476782&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F05%2F2023.09.23.23296017.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000243354200019&link_type=ISI) 26. 26.Guzzetta G, Mammone A, Ferraro F, Caraglia A, Rapiti A, Marziano V, et al. Early Estimates of Monkeypox Incubation Period, Generation Time, and Reproduction Number, Italy, May-June 2022. Emerging infectious diseases. 2022;28(10):2078–81. doi:10.3201/eid2810.221126 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2810.221126&link_type=DOI) 27. 27.Stan Development Team. Stan modeling language users guide and reference manual, 2.33.0. 2023 [cited 23 August 2023]; Available from: [https://mc-stan.org](https://mc-stan.org) 28. 28.Weber S, Gelman A, Lee D, Betancourt M, Vehtari A, Racine-Poon A. Bayesian aggregation of average data: An application in drug development. The Annals of Applied Statistics. 2018;12(3):1583–604, 22. doi:10.1214/17-AOAS1122 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/17-AOAS1122&link_type=DOI) 29. 29.Vehtari A, Gelman A, Simpson D, Carpenter B, Bürkner P-C. Rank-Normalization, Folding, and Localization: An Improved R-hat for Assessing Convergence of MCMC (with Discussion). Bayesian Analysis. 2021;16(2):667–718, 52. doi:10.1214/20-BA1221 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/20-BA1221&link_type=DOI) 30. 30.Borges V, Duque MP, Martins JV, Vasconcelos P, Ferreira R, Sobral D, et al. Viral genetic clustering and transmission dynamics of the 2022 mpox outbreak in Portugal. Nature Medicine. 2023 2023/09/11. doi:10.1038/s41591-023-02542-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-023-02542-x&link_type=DOI) 31. 31.Patel VM, Patel SV. Epidemiological Review on Monkeypox. Cureus. 2023 Feb;15(2):e34653. doi:10.7759/cureus.34653 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7759/cureus.34653&link_type=DOI) 32. 32.Mitjà O, Alemany A, Marks M, Lezama Mora JI, Rodríguez-Aldama JC, Torres Silva MS, et al. Mpox in people with advanced HIV infection: a global case series. The Lancet. 2023;401(10380):939–49. doi:10.1016/S0140-6736(23)00273-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(23)00273-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36828001&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F05%2F2023.09.23.23296017.atom) 33. 33.Lin K-Y, Sun H-Y, Chen Y-H, Lo Y-C, Hsieh S-M, Sheng W-H, et al. Effect of a Hepatitis A Vaccination Campaign During a Hepatitis A Outbreak in Taiwan, 2015– 2017: A Modeling Study. Clinical Infectious Diseases. 2019;70(8):1742–9. doi:10.1093/cid/ciz471 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciz471&link_type=DOI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/inline-graphic-1.gif [5]: /embed/inline-graphic-2.gif [6]: /embed/inline-graphic-3.gif [7]: /embed/graphic-4.gif [8]: /embed/inline-graphic-4.gif [9]: /embed/graphic-5.gif [10]: /embed/graphic-6.gif [11]: /embed/graphic-7.gif [12]: /embed/graphic-8.gif [13]: /embed/graphic-9.gif [14]: /embed/inline-graphic-5.gif [15]: /embed/graphic-10.gif [16]: /embed/inline-graphic-6.gif [17]: /embed/inline-graphic-7.gif [18]: /embed/graphic-11.gif [19]: /embed/inline-graphic-8.gif [20]: /embed/graphic-12.gif