Feb. 8th, 2020

Интересный вопрос, как считать летальность уханьской заразы, т.е. какой процент больных с диагностированным заболеванием скончается. Самый просто способ - поделить число умерших на общее число больных. Сейчас больных 34939, умерло - 725. Получаем смертность 2,075%. Ну, вроде неплохо, "смешная цифра!". В 40 раз выше, чем у обычного гриппа, правда, но кого это волнует. Однако есть один нюанс. В условиях быстро растущего числа заболевших мы постоянно загоняем в знаменатель новых людей, которым еще болеть и болеть, из-за чего получаем заведомо заниженные результаты. Самый наглядный пример в модельке (link):

Представим себе вирус, который ВСЕГДА (смертность 100%) убивает носителя на 8 день вечером. Каждое утро число зараженных удваивается.

Тогда на 8 день - 256 заболевших, 1 умерший. Смертность, подсчитанная так, как у многих аффторов - 0.4%. На 9 день у нас 511 заболевших и 2 умерших. Процент опять в районе 0.4%. Далее можно проверить - ситуация будет похожая. Мы хронически недооцениваем смертельность вируса.


Есть подход с другой стороны - если известно число умерших (D, death) и выздоровевших (R, recovered), то логично что летальность (CFR) будет равна CFR=D/(D+R)*100%. Так, например, посчитали в статье на хабре (link). Но тут тоже есть системная неточность. Дело в том, что выздоравливать могут дольше. Что ж, есть статья в Ланцете (link) где наблюдали группу больных в 41 человека. Статистика так себе, но пока я лучшего не нашел. Согласно статье, медианная смертность - на 15 день. При этом есть данные на 20 день: 68% выздоровели (discharged), 17% госпитализировано, 15% скончались.

Что ж, тогда берем умерших и выздоровевших с временным лагом в пять дней (20-15) ну, и домножаем число выздоровевших на (68+17)/68=1.25 - чтобы честно учесть группу, остающуюся в госпитале (считаем, что они все выздоровеют). Итоговая формула:
CFR=D(x)/(D(x)+R(x+t)*k)*100%
где
x - день, для которого берем число умерших D
t - число дней, на которое берем поправку чтобы учесть время выздоровления
k - коэффициент, дает поправку с учетом тех, кто не выздоровел и остается в госпитале.
Взяв данные с известного карто-сайта, я получил следующие данные (x=1 февраля, x+t=6 февраля)
D(x)=259
R(x+t)=1499

тогда CFR=259/(259+1499*1,25)*100%=12,14%
Вот как-то так пока получается, и это нижняя оценка.

А теперь хорошие новости! )

З.Ы. Идея рассчета с временным сдвигом не моя, а [livejournal.com profile] optimist2008, увидел в комментариях.

Profile

nicoljaus: (Default)
nicoljaus

March 2022

S M T W T F S
   1 2 345
678 9 101112
13 1415 16 17 1819
20 2122 23 24 2526
2728 29 30 31  

Most Popular Tags

Style Credit

Expand Cut Tags

No cut tags
Page generated Aug. 19th, 2025 02:43 pm
Powered by Dreamwidth Studios