D.A. Henderson, synek Franklina i racjonalność decyzji o szczepieniu

W tych dniach zmarł D.A. Henderson, epidemiolog, który walnie przyczynił się do zlikwidowania ospy na świecie. Był to wynik wieloletniej planowej pracy zespołu ludzi, którymi kierował najpierw w amerykańskiej CDC, a później w WHO. Fachowcy mówią, że to największy wymierny sukces w historii medycyny. Dramatem naszego świata jest fakt, że ludzie tacy jak on są niezbyt znani w przeciwieństwie do różnej maści celebrytów, skandalistów i kokainistów płci obojga.  OB-Henderson__13981471621450

Pisałem o epidemii w roku 1721 w Bostonie i tragicznym losie małego synka Benjamina Franklina. Stosując rachunek prawdopodobieństwa, nietrudno uzasadnić racjonalność decyzji o szczepieniu nawet przy niepełnych danych z XVIII wieku. Musimy pamiętać, że ówczesne szczepienie, tzw. inokulacja albo wariolizacja, różniły się od późniejszej metody. Zaszczepiano bowiem ludziom ospę ludzką, co w niektórych przypadkach kończyło się śmiercią. Dopiero pod koniec stulecia Edward Jenner odkrył, że bezpieczniejsze jest zaszczepianie ludziom ospy krowiej.

Zazwyczaj w podręcznikach matematyki mamy do czynienia z urnami, z których wyciąga się kule i w zależności od tego, co wyciągniemy, pojawiają się różne możliwości i budujemy drzewo rozmaitych ewentualności. Szczepienia są przykładem lepiej chyba przemawiającym do wyobraźni niż losowania białych i czarnych kul z urny.

Oto dane dla epidemii w Bostonie w roku 1721.

  • Liczba ludności miasta: 10 700
  • Poddanych inokulacji 281, z czego 6 zmarło
  • Spośród niepoddanych inokulacji 4917 zachorowało i przeżyło, 842 osoby zachorowały i zmarły, a 4654 osoby w ogóle nie zachorowały

Będziemy prawdopodobieństwa przybliżać częstościami, zazwyczaj nie mamy na to lepszego sposobu, należy pamiętać, że dane pochodzące z niewielkiej próby mogą się okazać niedokładne i dysponując większą statystyką, otrzymalibyśmy nieco inne wyniki. Mamy więc prawdopodobieństwo zgonu po inokulacji równe 6/281=0,021 i przeżycia inokulacji 1-0,021=0,979.

Prawdopodobieństwo zgonu wśród niepoddanych inokulacji oraz zarażonych jest równe 842/(842+4917)=0,146, a prawdopodobieństwo przeżycia w tej samej grupie równa się 1-0,146=0,854.

Prawdopodobieństwo zarażenia osoby niepoddanej inokulacji możemy próbować oszacować na podstawie naszych danych jako (4917+842)/(4654+4917+842)=0,553. Jest to szacowanie z dołu: musimy pamiętać, że część spośród 4654 osób, które nie zachorowały, przeszła już kiedyś ospę i była uodporniona na resztę życia. Jeśli prawdopodobieństwo zarażenia osoby, która nie przeszła ospy, oznaczymy przez x, mamy następujące drzewo możliwości.

qc23465.f1

Rysunek z pracy M Best, A Katamba, and D Neuhauser, Making the right decision: Benjamin Franklin’s son dies of smallpox in 1736.

Jeśli przyjmiemy x=0,553, to prawdopodobieństwo przeżycia bez inokulacji będzie równe (1-x)+x \cdot 0,854=0,919. Jak widać, wartość ta jest mniejsza od prawdopodobieństwa przeżycia inokulacji, zatem statystycznie biorąc, zabieg ten zwiększa szanse przeżycia. Gdybyśmy mieli więcej informacji, wartość x mogłaby się okazać jeszcze większa, a to by oznaczało, że prawdopodobieństwo przeżycia bez inokulacji jest jeszcze mniejsze (można zapisać to prawdopodobieństwo jako 1-x+0,854x=1-0,146x, jest to więc malejąca funkcja zmiennej x).

Można też się zastanowić, jaka musi być najmniejsza wartość x, żeby inokulacja była racjonalnym zabiegiem. Granicą racjonalności będą równe prawdopodobieństwa zgonu: x\cdot 0,146=0,021, skąd x> 0,144. Ponieważ dane wskazują, że prawie na pewno ostatni warunek jest spełniony, inokulacja jest racjonalnym zabiegiem.

Nie mamy, niestety, danych dla epidemii w 1736 roku w Filadelfii, gdzie mieszkał Benjamin Franklin z rodziną. Mamy jednak dane dla późniejszej epidemii w Bostonie w roku 1752.

  • Boston liczył wówczas 15 684 mieszkańców
  • 5998 osób przeszło już ospę i nie musiało się jej obawiać
  • 2124 osoby poddały się inokulacji (znacznie więcej niż w roku 1721), 30 z nich zmarło
  • 1843 osoby uciekły na wieś, by przeczekać epidemię, nie wiemy, jak wiele spośród nich zmarło.
  • 5719 osób nie poddało się inokulacji ani nie uciekło; 97% spośród nich zachorowało, a 539 zmarło

Prawdopodobieństwo zgonu po inokulacji równe jest 30/2124=0,014; prawdopodobieństwo przeżycia: 0,986. Wartości zbliżone są do tego, co otrzymaliśmy wyżej dla roku 1721.

Wśród niezaszczepionych i narażonych na zachorowanie śmiertelność była równa 539/(0,97\cdot 5719)=0,097, prawdopodobieństwo przeżycia choroby równało się 1-0,097=0,903. Oznaczało to, że nie robiąc nic, ma się prawdopodobieństwo przeżycia 0,03+0,97\cdot 0,903=0,906. Należy porównywać to z wartością 0,986 dla zaszczepionych. Inokulacja była więc znacznie lepszą decyzją.

Statystyka z roku 1752 obejmuje jeszcze możliwość ucieczki z miasta. Była to najprostsza metoda unikania chorób epidemicznych i kogo było na nią stać, ten ją stosował. Nie znamy prawdopodobieństwa zachorowania wśród tych, co uciekli. Oznaczmy je przez y. Mamy więc następujące drzewo możliwości.

qc23465.f2

(Rysunek z pracy jw.)

Można zadać pytanie, jakie powinno być y, aby ucieczka była lepszym wyjściem niż pozostanie w Bostonie i poddanie się inokulacji. Prawdopodobieństwo zgonu osoby uciekającej to 0,097y, należy je porównać z prawdopodobieństwem zgonu po inokulacji, równym 0,014. A zatem, jeśli y< 0,144, to ucieczka jest racjonalna. Trudno jest oczywiście oszacować wartość y, zależy ona np. od tego, czy uciekniemy, zanim jeszcze epidemia się rozwinie, czy w jej późniejszej fazie (choroba ma pewien okres inkubacji, możemy więc wyjeżdżając czuć się dobrze mimo zarażenia). W dodatku uciekając, nadal nie mamy odporności na ospę, a w Bostonie w ciągu osiemnastego wieku większe epidemie wystąpiły w latach 1721, 1730, 1752, 1764, 1776, 1778 oraz 1792. Można się było spodziewać, że za kilkanaście lat choroba znów się pojawi.

Skomentuj

Wprowadź swoje dane lub kliknij jedną z tych ikon, aby się zalogować:

Logo WordPress.com

Komentujesz korzystając z konta WordPress.com. Log Out / Zmień )

Zdjęcie z Twittera

Komentujesz korzystając z konta Twitter. Log Out / Zmień )

Facebook photo

Komentujesz korzystając z konta Facebook. Log Out / Zmień )

Google+ photo

Komentujesz korzystając z konta Google+. Log Out / Zmień )

Connecting to %s