Po době kladoucí důraz na co nejtvrdší data, co nejlepší soubory dat a co nejlepší a nejexaktnější hodnocení těchto dat přichází i čas na ohlédnutí. A tam má místo i slovesnost. Zachycení příhod, událostí, jednání hrdinů různých velikostí je možné a potřebné ve formě vyprávění může zpětně popsat to, co se do tabulek rozhodně nevejde. Pro poučení bude potřebné psát knihy, aby bylo co číst a porovnávat pro tvorbu historické paměti.
Takovým příspěvkem k reflexi je nová kniha Michala Pohludky, která naplňuje charakter literárního, nikoli faktografického záznamu, což ji činí zvláště cennou. Faktografie je oslabena i cíleně: neuvádí žádná jména, ani osob, ani měst či institucí. Získává tím svobodu napsat vše naplno, aniž by se musel ohlížet, zda někomu nešlápne na kuří oko.
Autor je přímou součástí děje, jde totiž o jeho zážitky a dojmy a nelze jej ze zachyceného vyjmout. Do textu se tudíž promítá autorova trojí kvalifikace: jednak specialista na laboratorní metodu, která zažila během pandemie opravdovou explozi, na polymerázovou řetězovou reakci. Za druhé je to zkušený produktový specialista s praxí v nadnárodních firmách, tudíž i školitel, rádce i obchodník, který pomohl v mnoha laboratořích ve zlomových chvílích. A nakonec je to i korporátem proškolený vyjednavač. V poslední roli působil hlavně v řadě jednání se státními orgány, úředníky a politiky.
Všechny tři roviny zachycuje v sérii vzpomínek, spíše fejetonového charakteru, občas připomíná báseň ve volném verši, ovšem každá kapitola má dramatický oblouk s různě silně vyjádřenou pointou.
V zobrazených motivech se rozpaky nad cílením plošného testování a zaváděných opatření, někdy ústící do frustrace z jednání na nejvyšších úrovních mísí se zápalem biochemika, zakládajícího firmu na dodávku řešení pro laboratoře, na něž jsou náhle kladeny nevídané požadavky, vyvíjejícího se spolupracovníky odběrový systém či pracovní roztok, v pozdních nočních hodinách řešícího zádrhele a spícího v prázdných hotelech.
Ve vrcholících vlnách pak převažuje adrenalin rostoucích počtů zvládnutých analýz, vyřešených problémů a neotřelých řešení (například když angažovaní biochemici odvážejí na dobré slovo půjčené výzkumné přístroje z akademických pracovišť, repasují dávno odložené stroje a hledají řešení na dálku pomocí osobních kontaktů). Další frustrace přicházejí, když laboratoře s vychytaným systémem přecházejí na čínskou levnější diagnostiku. Nebo když mediálně hyperaktivní odborníci v řadě případů hájí vlastní obchodní zájmy. Nebo černá beznaděj, když si vzpomene na populistická nařízení, která podle všeho znamenala mnoho zbytečných úmrtí.
Dojem z knihy však nekončí zklamáním. Pohludka si je vědom práce během epidemie poctivě odvedené a rovněž ředí dramatické chvíle humorem a paradoxními situacemi, které by bylo možno snadno dramatizovat jako veselohru.
Autor je sice slovesně nadaný, ale nikoli rutinér. Z jazyka je trochu cítit, že řadu let spíš četl a psal v angličtině. Textu by jistě prospěla důkladnější redakce. Nicméně poetická místa prozrazují výborný jazykový cit. To vše bohatě vyvažuje časnost uveřejnění textu a předložení čtenářům v době, kdy máme sami své vlastní zážitky v čerstvé paměti. Díky za to.
Před týdnem jsem vyhlásil malou prognostickou soutěž o velké pivo. Zájemci mohli predikovat / tipovatvýši meziroční změny počtu živě narozených dětí za rok 2022 oproti roku 2021. Známo bylo, že za období leden až září 2022 šlo o meziroční pokles o 9,5 %, takže bylo možno očekávat, že celoroční změna se bude pohybovat právě kolem této hodnoty.
Soutěžní predikce
Zaznamenal jsem celkem 22 predikcí (5 emailem, 17 na facebooku). Predikce se pohybovaly od 7,3 % do 12,6 %. Jen 5 účastníků predikovalo, že výsledný celoroční pokles bude menší než 9,5 %, naopak 16 účastníků predikovalo ještě větší propad. Kdybychom seřadili všechny predikce od nejmenší po největší a vzali hodnoty „uprostřed“, dostali bychom hodnoty 10,43 a 10,7, což jsou (jak se za chvíli ukáže) predikce velmi blízké skutečné hodnotě meziročního poklesu. Že by se projevila moudrost davu?
Správná odpověď byla 10,70 % Meziroční pokles v posledních třech měsících roku 2022 byl (bohužel) větší než „průměrných“ 9,5 % v měsících předešlých. Dnes zveřejněná data ukazují, že ve srovnání se stejnými měsíci roku 2021 došlo v posledních měsících roku 2022 meziročně k poklesu o 11,5, 14,1 a 14,2 procent. Celkový počet živě narozených za rok 2022 byl 99 834, což představuje pokles o 10,70 procent oproti roku 2021, kdy se narodilo 111 973 dětí.
Ve svých předešlých textech jsem ukázal, že na vrub změny četnosti a věkové struktury populace žen můžeme přičíst jen asi 2,6% pokles. Pokles o 10,7 % znamená velký pokles úhrnné plodnosti.
Kdo dobře predikoval
Nejpřesněji predikoval Aleš Dumek, který tipoval pokles o 10,7 % (v absolutních číslech 11 975, přičemž reálně byl pokles o 11 959). Za zmínku stojí skutečnost, že druhou nejlepší predikci zveřejnila na svém facebookovém účtu Zuzana Krátká (10,77 %) a vysvětluje v něm i to, jak by mohl pokles porodnosti souviset s vakcinací. Blízko se svými predikcemi byli i Ivo Petrouš (10,43 %), Vladimír Vladaxy (10,8 %) a Marek Hubert (10,9 %). Vítěz nechť se přihlásí a domluvíme se, kdy a kde si vyzvedne slíbené velké pivo. Já vzhledem k tomu, že meziroční pokles je nakonec větší než 10 %, pociťuji hořkost už nyní. Až nyní mi dochází, že hlavní cena byla už od počátku hořká. Souvislost mezi výší poklesu a stupni piva je nezamýšlená.
Imunita získaná po infekci spolehlivě chrání před těžkým průběhem reinfekce. Výsledky metaanalýzy, která sledovala efekt imunity získané po prodělání infekce, byly publikované nedávno v časopise Lancet. Souhrn 65 publikací z 19 různých zemí ukazuje, že ochrana proti těžkému průběhu případné reinfekce dosahuje přibližně 90 procent a přetrvává téměř rok. Tomáš Fürst o ní píše v tomto článku. Proč ale vydání metanalýzy trvalo tak dlouho? Proč se odborníci tak dlouho bránili uznat zcela zjevnou skutečnost, že lidé, kteří prodělali onemocnění, nepotřebují očkování?
Diskuse o tom, koho v sezóně 2022/2023 očkovat, proběhla na půdě Národního institutu pro zvládnutí pandemie již v dubnu 2022 . Kolektiv odborníků spolupracujících se SMIS doporučil jen očkování rizikových skupin, ale vedení NIZP doporučilo očkovat třetími a čtvrtými dávkami všechny skupiny bez omezení. Ministr Válek poté vyhlásil výběrové řízení na reklamní kampaň a při zadávání zakázky zveřejnil vskutku ambiciózní cíle počtu očkovaných osob.
V grafu porovnáváme jak cílové počty dávek podle přání pana ministra resp. jeho poradců, tak počty skutečně aplikovaných dávek evidovaných na stránkách MZČR (počty do 9.3.2023). V zadávací dokumentaci se uvádí: „Zadavatel dále očekává, že kampaň přispěje ke změně chování občanů ČR, jejímž výsledkem budou zvýšená aplikace očkování proti covid-19 a to v následujícím očekávaném rozsahu…“ a zveřejňují se počty pro každou kategorii zvlášť ve dvou variantách – jako ambiciózní „cíl“, tak méně ambiciózní „úspěch“ . Např. u vakcinací dosud neočkovaných osob bylo definovaným „cílem“ více jak 187 tisíc osob, „úspěchem“ téměř 75 tisíc osob. Realitou, které se podařilo dosáhnout, bylo 8 tisíc osob očkovaných od srpna 2022 do března 2023, resp. 12,4 tisíc osob, započteme-li všechny nově očkované od července 2022, kdy se spustilo očkování na novou sezónu.
Druhá dávka – https://onemocneni-aktualne.mzcr.cz/vakcinace-cr
Zajímalo by mne, jak je pan ministr spokojený s úspěšností reklamní kampaně. Připomeňme si, že byla spuštěna 15. srpna 2022 a její hlavní část startovala nakonec až v říjnu 2022. Výsledek? Od půlky srpna 2022 do počátku března 2023 bylo podáno 800 tisíc dávek, tedy přibližně čtvrtina plánovaného počtu – „úspěchem“ by byly aplikované 3 milióny dávek a „cílem“ bylo aplikovat více než 4 milióny dávek. Pokud započítáme i dobu před kampaní, tak bylo podáno cca 952 tisíce dávek. Novou bivalentní vakcínou působící i proti omikronu, která byla k dispozici od září 2022, bylo naočkováno pouze 551 tisíc osob, tedy necelých 8 % z těch, kteří na ní měli nárok – tedy z osob, které obdržely základní dávku původními vakcínami.
Všechny dávky – https://onemocneni-aktualne.mzcr.cz/vakcinace-cr
Ukázalo se, že kdo nebyl očkovaný v sezóně 2021/2022, nově se očkovat nenechal. Buď onemocnění prodělal nebo očkování mRNA vakcínou odmítá a odmítat bude prostě z principu. A potvrdilo se, že ochota lidí nechat se opakovaně očkovat s každou další dávkou klesala. Zda polevil strach z onemocnění či zda převládla skepse z účinnosti a lidé měli obavy z bezpečnosti vakcíny, to už se nedozvíme. Nicméně pro ministerstvo zdravotnictví jde o jasný signál, že plány na další očkování je třeba dramaticky seškrtat, jelikož jsou zcela odtržené od reality.
Vážení čtenáři, dovoluji si nabídnout Vám malou prognostickou soutěž.
K jejímu vypsání mě inspiroval nejmenovaný kolega (fyzik), který rád pátrá po tom, jak funguje svět kolem nás, a pro jeho lepší pochopení se nebojí činit různé predikce, o jejichž splnění je ochoten se vsadit.
V nedávné době mě zaujala problematika poklesu počtu narozených dětí v roce 2022. Zatím víme, že v období ledna až září roku 2022 se v ČR narodilo 68 505 živě narozených dětí, což představuje meziroční pokles o 9,5 procent, neboť ve stejném období roku 2021 se narodilo 75 706 živě narozených dětí. Počty živě narozených za poslední tři měsíce roku 2022 by měly být zveřejněny 21. března 2023, tedy příští úterý.
Cílem soutěže je odhadnout co nejpřesněji meziroční změnu počtu živě narozených dětí za celý rok 2022 oproti roku 2021.
Potřebná data k této malé prognostické soutěži najdete na stránkách Českého statistického úřadu, konkrétně na konci sloupce E tabulky 1 (řádky 583-591 obsahují údaje o počtu živě narozených za měsíce leden až září 2022). Grafickou nápovědu naleznete na konci tohoto článku.
Doplním zde pár svých úvah k tomuto tématu:
Ve veřejném prostoru se objevilo několik hypotéz vysvětlujících meziroční pokles o 9,5 procent v prvních devíti měsících roku 2022. I kdyby se v posledním čtvrtletí roku 2022 narodilo stejně dětí jako za poslední tři měsíce roku 2021 (tedy skokově by se po propadu o 9,5 procent hodnoty vrátily na 0 procentní změnu), byl by celkový meziroční pokles o 8,1 procent. Ve skutečnosti však zřejmě bude ještě větší.
Faktor 1: Četnost a struktura populace žen v reprodukčním věku
První hypotézou je, že k poklesu došlo výhradně z důvodů demografických, že totiž vyplývá ze změny četnosti a struktury populace žen v reprodukčním věku. Tuto hypotézu prezentoval Petr Koubský v Deníku N, když napsal: „Pokles počtu nově narozených dětí je nesporný, ale k jeho vysvětlení stačí vzít v úvahu, jak klesá počet žen v reprodukčním věku, další mechanismy nejsou nutné.“
Pokud by tomu tak bylo a zůstala zachována tzv. úhrnná plodnost, tedy byla by zachována „intenzita procesu rození“, měl by dle mých výpočtů meziroční pokles počtu narozených dětí činit 2,6 procent. Aby však bylo této hodnoty za celý rok dosaženo, musely by být poslední tři měsíce výrazně v plusu oproti roku 2021. To ale asi očekávat nelze. Koubského hypotézu tudíž odmítám.
Výše uvedený pokles počtu narozených o 2,6 procent předpokládal zachování úhrnné plodnosti, tedy intenzity procesu rození dětí. Jelikož skutečný pokles bude zřejmě větší, je už nyní evidentní, že v roce 2022 došlo k snížení úhrnné plodnosti. Otázkou je, proč k tomu došlo? Připomeňme, že změnu úhrnné plodnosti není možno zdůvodňovat změnou četnosti a struktury populace žen v reprodukčním věku.
Příčiny snížení úhrnné plodnosti
Hypotézu, která pokles úhrnné plodnosti vysvětluje, prezentovala demografka Klára Hulíková Tesárková. Ta uvádí, že „po roce mimořádného zvýšení plodnosti přijde spíše rok poklesu, neboť tento výkyv mohl být změnou časování. Ženy, které by možná měly dítě až v roce 2022, ho měly již v roce 2021. Tím došlo k růstu intenzity plodnosti v roce 2021, ale logicky na úkor roku 2022.“
Úhrnná plodnost opravdu dosáhla v roce 2021 rekordní úrovně (v kontextu uplynulého desetiletí). Po třech letech, kdy docházelo k mírnému poklesu v počtu narozených dětí (meziročně o 0,3; 1,6 a 1,8 procent) při mírně rostoucí úhrnné plodnosti, došlo v roce 2021 k meziročnímu nárůstu počtu živě narozených o 1,4 % a úhrnná plodnost tak zaznamenala zřetelný růst. Výše uvedená hypotéza tedy má opodstatnění.
Otázkou je, jak velký pokles úhrnné plodnosti bylo možno na základě této hypotézy očekávat?
Nápad 1:
Někdo by mohl očekávat návrat k hodnotám z roku 2020 (rok před „mimořádným“ rokem 2021). Pokud bychom se přidrželi oficiálně publikovaných odhadů, dospěli bychom k tomu, že očekávaný meziroční pokles počtu narozených mezi roky 2021 a 2022 je 9,5 procent (při předpokládané úrovni plodnosti z roku 2020 dle oficiálních statistik). Z tohoto pohledu by tedy bylo rozumné očekávat, že také v posledním čtvrtletí roku 2022 bude pokles oproti stejným měsícům roku 2021 přibližně 9,5 procent, neboť bude-li za měsíce říjen až prosinec pokles oproti stejným měsícům roku 2021 opět o 9,5 procent, bude i celkový meziroční pokles o 9,5 procent. Jenže skutečná úhrnná plodnost byla zřejmě vyšší, než bylo v roce 2020 odhadnuto, jak vysvětluji ve svém článku o opomenutých aspektech této problematiky.
Při aktuální četnosti a věkové struktuře populace žen v reprodukčním věku by podle mých hrubých propočtů návrat k hodnotám úhrnné plodnosti z roku 2020 či dokonce z roku 2019 znamenal v roce 2022 meziroční pokles o 5 až 6 procent. Toto je odhad meziročního poklesu při předpokladu současné velikosti a struktury populace žen a návratu úhrnné plodnosti na hodnoty z let 2019 či 2020.
Nápad 2:
Vezmeme-li výše uvedenou hypotézu doslovně, dostaneme následující model: předpokládejme, že děti narozené v roce 2021 se dělí do dvou skupin: na „základní“ a „předrozené“ (při běžných podmínkách by se narodily až v roce 2022). Kolik je těch „předrozených“? Pokud by úhrnná plodnost byla i v roce 2021 na úrovni roku předešlého (2020), narodilo by se v něm přibližně 107 650 dětí, tedy o přibližně 4150 méně než ve skutečnosti. To je tedy počet „předrozených“ dětí. Pokud by se tyto děti narodily až v roce 2022 a přitom byla úroveň úhrnné plodnosti i v roce 2022 stejná jako v roce 2020, očekávali bychom (vzhledem k současné populaci žen v reprodukčním věku) narození celkem asi 105 850 (živě narozených) dětí za rok 2022. To odpovídá výše zmíněnému meziročnímu poklesu o cca 5,3 %. Pokud odpočítáme „předrozených“ cca 4150 dětí, dostáváme se k počtu cca 101 700 živě narozených, což by představovalo meziroční pokles o 9,0 % (pokud místo úhrnné plodnosti z roku 2020 vezmeme tu z roku 2019, dojdeme k odhadu velikosti meziročního poklesu 9,9 %.) Připomeňme, že v tomto scénáři bude po „nadprůměrném“ roce 2021 následovat rok nikoliv průměrný, ale „podprůměrný“, protože dojde k bezprostřednímu vyrovnání „schodku“ z roku 2021, načež by se v roce 2023 měla úhrnná plodnost opět zvýšit. Právě bezprostřednost vyrovnání zmíněného „schodku“ (jde o přesun pouze mezi roky 2021 a 2022) je dosti silným předpokladem tohoto modelu.
Závěry:
Existuje poměrně věrohodné vysvětlení skokového meziročního propadu počtu narozených dětí o 9-10 procent, a to v kombinaci dvou faktorů – změny v četnosti a struktuře populace žen v reprodukčním věku (odpovídající meziročnímu poklesu o 2,6 %) a změny časování těhotenství ve prospěch roku 2021 na úkor roku 2022, který se projevil nárůstem úhrnné plodnosti v roce 2021 a výrazným poklesem této míry v roce 2022.
Samotné změny v četnosti a struktuře populace žen v reprodukčním věku k vysvětlení poklesu o 9-10 procent v žádném případě nestačí. Dominantní faktor určující tuto změnu je třeba hledat jinde.
V roce 2022 došlo nejen k výraznému poklesu počtu živě narozených dětí (který je zčásti dán změnami v četnosti a struktuře populace žen v reprodukčním věku), ale také k výraznému poklesu úhrnné plodnosti. Existence jednoho možného vysvětlení pozorovaného poklesu ani zdaleka nevylučuje působení dalších faktorů.
Pokud v pozorovaném 9,5% poklesu počtu živě narozených dětí hrálo nějakou roli očkování probíhající v roce 2021 (a nemusí jít přitom o efekt biologický, ale např. o efekt psychologický), šlo jen o jeden z působících faktorů. Celý pokles rozhodně nemůže být tomuto faktoru přičítán.
Slíbená grafická nápověda:
Následující graf znázorňuje vývoj meziroční změny počtu živě narozených. Porovnáván je vždy daný měsíc se stejným měsícem předchozího roku a toto srovnání je znázorněno modrými body. Červené vodorovné úsečky pak znázorňují meziroční srovnání celkových počtů – daný rok je vždy porovnán s rokem předchozím.
Cílem soutěže je odhadnout hodnotu pro červenou úsečku za rok 2022, jehož údaje jsou zatím nekompletní. Své tipy můžete zasílat na email ondrej.vencalek@upol.cz nebo zveřejnit v diskusi u našeho FB příspěvku k tomuto článku.
Velmi nás těší, že naprostá většina diskutujících pochopila, o co jde v jednoduchém modelu od Tomáše Fürsta popsaném v článcích “Znamení slunce aneb observační studie a selection bias” a “Znamení slunce útočí”. Předtím, než se pustíme do dalšího vysvětlení pro ty, co mají z nějakého důvodu problém nezaujatě pochopit podstatu popsaného efektu, popřípadě skriptu v R, si přečtěte následující:
CO TO JE: Grafy jsou realizace jednoduchého modelu, na kterém je ukázán efekt výběrového vychýlení a efekt špatné klasifikace. Model nepotřebuje ŽÁDNÁ řádková data, časové řady ani nic podobného.
CO TO NENÍ: V žádném případě se nejedná o výpočet efektivity vakcinace na základě reálných dat, skript nepotřebuje žádný externí soubor, pouze pár knihoven jazyka R.
Model má několik vstupních PARAMETRŮ, které si můžete nastavit (aktuální skript v jazyce R je ke stažení zde https://github.com/strakrob/znameniSLUNCE, R a Rstudio lze stáhnout zde https://posit.co/download/rstudio-desktop/) dle libosti. Potřebujete 5 – 10 parametrů (záleží, jaký efekt chcete modelovat), všechny jsou okomentovány spolu s celým zdrojem pro usnadnění práce se skriptem a případně jeho další úpravou pro zájemce. Nastavení parametrů pro jednotlivé grafy z výše uvedených článků je v nich důsledně popsáno. Skript nedělá nic jiného, než co je popsáno v článcích, tj., v prvním z článků:
Vybere náhodně části populace, které zemřou a které dostanou Znamení slunce. Z konstrukce modelu i jeho implementace v R je zřejmé, že Znamení nemá v tomto případě žádný vliv na úmrtí, což je vidět i v grafech.
U efektu výběrového vychýlení odstraníme udělení Znamení slunce (s určitou pravděpodobností různou od 0), ale pouze pokud k němu dojde 14 dní před úmrtím, na grafech se objeví virtuální účinnost Znamení – tím větší, čím větší je zmíněná pravděpodobnost.
V druhém z článků, kde určitá malá část Znamení slunce zabíjí:
Vybere náhodně určitou část populace se Znamením a nechá jí náhodně umřít mezi 0. a 14. dnem od udělení Znamení, na grafech je tento efekt patrný od samého začátku.
Totéž, co v 1), ale špatně klasifikujeme VŠECHNY jedince se Znamením tak, že je počítáme jako se znamením až od 14. dne po jeho udělení. Zvýšené riziko nešťastníků z bodu 1) zmizí díky tomu, že oficiálně umírají jako jedinci bez Znamení (tady si můžete pohrát s nastavením parametrů a sledovat výsledek).
Totéž, co v 1), ale špatně klasifikujeme POUZE jedince, kteří zemřeli do 14 dnů od udělení Znamení (prohlásíme je za jedince bez Znamení) a rázem se nám objeví úplně opačný výsledek než v bodě 1).
Tento jednoduchý model poukazuje na úskalí interpretace výsledků na základě observačních dat, kde může lehce dojít k různým typům vychýlení až k nepoctivé klasifikaci, přičemž v tomto modelu není důležité zda při “prevenci” jde o Znamení slunce, vakcíny, léky, nebo červené svetry. Na druhou stranu nás trochu mrzí, že členové a skalní fanoušci ze sdružení Sníh (včetně těch, kteří vystupují jako matematici) mají problém s kontrolou tak jednoduchého modelu a pochopení skriptu na pár desítek řádků. Nic však pro ně ještě není ztraceno, pokud mají zájem (i kdokoli další), každou středu od 9:30 do 11:30 mám konzultace a rád vše zájemcům vysvětlím (online, pokud zrovna nebudete cestovat do Krakowa), email lze dohledat na stránkách AGH (https://skos.agh.edu.pl).
Vypadá to, že naše země není jediná, kde se ke zdůvodnění předem rozhodnutých protiepidemických opatření používala chybná či zmanipulovaná data. Přečtěte si v The Telegraph frustrující článek od Carla Heneghana, který líčí, jak v Downing Street marně vysvětloval, že data, kterými chtějí britští politici zdůvodnit další lockdown, jsou chybná a zastaralá.
Pak si vzpomeňte, že u nás používali politici přesně stejnou strategii. Lockdown zdůvodňovali nesmyslným modelem z dílny ÚZISu, který byl tak absurdní, že jsme mu v dubnu 2021 věnovali samostatné pojednání s názvem Úsvit šamanů. Později na ÚZISu vykouzlili naprosto zcestnou „evidenci“ k uzavření škol, za kterou jsme Duškův ústav překřtili na Ústav pro zkreslování informací a statistiky.
Zatímco obyvatelé Britského království se své covidistické věrchušky dávno zbavili, u nás se nezměnilo vůbec nic. Ve vládě vystřídali jedni šílenci druhé a ÚZIS pokračuje v manipulaci s daty jako by se nechumelilo.
Holt každý národ má takové instituce, jaké si zaslouží.
Znamení slunce už nám jednou posloužilo k ilustraci mimořádně zajímavého fenoménu, který je samozřejmě zcela fiktivní a jakákoliv jeho podobnost s reálnými daty je čistě náhodná. Dnes Znamení probudíme znova a ukážeme další zajímavost, která je ještě více fiktivní a jakákoliv její podobnost s realitou – a to zejména realitou popsanou zde – je ještě více náhodná.
Představme si populaci deseti milionů lidí, kterou budeme sledovat po 365 po sobě jdoucích dní. Náhodně vybereme setinu populace, která během sledované doby zemře. Čas úmrtí každého nešťastníka je vybrán náhodně s rovnoměrným rozdělením přes celý sledovaný rok. Dále náhodně vybereme asi 65 procent populace, která si někdy během roku dojde na Hlavní nádraží pro Znamení slunce (to je malý obrázek sluníčka, který je zájemcům nakreslen henou na rameno). Čas udělení Znamení slunce je také náhodně vybrán s rovnoměrným rozdělením přes sto sledovaných dnů.
Dále předpokládejme, že jeden z tisíce sluníčkových do 14 dnů od obdržení Znamení zemře (jedno z tisíce Znamení je totiž ve skutečnosti Znamením smrti). Vybereme tedy náhodně tisícinu sluníčkových a necháme je někdy během 14 dnů od udělení Znamení zemřít. K přirozeným 100 tisícům úmrtí tedy přibyde během sledovaného období ještě asi 6 500 úmrtí po obdržení Znamení slunce (smrti).
Pro každý den spočteme úmrtnost sluníčkových (tedy procento lidí se Znamením, kteří daný den zemřeli) a úmrtnost nesluníčkových (procento lidí bez Znamení, kteří daný den zemřeli). Není nijak překvapivé, že sluníčkoví umírají více a efekt je nejvíce patrný na začátku sledovaného období, kdy je sluníčkových ještě málo. Graf ze sta opakování výše popsaného modelu je níže. Je zobrazen sedmidenní centrální klouzavý průměr denní úmrtnosti, abychom mírně potlačili šum.
A teď se podívejme, co se stane, když se někdo rozhodne, že Znamení slunce začíná platit až 14 dnů od udělení. Všechny, kteří mají Znamení slunce kratší dobu než 14 dnů, budeme tedy považovat za neoznačené. Pokud provedeme tuto jednoduchou účetní operaci, dopadne vyhodnocení výše uvedeného modelu takto:
Sluníčková mortalita je rázem pryč, a ještě k tomu to v druhé polovině roku vypadá, že Znamení slunce před úmrtím mírně chrání!
A když už jsme u toho účetnictví, můžeme se podívat, jak by se projevila ještě jedna změna. Za zemřelé se znamením budeme opět počítat až ty, kteří jsou alespoň 14 dnů od znamení. Ovšem celkové množství lidí se znamením je počítáno dle původní metodiky, to jest ihned od udělení znamení. Tento typ „kreativního účetnictví“ skončí následovně:
Dne 16. ledna 2023 jsme se Ústavu zdravotnických informací a statistiky (ÚZIS) dotázali, proč došlo v reportech, které ÚZIS zveřejňuje, k dramatickým změnám se zpětnou platností. Tyto změny podrobně popisujeme zde. ÚZIS nám odpověděl 8. února 2023 velmi stručným textem, kde sděluje následující:
1. ÚZIS potvrzuje, že k námi popsaným změnám skutečně došlo.
2. ÚZIS tvrdí, že v databázích se nic nezměnilo, změna se udála pouze v tom, jak se z databáze generuje každodenní report, který je k nalezení zde jako příloha 4. Zde je potřeba dodat, že právě příloha 4 byla celou dobu tím zásadním dokumentem, ze kterého veřejnost, novináři, politici a jiní čerpali podklady pro svá tvrzení o účinnosti vakcín či obsazenosti nemocnic.
3. ÚZIS sděluje, že nově (tedy po úpravě generování reportů) začal lidi do 14 dnů po první dávce považovat za neočkované. Tím vysvětluje naše pozorování, že 22. října 2022 najednou z přílohy 4 zmizelo mnoho covidových úmrtí u jedinců očkovaných první dávkou. Abychom byli konkrétní, před úpravou dat byly hodnoty v kategorii „Zemřelí po první dávce“ (myšlena jsou jen úmrtí s covidem) 72 lidí v lednu 2021, 210 lidí v únoru a 565 lidí v březnu. Po úpravě dat se tato čísla změnila na 13, 104 a 236. Pokud by byla pravda to, co ÚZIS tvrdí – tedy že lidé, kteří jsou méně než 14 dnů od první dávky, byli přeřazeni do kategorie neočkovaných – musela by adekvátně poklesnout i celková velikost kohorty očkovaných prvních dávkou. To se ovšem nestalo. Před úpravou dat byla čísla v kategorii „Průměrný denní počet osob v kategorii s 1 dávkou v daném intervalu“ za leden, únor a březen 2021 tato: [121 912, 252 326 a 605 128]. Po úpravě dat se tato čísla změnila na [121 817, 252 126 a 604 616]. Je tedy zřejmé, že v této kategorii k odpovídající úpravě nedošlo. Jinými slovy, z této kategorie byli odstraněni v podstatě pouze zemřelí (potud je informace ÚZISu pravdivá), ale téměř nikdo další. Rozhodně nebyli odstraněni všichni, kteří byli do 14 dnů po první dávce. Vysvětlení ÚZISu je tedy prokazatelně nepravdivé.
Pro lepší porozumění nabízíme následující příklad. Řekněmě, že 1. ledna 2021 máme sto tisíc neočkovaných lidí. Dne 18. ledna jich tisíc naočkujeme. Deset z nich v následujících 14 dnech zemře. Úmrtnost očkovaných v lednu by se tedy logicky měla spočítat jako 10 / 1000. ÚZIS ve své odpovědi tvrdí, že dle nové metodiky by v tomto případě všech 1000 očkovaných považoval až do konce ledna za neočkované, tedy úmrtnost očkovaných v lednu by byla 0 / 0 neboli nedefinovaná. Ovšem podle změn, které ÚZIS v datech ve skutečnosti provedl, to vypadá, že za neočkované považoval jen těch 10 mrtvých, ale všech 990 přeživších nechal v kategorii očkovaných. Úmrtnost by tedy vyšla 0 / 990. To jsou dramatické rozdíly.
V mezičase (mezi 8. únorem 2023 a dneškem) se data v příslušném reportu opět trochu změnila (důvod nám nebyl sdělen), ani současné hodnoty však neodpovídají ÚZISem popsanému stavu.
4. Dále ÚZIS píše: “zároveň bohužel došlo ve výstupu k chybě, kdy blok výsledku pro JIP ve věku 75+ byl zobrazen i v sekci 65+, to je již po předchozím upozornění uživatelů, za něž děkujeme, opraveno.” V tomto bodě souhlasíme s vysvětlením ÚZISu, považujeme ovšem takové chyby za mimořádně nešťastné. Jde o tisíce pacientů, kteří se přesouvají tam a zpět, a to právě v kolonkách, podle kterých se v této zemi řídila nejdrastičtější opatření v jejích dějinách.
5. Zbylý text odpovědi ÚZISu sděluje, že popsané rozdíly nehrají roli. ÚZIS píše „U počtů po dokončeném očkovaní, ze kterých se počítají ochranné efekty, jsou rozdíly v počtech eventů minimální, např. v únoru 2021 […] u pozitivity to je cca 5 % eventů, což se do ochranného efektu promítá jako 89,3 % vs. 89,7 %, což je minimální rozdíl.“ Toto sdělení je nicméně irelevantní, neboť – jak vysvětlujeme výše – pokud nejsou správně spočtené velikosti kohort oba „ochranné efekty“ citované ÚZISem výše nemají nic společného s realitou. To, že se tyto dva odhady vzájemně liší jen málo, je tedy zcela irelevantní.
Dále konstatujeme, že data o počtech očkovaných v příloze 4 neodpovídají datům na stránce „Přehled vykázaných očkování“, a to ani po zahrnutí „vakcinačního okna v délce 14 dnů“. Je také pozoruhodné, že dle Českého statistického úřadu kolísal celkový počet obyvatel od ledna 2021 do září 2022 mezi 10 492 910 a 10 527 732. ÚZIS ovšem vykazuje kolísání mezi 10 701 749 a 10 701 760.
Na náš formální dotaz dle zákona č. 106/1999 Sb., o svobodném přístupu k informacím nám tedy ÚZIS poskytl odpověď, která je prokazatelně nepravdivá. Pokud byla mezi neočkované přesunuta část očkovaných, která zemřela (s covidem), ale ne odpovídající část očkovaných, která nezemřela, jedná se o manipulaci s daty, která je nepřípustná a vede k chybným odhadům účinnosti (a bezpečnosti) vakcín proti covidu. Nejde o první selhání tohoto úřadu. Již dříve jsme upozorňovali na zásadní nesrovnalosti v činnosti ÚZIS (viz například zde, zde, zde, zde nebo zde).
Odpověď ÚZISu tak jen zdůraznila důležitost našeho opakovaného požadavku, aby byl nezávislým datovým analytikům umožněn přístup k anonymizované datové sadě na úrovni jednotlivců, včetně kalendářního data všech covidových očkování, data a výsledku všech PCR testů, data hospitalizací a data úmrtí z jakýchkoliv příčin. Nezávislí analytici by tak mohli důsledně ověřovat tvrzení ÚZISu a ministerstva, která se v minulosti ukázala tak často být nepravdivá či zavádějící.
Tři roky od začátku epidemie je nejvyšší čas začít krok po kroku ze skutečných dat rekonstruovat, co se v této zemi během pandemie dělo. Momentálně máme jedinou jistotu – obraz podávaný analytiky ÚZISu neodpovídá realitě.
V uplynulých přibližně dvou týdnech proběhla čilá debata o vývoji počtu narozených dětí a plodnosti v České republice v posledních letech. Zkusím nejprve tuto debatu krátce rekapitulovat a poté poukázat na některé zatím opomíjené aspekty této problematiky. Zejména se zaměřím na hodnocení roku 2021 jako roku s (výrazně) zvýšenou plodností, po níž by bylo možné očekávat pokles („návrat k průměru“). Zjistíme totiž, že plodnost nebyla zvýšená až tak výrazně, jak se může na první pohled zdát. Překvapivě důležitou roli totiž hraje skutečnost, že v roce 2021 proběhlo sčítání lidu.
Rekapitulace únorové debaty o vývoji počtu narozených dětí a plodnosti
Debatu o vývoji počtu narozených dětí (a potažmo plodnosti) v ČR odstartoval článek Tomáše Fürsta z 6. února 2023, v němž se autor pozastavuje nad „skokovým desetiprocentním poklesem počtu narozených dětí“ v roce 2022 oproti roku 2021. Dále upozorňuje na „velmi dramatický pokles počtu narozených dětí i v mnoha jiných zemích“ a zamýšlí se nad možnými příčinami poklesu. Ty je podle Fürsta „třeba hledat v něčem, co nastalo na jaře roku 2021“. Možnou souvislost nachází v časové koincidenci s počátkem očkování: „Data z ostatních zemí ukazují přesně stejnou časovou souvislost – pokles porodnosti nastává asi devět měsíců po začátku očkování věkové kohorty rodičů.“ Zcela správně upozorňuje, že může jít o souběh čistě náhodný, zároveň však nenachází evidenci pro jiné vysvětlení pozorovaných změn a vybízí proto k důkladnému zkoumání příčin poklesu počtu narozených dětí.
Reakcí na výše uvedený článek byl text Petra Koubského z 15. února 2023 v Deníku N. Koubský v něm zcela správně poukazuje na zjevnou slabinu Fürstovy úvahy, která spočívá v práci s absolutními počty narozených dětí bez zohlednění četnosti a věkové struktury populace žen v reprodukčním věku (15 až 49 let). Tvrdí, že „pokles počtu nově narozených dětí je nesporný, ale k jeho vysvětlení stačí vzít v úvahu, jak klesá počet žen v reprodukčním věku, další mechanismy nejsou nutné.“ Článek je shrnut tvrzením, že „K vysvětlení [současných] populačních trendů postačují známé a prozkoumané faktory, žádný nový se neprojevuje“. Svá tvrzení dokládá vývojem jednoho z demografy používaných ukazatelů plodnosti („intenzity rození dětí“), tzv. úhrnné plodnosti. K této míře se v následujícím textu budeme opakovaně vracet.
Třetí text publikovaný na stránkách skupiny SMIS je datován 17. února 2023 a má formu rozhovoru imunoložky Zuzany Krátké s demografkou Klárou Hulíkovou Tesárkovou. Ta zopakovala potřebu pracovat nikoliv přímo s počty narozených dětí, ale s mírami plodnosti, které zohledňují počet a věkové složení populace žen v reprodukčním věku: „Pro možné analýzy souvislostí mezi plodností a jakýmikoli možnými faktory je zcela nutné pracovat s intenzitami očištěnými od strukturálních vlivů, protože struktura populace a počty žen jsou zde jednoznačně rušivým faktorem zkreslujícím celé vyznění tématu.“ Pokles počtu narozených dětí Hulíková Tesárková vysvětluje strukturální změnou populace, která je „logickou, plně předvídatelnou a hlavní příčinou pozorovaného poklesu.“ O roku 2021, kdy „došlo nejen v ČR, ale i v jiných evropských státech ke zvýšení úrovně plodnosti, tedy kdy se zvýšil průměrný počet živě narozených dětí na ženu“, hovoří jako o výjimečném. Uvádí: „Pokud by nedošlo k výjimečnosti roku 2021, tak i v něm by nastal další pokles počtu živě narozených dětí o několik tisíc (mezi roky 2019 a 2020 byl pokles o více než 2 tisíce) a pokračující pokles v roce 2022 by nebyl tak mimořádně působící.“
Zatím poslední text zveřejnil na stránkách skupiny SMIS 19. února 2023 Arnošt Komárek. V tomto textu vypočítává meziroční změnu počtu narozených pro jednotlivé měsíce a dochází k závěru, že „S ohledem na bezprecedentní skokovou změnu v porodnosti pozorovanou z měsíce na měsíc (po očištění o sezónnost) se však důvodně domnívám, že ji nelze vysvětlit běžným demografickým vývojem, ale musí jím být nějaká téměř jednorázová událost.“
Proč je rozumné pracovat spíše s mírami plodnosti než přímo s počty narozených
Počet narozených dětí za určité období (např. daný kalendářní rok) na určitém území (např. v daném státě) je zásadním způsobem ovlivňován počtem obyvatel žijícím v daném čase na tomto území. Podstatné je rovněž věkové složení obyvatelstva – počet narozených dětí je těsněji spjat s počtem žen v reprodukčním věku (15-49 let) než s celkovým počtem žen, respektive celkovým počtem obyvatel. Vztáhneme-li počet narozených dětí k počtu žen v reprodukčním věku, získáme míru plodnosti, tedy ukazatel, který vypovídá o „intenzitě“, s jakou ženy v dané (časově a prostorově určené) populaci rodí děti. Demografové používají různé míry plodnosti, přičemž společné těmto mírám je to, že počty narozených dětí vztahují k počtům obyvatel, typicky při zohlednění věku těchto obyvatel.
Jednou z nejčastěji používaných měr plodnosti je tzv. úhrnná plodnost, která představuje „průměrný počet živě narozených dětí, které by se narodily jedné ženě za předpokladu, že by míry plodnosti podle věku zaznamenané ve sledovaném kalendářním roce zůstaly během jejího reprodukčního věku (15–49 let) neměnné.“, viz Metodické poznámky k Demografické příručce ČR 2021. Úhrnná plodnost se vypočítává jako součet věkově specifických měr plodnosti.
Co má vliv na to, s jakou intenzitou se děti rodí? Co ovlivňuje plodnost? Odpovědi hledáme jednak v prostorovém srovnání měr plodnosti (v daném čase), ale také hodnocením vývoje měr plodnosti na daném území v čase.
Vývoj úhrnné plodnosti v ČR od 90. let 20. století
Zaměříme-li se vývoj plodnosti v České republice od 90. let 20. století, budeme pozorovat její nárůst až k rekordně vysokým hodnotám z roku 2021. V tiskové zprávě Českého statistického úřadu (ČSÚ) ze dne 15. prosince 2022 se uvádí: „Díky nárůstu intenzity plodnosti v posledních jedenácti letech se Česko v roce 2021 zařadilo k zemím s nejvyšší úrovní úhrnné plodnosti v Evropě. Průměrný počet dětí na jednu ženu u nás dosáhl 1,83 a byl nejvyšší od roku 1992.“ Tuto zprávu převzala ČTK [1] a nejsledovanější veřejnoprávní [2, 3, 4] i komerční média [5, 6, 7]. Česká televize ke zprávě připojila také graf vývoje úhrnné plodnosti.
Obrázek 1: Vývoj úhrnné plodnosti v ČR v letech 2011-2021. Data: ČSÚ; graf formátován dle [2].
Z tohoto grafu je dobře patrné, že v roce 2021 došlo k značně velkému zvýšení úhrnné plodnosti. Zatímco ve čtyřech předešlých letech byly hodnoty stabilní (1,69 v roce 2017 a 1,71 v následujících třech letech), v roce 2021 byla hodnota úhrnné plodnosti 1,83. Ponechme nyní stranou skutečnost, že díky omezení hodnot vykreslených na svislé ose (1,3 až 1,9 dětí na ženu) se změny hodnoty zdají „dramatičtější“, než kdyby svislá osa začínala hodnotou nula. Meziroční nárůst z hodnoty 1,71 na hodnotu 1,83, tj. nárůst o cca 7 %, je v posledních 10 letech největší a budí tudíž pozornost. Co je jeho příčinou? Než se pustíme do spekulací, podívejme se do dat trochu blíže.
Vývoj počtu narozených dětí a počtu žen v reprodukčním věku („záhada mizejících žen“)
Pro výpočet úhrnné plodnosti jsou zásadní údaje o počtu živě narozených dětí, věku jejich matek a počtu žen v daném věku. Potřebné údaje najdeme v Demografické příručce, kterou každoročně publikuje ČSÚ. V současné době je k dispozici demografická příručka za rok 2021, která byla publikována 31.1.2023.
Jak se dá očekávat, k výpočtu „intenzity rození dětí“ je třeba znát především počty narozených dětí. Celkové počty narozených dětí můžeme zjistit hned z několika různých tabulek demografické příručky, např. z tabulky 1-5a Přirozený pohyb obyvatel na území Česka v letech 1919–2021 (absolutní údaje). Tyto počty jsou znázorněny na obrázku 2.
Počty (živě) narozených se v posledních několika letech (od roku 2015) pohybují v rozmezí 110 až 115 tisíc. Meziroční změny počtu narozených v těchto letech nepřekračují 2 procenta (v letech 2015 až 2021 představují postupně +0,8, +1,7, +1,5, -0,3, -1,6, -1,8 a +1,4 procent).
Obrázek 2: Vývoj počtu živě narozených dětí a počtu žen v reprodukčním věku v letech 1951 – 2021. Data: ČSÚ.
Zabývejme se nyní počtem žen v reprodukčním věku. Vyjdeme z tabulky 1-11 Ženy ve věku 15–49 let podle jednotek věku v letech 1920–2021 (stav k 1. 7.), kde pro jednotlivé roky nalezneme i celkový počet žen ve věku 15-49 let. A zde nás čeká možná velké překvapení. Jak je patrné z obrázku 2 (povšimněme si prudkého poklesu červené čáry na jejím konci), mezi roky 2020 a 2021 došlo k největšímu poklesu počtu žen v reprodukčním věku za období od roku 1951.
Počet žen v reprodukčním věku uváděný v datech ČSÚ poklesl mezi roky 2020 a 2021 o 83 471, tj. přibližně o 3,5 %. Nejde přitom jen o rozdíl mezi nejstarší generací 49letých, která v roce 2020 čítala 73 617 (a která už se v roce 2021 nezapočítává, protože přesáhla vymezené věkové rozmezí), a nejmladší generací 15letých, která v roce 2021 čítala 51 023 žen. Pokud by byly všechny ostatní věkové skupiny stabilní, představovala by tato generační obměna úbytek o velikosti 73 617 – 51 023 = 22 594. Navíc by šlo o úbytek v „okrajových“ věkových skupinách, kde je plodnost dosti nízká, a na úhrnnou porodnost by tato samotná změna měla jen omezený vliv. Jenže došlo (alespoň formálně) k výraznému „úbytku“ žen i ve většině dalších věkových kohort. V každé z věkových skupin 31 až 39 let šlo o pokles o 4 až 5 procent. Opět zdůrazněme, že nejde o efekt různě populačně silných ročníků, neboť zde srovnáváme ve dvou po sobě jdoucích letech velikost „stejné“ skupiny. K výrazné změně došlo například ve skupině 30letých žen (k 1.7.2020) čítající 68 414 žen. Tato skupina tvoří „základ“ skupiny 31letých žen v roce následujícím (k 31.7.2021), v něm však skupina 31letých čítá jen 65 419 žen. Za rok se tedy někde „ztratilo“ jen v této věkové skupině téměř 3 tisíce žen, tj. skoro 4,4 %. Celkově se skupina 15-48 letých z roku 2020 za rok zmenšila o 60 877 žen, tj. o téměř 2,6 %.
Počet žen v reprodukčním věku a sčítání lidu („rozluštění záhady mizejících žen“)
Výše uvedený „úbytek“ počtu žen v reprodukčním věku samozřejmě v žádném případě nepředstavuje skutečný úbytek, ale je dán tím, že v roce 2021 proběhlo v ČR sčítání lidu, domů a bytů. Počty žen v jednotlivých věkových skupinách jsou totiž ve skutečnosti jen odhadovány. V roce 2021 došlo k bezprecedentní úpravě odhadu počtu žen v reprodukčním věku, jak je vidět z obrázků 3 a 4.
Obrázek 3: Meziroční procentuální změna počtu žen v reprodukčním věku. Data: ČSÚ.
Z obrázku 3 je patrné, že k takto dramatickému „poklesu“ počtu žen v reprodukčním věku předtím nikdy nedošlo. Pokles o více než jedno procento byl zaznamenán ještě v letech dvou předešlých sčítání lidu, tj. v letech 2001 (o 1,0 %) a 2011 (o 1,4 %). O to více vynikne mimořádnost poklesu v roce 2021 (o 3,5 %).
Meziroční změna počtu žen v reprodukčním věku pozorovatelná na obrázku 3 je z velké části dána nestejnou velikostí nastupující generace (15letých žen) a skupiny, která reprodukční věk opouští (49 let). Nejde však jen o tuto změnu. Vlivem přeshraniční migrace a (v malé míře) také úmrtnosti se mění také velikost skupiny žen, které žijí v České republice a zůstávají v reprodukčním věku (věk 15-48 let, v následném roce 16-49 let). Velikost těchto změn je zachycen na obrázku 4.
Obrázek 4: Meziroční procentuální změna v četnosti skupiny žen, které zůstávají v reprodukčním věku. Data: ČSÚ.
Z obrázku 4 je patrné, že odhady počtu žen v reprodukčním věku se výrazně upravovaly prakticky při každém sčítání lidu, a to směrem dolů (počet zjištěný sčítáním byl vždy menší, než byl odhad v předchozím roce). Největší korekce byla provedena při sčítání 2021. Po odhlédnutí od přirozeného stárnutí populace šlo o přibližně 2,6% „pokles“ (opravu). Druhá největší oprava (z roku 1970) byla poloviční (o 1,3 %).
Za povšimnutí stojí také skutečnost, že v období mezi posledními sčítáními (2012 až 2020) byly prováděny „aktualizace“ odhadu počtu žen v reprodukčním věku vždy směrem nahoru, a to až o 0,5 % ve všech třech letech před sčítáním. Celkově se tak za roky 2012 až 2020 skupina žen v reprodukčním věku „zvětšila“ o 56 079 osob, aby byla následně v roce 2021 o 60 877 osob „zmenšena“.
Změna plodnosti mezi roky 2020 a 2021
Z výše uvedeného vyplývá, že k velkému zvýšení úhrnné plodnosti mezi roky 2020 a 2021 došlo zřejmě do značné míry díky „aktualizaci“ (snížení) odhadu počtu žen v reprodukčním věku.
Jak moc se tato skoková změna promítla do výpočtu úhrnné plodnosti? Záleží na tom, jak se postavíme k velké diskrepanci v odhadech velikosti věkových kohort před a po sčítání lidu 2021. Skokovou změnu (emigrace či úmrtí desítek tisíc žen v této skupině během jediného roku) nepovažuji za věrohodné vysvětlení.
Věřím, že výsledky ze sčítání lidu jsou nejlepšími dostupnými odhady pro rok 2021 a skutečná hodnota úhrnné plodnosti se tedy v roce 2021 opravdu pohybovala kolem hodnoty 1,83. Při vyloučení výše uvedené reálné skokové změny to ale znamená, že v letech bezprostředně předcházejících byl odhad počtu žen v reprodukčním věku mírně nadhodnocen, a úhrnná plodnost tedy byla mírně podhodnocena.
Vzhledem k výše uvedené „bilanci“ za roky 2012 až 2020, kdy bylo ke skupině žen v reprodukčním věku postupně přičteno přibližně 56 tisíc osob a následně bylo v roce 2021 přibližně 61 tisíc osob odečteno, se jeví jako rozumné počítat s nulovou bilancí za roky 2012 až 2020 a jednotlivé skupiny nechat „stárnout“, ale nijak je nezvětšovat ani nezmenšovat. Skok v roce 2021 v počtu žen v reprodukčním věku pak nebude tak dramatický.
Takto odhadnuté počty žen v reprodukčním věku následně použijeme k standardnímu výpočtu úhrnné plodnosti za roky 2011 až 2021. K výpočtu dále použijeme údaje o počtu živě narozených dětí podle věku matky, které v Demografických ročenkách ČR publikuje pravidelně ČSÚ. Konkrétně vycházíme z tabulek D.04 Živě narození podle věku matky, legitimity a pořadí narození.
Výsledný odhad vývoje úhrnné plodnosti v ČR v letech 2011 až 2021 je zobrazen na obrázku 5 zeleně, červená linka ukazuje původní, oficiálně publikované hodnoty.
Obrázek 5: Vývoj úhrnné plodnosti v ČR v letech 2011 až 2021 dle oficiálních údajů ČSÚ (červeně) a po přepočtu na základě údajů ze sčítání lidu 2021 (zeleně). Data: ČSÚ.
Je třeba si uvědomit, že oficiálně uváděné hodnoty úhrnné plodnosti byly publikovány postupně v jednotlivých letech mezi sčítáními lidu, kdy statistici upravovali odhady počtu žen v reprodukčním věku. Při posledním sčítání lidu (2021) se však ukázalo, že tyto odhady byly značně nadhodnocené a skutečná úhrnná plodnost byla tudíž zřejmě větší, než bylo oficiálně publikováno.
To znamená, že oproti roku 2020, v němž byla oficiálně uváděna úhrnná plodnost 1,71 dítěte na ženu, opravdu došlo v roce 2021 ke zvýšení úhrnné plodnosti, ale změna ve skutečnosti nebyla tak dramatická, jak by se z oficiálních dat (uvádějících hodnotu 1,83 za rok 2021) mohlo zdát.
Jak postupovat při hodnocení dat za rok 2022 – „průměrná“ hodnota může být dosti matoucí pojem
Při hodnocení nejnovějších dat je třeba si uvědomit, že oficiálně publikované hodnoty úhrnné plodnosti před rokem 2021 na jedné straně a její hodnoty z roku 2021 a dalších let na straně druhé jsou (resp. budou) obtížně srovnatelné, protože v roce 2021 došlo k bezprecedentně velké úpravě odhadu počtu žen v reprodukčním věku, přičemž z těchto aktuálních hodnot budou zřejmě vycházet i odhady v letech následujících.
Data publikovaná ČSÚ jsou přebírána do mezinárodní databáze Human Fertility Database. Nejnovější měsíční počty narozených v ní najdeme v sekci Short-Term Fertility Fluctuations v tabulce pro ČR obsahující neupravená zdrojová data Crude data. Počty živě narozených dětí pak slouží k odhadu úhrnné plodnosti. Jelikož počet narozených dětí závisí mj. na ročním období (jsou měsíce, kdy se pravidelně rodí dětí více), jsou k dispozici i odhady měsíčních hodnot úhrnné plodnosti očištěné od tohoto vlivu, které lépe slouží k posuzování dlouhodobých trendů. Pokud by se naplnil odhad, který na základě měsíčních odhadů úhrnné plodnosti publikovaných na stránkách Human Fertility Database uvádí ve svém článku Koubský, tedy pokud by hodnota úhrnné plodnosti klesla v roce 2022 na 1,65, šlo by pokles, který je rozhodně hoden naší pozornosti. Proč si to myslím? Ať už vyjdeme z oficiálně publikovaných (neaktualizovaných) odhadů úhrnné plodnosti, nebo z aktualizovaných hodnot, jde o poměrně zásadní změnu trendu ve vývoji úhrnné plodnosti. Ta od roku 2011 až do roku 2021 rostla. V roce 2022 pak zřejmě došlo k poklesu, a to ne zcela zanedbatelnému. Ilustrujeme tuto změnu na obrázku 6.
Obrázek 6: Úhrnná plodnost v ČR 2011-2021 (přepočítané hodnoty) a predikce pro rok 2022. Hodnoty na základě dat ČSÚ za roky 2011-2021 spočítal O. Vencálek, predikce na rok 2022 je odhadem na základě měsíčních dat Human Fertility Database.
V diskusi o vývoji úhrnné plodnosti po roce 2021 je klíčové hodnocení hodnoty z roku 2021. Pokud bychom připustili, že tato hodnota byla mimořádně vysoká, bylo by možno v následujícím roce (2022) očekávat mírný pokles (princip „návratu k průměru“). Zde jsme však ukázali, že nárůst úhrnné plodnosti v roce 2021 extrémní nebyl. Možná, že expertně můžeme mírný pokles v roce 2022 předjímat (zejména, když už průběžná data tuto „predikci“ potvrzují), těžko by však někdo z dat do roku 2021 predikoval pokles o více než několik setin (modely, které po sérii rostoucích hodnot „predikují“ výraznější pokles vznikají často až poté, co k poklesu opravdu došlo).
Pro úplnost dodejme, že hodnota 1,65 je pouze „predikcí“ a není zcela zřejmé, nakolik je tato predikce ovlivněna skokovou změnou odhadu počtu žen v reprodukčním věku v roce 2021. Dokonce i hodnoty měsíční úhrnné plodnosti za měsíce, pro něž už jsou počty narozených známy (byť jde o předběžné údaje), jsou totiž netriviálními odhady. Jejich autoři při jejich výpočtu nepočítají s detailními daty propojujícími počty narozených dětí s věkem jejich matek, jak je uvedeno na str. 7 v metodologické poznámce k těmto datům.
Za pozornost stojí ještě jedna poněkud nečekaná nesnáz. Jde o interpretaci pojmu „průměrná hodnota“. Ve svém výše citovaném článku Petr Koubský uvádí: „Každopádně lze konstatovat, že rok 2022 není z hlediska úhrnné plodnosti nijak výjimečný, v rámci uplynulého desetiletí se jeví jako průměrný.“ Pokud se na zjištěné hodnoty úhrnné plodnosti pohybující se v rozmezí 1,43 až 1,83 budeme dívat bez jejich časového uspořádání, pak dojdeme k závěru, že hodnota 1,65 je „průměrná“. Zde ale hodnotíme vývoj úhrnné plodnosti v čase. Pokud by šlo o ustálený proces, který osciluje kolem jakési „průměrné hodnoty“, i v tom případě by pojmenování hodnoty 1,65 jako „průměrná hodnota“ bylo v pořádku. Zde však o ustálený proces nejde. Těžko budeme hledat někoho, kdo by trend v letech 2011 až 2021 označil za „konstantní“. Spíše se shodneme na tom, že byl rostoucí. V situaci, kdy pozorujeme rostoucí trend a následný pokles z nejvyšší dosažené hodnoty na hodnotu „průměrnou“, je závěr „Rok 2022 […] byl z hlediska úhrnné plodnosti průměrný.“ zavádějící.
Ze stejného důvodu považuji v Koubského článku za matoucí graf znázorňující krátkodobé fluktuace plodnosti v Česku. I zde je totiž opominuta časová posloupnost zjištěných hodnot. Jen velmi pozorný čtenář si všimne, že „podprůměrné hodnoty“ byly zjištěny převážně v první části sledovaného období, zatímco ty nadprůměrné v části druhé, po níž následuje již zmiňovaný „průměrný rok 2022“.
Koubský na základě toho, že hodnoty úhrnné plodnosti v roce 2022 byly průměrné, vyvozuje tyto závěry: „Nic nenasvědčuje, že by se během něj [roku 2022] projevovaly atypické efekty snižující plodnost. Pokles počtu nově narozených dětí je nesporný, ale k jeho vysvětlení stačí vzít v úvahu, jak klesá počet žen v reprodukčním věku, další mechanismy nejsou nutné.“ Ve své analýze jsem ukázal, že pokud v roce 2022 došlo k poklesu úhrnné plodnosti na hodnotu 1,65 (či podobně nízkou), je třeba nebrat tento pokles na lehkou váhu, ale pátrat po jeho příčinách.
Stojí za to zopakovat, že úhrnná plodnost je charakteristika intenzity rození dětí a je konstruována tak, aby se na ní neprojevovaly strukturní změny v populaci, tedy aby se „odfiltroval“ vliv měnícího se počtu a věkové struktury žen v reprodukčním věku, kteréžto faktory jsou bezesporu hlavními determinantami počtu narozených dětí. Jestliže tedy dochází ke změně trendu ve vývoji tohoto ukazatele (např. mění-li se trend z rostoucího na klesající), je potřeba takovémuto jevu věnovat pozornost. To platí zejména, pokud se potvrdí, že tato změna je výrazná.
Závěr
Svým příspěvkem jsem měl ambici upozornit zejména na následující aspekty hodnocení nejnovějších údajů o počtu narozených dětí a jeho determinantách:
V diskusi o vývoji úhrnné plodnosti po roce 2021 je klíčové hodnocení hodnoty z roku 2021. Ukázali jsme, že v roce 2021 došlo k méně extrémní odchylce ve vývoji úhrnné plodnosti, než by se mohlo z oficiálních dat zdát. Změna, ke které došlo v oficiálně publikovaných odhadech úhrnné plodnosti mezi lety 2020 a 2021, je do značné míry dána bezprecedentně velkou opravou odhadu počtu žen v reprodukčním věku po proběhlém sčítání lidu.
Pokud by hodnota úhrnné plodnosti klesla v roce 2022 na 1,65, jak bylo predikováno z měsíčních dat Human Fertility Database, šlo by o pokles, který je rozhodně hoden naší pozornosti. Jde o výraznou změnu trendu, jejíž příčiny by měly být odhaleny.
K hodnocení vývoje úhrnné plodnosti je třeba přistupovat jako k hodnocení trendu a jeho změn, nikoliv jako k hodnocení ustáleného procesu. Označení hodnoty jako „průměrné“ může být při takovéto analýze silně matoucí.
Před několika dny se sešly dva poměrně důležité výsledky. Ukazují, že my – odborníci sdružení ve SMIS – jsme měli po celou dobu epidemie pravdu ve dvou důležitých otázkách: V problému dlouhodobé imunity získané předchozím proděláním nemoci a v problému protektivity roušek a respirátorů.
Článek v Lancetu přináší dosti definitivní evidenci, že předchozí prodělání infekce koronavirem chrání před opakovanou infekcí příbuznou variantou, zejména před těžkým průběhem covidu, a to velmi dobře a dlouhodobě. Souhrn 65 publikací z 19 různých zemí ukazuje, že ochrana proti těžkému průběhu případné reinfekce dosahuje přibližně 90 procent a přetrvává téměř rok. Po prodělání nemoci byla po 40 týdnech ochrana před reinfekcí deltou či alfou 78,6% (před reinfekcí, nikoliv hospitalizací, což byl v porovnání s očkováním fantastický výsledek). Pro nás, kteří se na stránkách SMIS imunitou dlouhodobě zabýváme, není tento výsledek nijak překvapivý. Je ve shodě s učebnicovou teorií – měl tedy být očekáván – a veškerá data, která byla od začátku epidemie k dispozici, jej potvrzovala. Poprvé jsme o dlouhodobě přetrvávající ochraně psali již v prosinci roku 2020, v létě 2021 jsme dokonce publikovali článek s názvem „Protilátky po prodělané infekci poskytují dlouhodobou ochranu proti onemocnění COVID-19“. Později jsme byli jedni z mála, kteří již v polovině roku 2021 upozorňovali na to, že očkování proti koronaviru v kohortách mladých rekonvalescentů nedává smysl. Za veškeré tyto publikace jsme byli veřejně uráženi, mazáni ze sociálních sítí a označováni za dezinformátory ze strany vlády, různých samozvaných fakt-čekistů ale i českých akademických „odborníků“. Výše citovaná metastudie otištěná v jednom z nejprestižnějších lékařských časopisů na světě nám dává plně za pravdu.
Podobně dopadl spor o protektivitu roušek. Dlouhodobě a opakovaně jsme kritizovali nesmyslná (a často protizákonná) vládní nařízení o nošení roušek ve školách, opakovaně se vysmívali zfalšovaným studiím, které účinnosti roušek „dokazovaly“, psali odborná stanoviska k tomuto tématu, vyvraceli oficiální dezinformace, a dokonce o tom povídali senátorům na půdě Parlamentu České republiky. Koncem ledna 2023 tuto debatu ukončila rozsáhlá metastudie z dílny Cochranovské společnosti. Tato společnost byla – do příchodu globální epidemie ztráty soudnosti – všeobecně považována za klenot v koruně Evidence Based Medicine. Je až přízračné, jak nepřijatelný nátlak byl na autory této studie vyvíjen, aby ji nevydali či alespoň pozdrželi. Studie shrnuje evidenci z jedenácti prospektivních randomizovaných klinických studií, které sledovaly účinnost nošení roušek či respirátorů vzhledem k šíření infekčních respiračních onemocnění. V případě roušek (chirurgických či jiných) dochází k závěru, že jejich nošení má pravděpodobně malý či žádný vliv na šíření infekcí. V případě respirátorů je závěr podobný a výsledek se nezměnil, ani když autoři vybrali jen studie, které se týkaly zdravotnického personálu.
Obě výše uvedené otázky jsou tedy z hlediska standardní Evidence Based Medicine uzavřeny, a to zcela ve shodě se stanoviskem, které sdružení SMIS dlouhodobě zastávalo a pečlivě zdůvodňovalo dostupnými daty.
