-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathindex.tmpl
80 lines (62 loc) · 3.53 KB
/
index.tmpl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
<html>
<head>
<meta http-equiv="CONTENT-TYPE" content="text/html; charset=iso-8859-1"/>
<title>COVID-19 forecasts</title>
</head>
<body dir="ltr" lang="en-US" bgcolor="White">
<div style="margin: 10px 10px 10px 10px; width: 1200px">
<center>
<br/>
<big>
<big>
<b><u>COVID-19 statistical forecasts</u></b>
</big>
<br/>
<font size="1" > If the page does not update, press shift + refresh </font>
<br/>
<font size="3" > Last update <DATE> </font>
<br/><br/>
<b>Disclaimer: this is a scientific exercise, do not take it too seriously.
<br/>
<br> <i> Plots on the left report the cumulated number of deaths, plots on the right the daily number of deaths.
<br/> FORECASTS SHOWING DASHED LINES AND "NA" IN THEIR ESTIMATE ARE NOT RELIABLE SO FAR.
<br/> Since April 3rd, a 5 day running mean is applied to all data. </i></b>
<br/> Since April 9th, with Davide Proment we implemented the same exercise for England regions on his <a href=https://github.com/davideproment/covid-19-analysis >GitHub page.</a>
<br/>
<br/>
<br/>
<img style="border: 0px solid; width: 1024px" alt="" src="covid-status-today.svg"/>
<br/><br/><br/>
<i>Italian data from <a href=https://github.com/pcm-dpc/COVID-19 >GitHub repository of Italian Civil Protection.</a></i>
<br>
<i>World data from <a href=https://github.com/CSSEGISandData/COVID-19 >GitHub repository of Johns Hopkins University.</a></i>
<br>
<br/>
<div style="margin: 10px 10px 10px 10px; width: 1000px">
Methodology note </br>
<font size="3" > <i> Forecasts are estimated with Gompertz functions computed via the nls.lm R function. Ribbons show the 95% confidence interval computed with the confint2 R function.
<br> <font color="red"> Since April 3rd, a 5 day running mean is applied to all data</font>
</br> The non-linear model to forecast the evolution is reliable only after one third of its path (Winsor, PNAS, 1932). Hence errors are likely larger than currently estimated in situation where the outbreak is just started.
</br> Saturation values are the alpha parameter of the Gompertz fit. Errors are estimated using the half of the distance between the 2.5% and the 97.5% percentiles. When forecasts error is larger than half of the forecast value, i.e. the forecast is completely unreliable, NA is shown.
</br> Forecasts are updated at least twice a day (18.30 and 05:30).
</br> Some inconsistencies in Italian prediction may arise from slight data differences between the Civil Protection and the Johns Hopkins archives. </br> The code can be found on <a href=https://github.com/oloapinivad/Cointreau >my GitHub repository</a>, althought it is really badly written. </i> </font>
</div>
<br/><br/><br/>
<b> Italian Time Evolutions of the Prediction </b>
<br/>
<img style="border: 0px solid; width: 800px" alt="" src="forecast_predict_evolution.svg"/>
<br/>
<b> World Time Evolutions of the Prediction </b>
<br/>
<img style="border: 0px solid; width: 800px" alt="" src="forecast_predict_world_evolution.svg"/>
<br/>
<b> Italian Forecasts Evolutions </b>
<br/>
<img style="border: 0px solid; width: 800px" alt="" src="italyGIF.gif"/>
<br/>
<b> Forecasts from the past 10 days </b>
<br/><br/><br/>
<font size="3" > P. Davini (CNR-ISAC), March 2020 </font>
</div>
</body>
</html>