Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Mathematical prediction of the time evolution of the COVID-19 pandemic in Italy by a Gauss error function and Monte Carlo simulations</title>
<author><name sortKey="Ciufolini, Ignazio" sort="Ciufolini, Ignazio" uniqKey="Ciufolini I" first="Ignazio" last="Ciufolini">Ignazio Ciufolini</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="GRID">grid.9906.6</institution-id>
<institution-id institution-id-type="ISNI">0000 0001 2289 7785</institution-id>
<institution>Dipartimento di Ingegneria dell’Innovazione,</institution>
<institution>University of Salento, Lecce, and Centro Fermi,</institution>
</institution-wrap>
Rome, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Paolozzi, Antonio" sort="Paolozzi, Antonio" uniqKey="Paolozzi A" first="Antonio" last="Paolozzi">Antonio Paolozzi</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.7841.a</institution-id>
<institution>Scuola di Ingegneria Aerospaziale, Sapienza,</institution>
<institution>University of Rome,</institution>
</institution-wrap>
Rome, Italy</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">32309108</idno>
<idno type="pmc">7156796</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC7156796</idno>
<idno type="RBID">PMC:7156796</idno>
<idno type="doi">10.1140/epjp/s13360-020-00383-y</idno>
<date when="2020">2020</date>
<idno type="wicri:Area/Pmc/Corpus">000115</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000115</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Mathematical prediction of the time evolution of the COVID-19 pandemic in Italy by a Gauss error function and Monte Carlo simulations</title>
<author><name sortKey="Ciufolini, Ignazio" sort="Ciufolini, Ignazio" uniqKey="Ciufolini I" first="Ignazio" last="Ciufolini">Ignazio Ciufolini</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="GRID">grid.9906.6</institution-id>
<institution-id institution-id-type="ISNI">0000 0001 2289 7785</institution-id>
<institution>Dipartimento di Ingegneria dell’Innovazione,</institution>
<institution>University of Salento, Lecce, and Centro Fermi,</institution>
</institution-wrap>
Rome, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Paolozzi, Antonio" sort="Paolozzi, Antonio" uniqKey="Paolozzi A" first="Antonio" last="Paolozzi">Antonio Paolozzi</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.7841.a</institution-id>
<institution>Scuola di Ingegneria Aerospaziale, Sapienza,</institution>
<institution>University of Rome,</institution>
</institution-wrap>
Rome, Italy</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">European Physical Journal plus</title>
<idno type="eISSN">2190-5444</idno>
<imprint><date when="2020">2020</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p id="Par1">In this paper are presented mathematical predictions on the evolution in time of the number of positive cases in Italy of the COVID-19 pandemic based on official data and on the use of a function of the type of a Gauss error function, with four parameters, as a cumulative distribution function. We have analyzed the available data for China and Italy. The evolution in time of the number of cumulative diagnosed positive cases of COVID-19 in China very well approximates a distribution of the type of the error function, that is, the integral of a normal, Gaussian distribution. We have then used such a function to study the potential evolution in time of the number of positive cases in Italy by performing a number of fits of the official data so far available. We then found a statistical prediction for the day in which the peak of the number of daily positive cases in Italy occurs, corresponding to the flex of the fit, that is, to the change in sign of its second derivative (i.e., the change from acceleration to deceleration), as well as of the day in which a substantial attenuation of such number of daily cases is reached. We have also analyzed the predictions of the cumulative number of fatalities in both China and Italy, obtaining consistent results. We have then performed 150 Monte Carlo simulations to have a more robust prediction of the day of the above-mentioned peak and of the day of the substantial decrease in the number of daily positive cases and fatalities. Although official data have been used, those predictions are obtained with a heuristic approach since they are based on a statistical approach and do not take into account either a number of relevant issues (such as number of daily nasopharyngeal swabs, medical, social distancing, virological and epidemiological) or models of contamination diffusion.
</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Li, R" uniqKey="Li R">R Li</name>
</author>
<author><name sortKey="Pei, S" uniqKey="Pei S">S Pei</name>
</author>
<author><name sortKey="Chen, B" uniqKey="Chen B">B Chen</name>
</author>
<author><name sortKey="Song, Y" uniqKey="Song Y">Y Song</name>
</author>
<author><name sortKey="Zhang, T" uniqKey="Zhang T">T Zhang</name>
</author>
<author><name sortKey="Yang, W" uniqKey="Yang W">W Yang</name>
</author>
<author><name sortKey="Shaman, J" uniqKey="Shaman J">J Shaman</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Lessler, J" uniqKey="Lessler J">J Lessler</name>
</author>
<author><name sortKey="Reich, Ng" uniqKey="Reich N">NG Reich</name>
</author>
<author><name sortKey="Brookmeyer, R" uniqKey="Brookmeyer R">R Brookmeyer</name>
</author>
<author><name sortKey="Perl, Tm" uniqKey="Perl T">TM Perl</name>
</author>
<author><name sortKey="Nelson, Ke" uniqKey="Nelson K">KE Nelson</name>
</author>
<author><name sortKey="Cummings, Da" uniqKey="Cummings D">DA Cummings</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Press, Wh" uniqKey="Press W">WH Press</name>
</author>
<author><name sortKey="Flannery, Bp" uniqKey="Flannery B">BP Flannery</name>
</author>
<author><name sortKey="Teukolsky, Sa" uniqKey="Teukolsky S">SA Teukolsky</name>
</author>
<author><name sortKey="Vetterling, Wt" uniqKey="Vetterling W">WT Vetterling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bevington, Pr" uniqKey="Bevington P">PR Bevington</name>
</author>
<author><name sortKey="Robinson, D Keith" uniqKey="Robinson D">D Keith Robinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ciufolini, I" uniqKey="Ciufolini I">I Ciufolini</name>
</author>
<author><name sortKey="Monge, Bm" uniqKey="Monge B">BM Monge</name>
</author>
<author><name sortKey="Paolozzi, A" uniqKey="Paolozzi A">A Paolozzi</name>
</author>
<author><name sortKey="Konig, R" uniqKey="Konig R">R König</name>
</author>
<author><name sortKey="Sindoni, G" uniqKey="Sindoni G">G Sindoni</name>
</author>
<author><name sortKey="Michalak, G" uniqKey="Michalak G">G Michalak</name>
</author>
<author><name sortKey="Pavlis, Ec" uniqKey="Pavlis E">EC Pavlis</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Eur Phys J Plus</journal-id>
<journal-id journal-id-type="iso-abbrev">Eur Phys J Plus</journal-id>
<journal-title-group><journal-title>European Physical Journal plus</journal-title>
</journal-title-group>
<issn pub-type="epub">2190-5444</issn>
<publisher><publisher-name>Springer Berlin Heidelberg</publisher-name>
<publisher-loc>Berlin/Heidelberg</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">32309108</article-id>
<article-id pub-id-type="pmc">7156796</article-id>
<article-id pub-id-type="publisher-id">383</article-id>
<article-id pub-id-type="doi">10.1140/epjp/s13360-020-00383-y</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Regular Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Mathematical prediction of the time evolution of the COVID-19 pandemic in Italy by a Gauss error function and Monte Carlo simulations</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><name><surname>Ciufolini</surname>
<given-names>Ignazio</given-names>
</name>
<address><email>ignazio.ciufolini@gmail.com</email>
</address>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Paolozzi</surname>
<given-names>Antonio</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.9906.6</institution-id>
<institution-id institution-id-type="ISNI">0000 0001 2289 7785</institution-id>
<institution>Dipartimento di Ingegneria dell’Innovazione,</institution>
<institution>University of Salento, Lecce, and Centro Fermi,</institution>
</institution-wrap>
Rome, Italy</aff>
<aff id="Aff2"><label>2</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.7841.a</institution-id>
<institution>Scuola di Ingegneria Aerospaziale, Sapienza,</institution>
<institution>University of Rome,</institution>
</institution-wrap>
Rome, Italy</aff>
</contrib-group>
<pub-date pub-type="epub"><day>15</day>
<month>4</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>15</day>
<month>4</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="ppub"><year>2020</year>
</pub-date>
<volume>135</volume>
<issue>4</issue>
<elocation-id>355</elocation-id>
<history><date date-type="received"><day>30</day>
<month>3</month>
<year>2020</year>
</date>
<date date-type="accepted"><day>3</day>
<month>4</month>
<year>2020</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2020</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
.</license-p>
</license>
</permissions>
<abstract id="Abs1"><p id="Par1">In this paper are presented mathematical predictions on the evolution in time of the number of positive cases in Italy of the COVID-19 pandemic based on official data and on the use of a function of the type of a Gauss error function, with four parameters, as a cumulative distribution function. We have analyzed the available data for China and Italy. The evolution in time of the number of cumulative diagnosed positive cases of COVID-19 in China very well approximates a distribution of the type of the error function, that is, the integral of a normal, Gaussian distribution. We have then used such a function to study the potential evolution in time of the number of positive cases in Italy by performing a number of fits of the official data so far available. We then found a statistical prediction for the day in which the peak of the number of daily positive cases in Italy occurs, corresponding to the flex of the fit, that is, to the change in sign of its second derivative (i.e., the change from acceleration to deceleration), as well as of the day in which a substantial attenuation of such number of daily cases is reached. We have also analyzed the predictions of the cumulative number of fatalities in both China and Italy, obtaining consistent results. We have then performed 150 Monte Carlo simulations to have a more robust prediction of the day of the above-mentioned peak and of the day of the substantial decrease in the number of daily positive cases and fatalities. Although official data have been used, those predictions are obtained with a heuristic approach since they are based on a statistical approach and do not take into account either a number of relevant issues (such as number of daily nasopharyngeal swabs, medical, social distancing, virological and epidemiological) or models of contamination diffusion.
</p>
</abstract>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© Società Italiana di Fisica and Springer-Verlag GmbH Germany, part of Springer Nature 2020</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1"><title>Introduction</title>
<p id="Par2">By considering the cumulative diagnosed positive cases of COVID-19 infections and fatalities available in the Web site of the Italian “Ministero della Salute” [<xref ref-type="bibr" rid="CR1">1</xref>
], Worldometer [<xref ref-type="bibr" rid="CR2">2</xref>
] and World Health Organization [<xref ref-type="bibr" rid="CR3">3</xref>
], we found that they can be well approximated by a cumulative distribution function (CDF) with four parameters of the type of the Gauss error function, that is, the integral of a normal, Gaussian distribution (see Sect. <xref rid="Sec2" ref-type="sec">2</xref>
). By positive cases, we mean the positive cases actually diagnosed plus, for future days, the positive cases that we expect to be diagnosed. Indeed, it is well known among the virologists that the actual number of positive cases is much higher than the diagnosed ones [<xref ref-type="bibr" rid="CR4">4</xref>
]. However, it is assumed that the diagnosed cases are a good statistical representation of the entire population of the positive cases. In Fig. <xref rid="Fig1" ref-type="fig">1</xref>
, we report the result of our fit of the cumulative diagnosed positive cases in China. We have then applied such a CDF to study the evolution in time of the number of positive cases in Italy, in the attempt to possibly, statistically, predict the peak in the number of daily positive cases and the possible date of a substantial decrease in the number of daily positive cases. Furthermore, we have also applied our method to the CDF of the number of fatalities in both China and Italy, confirming our predictions obtained with the cumulative positive cases. Finally, in Sect. <xref rid="Sec4" ref-type="sec">4</xref>
, we have performed a number of Monte Carlo simulations to get a more robust prediction of both the date of the peak in the number of positive cases, diagnosed each day, and in the date after which the number of new positive cases will be below a certain threshold [<xref ref-type="bibr" rid="CR5">5</xref>
].</p>
</sec>
<sec id="Sec2"><title>Fit of cumulative diagnosed positive cases and fatalities of COVID-19 in China</title>
<p id="Par3">Based on the number of diagnosed positive cases and fatalities of COVID-19 in China, we have fitted the cumulative numbers with a function of the type of the Gauss error function<disp-formula id="Equ1"><label>1</label>
<alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} a+b\,\, \text{ erf }\left( {cx-d}\right) \end{aligned}$$\end{document}</tex-math>
<mml:math id="M2" display="block"><mml:mrow><mml:mtable><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:mi>a</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>b</mml:mi>
<mml:mspace width="0.166667em"></mml:mspace>
<mml:mspace width="0.166667em"></mml:mspace>
<mml:mspace width="0.333333em"></mml:mspace>
<mml:mtext>erf</mml:mtext>
<mml:mspace width="0.333333em"></mml:mspace>
<mml:mfenced close=")" open="("><mml:mrow><mml:mi>c</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<graphic xlink:href="13360_2020_383_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
containing the four parameters <italic>a</italic>
, <italic>b</italic>
, <italic>c</italic>
, <italic>d</italic>
, that we have fitted using the available official data. A similar distribution is also observed in other studies of seasonal influenza and pandemic [<xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR7">7</xref>
]. The result of the fit is reported in Fig. <xref rid="Fig1" ref-type="fig">1</xref>
, which shows the good level of the fit using function (<xref rid="Equ1" ref-type="">1</xref>
) with those four parameters. In Fig. <xref rid="Fig2" ref-type="fig">2</xref>
, we report the fit of the cumulative number of fatalities due to COVID-19 in China, which also shows the good level of the fit using function (<xref rid="Equ1" ref-type="">1</xref>
). In Sect. <xref rid="Sec3" ref-type="sec">3</xref>
, the same procedure was applied to the Italian data.<fig id="Fig1"><label>Fig. 1</label>
<caption><p>Fit of the cumulative number of diagnosed positive cases of COVID-19 in China (red dots) from January 22, 2020 (included), to March 27, 2020 (included), and the fitting function of the type of a Gauss error function with four parameters (solid line). The horizontal axis reports the days from January 22, 2020; the vertical axis reports the cumulative number of diagnosed positive cases</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig1_HTML" id="MO9"></graphic>
</fig>
<fig id="Fig2"><label>Fig. 2</label>
<caption><p>Fit of the cumulative number of fatalities due to COVID-19 in China (red dots) from January 22, 2020 (included), to March 28, 2020 (included), and the fitting function of the type of a Gauss error function with four parameters (solid line). The horizontal axis reports the days from January 22, 2020; the vertical axis reports the cumulative number of fatalities in China</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig2_HTML" id="MO10"></graphic>
</fig>
</p>
</sec>
<sec id="Sec3"><title>Fit and predictions of cumulative positive cases and fatalities of COVID-19 in Italy</title>
<p id="Par4">In this section, we report the results of the fit of the cumulative diagnosed positive cases of COVID-19 in Italy using a function of the type of the Gauss error function, given in Sect. <xref rid="Sec2" ref-type="sec">2</xref>
. We obtained very similar results using a logistic function with four parameters, which are not reported here for the sake of brevity. Figure <xref rid="Fig3" ref-type="fig">3</xref>
shows the fit of the data from February 15, 2020, to March 29, 2020. According to this fit, the flex, i.e., the point where the second derivative of the fit is becoming negative, that is, the difference between two successive daily cases becomes negative, or in other words the point where there is a deceleration in the number of positive cases is reached at March 25, 2020, plus or minus 2 days (2-sigma), this uncertainty is derived in Sect. <xref rid="Sec4" ref-type="sec">4</xref>
. According to this fit, the date of a substantial reduction in the number of cumulative positive cases in Italy (about 100 cases), is April 22, 2020, plus or minus 4.6 days (2-sigma) (see Sect. <xref rid="Sec4" ref-type="sec">4</xref>
).</p>
<p id="Par5">To check whether we obtain the same predictions using the cumulative number of fatalities of COVID-19 in Italy, instead of using the number of positive cases of COVID-19 in Italy, we repeated the previous analysis using the cumulative fatalities enumerated in the Web sites of the Italian “Ministero della Salute” [<xref ref-type="bibr" rid="CR1">1</xref>
], Worldometer [<xref ref-type="bibr" rid="CR2">2</xref>
] and World Health Organization [<xref ref-type="bibr" rid="CR3">3</xref>
]. The results are shown in Fig. <xref rid="Fig4" ref-type="fig">4</xref>
. The analysis of the curve shows that the flex and the number of a substantial reduction in the number of fatalities (about 11 fatalities, i.e., approximately proportional to the threshold of 100 positive cases, chosen in the previous analysis, times the ratio of fatalities to positive cases) are both only four days after the previous predictions, thus confirming the consistency of our results using the positive cases.<fig id="Fig3"><label>Fig. 3</label>
<caption><p>Fit of the cumulative number of diagnosed positive cases of COVID-19 in Italy (red dots) from February 15, 2020 (included), to March 29, 2020 (included), using a fitting function of the type of the Gauss error function with four parameters (solid line). A similar result was obtained with a fitting function of the type of the logistic function</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig3_HTML" id="MO11"></graphic>
</fig>
<fig id="Fig4"><label>Fig. 4</label>
<caption><p>Fit of the cumulative number of fatalities due to COVID-19 in Italy (red dots) from February 15, 2020 (included), to March 29, 2020 (included), using a fitting function of the type of the Gauss error function with four parameters (solid line). The analysis of the curve shows that the date of the flex and the date of a substantial reduction in the number of fatalities are four days after the predictions based on the cumulative number of diagnosed positive cases of COVID-19 in Italy (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
), thus confirming the consistency of our results</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig4_HTML" id="MO12"></graphic>
</fig>
</p>
<p id="Par6">To evaluate the standard deviation relative to the date of the flex, we have used two methods. The first one is to fit the cumulative diagnosed positive cases in Italy using the data from the beginning, February 15, 2020 (included), to March 16, 2020 (included), then from February 15, 2020, to March 17, 2020, and so on, until March 29, 2020 (included), thus getting 14 evaluations of the date of the flex. We then evaluated the standard deviation of these 14 points, and we obtained a 1-sigma standard deviation for <italic>the flex</italic>
of 3 days. The 2-sigma (95.5% probability) and 3-sigma (99.7% probability) uncertainties in the day of the flex are then 6 days and 9 days, respectively. In regard to <italic>the date of a substantial reduction in the number of cumulative positive cases</italic>
in Italy (about 100 cases), the standard deviation is 6 days (68.2% probability), 12 days (95.5% probability) and 18 days (99.7% probability). Both the date of the flex and the date of a substantial reduction in the number of cases have a quasi oscillating behavior whose amplitude is decreasing by increasing the final date of the analysis, i.e., increasing the number of points used in the fit. Indeed, using the last seven evaluations out of the 14 ones, the previous standard deviations reduce by a factor of about two. However, using a second more robust method, by a Monte Carlo analysis [<xref ref-type="bibr" rid="CR8">8</xref>
–<xref ref-type="bibr" rid="CR10">10</xref>
], we evaluated narrower uncertainties for the day of the flex and for the day of a substantial reduction in the number of diagnosed positive cases, as described in Sect. <xref rid="Sec4" ref-type="sec">4</xref>
.</p>
</sec>
<sec id="Sec4"><title>Monte Carlo simulations of cumulative positive cases of COVID-19 in Italy</title>
<p id="Par7">The Monte Carlo simulations [<xref ref-type="bibr" rid="CR8">8</xref>
–<xref ref-type="bibr" rid="CR10">10</xref>
] have been designed to possibly take into account the measurement error in each daily number of the cumulative positive cases of COVID-19 in Italy. This error should describe the uncertainty in the process of measuring the daily number of positive cases due to fluctuations in the measurement procedures (such as a different number of performed daily nasopharyngeal swabs of one day with respect to another day); of course, this error does <underline>not</underline>
describe the difference between the actual total positive cases and the diagnosed ones which can be very large [<xref ref-type="bibr" rid="CR4">4</xref>
]. However, the diagnosed cases are hypothesized to be a representative sample of the actual population (i.e., of the total number of positive cases, which is unknown). To get an estimate of the uncertainty in each daily number, we applied the following heuristic approach. We have assumed a measurement uncertainty in the total number of positive cases equal to 10% of each daily number (Gaussian distributed).</p>
<p id="Par8">The second step was to generate a random matrix <inline-formula id="IEq1"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left( {m\times n} \right) $$\end{document}</tex-math>
<mml:math id="M4"><mml:mfenced close=")" open="("><mml:mrow><mml:mi>m</mml:mi>
<mml:mo>×</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:math>
<inline-graphic xlink:href="13360_2020_383_Article_IEq1.gif"></inline-graphic>
</alternatives>
</inline-formula>
, where <italic>n</italic>
(columns) is the number of observed days and <italic>m</italic>
(rows) is the number of random outcomes, which we have chosen to be 150. Each number in the matrix is part of a Gaussian distribution with mean equal to 1 and sigma equal to 0.1 (i.e., 10% of 1), both row-wise and column-wise. In such a way, each day will be characterized by 150 simulated outcomes that allow to apply a statistical approach. The 150 outcomes represent a reasonably large number of simulated deviations from the official data.</p>
<p id="Par9">We then multiplied the nominal value of the diagnosed positive cases of the <italic>j</italic>
th day for the 150 numbers of the <italic>j</italic>
th column of the random matrix mentioned above. In such a way, each day will be associated with a series of 150 numbers (Gaussian distributed with a 10% standard deviation), which simulate the statistical nature of a single datum. The index <italic>j</italic>
will run from February 15, 2020, to March 26, 2020. Finally, we integrated the daily data to obtain 150 series of the cumulative diagnosed positive cases that allowed to perform the statistical analysis.</p>
<p id="Par10">In summary, starting from the <italic>n</italic>
nominal values of the daily data, we generated <italic>n</italic>
Gaussian distributions with 150 outcomes, with mean equal to the <italic>n</italic>
nominal values and with 10% standard deviation. Then, for each of the 150 simulations, these <italic>n</italic>
values (corresponding to the cumulative positive cases of <italic>n</italic>
days) were fitted with a four-parameter function of the type of the Gauss error function, and we then determined the date of the flex with such fitted function for each simulation. Using the fitted function, we also determined the date at which the number of daily positive cases will be at a certain threshold that, for example, we have chosen to be 100. Finally, we calculated the standard deviation of these 150 simulations. In Fig. <xref rid="Fig5" ref-type="fig">5</xref>
and Fig. <xref rid="Fig6" ref-type="fig">6</xref>
, we report the values (red dots) and the mean (horizontal solid line) of the Monte Carlo simulations, respectively, for the date of the flex and for the date of a substantial reduction in the number of daily positive cases.<fig id="Fig5"><label>Fig. 5</label>
<caption><p>Monte Carlo simulations: each red dot corresponds to the day of occurrence of the flex (after which there is a reduction in the number of daily cases, i.e., there is a deceleration in the number of daily cases), obtained with each of the 150 Monte Carlo simulations. The vertical axis reports the number of days from February 15, 2020</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig5_HTML" id="MO13"></graphic>
</fig>
<fig id="Fig6"><label>Fig. 6</label>
<caption><p>Monte Carlo simulations: each red dot corresponds to the day in which a substantial reduction in the number of daily cases (about 100) occurs. Each red dot is obtained with each of the 150 Monte Carlo simulations. The vertical axis reports the number of days from February 15, 2020</p>
</caption>
<graphic xlink:href="13360_2020_383_Fig6_HTML" id="MO14"></graphic>
</fig>
</p>
<p id="Par11">Using <inline-formula id="IEq2"><alternatives><tex-math id="M5">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$n = 41\hbox { days}$$\end{document}</tex-math>
<mml:math id="M6"><mml:mrow><mml:mi>n</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>41</mml:mn>
<mml:mspace width="0.333333em"></mml:mspace>
<mml:mtext>days</mml:mtext>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="13360_2020_383_Article_IEq2.gif"></inline-graphic>
</alternatives>
</inline-formula>
(i.e., the number of daily diagnosed positive cases up to March 26, 2020), the mean of the 150 Monte Carlo simulations gives the expected dates of March 25, 2020, and April 22, 2020, for the flex and the day of a substantial reduction in the number of daily cases (about 100), respectively. We then obtained a standard deviation (1-sigma) of 1 day for the date of the flex and of 2.3 days for the date in which a substantial reduction in the number of daily cases would be about 100. This result corresponds to a probability of 68.2% that the date of the flex will be at a certain date plus or minus 1 day and that the date of a substantial reduction in the number of cases will be at a certain date plus or minus 2.3 days. A 2-sigma standard deviation will give a more robust probability of 95.4% of the day of the flex and of the day of a substantial reduction in the number of cases. The 2-sigma values correspond to plus or minus 2 days for the day of the flex and plus or minus 4.6 days for the day of a substantial reduction in the number of cases. A similar uncertainty was obtained for the day of the flex and the day of a substantial reduction in the number of fatalities.</p>
</sec>
<sec id="Sec5"><title>Conclusions</title>
<p id="Par12">In this paper, we considered the first 44 days for fitting the cumulative diagnosed positive cases of COVID-19 in Italy (i.e., from February 15, 2020, to March 29, 2020) and the first 41 days for the Monte Carlo simulations. The function used for the fitting is of the type of the Gauss error function with four parameters (a distribution function which well fits the corresponding cumulative diagnosed positive cases in China). We obtained that the day of the flex (i.e., the day of deceleration in the number of daily positive cases) is in Italy, with a 95.4% probability, between March 23, 2020, and March 27, 2020; the 2-sigma uncertainty of <inline-formula id="IEq3"><alternatives><tex-math id="M7">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$+/-$$\end{document}</tex-math>
<mml:math id="M8"><mml:mrow><mml:mo>+</mml:mo>
<mml:mo stretchy="false">/</mml:mo>
<mml:mo>-</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="13360_2020_383_Article_IEq3.gif"></inline-graphic>
</alternatives>
</inline-formula>
2 days was obtained with Monte Carlo simulations. In regard to the day of a substantial reduction in the number of daily positive cases (which, for example, we took to be 100), this day will be in Italy, with a 95.4% probability, between April 17, 2020, and April 27, 2020; the 2-sigma uncertainty of <inline-formula id="IEq4"><alternatives><tex-math id="M9">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$+/-$$\end{document}</tex-math>
<mml:math id="M10"><mml:mrow><mml:mo>+</mml:mo>
<mml:mo stretchy="false">/</mml:mo>
<mml:mo>-</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="13360_2020_383_Article_IEq4.gif"></inline-graphic>
</alternatives>
</inline-formula>
4.6 days was also obtained with the Monte Carlo simulations. In regard to the day of the flex and the day of a substantial reduction in the number of fatalities, we obtained a four-day forward shift with respect to the corresponding dates of the positive cases. The predictions discussed in this paper are statistical in nature and do not explicitly take into account the relevant factors of the daily number of nasopharyngeal swabs, social distancing, and epidemiological and virology studies, which are outside the analysis of the present paper.</p>
<p id="Par13">It is important to stress that the mathematical predictions reported in this paper provide an indication of the approximate period when a substantial reduction of positive cases and fatalities is expected and not the end of the extremely important mitigation measures.</p>
</sec>
</body>
<back><ack><title>Acknowledgements</title>
<p>We gratefully thank Richard Matzner of the University of Texas at Austin for helpful suggestions, Alessandro Paolozzi, and Claudio Paris of Centro Fermi.</p>
</ack>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<mixed-citation publication-type="other"><ext-link ext-link-type="uri" xlink:href="http://www.salute.gov.it/portale/home.html">http://www.salute.gov.it/portale/home.html</ext-link>
</mixed-citation>
</ref>
<ref id="CR2"><label>2.</label>
<mixed-citation publication-type="other"><ext-link ext-link-type="uri" xlink:href="https://www.worldometers.info/coronavirus/country/italy/">https://www.worldometers.info/coronavirus/country/italy/</ext-link>
</mixed-citation>
</ref>
<ref id="CR3"><label>3.</label>
<mixed-citation publication-type="other"><ext-link ext-link-type="uri" xlink:href="https://www.who.int/emergencies/diseases/novel-coronavirus-2019">https://www.who.int/emergencies/diseases/novel-coronavirus-2019</ext-link>
</mixed-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>R</given-names>
</name>
<name><surname>Pei</surname>
<given-names>S</given-names>
</name>
<name><surname>Chen</surname>
<given-names>B</given-names>
</name>
<name><surname>Song</surname>
<given-names>Y</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>T</given-names>
</name>
<name><surname>Yang</surname>
<given-names>W</given-names>
</name>
<name><surname>Shaman</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2)</article-title>
<source>Science</source>
<year>2020</year>
<volume>1</volume>
<fpage>eabb3221</fpage>
<pub-id pub-id-type="doi">10.1126/science.abb32214</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<mixed-citation publication-type="other">I. Ciufolini, A. Paolozzi, Prediction of the time evolution of the Covid-19 Pandemic in Italy by a Gauss Error Function and Monte Carlo simulations. Submitted to BioRxiv on 03.26.2020 and transferred on 03.27.2020 to MedRxiv. 10.1101/2020.03.27.20045104</mixed-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lessler</surname>
<given-names>J</given-names>
</name>
<name><surname>Reich</surname>
<given-names>NG</given-names>
</name>
<name><surname>Brookmeyer</surname>
<given-names>R</given-names>
</name>
<name><surname>Perl</surname>
<given-names>TM</given-names>
</name>
<name><surname>Nelson</surname>
<given-names>KE</given-names>
</name>
<name><surname>Cummings</surname>
<given-names>DA</given-names>
</name>
</person-group>
<article-title>Incubation periods of acute respiratory viral infections: a systematic review</article-title>
<source>Lancet Infect Dis</source>
<year>2009</year>
<volume>9</volume>
<issue>5</issue>
<fpage>291</fpage>
<lpage>300</lpage>
<pub-id pub-id-type="doi">10.1016/S1473-3099(09)70069-6</pub-id>
<pub-id pub-id-type="pmid">19393959</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<mixed-citation publication-type="other">IHME COVID-19 Health Service Utilization Forecasting Team, University of Washington, Forecasting COVID-19 impact on hospital bed-days, ICU-days, ventilator- days and deaths by US state in the next 4 months. Preprint submitted to MedRxiv. 10.1101/2020.03.27.20043752</mixed-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Press</surname>
<given-names>WH</given-names>
</name>
<name><surname>Flannery</surname>
<given-names>BP</given-names>
</name>
<name><surname>Teukolsky</surname>
<given-names>SA</given-names>
</name>
<name><surname>Vetterling</surname>
<given-names>WT</given-names>
</name>
</person-group>
<source>Numerical Recipes</source>
<year>1989</year>
<publisher-loc>Cambridge</publisher-loc>
<publisher-name>Cambridge University Press</publisher-name>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Bevington</surname>
<given-names>PR</given-names>
</name>
<name><surname>Robinson</surname>
<given-names>D Keith</given-names>
</name>
</person-group>
<source>Data Reduction and Error Analysis for the Physical Sciences</source>
<year>1969</year>
<publisher-loc>New York</publisher-loc>
<publisher-name>McGraw-Hill</publisher-name>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ciufolini</surname>
<given-names>I</given-names>
</name>
<name><surname>Monge</surname>
<given-names>BM</given-names>
</name>
<name><surname>Paolozzi</surname>
<given-names>A</given-names>
</name>
<name><surname>König</surname>
<given-names>R</given-names>
</name>
<name><surname>Sindoni</surname>
<given-names>G</given-names>
</name>
<name><surname>Michalak</surname>
<given-names>G</given-names>
</name>
<name><surname>Pavlis</surname>
<given-names>EC</given-names>
</name>
</person-group>
<article-title>Monte Carlo simulations of the LARES space experiment to test General Relativity and fundamental physics</article-title>
<source>Class. Quant. Grav.</source>
<year>2013</year>
<volume>30</volume>
<issue>23</issue>
<fpage>235009</fpage>
<pub-id pub-id-type="doi">10.1088/0264-9381/30/23/235009</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/StressCovidV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000115 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000115 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= StressCovidV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |