По состоянию на 30 марта 2020 г. на основе официальных данных была построена зависимость зарегистрированных случаев COVID-19 в Москве. Первая точка на графике соответствует 17 марта (55 ± 7). По горизонтальной оси отложена дата в днях. По вертикальной - число зарегистрированных случаев. Масштаб вертикальной оси логарифмический. Видно, что точки хорошо аппроксимируются прямой линией, что говорит о том, что на данном этапе наблюдается геометрический рост случаев, где знаменатель геометрической прогрессии равен 1,25. Это означает, что количество случаев удваивается приблизительно через каждые три дня (3,08).
Точки на графике аппроксимируются функцией:
Lg(N) = kD + c,
где N - число случаев, k = 0,0976 ± 0,002, c = 2, D - число целых дней (100 случаев соответствует 18 марта + 0,8 суток).
Также была сделана оценка относительно количества случаев, если приведенная зависимость будет сохраняться. 8 апреля 2020 г. в Москве ожидается 10 000 случаев.
Рис. 1. Данные по официальному количеству случаев COVID-19 в г. Москве по состоянию на 30 марта 2020 г.
Рис. 2. Аппроксимация данных по официальному количеству случаев COVID-19 в г. Москве по состоянию на 14 апреля 2020 г.
Рис. 3. Данные по официальному количеству случаев COVID-19 в г. Москве по состоянию на 18 апреля 2020 г.
Рис. 4. Данные по официальному количеству случаев COVID-19 в г. Москве по состоянию на 26 апреля 2020 г. Красная линия - результат аппроксимации МНК с помощью функции, указанной в заголовке рисунка. Данные обработаны в среде ROOT [6]
Рис. 5. Поиск даты "нулевого" распространителя в Москве: 28,4 февраля 2020 г.
Рис. 6. Данные в Москве: 17 марта - 28 апреля 2020 г. Черная линия - подгонка экспонентой. Красная - подгонка функцией e в степени полином третьего порядка
Рис. 7. Данные в Москве: 2 апреля - 30 апреля 2020 г. Черная линия - подгонка экспонентой данных на начальном этапе (рис. 5).
Рис. 8. Данные в Москве: 24 апреля - 30 апреля 2020 г. Красная линия - подгонка экспонентой данных 24-28 апреля 2020 г.
Рис. 9. Данные в Москве: 17 марта - 2 мая 2020 г. Аппроксимация тремя экспонентами
Рис. 10. Количество выздоровевших за один день (11 апреля - 4 мая). Аппроксимация экспонентой
Рис. 11. Количество заболевших за один день (11 апреля - 4 мая). Аппроксимация экспонентой
Рис. 12. Решение задачи о том, когда количество выздоровевших станет равно числу заболевших при скорости роста на 4 мая. Ответ получился порядка 140 дней после 1 мая при количестве 10 млн в день. Это означает, что на 4 мая конец не виден
Рис. 13. Количество подтвержденных случаев (18 апреля - 5 мая). Аппроксимация экспонентой, см. рис. 9
Рис. 14. Количество подтвержденных случаев (18 апреля - 7 мая). Аппроксимация экспонентой, см. рис. 9
Рис. 15. Количество подтвержденных случаев (18 апреля - 12 мая). Аппроксимация экспонентой, см. рис. 9
Рис. 16. Количество подтвержденных случаев (18 апреля - 14 мая). Аппроксимация экспонентой, см. рис. 9
Рис. 17. Количество подтвержденных летальных случаев (1 апреля - 14 мая)
Рис. 18. Количество подтвержденных случаев (17 мая - 24 мая). Аппроксимация экспонентой, см. рис. 9
Рис. 19. Полное количество летальных случаев за период 1 апреля - 29 мая (слева). Количество летальных случаев в день (справа). Аппроксимация прямой линией. Средний рост составляет 1,38 ± 0,04 летальных случаев в день
Рис. 20. Полное количество подтвержденных случаев за период 17 марта - 2 июня (слева). Количество случаев в день (справа). Аппроксимация на правом рисунке гауссианом
Рис. 21. Количество летальных случаев в Москве в день за период 1 апреля - 7 июня (слева). Распределение по последней цифре (справа). Количество событий в гистограмме - 67. Цифры 0, 2, 6, 9 встретились 4 раза. 1 - 12. 7 и 8 - 11. Наиболее вероятное значение - 7. Одно стандартное отклонение - 2,6
Рис. 22. Количество летальных случаев в Москве за период 1 апреля - 12 июня (слева). Количество летальных случаев в день (справа). Аппроксимация на правом рисунке прямыми линиями. Красная: 1 апреля - 28 мая. Синяя: 1 апреля - 12 июня. Начиная с 28 мая, наблюдается резкое изменение наклона прямой. Критерий согласия Пирсона (χ2) на количество степеней свободы указывает на изменение механизма формирования летальных данных и невозможность описания всего набора данных одной прямой линией
Рис. 23. Количество летальных случаев в день в Москве за период 1 апреля - 12 июня (слева). Распределение по последней цифре (справа). Синий линией выполнена подгонка распределения горизонтальной линий - равномерное распределение. Полученное значение 6,1 ± 0,8 отличается от математического ожидания 7,2 больше одного стандартного отклонения. Значение 0 и 2 встретились всего четыре раза. Значение 1 встретилось 13 раз
Рис. 24. Функция вероятности биномиального распределения для n = 72 и вероятности p = 0,1. Для к = 13 вероятность (вероятность, что одна и та же цифра в выборке встретилась к раз) составляет 0,014. Для к = 4 - 0,080
Рис. 25. Количество летальных случаев в день в Москве за период 1 апреля - 3 июля (слева). Распределение по последней цифре (справа). Синий линией выполнена подгонка распределения горизонтальной линий - равномерное распределение. Значение, заканчивающееся на 1, не появлялось 25 дней подряд
Рис. 26. Количество летальных случаев в день в Москве за период 1 апреля - 15 июля (слева). Распределение по последней цифре (справа). Синей линией показана подгонка распределения горизонтальной линией - равномерное распределение.
Рис. 27. Количество летальных случаев в день в Москве за 23 дня на 14 августа 2020 (слева). Красная линия - подгонка прямой линией МНК. Ошибка показана статистическая. Распределение по количеству летальных случаев в день в Москве за 23 дня на 14 августа (справа). Синей линией показана подгонка распределения функцией Гаусса. Ожидаемое (статистическое) значение параметра Sigma: 3,6. Аномально малое значение Sigma противоречит гипотезе о нормальном распределении данных
Рис. 28. Количество летальных случаев в день в Москве за 37 дней на 28 августа 2020 (слева). Красная линия - подгонка прямой линией МНК. Ошибка показана статистическая. Распределение по количеству летальных случаев в день в Москве (справа). Синей линией показана подгонка распределения функцией Гаусса
Рис. 29. Рисунок официальной страницы штаба за 10 сентября 2020 г. по Москве. Красной стрелкой на левом нижнем распределении показано изменение дисперсии случайной величины 22 августа, заметное органолептически
Рис. 30. Официальные данные штаба (см. рис. 29) случаев коронавируса в день в Москве за период 18 марта - 10 сентября 2020 г. Cтрелками отмечены аномалии в среднеквадратичном отклонении точек от среднего значения
Рис. 31. Официальные данные штаба (см. рис. 30) случаев коронавируса в день в Москве за период 26 июля - 22 августа 2020 г. На левом рисунке в результате аппроксимации прямой линией с помощью МНК получено аномально низкое значение параметра согласия хи-квадрат на число степеней свободы (существенно меньше единицы). На правом рисунке представлено распределение случайной величины зарегистрированных случаев в день за тот же период и результат аппроксимации функцией Гаусса. RMS - среднеквадратичное отклонение. Sigma - параметр функции Гаусса (RMS), полученный в результате аппроксимации. Аномально малое значение Sigma противоречит гипотезе о нормальном распределении данных
Рис. 32. Официальные данные штаба (см. рис. 30) случаев коронавируса в день в Москве за период 22 августа - 12 сентября 2020 г. На левом рисунке показан результат аппроксимации данных прямой линией с помощью МНК. На правом рисунке представлено распределение случайной величины зарегистрированных случаев в день за тот же период и результат аппроксимации функцией Гаусса. RMS - среднеквадратичное отклонение. Sigma - параметр функции Гаусса (RMS), полученный в результате аппроксимации
Рис. 33. Официальные данные штаба (см. рис. 30) случаев коронавируса за сутки в Москве за период 3 апреля - 19 сентября 2020 г.
Рис. 34. Официальные данные штаба (см. рис. 30) случаев коронавируса за сутки в Москве за период 17 марта - 20 сентября 2020 г.
Рис. 35. Официальные данные штаба (см. рис. 30) случаев коронавируса за сутки в Москве за период 17 марта - 25 сентября 2020 г. Справа в правой колонке показано фактическое количество случаев на день эпидемии, который отсчитывается от 17-го марта. Видно, что за последние 11 дней данные округлены до числа, кратного 5
Рис. 36. Официальные данные штаба (см. рис. 30) случаев коронавируса за сутки в Москве за период 17 марта - 4 октября 2020 г. Видно, что к 4-му октября значение достигло приблизительно половины максимального значения в мае.
Рис. 37. Официальные данные штаба (см. рис. 30) летальных случаев коронавируса за сутки в Москве за период 18 марта - 4 октября 2020 г. Видно, что к 4-му октября значение достигло приблизительно половины максимального значения в мае.
Рис. 38. Официальные данные штаба (см. рис. 30) зарегистрированных и летальных случаев коронавируса за сутки в Москве за период 18 марта - 16 октября 2020 г. Видно, что ширина на полувысоте для зарегистрированных случаев составляет порядка 27 дней. Сдвиг пика летальности в мае относительно зарегистрированных случаев составил порядка 17 дней.
Рис. 39. Официальные данные штаба (см. рис. 30) летальных и зарегистрированных случаев коронавируса за сутки в Москве за период 18 марта - 23 октября 2020 г.
Рис. 40. Официальные данные штаба (см. рис. 30) случаев смерти (слева) и количество смертей за сутки (справа) в Москве за период 18 марта - 3 декабря 2020 г.
Рис. 41. Аппроксимация количества смертей за сутки (слева) в Москве прямой линией за период 5 ноября - 3 декабря 2020 г. Ошибка отложена статистическая. Наблюдается аномально малое значение критерия согласия Пирсона (χ2) 1,98 на 27 степеней свободы (ожидаемое значение около 27). Справа показано распределение количества смертей за сутки в Москве за период 5 ноября - 3 декабря 2020 г. Подгонка выполнена функцией Гаусса. Критерий согласия Пирсона (χ2) говорит о хорошем качестве подгонки. Значение полуширины на полувысоте составило 2,57 ± 0,8, что существенно меньше значения, ожидаемого для нормального распределения случайной величины (8,6). Наблюдаемое значение полуширины на половине высоты распределения отклонилось от ожидаемого значения на 7,5 стандартных отклонений. Случайное отклонение на такую величину невозможно
Рис. 42. Аппроксимация количества смертей за сутки (слева) в Москве прямой линией за период 5 ноября - 18 января 2020-2021 гг. Ошибка отложена статистическая. Наблюдается аномально малое значение критерия согласия Пирсона (χ2) 7,7 на 73 степени свободы (ожидаемое значение около 73). Коэффициент наклона аппроксимирующей прямой за указанный период в пределах одного стандартного отклонения равен нулю (0,025 ± 0,046). Все 75 точек отклоняются от аппроксимирующей прямой на величину меньше одного стандартного отклонения. Справа показано распределение количества смертей за сутки в Москве за период 5 ноября - 18 января 2020-2021 гг. Подгонка выполнена функцией Гаусса. Критерий согласия Пирсона (χ2) говорит о хорошем качестве подгонки. Значение полуширины на полувысоте составило 2,35 ± 0,26, что существенно меньше значения, ожидаемого для нормального распределения случайной величины (8,6). Наблюдаемое значение полуширины на половине высоты распределения отклонилось от ожидаемого значения на 24 стандартных отклонения. Случайное отклонение измеренного значения на величину больше трех стандартных отклонений от ожидаемого значения невозможно
Рис. 43. Количество смертей (слева), вызванных коронавирусом, в Москве за период 17 марта 2020 г. - 19 марта 2021 г. Ошибка отложена статистическая. Масштаб полулогарифмический. Количество умерших за сутки из-за коронавируса (справа) в Москве за тот же период. Ошибка отложена статистическая (квадратный корень из значения). На правом плато отчетливо заметна аномалия (см. также рис. 42), когда все точки находятся в пределах одного стандартного отклонения (внутри коридора ошибок)
Рис. 44. Количество смертей (слева), вызванных коронавирусом, в Москве за период 17 марта 2020 г. - 7 июля 2021 г. Ошибка отложена статистическая. Масштаб полулогарифмический. Количество умерших за сутки из-за коронавируса (справа) в Москве за тот же период. Ошибка отложена статистическая (квадратный корень из значения). Стрелками отмечены даты, показывающие характерные точки
Рис. 45. Количество смертей за сутки (слева), вызванных коронавирусом, в Москве за период 24 февраля - 12 июня 2021 г. (период между двумя так называемыми волнами). Ошибка отложена статистическая. Красная линия - фитирование МНК прямой линией. Получен статистически значимый рост по смертям в день в Москве: 0,12 ± 0,02. На втором графике показано распределение количества смертей за сутки в Москве за тот же период. Подгонка выполнена функцией Гаусса. Критерий согласия Пирсона (χ2) говорит о хорошем качестве подгонки. На третьем графике показано распределение отклонений точек от фитирующей прямой, см. левый график. Подгонка выполнена функцией Гаусса. Значение полуширины на полувысоте составило 2,5 ± 0,2, что существенно меньше значения, ожидаемого для нормального распределения случайной величины (7,4). Наблюдаемое значение полуширины на половине высоты распределения отклонилось от ожидаемого значения на 24 стандартных отклонения. Случайное отклонение измеренного значения на величину больше трех стандартных отклонений от ожидаемого значения невозможно
Рис. 46. Распределение случаев коронавируса за сутки в Москве по состоянию на 23 июля 2021 г. (слева). Распределение смертей за сутки, вызванных коронавирусом, в Москве по состоянию на тоже число (в середине). Данные смертей смещены к началу эпидемии на 13 дней. Двумерное распределение случаев коронавируса за сутки (по вертикали) и смертей от коронавируса за сутки (по горизонтали). Красной линией показана корреляция событий, если не учитывать аномальные выбросы вверх при постоянном количестве смертей за сутки. На втором и третьем рис. видно аномальное отсутствие событий по количеству умерших за сутки для интервала значений 80 - 100
Рис. 47. Количество смертей от коронавируса за сутки в Москве за период 1 апреля - 1 июня 2020 г. (слева). Красная линия - фит прямой линией (МНК). Распределение смертей за сутки, вызванных коронавирусом, в Москве за тот же интервал (в середине). Распределение отклонений точек от фитирующей прямой, нормированных на статистическую ошибку,- справа. Синяя линия - фит функцией Гаусса (МНК). Параметр Sigma соответствует ожидаемому значению 1.
Рис. 48. Количество смертей от коронавируса за сутки в Москве за период 28 июня - 27 июля 2021 г. (слева). Красная линия - фит прямой линией (МНК). Распределение смертей за сутки, вызванных коронавирусом, в Москве за тот же интервал (в середине). Распределение отклонений точек от фитирующей прямой, нормированных на статистическую ошибку,- справа. Синяя линия - фит функцией Гаусса (МНК). Параметр Sigma отклонился от ожидаемого значения 1 на 6,5 стандартных отклонений
Рис. 49. Количество зарегистрированных случаев коронавируса за сутки в Москве за период 17 марта 2020 г.- 11 августа 2021 г., умноженное на 0,0125 (синие точки). Количество смертей от коронавируса за сутки в Москве за тот же период (красные точки). График сдвинут к левой границе первого пика для данных по зарегистрированным случаям за сутки. С 220 по 300 день эпидемии наблюдается аномально малый разброс данных по количеству смертей за сутки по сравнению с количеством зарегистрированных случаев. С 330 по 440 день эпидемии наблюдается значительное (в 1,7 - 2,5 раза) превышение смертности по сравнению с интервалом 1 - 330
Таблица 1. Официальные данные по количеству случаев
COVID-19 в г. Москве по состоянию на март-июнь 2020 г. [7]
------------------------------
Дата N Ошибка
(дни) (корень из N)
------------------------------
-14 55 7 Значение в первой колонке соответствует 17 марта 2020
-13 86 9
-12 98 10
-11 131 11
-10 137 12
-9 191 14
-8 262 16
-7 290 17
-6 410 20
-5 546 23
-4 703 27
-3 817 29
-2 1014 32
-1 1226 35
0 1613 40
1 1880 43 Значение в первой колонке соответствует 1 апреля
2 2475 50
3 2923 54
4 3357 58
5 3893 62
6 4484 67
7 5181 72
8 5841 76
9 6698 82
10 7822 88
11 8852 94
12 10158 101
13 11513 107
14 13002 114
15 14776 122
16 16146 127
17 18105 135
18 20754 144
19 24324 156
20 26350 162
21 29433 172
22 31981 179
23 33940 184
24 36897 192
25 39509 199
26 42480 206
27 45351 213
28 48426 220
29 50646 225
30 53739 232
1 57300 240
2 62658 250
3 68606 262
4 74401 273
5 80115 283
6 85973 293
7 92676 304
8 98522 314
9 104189 323
10 109740 331
11 115909 340
12 121301 348
13 126004 355
14 130716 362
15 135464 368
16 138969 373
17 142824 378
18 146062 382
19 149607 387
20 152306 390
21 155219 394
22 158207 398
23 161397 402
24 163913 405
25 166473 408
26 169303 411
27 171443 414
28 173497 417
29 175829 419
30 178196 422
31 180791 425
1 183088 428
2 185374 431
3 187216 433
4 189214 435
5 191069 437
6 193061 439
7 195017 442
8 197018 446
9 198590 446
10 199785 447
11 201221 449
12 202935 450
------------------------------
Таблица 2. Официальные данные по количеству летальных случаев
COVID-19 в г. Москве по состоянию на апрель-июнь 2020 г. [7]
------------------------------
Дата N Ошибка
апрель (корень из N)
------------------------------
1 16 4
2 19 4.358898944
3 20 4.472135955
4 27 5.196152423
5 29 5.385164807
6 29 5.385164807
7 31 5.567764363
8 31 5.567764363
9 38 6.164414003
10 50 7.071067812
11 58 7.615773106
12 72 8.485281374
13 82 9.055385138
14 95 9.746794345
15 106 10.29563014
16 113 10.63014581
17 127 11.26942767
18 148 12.16552506
19 176 13.26649916
20 204 14.28285686
21 233 15.26433752
22 261 16.15549442
23 288 16.97056275
24 325 18.02775638
25 366 19.13112647
26 404 20.09975124
27 435 20.85665361
28 479 21.88606863
29 546 23.36664289
30 611 24.71841419
31 658 25.65151068
32 695 26.36285265
33 729 27
34 764 27.64054992
35 816 28.56571371
36 866 29.42787794
37 905 30.08321791
38 956 30.91924967
39 1010 31.78049716
40 1068 32.68026928
41 1124 33.52610923
42 1179 34.33656943
43 1232 35.09985755
44 1290 35.91656999
45 1358 37
46 1432 38
47 1503 39
48 1580 40
49 1651 41
50 1726 42
51 1794 42
52 1867 43
53 1934 44
54 1993 45
55 2034 45
56 2110 46
57 2183 47
58 2254 47
59 2330 48
60 2408 49
61 2477 50
62 2553 51
63 2624 51
64 2685 52
65 2749 52
66 2806 53
67 2864 54
68 2919 54
69 2970 55
70 3029 55
71 3085 56
72 3138 56
73 3187 56
------------------------------
Таблица 3. Официальные данные по количеству различных случаев
COVID-19 в РФ и г. Москве по состоянию на май 2020 г. - январь 2021 г. [7]
----------------------------------------------------------------------
дата РФ Москва
всего выздоровело умерло всего выздоровело умерло
----------------------------------------------------------------------
8 (май) 187859 26608 1723 98522 10259 956
9 198676 31916 1827 104189 12779 1010
10 209688 34306 1915 109740 13790 1068
11 221344 39801 2009 115909 17822 1124
12 232243 43512 2116 121301 19642 1179
13 242271 48003 2212 126004 21506 1232
14 252245 53530 2305 130716 23327 1290
15 262843 58226 2418 135464 24562 1358
16 272043 63166 2537 138969 26032 1432
17 281752 67373 2631 142824 27490 1503
18 290678 70209 2722 146062 28913 1580
19 299941 76130 2837 149607 31496 1651
20 308705 85392 2972 152306 36936 1726
21 317554 92681 3099 155219 40682 1794
22 326448 99825 3249 158207 43582 1867
23 335882 107936 3388 161397 47413 1934
24 344481 113299 3541 163913 49840 1993
25 353427 118798 3633 166473 53586 2034
26 362342 131129 3807 169303 61619 2110
27 370680 142208 3968 171443 67458 2183
28 379051 150993 4142 173497 71251 2254
29 387623 159257 4374 175829 74725 2330
30 396575 167469 4555 178196 78324 2408
31 405843 171883 4693 180791 80179 2477
1 414878 175877 4855 183088 82239 2553
2 423741 186985 5037 185374 88269 2624
3 432277 195957 5215 187216 91654 2685
4 441108 204623 5384 189214 94715 2749
5 449834 212680 5528 191069 97358 2806
6 458689 221388 5725 193061 100431 2864
7 467673 226731 5859 195017 102714 2919
8 476658 230688 5971 197018 104347 2970
9 485253 242397 6142 198590 109737 3029
10 493657 252783 6358 199785 113533 3085
11 502436 261150 6532 201221 115769 3138
12 511423 269370 6715 202935 118024 3187
13 520129 274641 6829 204428 119558 3231
14 528964 280050 6948 205905 121691 3281
15 537210 284539 7091 207264 123240 3334
16 545458 545458 7284 208680 126114 3386
17 553301 304342 7478 209745 128385 3434
18 561091 313963 7660 210785 130431 3483
19 569063 324406 7841 211921 132434 3531
20 576952 334592 8002 212978 134549 3565
21 584680 339711 8111 213946 135556 3597
22 592280 344416 8206 215014 136670 3617
23 599705 356429 8359 216095 139248 3643
24 606881 368822 8513 216906 141268 3657
25 613994 375164 8605 217791 142194 3669
26 620794 384152 8781 218604 144191 3694
27 627646 393352 8969 219354 146043 3714
28 634437 399087 9073 220071 147095 3738
29 641156 403430 9166 220853 148102 3761
30 647849 412650 9320 221598 150397 3796
1 654405 422931 9536 222209 152193 3831
2 661165 428978 9683 222871 153248 3870
3 667883 437893 9859 223530 154839 3904
4 674515 446879 10027 224210 156459 3929
5 681251 450750 10161 224860 157020 3953
6 687862 454329 10296 225545 157652 3975
7 694230 463880 10494 226174 169718 3999
8 700792 472511 10667 226795 160913 4027
9 707301 481316 10843 227363 162086 4059
10 713936 489068 11017 228000 163093 4087
11 720547 497446 11205 228678 164095 4116
12 727162 501061 11335 229357 164560 4143
13 733699 504021 11439 230029 164954 4168
14 739947 512825 11614 230642 165966 4205
15 746369 523249 11770 231270 167810 4234
16 752797 531692 11937 231801 168955 4258
17 759203 539373 12123 232376 169637 4271
18 765437 546863 12247 232954 170180 4285
19 771546 550344 12342 233545 170539 4299
20 777486 553602 12427 234123 170918 4314
21 783328 562384 12580 234725 171771 4331
22 789190 572053 12745 235363 172939 4350
23 795038 580330 12892 235971 173392 4364
24 800849 588774 13046 236616 173854 4375
25 806720 597140 13192 237264 174415 4389
26 812485 600250 13269 237947 174721 4398
27 818120 603329 13354 238641 174985 4411
28 823515 612217 13504 239315 175812 4421
29 828990 620333 13673 239986 176437 4434
30 834499 629655 13802 240664 177908 4446
31 839981 638410 13963 241359 179168 4460
1 845443 646524 14058 242049 180581 4473
2 850870 650173 14128 242713 181342 4485
3 856264 653593 14207 243406 182184 4498
4 861426 661471 14351 244097 183640 4510
5 866627 669026 14490 244784 184969 4521
6 871894 676357 14606 245468 186223 4534
7 877135 683592 14725 246154 187576 4546
8 882347 690207 14854 246845 188854 4560
9 887536 693422 14931 247534 189657 4572
10 892654 696681 15001 248228 190555 4585
11 897599 703175 15131 248922 191598 4599
12 902701 710298 15260 249611 192686 4611
13 907758 716396 15384 250303 193473 4622
14 912823 722964 15498 250991 194823 4633
15 917884 729411 15617 251686 196098 4645
16 922853 732968 15685 252374 197143 4656
17 927745 736101 15740 253064 197962 4666
18 932493 742628 15872 253757 199117 4677
19 937321 749423 15989 254448 200359 4687
20 942106 755713 16099 255136 201567 4698
21 946976 761330 16189 255826 202727 4710
22 951897 767477 16310 256513 203821 4720
23 956749 770639 16383 257124 204786 4731
24 961493 773095 16448 257749 205649 4741
25 966189 779747 16568 258430 206804 4753
26 970865 786150 16683 259070 208193 4764
27 975576 792561 16804 259707 209628 4776
28 980405 798466 16914 260361 210908 4786
29 985346 804383 17025 261038 212301 4798
30 990326 806982 17093 261733 213205 4809
31 995319 809387 17176 262418 214121 4821
1 1000048 815705 17299 263059 215383 4832
2 1005000 821169 17414 264684 216742 4844
3 1009995 826935 17528 264374 217899 4857
4 1015105 832747 17649 265066 219095 4867
5 1020310 838126 17759 265737 220407 4878
6 1025505 840949 17820 266357 221370 4891
7 1030690 843277 17871 267047 222252 4905
8 1035789 850049 17993 267742 223502 4921
9 1041007 856458 18135 268384 224840 4933
10 1046370 862373 18263 269079 226274 4947
11 1051874 868107 18365 269777 227530 4956
12 1057362 873535 18484 270447 228673 4968
13 1062811 876225 18578 271097 229573 4982
14 1068320 878700 18635 271793 230397 4993
15 1073849 884305 18785 272523 231483 5006
16 1079519 890114 18917 273273 232706 5016
17 1085281 895868 19061 273003 234029 5025
18 1091186 901207 19195 274808 235136 5033
19 1097251 906462 19339 275633 236055 5044
20 1103399 909357 19418 276493 236947 5057
21 1109595 911973 19489 277408 237935 5069
22(175) 1115810 917949 19649 278388 239126 5084
23 1122241 923699 19799 279358 240352 5100
24 1128836 929829 19948 280408 241656 5115
25(178) 1136048 934146 20056 281968 242883 5129
26 1143571 940150 20225 283760 244006 5146
27 1151438 943218 20324 285776 245132 5164
28(181) 1159573 945920 20385 287993 246074 5180
29 1167805 952399 20545 290293 247301 5203
30 1176286 958257 20722 292601 248650 5230
1 (184) 1185231 964242 20891 295025 249861 5254
2 1194643 970296 21077 297729 251018 5282
3 1204502 975859 21251 300613 252087 5314
4 (187) 1215001 979143 21358 303940 253031 5343
5 1225889 982324 21475 307477 253944 5370
6(189) 1237504 988576 21663 311559 255100 5401
7 1248619 995275 21865 314788 256362 5442
8(191) 1260112 1002329 22056 318111 257703 5497
9 1272238 1009421 22257 321812 259127 5530
10(193) 1285084 1016202 22454 325917 260433 5560
11 1298719 1020442 22597 330418 261619 5595
12 1312310 1024235 22722 334813 262750 5629
13 1326178 1031785 22966 339431 264500 5687
14 1340409 1039705 23205 344004 266470 5739
15 1354163 1048097 23491 347946 268423 5796
16(199) 1369313 1056582 23723 352995 270226 5850
17 1384235 1065199 24002 357643 272227 5906
18 1339334 1070576 24187 362253 273881 5958
19 1415316 1075904 24366 367629 275648 6009
20 1431635 1085608 24635 372628 277608 6058
21 1447335 1096560 24952 377017 279973 6121
22 1463306 1107988 25242 381430 282492 6187
23 1480646 1119251 25525 386908 285271 6249
24 1497167 1130818 25821 391361 287910 6312
25 1513877 1138522 26050 395816 290524 6380
26 1531224 1146096 26269 401040 293205 6442
27 1547774 1158940 26589 405352 296474 6503
28 1563976 1171301 26935 409022 300051 6578
29 1581693 1186041 27301 413928 303831 6644
30 1599976 1200560 27656 419196 307816 6713
31 1618116 1215414 27990 424148 311450 6768
1 1636781 1225673 28235 429409 314791 6820
2 1655038 1236033 28473 434205 318032 6873
3 1673686 1251364 28828 439355 321679 6936
4 1693454 1266931 29217 445181 325531 7004
5 1712858 1279169 29509 450436 329201 7071
6(220) 1733440 1296124 29887 456689 332984 7140
7 1753836 1312927 30251 462518 336659 7215
8 1774334 1324419 30537 468269 339913 7289
9 1796132 1335141 30793 475166 343226 7361
10 1817109 1350741 31161 481068 346889 7429
11 1836960 1369357 31593 485545 350746 7502
12 1858568 1388168 32032 491542 354705 7573
13 1880551 1406903 32443 497516 358924 7643
14 1903253 1425529 32834 503943 363347 7712
15 1925825 1439985 33186 510214 367613 7787
16 1948603 1453849 33489 516574 371857 7859
17 1971013 1475904 33931 522456 377131 7933
18 1991998 1501083 34387 526630 382634 8009
19 2015608 1526656 34850 533068 388416 8082
20 2039926 1551414 35311 539970 393252 8159
21 2064748 1577435 35778 547138 400332 8233
22 2089329 1595443 36179 553713 405577 8308
23 2114502 1611445 36540 560579 410641 8379
24 2138828 1634671 37031 566417 416242 8455
25 2162503 1660419 37538 571102 422189 8530
26 2187990 1685492 38062 577177 428418 8603
27 2215533 1712174 38558 585095 434902 8680
28 2242633 1739470 39068 592415 441535 8756
29 2269316 1761457 39527 599213 447111 8828
30 2295654 1778704 39895 605724 452184 8902
1 2322056 1803467 40464 612248 458409 8978
2 (246) 2347401 1830349 41053 617439 464984 9053
3 (247) 2375546 1859851 41607 625189 471833 9126
4 2402949 1888752 42176 632057 478424 9203
5 2431731 1916396 42684 640050 485308 9277
6 2460770 1937738 43141 647562 490844 9349
7 2488912 1956588 43597 654841 495871 9425
8 2515009 1981526 44159 660073 501403 9496
9 2541199 2007792 44718 665218 507149 9571
10 2569126 2033669 45280 671948 512969 9645
11 2597711 2059840 45893 679163 518889 9722
12 2625848 2085958 46453 685785 524779 9798
13 2653928 2106235 46941 692210 529620 9870
14 2681256 2124797 47391 698084 534552 9945
15 2707945 2149610 47968 703502 539859 10022
16 2734454 2176100 48564 708530 545430 10095
17 2762668 2202540 49151 715241 551207 10171
18 2791220 2228633 49762 722178 557028 10243
19 2819429 2254742 50347 728637 562953 10317
20 2848377 2275657 50858 735900 567743 10394
21 2877727 2295362 51351 743697 572541 10469
22 2906503 2319520 51912 750934 577544 10540
23 2933753 2343967 52461 756586 582772 10613
24 2963688 2370857 53096 764789 588280 10689
25 2992706 2398254 53659 772104 593939 10763
26 3021964 2426439 54226 779584 599541 10840
27 3050248 2450828 54778 786204 604604 10913
28 3078035 2471309 55265 792457 609770 10989
29 3105037 2496183 55827 798098 615396 11064
30 3131550 2525418 56426 803203 621292 11136
31 3159297 2554340 57019 809769 627368 11209
1 3186336 2580138 57555 815676 633702 11280
2 3212637 2559035 58002 821128 639809 11356
3 3236787 2618882 58506 824979 645923 11433
4 3260138 2640036 58988 828570 652324 11504
5 3284384 2662668 59506 833412 659076 11582
6 3308601 2685723 59951 838554 666282 11652
7 3332142 2709452 60457 843425 673692 11728
8 3335794 2731129 60911 848452 680759 11799
9 3379103 2754809 61381 853352 688061 11873
10 3401954 2778889 61837 857568 695106 11950
11 3425269 2800675 62273 862214 701790 12025
12 3448203 2825430 62804 867215 708568 12097
13 3471053 2854088 63370 871535 715438 12173
14 3495816 2882044 63940 877428 722487 12241
15 3520531 2909680 64495 882962 729089 12322
16 3544623 2936991 65085 887636 734728 12401
17 3569209 2960431 65566 891648 739340 12468
18 3591066 2978764 66037 895327 743349 12541
----------------------------------------------------------------------
Приложение A
------------
Код C++ для построения двух графиков из одного набора данных в файле в среде ROOT. «Фитирование» функцией закомментировано. Файл с кодом должен называться deathsday.c. Файлы deathsday.c и 18-deaths-now.txt должны резмещаться в одной папке (директории). Строка для запуска функции deathsday в ROOT:
root [0] .x deathsday.c
// Plot data from file and draw a new plot
void deathsday ()
{
TCanvas *c1 = new TCanvas("c1","Deaths, COVID-19, Moscow",0,0,700,512);
c1->Divide(1,2);
c1->cd(1);
// Read data from file
TGraphErrors *g = new TGraphErrors("18-deaths-now.txt", "%lg %lg %lg", "");
g->SetMaximum(30000);
g->SetMinimum(10);
// Draw data from file
g->SetTitle("Deaths, COVID-19, Moscow, from 17 March 2020;");
g->SetMarkerStyle(7);
g->SetLineColor(kBlack);
g->Draw("AP");
gPad->SetLogy();
c1->cd(2);
// Read data from file
TGraphErrors *g1 = new TGraphErrors("18-deaths-now.txt", "%lg %lg %lg", "");
g1->SetTitle("Deaths per day, COVID-19, Moscow, from 17 March 2020");
g1->SetMarkerStyle(7);
g1->SetMarkerColor(kRed);
int msizeFromGet, msize, val;
double x, y, yprev, deathday, errorY;
// Get number of lines in the file (points in TGraf g)
msize = g->GetN();
yprev = 0;
for(int i = 0; iGetPoint(i, x, y);
cout
Fit(f);
g1->Draw("AP");
TF1 *f1 = new TF1("f1", "1.376*x-4.053", 2, 70);
f1->SetLineColor(kRed);
f1->SetLineWidth(1);
// f1->Draw("AP SAME");
}
Содержимое файла "18-deaths-now.txt"
------------------------------------
-14 0 0
-13 0 0
-12 0 0
-11 0 0
-10 0 0
-9 0 0
-8 0 0
-7 0 0
-6 0 0
-5 0 0
-4 0 0
-3 0 0
-2 0 0
-1 0 0
0 0 0
1 16 4
2 19 4.358898944
3 20 4.472135955
4 27 5.196152423
5 29 5.385164807
6 29 5.385164807
7 31 5.567764363
8 31 5.567764363
9 38 6.164414003
10 50 7.071067812
11 58 7.615773106
12 72 8.485281374
13 82 9.055385138
14 95 9.746794345
15 106 10.29563014
16 113 10.63014581
17 127 11.26942767
18 148 12.16552506
19 176 13.26649916
20 204 14.28285686
21 233 15.26433752
22 261 16.15549442
23 288 16.97056275
24 325 18.02775638
25 366 19.13112647
26 404 20.09975124
27 435 20.85665361
28 479 21.88606863
29 546 23.36664289
30 611 24.71841419
31 658 25.65151068
32 695 26.36285265
33 729 27
34 764 27.64054992
35 816 28.56571371
36 866 29.42787794
37 905 30.08321791
38 956 30.91924967
39 1010 31.78049716
40 1068 32.68026928
41 1124 33.52610923
42 1179 34.33656943
43 1232 35.09985755
44 1290 35.91656999
45 1358 37
46 1432 38
47 1503 39
48 1580 40
49 1651 41
50 1726 42
51 1794 42
52 1867 43
53 1934 44
54 1993 45
55 2034 45
56 2110 46
57 2183 47
58 2254 47
59 2330 48
60 2408 49
61 2477 50
62 2553 51
63 2624 51
64 2685 52
65 2749 52
66 2806 53
67 2864 54
68 2919 54
69 2970 55
70 3029 55
71 3085 56
72 3138 56
73 3187 56
------------------------------------
Литература
[1] Menachery, V., Yount, B., Debbink, K. et al. A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence. Nat Med 21, 1508-1513 (2015). https://www.nature.com/articles/nm.3985
[2] Coronavirus: the first three months as it happened. 17 April 16:00 bst - Wuhan death toll jumps by 50%. Nature News 2020
[3] Собянин сообщил о сокращении тестирования на коронавирус из-за дефицита пипеток. «Коммерсантъ» 23 апреля 2020
[4] Мэрия Москвы пообещала «миллион призов» для участников голосования. «РБК» 11 июня 2020
[5] Binomial distribution and probability mass function. WikipediA
[6] ROOT. Data Analysis Framework. https://root.cern/
[7] Официальная информация о коронавирусе в России. https://стопкоронавирус.рф/
[8] Enrico Bucci. Note of concern https://cattiviscienziati.com/2020/09/07/note-of-concern/
[9] РБК: Число суточных заражений COVID в России достигло максимума с 26 июля https://www.rbc.ru/society/17/09/2020/5f63118a9a7947d44b488efd?
[10] S. Mallapaty, Nature, January 22, 2021: Are COVID vaccination programmes working? Scientists seek first clues
[11] D. Y. Logunov et al., The Lancet, February 02, 2021: Safety and efficacy of an rAd26 and rAd5 vector-based heterologous prime-boost COVID-19 vaccine: an interim analysis of a randomised controlled phase 3 trial in Russia
[12] «Коммерсантъ». Онлайн-статистика, графики и карты распространения COVID-19: Последние данные по коронавирусу в России и мире
[13] B. Nogrady, Nature, July 06, 2021: Mounting evidence suggests Sputnik COVID vaccine is safe and effective