Application of the Logistic Model to the COVID-19 Pandemic in South Africa and the United States: Correlations and Predictions ============================================================================================================================== * David H. Roberts ## ABSTRACT We apply the simple logistic model to the four waves of COVID-19 taking place in South Africa over the period 2020 January 1 through 2022 January 11. We show that this model provides an excellent fit to the time history of three of the four waves. We then derive a theoretical correlation between the growth rate of each wave and its duration, and demonstrate that it is well obeyed by the South Africa data. We then turn to the data for the United States. As shown by Roberts (2020a, 2020b), the basic logistic model provides only a marginal fit to the early data. Here we break the data into six “waves,” and treat each one separately. For four of the six the logistic model is useful, and we present full results. We then ask if these data provide a way to predict the length of the ongoing Omicron wave in the US (commonly called “wave 4,” but the sixth wave as we have broken the data up). Comparison of these data to those from South Arica, and internal comparison of the US data, *suggest* that this last wave will die out by about 2022-January-20. ## 1. Introduction Since the spring of 2020 a pandemic of infection of a novel coronavirus (SARS-CoV-2) has overspread the world. Epidemiologists have struggled to describe adequately this pandemic and to predict its future course. This is a complex undertaking, involving the mathematics of epidemiology and a huge amount of uncertain data that serve as input to the models. In this paper we examine the ability of a simple epidemiological model, the *logistic model*, to describe the course of the pandemic in the Republic of South Africa and the United States. After showing that it is adequate to the task for South Africa, we derive a theoretical correlation between the growth rate of each wave of the pandemic and its duration, and show that it is well obeyed for the South Africa data. ## 2. The Logistic Model ### 2.1. Derivation The exposition in this section is from Roberts (2020). A simple model for the evolution of an pandemic is based on the logistic differential equation. This describes an pandemic that begins with a small number *f* of infected individuals, and subsequently spreads through a population. The motivation for this model is as follows. If the population of infected individuals as a function of time is *f*(*t*), simple exponential growth with growth rate *r* is determined by the differential equation ![Formula][1] with the growing exponential solution ![Formula][2] where *f* = *f*(0). This is what happens with an infinite pool of subjects. However, for a finite pool of subjects, as the population of infected individuals grows the number of subjects available to be infected gets smaller. This is taken into account by modifying the exponential differential equation to become ![Formula][3] where *K* is the total available pool of individuals. The solution of this equation2 is ![Formula][4] which satisfies the required limits *f*(0) = *f* and *f*(∞) = *K*. The time course described by Equation 2 is the familiar “S curve” used to describe bacterial growth and other phenomena (see Figure 3). An analysis of the total number of cases as a function of time *f*(*t*) is just one way to compare model and data. Instead we can examine the number of new cases per day as a function of time; this is the time derivative of *f*(*t*).3 This is easily found by substituting Equation 2 into Equation 1, ![Formula][5] In either case, for the logistic model the three parameters to be adjusted are *f*, *K*, and *r*. ### 2.2. COVID data for the United States and the Republic of South Africa ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F1) Fig. 1. Cases data for the entire pandemic for the United States. ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F2) Fig. 2. Cases data for the entire pandemic for the Republic of South Africa.. ### 2.3. Fit to the COVID-19 Pandemic in the Republic of South Africa When applied to the ongoing COVID-19 pandemic in the United States the logistic model does only a fair job of accounting for the actual history of the total number of infected individuals as a function of time in the first wave of COVID-19, which took place in the first half of 2020 (see Roberts 2020a). It fails even more spectacularly for later waves (Roberts 2020b). ### 2.4. RSA Wave 1 ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F3) Fig. 3. Fits of the logistic model to RSA cases of COVID-19 in wave 1. ### 2.5. RSA Wave 2 ![Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F4.medium.gif) [Fig. 4.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F4) Fig. 4. Fits of the logistic model to RSA cases of COVID-19 in wave 2. ### 2.6. RSA Wave 4 ![Fig. 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F5.medium.gif) [Fig. 5.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F5) Fig. 5. Fits of the logistic model to RSA cases of COVID-19 in wave 4. ## 3. Correlation of the Growth Rate and Duration for the COVID-19 Waves in the Republic of South Africa It is natural that the growth rate and the duration of the waves of a pandemic should be correlated, higher growth rates leading to shorter durations. It can be shown (see the appendix) that for a logistic model the daily cases distribution has a full-with-half-maximum given by ![Formula][6] Note that this result is independent of *f* and *K*. The data and the prediction of Eq.4 are compared in Figure 6, where we see that the agreement is excellent. ![Fig. 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F6.medium.gif) [Fig. 6.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F6) Fig. 6. Correlation of the growth rate *r* and the FWHM of each wave for the COVID-19 pandemic in the Republic of South Africa. The curve is Eq. 4. The two pairs of points at the ends are trom the LSQ and L1 fits to waves 1 and 4. ## 4. Six Waves in the United States ### 4.1. USA Wave 1 ![Fig. 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F7.medium.gif) [Fig. 7.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F7) Fig. 7. L1 fits of the logistic model to USA cases of COVID-19 in wave 1. ### 4.2. USA Wave 3 ![Fig. 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F8.medium.gif) [Fig. 8.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F8) Fig. 8. LSQ fits of the logistic model to total USA cases of COVID-19 in wave 3. ### 4.3. USA Wave 5 ![Fig. 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F9.medium.gif) [Fig. 9.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F9) Fig. 9. L1 fits of the logistic model to USA cases of COVID-19 in wave 5 ### 4.4. USA Wave 6 ![Fig. 10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F10.medium.gif) [Fig. 10.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F10) Fig. 10. L1 fits of the logistic model to total USA cases of COVID-19 in wave 6 ![Fig. 11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F11.medium.gif) [Fig. 11.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F11) Fig. 11. LSQ fits of the logistic model to USA cases of COVID-19 in wave 6 ## 5. Predictions for United States Wave 6 Fig.. 12 Shows the correlation between growth and duration of the various waves in the United States. Because the various waves tend to overlap, for all but wave 6 we can find only lower limits to the widths. The curve is Eq. 4, identical to the one in Fig. 6. ![Fig. 12.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F12.medium.gif) [Fig. 12.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F12) Fig. 12. Correlation of the growth rate *r* and the FWHM of each wave for the COVID-19 pandemic in the United states. The curve is Eq. 4. The black point is the measurement for wave 6 using the LSQ solution, the red points are lower limits for the widths of the other waves. In Figs. 13 & 14 we show the predictions for wave 6 that follow from the L1 and LSQ solutions. The latter provides a better fit to the data (see Figs. 10 & 11), so we prefer the predictions made from those solutions. Thus we expect the Omicron wave in the US to be dramatically reduced within about ten days from now, roughly 22 January 2022. Needless to say, this would be very welcome. ![Fig. 13.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F13.medium.gif) [Fig. 13.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F13) Fig. 13. Predictions of the logistic model for the total USA cases of COVID-19 in wave 6. ![Fig. 14.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/13/2022.01.12.22269193/F14.medium.gif) [Fig. 14.](http://medrxiv.org/content/early/2022/01/13/2022.01.12.22269193/F14) Fig. 14. Predictions of the logistic model for the daily USA cases of COVID-19 in wave 6. ## 6. Conclusions In this paper we have used the logistic model for epidemics to describe the COVID-19 outbreaks in the Republic of South Africa and in the United States. We find a universal analytic relationship between the growth rates and the durations of each wave, and this is closely followed by the (very clean) data from South Africa. When applied to the United States, the data are messier, but a tentative prediction is possible for the expected duration of the current Omicron wave in the US – it should be mostly over by about ten days from now (roughly 22-January-2022). ## Data Availability WHO COVID-19 Website [https://covid19.who.int/WHO-COVID-19-global-data.csv](https://covid19.who.int/WHO-COVID-19-global-data.csv) ## 7. Acknowledgements We thank Brian Boyle, Mary Roberts, and Bob Sauer for their helpful comments and suggestions. ## 9. Appendix ### Derivation of Equation 4 The differential distribution for the logistic function is4 ![Formula][7] and is derivative is ![Formula][8] This is zero at ![Formula][9] and *df /dt* at this maximum is ![Formula][10] Setting *df /dt* equal to half of this value and solving for the locations yields a complicated expression that reduces to ![Formula][11] Note that neither *f* nor *K* enters this expression. ## Footnotes * 2 This is a version of the Bernoulli differential equation *f**′*= *f*(1 − *f*). * 3 Fitting the daily numbers is statistically the preferred procedure as the data points are independent, unlike those for daily totals. * 4 All the calculations were done with Mathematica * Received January 12, 2022. * Revision received January 12, 2022. * Accepted January 13, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder has placed this preprint in the Public Domain. It is no longer restricted by copyright. Anyone can legally share, reuse, remix, or adapt this material for any purpose without crediting the original authors. ## 8. References 1. Roberts, D. H. (2020a), “Two New Models for Epidemics with Application to the COVID-19 Pandemic in the United States, Italy, and the United Kingdom,” MedRxiv:2020.07.13.20152686v1 2. Roberts, D. H. (2020b), “A New Adaptive Logistic Model for Epidemics and the Resurgence of COVID-19 in the United States,” MedRxiv:2020.07.17.20156109v1 3. World Health Organization (2022), world-wide COVID histories, version of 11 January 2022. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-11.gif [7]: /embed/graphic-21.gif [8]: /embed/graphic-22.gif [9]: /embed/graphic-23.gif [10]: /embed/graphic-24.gif [11]: /embed/graphic-25.gif