Изменил немного методику анализа, и добавил статистику за 06 апреля (данные на 07 апреля).
Что бы правильно анализировать данные стал использовать суммарное число больных на данный день, т.е. из числа заболевших всего вычел число выздоровевших всего и число умерших всего. Теперь "идеологически" можно фитировать такие данные Гауссианом, так же как и данные за день.
xc - позиция пика
w - "ширина пика"
На графиках так же вставки с результатами фитирования. Если верить результатам фитирования, то пик по больным будет через 3 недели, а пик по заболеваниям за день через 5 недель плюс минус "большой лапоть".
По количеству умерших картина выглядит более "оптимистичной", согласно результатам фитирования, пик пройден.
Все это игра цифр и не более того, пока что. Но некоторые оценки и предположения можно делать.
Доля летальных исходов к числу зарегистрированным заболевшим по суточным данным держится в среднем на уровне 1%.
Комментарии
17 марта статистика по Германии
Заражённых - 7257
Смертей - 17
Выздоровлений - 67
Обращаю внимание на Германию по причине похожего количества умерших к выздоровлениям.
А у них есть такая статистика в открытом доступе?
Думаю, начало вспышки эпидемии в России, октябрь-ноябрь 2019 года. Мои родные в это время переболели. Сильный кашель, температура. Обращались в клинику и делали по собственной инициативе флюорографию, подозревали пневмонию.
https://www.1tv.ru/news/issue/2019-11-19/09:00#1
Меня в декабре 2019 подкосило. Три недели промучился.
Но тут нюанс. За два дня до первых симптомов меня уговорили сделать прививку от гриппа. Поэтому все ухудшения я списывал на эту прививку. Так три недели и проболел ходя на работу и ожидая, что вот-вот полегчает.
Насколько знаю, то современная российская прививка от гриппа не предполагает немедленную реакцию иммунной системы организма на неё. Поэтому в интернете "нахваливают" какие-то забугорные, которые якобы "живые". Вам какую делали?
Точно. Из-за взрыва на "Векторе"
А в штатах летом был взрыв пневмонии.
Я и мой муж переболели чем-то подобным в 2019 году. А сын и свекровь - в январе. Из-за сильных буранов занятий в школах не было, получился что-то типа карантина. Помню у всех нас был сильный сухой кашель, боль в груди, очень болело горло, как при ангине. Температура высокая, а вот насморка не было абсолютно.Также в январе слёг почти весь офис, где муж работает. Он один остался не заболевшим! И ещё до кучи знакомых и их детей переболело так же. Я не утверждаю что это коронавирус, но удивительно похоже.
У знакомых в Подмосковье в январе 2020 г. тяжело переболели последовательно их внук и несколько детей в его классе, а потом учительница. Первоначальные симптомы ОРВИ: кашель, температура, потом воспаление лёгких.
у вас был банальный H1N1 наверное
А я в 2017 болел, думаю, это коронавирус.
Разработали же тест на антитела к covid-19, обещали запустить исследования. Можно будет выяснить, кто на самом деле им переболел. Не исключено, что за плату, человек ради интереса сможет сделать этот тест и сам. Дали же возможность самому сдать анализ на covid. 1250 рублей за тест + рублей 600 за вызов - сравните с ценами в США.
https://www.mk.ru/science/2020/04/05/gruppa-rossiyskikh-biologov-sozdala...
Заголовок нужно поправить...
А чего с вирусом в Байконуре(территория РФ) и Антарктиде на российских станциях?
Спасибо! Хорошо делано. Откуда взяли данные по ежедневным тестам?
Спасибо. Там ссылка на сайт стопкоронавирус. Для постройки графиков использовал http://scidavis.sourceforge.net/ -- софт аналогичный Origin, но бесплатный.
Маловероятно, что будущее это не опровергнет.
Не отрицаю. Достаточно двух "точек" с большим значением и все сдвинется. Но не на сильно. Что бы сильно сдвинуть должно что то реально массовое произойти.
фиттинг? что-то типа метода наименьших квадратов?
так этот способ никак не учитывает то, что:
1. сначала вводится карантин
2. через некоторое время (~5-20) достигается пик ежедневных новых выявленных
3. ещё через некоторое время (~5-20) достигается пик активных больных (выявленных живых не выздоровевших)
4. а количество смертей в день уже скорее некоторая функция от п.3
всё имхо
Пока что это то что можно сделать быстро и в одиночку. А что бы учесть все Ваши замечания "достаточно точно", целый институт и все данные по предыдущим эпидемиям понадобятся.
и десяток лучших суперкомпьютеров мира
или просто посмотреть в Вики, что пики смертности в Италии, Испании, Франции были ~~ через неделю после пиков еженедневных новых случаев
Сейчас ситуация немного отличается. В РФ основная масса зарегистрированных больных это прибывшие из-за рубежа и их ближайшее окружение. Пока что не видно больших очагов "внутреннего" происхождения.
я откланяюсь
Картинки не показывает без авторизации.
Какая-то скрытая перепись аш-евцев.
Большой брат уже на АШ из каждого утюга.
Добавил бы запрет на такие приколы.
Картинки не показывает Опера. Хоть с регистрацией, хоть без
у меня кажет
Перестаньте натягивать сову на глобус, ей больно. Россия только начала сходить с чёткой прямой "удвоение за три дня", и после шашлычников, видимо пойдёт по прямой "удвоение за неделю", как Японии. Пика пока ещё не видно даже на горизонте. Может числу к 15.
Я специально написал "если верить результатам фитирования". Там ошибка по xc большая. Да и модель в форме Гауссиана не обязательно должна быть самой правильной. Ее я использовал интуитивно.
Снова неверная интерпретация данных. исходными являются не число заболевших в день, а число озвученых новых подтверждений по заразившимся в день.
То есть фактически, количество заразившихся в день - неизвестно, количество изолированных заразившихся - не доступно, а видим лишь количество объявленных подтверждений о ранее заразившихся. Поэтому вся дальнейшая "аналитика" просто не имеет смысла.
Поэтому, правильнее смотреть на количество умерших, эта цифра ближе к истинному значению будет. Я не стал длинно писать в подписях, но везде в данных указывается что это число зарегистрированных. Понятно что не всех регистрируют. Но если считать что процент регистрируемых не меняется, то от этого форма кривой не поменяется.
Не знаю, как где, но вот в Екатеринбурге на сегодня лаборатория, которая проводит первичный анализ на коронавирус на сегодня просто не справляется с объёмом работы, и, соответственно, статистика "плывет". Задержка составляет уже несколько дней и далее нарастает. Так что объективной информации пока нет. Есть два тренда: рост выявленных носителей с ростом проверенного народа и замедление роста выявленных, пока лаборатории захлебываются. Понятно, что сейчас в них вкинут деньги, поставят дополнительное оборудование (неделя, как мне кажется) и разгребут завалы анализов. И мы получим всплеск числа выявленных..... Так что рановато за статистику ...
Это понятно. А теперь сделайте поправку на коэффициент выявления с учетом количества тестов, и тогда, получится, что сейчас регистрируют больше (частично) не потому что их реально больше а потому что больше охват и регистрируют тех кого не регистрировали раньше, когда тестов было меньше. Поэтому приведена статистика по тестам.
График "умерших за день / заболевших за день" не имеет смысла. Чтобы видеть оперативную летальность нужно делить умерших сегодня на заболевших, но не сегодня, а со сдвигом назад на количество дней, соответствующих количеству дней болезни этих умерших.
По данным из Ухани, средний период времени от госпитализации до смерти 9.4 дня. Я бы прибавил дня 3-4 на время от обнаружения симптомов до госпитализации. И получил бы искомый период в 12-14 дней от обнаружения симптомов до смерти. А то и всё 15 - с учётом времени от обнаружения симптомов до регистрации случая как заболевшего.
Вот со сдвигом на этот период назад и нужно смотреть количество заболевших для вычисления летальности.
Тогда у вас показатели летальности перестанут быть такими бодрыми.
Это понятно, но ведь могут регистрировать не сразу, а как раз через две недели, когда человеку стало плохо, и на следующий день человек умер. Поэтому как это правильно учесть не знаю, но от этого форма кривой не поменяется сильно, только сдвинется. Но пока что ожидаемого всплеска (исходя из ваших рекомендаций) не видно. Подождем.
Вам не нужно гадать. У вас есть цифры. Это какая-никакая, но статистика. Отнимаете от текущей даты 12 дней - это будет минимальная оценка летальности. Отнимаете 15 - получаете максимальную оценку. Выводите средние показатели. Или делаете еще проще. Берете умерших за последние 4 дня и делите на сумму заболевших в диапазоне 12-15 дней назад. И так двигаетесь до даты первого умершего.
Такая оценка будет точнее 99% всех оценок, которые сейчас есть в рунете. Я гарантирую это.
Жаль, что вам не видно. То, что видно мне говорит о разнице в летальности на пол-порядка-порядок. Просто попробуйте.
То что на порядок это ясно, можно просто визуально графики сдвинуть. Но..., почему рост количества умерших за день сошел на нет? Ведь, если этот параметр запаздывает, то он должен был сейчас как раз "взлетать"?
Больше смущает следующее:
Выздоровели: Мир - 277 280 (+20 059), РФ - 494 (+88)
Умерли: Мир - 73 608 (+4 183), РФ - 58 (+11)
Процент умерших по отношению к выздоровевшим???
Мир - 25%, РФ - 12%
Для меня слишком большой процент в обоих случаях. Т.о. полагаю, что умирают в среднем к примеру если через ~10 дней, то выздоравливают тогда ~ через 14 дней (если не умрут).
Вчера сделал тест. жду результатов…
Как то надо разделить тесты которые выполнены на первичные и повторные. Много людей тестируются регулярно: доктора, аппарат правительства, президента. Думаю среди выполненных тестов за сутки первичных не больше половины.
Вывод один: пока что прём по экспоненте, переходом на плато и не пахнет
Твои глаза тебя не обманывают, потому что автор натянул одну функцию на всю статистику. Если бы апроксимация была с плавающим окном, например, 5(6-7-8-) последних дней, то перегибы были бы заметны.
Вообще то аппроксимация была свободной, без фиксации параметров. Там в указаны значения параметров и диапазон возможных значений.
И вот, специально написал: "Все это игра цифр и не более того, пока что. Но некоторые оценки и предположения можно делать."
Вот так получится:
Всегда лучше работать с полным набором данных, чем с частью.
Сравнение фитирования Гауссианом и Экспонентой. По всем "статьям" Гауссиан лучше подходит.
И по качеству R^2 -- чемближе к единице тем точнее подгонка, и по Chi^2 -- сумме квадратов отклонения, и визуально.
Да какую аппроксимацию не применяй, кривая все равно выглядит "вогнутой":)))
Понятно, что она перейдет в плато. Но не не завтра
Не данные взял. Нужно было взять количество зарегистрированных больных (т.е. заболевшие минус выздоровевшие и минус умершие)
Мне кажется, что метод фитирования неверный. Слишком много зависит от социального восприятия эпидемии, и следующим из него социальным типом поведения.
Переломом в сознании будет момент, когда начнется реальная паника, что больницы не справляются, пациентов сортируют по перспективности и ставят в очередь на ИВЛ, вместо временных больниц завозят мобильные морги. Только тогда народ опомнится и начнет соблюдать меры самоизоляции.
Поэтому сейчас важнее оценить, как медики и власти справляются со своими обязанностями, раз уж граждане свои обязанности засаботировали.
Поэтому желательно просчитать именно экспоненту, но с плавающим окном. И посмотреть, как меняется ежедневный прирост, но уже по апроксимированной кривой.
У меня такое ощущение, что экспонентциальный рост пока есть, но он немножко замедляется.
Как бы то ни было, но форма кривой будет в виде волны. Не важно по какой причине будет спад больных, по причине повальной смерти (и просто некому уже будет заболеть) или по причине остановки пандемии. Будут новые данные, будут новые кривые. Все эти фитирования верны, если ситуация кардинально не измениться по причинам, которых еще нет.