Včera jsem narazil na zajímavou studii, ve které je jedna tabulka, která je na první pohled poněkud matoucí, ovšem po překopírování, opravě jednoho špatně opsaného údaje a přepočtu jsem se dostal k velmi zajímavému zjištění ohledně diabetes, které následně potvrzuje další studie.
https://onlinelibrary.wiley.com/doi/full/10.1002/ijc.32105
Kdo prováděl výzkum:
"The Copenhagen General Practitioners Laboratory serviced general practitioners from the Greater Copenhagen area with approximately 1.1 million inhabitants (first quarter of 2010). A database (CopD) included measurements of vitamin D from this laboratory in the period from 2004 to 2010. A reverse J‐shaped association between serum levels of vitamin D and all‐cause mortality has been reported based on the CopD database indicating an upper beneficial limit of vitamin D8."
Obyvatelům tedy měřili vitamín D (kalcifediol). Nejsem sice statistik, ale statistiky obecně a analýza dat mě baví a tak jsem si udělal několik tabulek. Původní tabulka s nemocma, tabulka 1 totiž obsahuje pouze počty pacientů a jakési procenta.
Co jsem udělal. Nejprve jsem jsem použil čísla, kolik bylo vzorků (N) v každém sloupci reprezentující hladinu vitamínu D (suma vzorků u dané nemoci). Pak jsem pro každý sloupec spočítal podíl kohorty v daném sloupci, tj. podíl vzhledem k sumě vzorků pro danou nemoc. Tyto údaje pro mě jsou dále podstatné pro přepočet:
nmol/L | <12.5 | 12.5–25 | 25–50 | 50–75 | 75–100 | 100–125 | ≥125 |
μg/L | 5 | 5-10 | 10-20 | 20-30 | 30-40 | 40-50 | 50> |
podíl | 17,71 | 5,91 | 3,03 | 3,74 | 8,32 | 26,57 | 53,15 |
Ještě než vysvětlím k čemu tyto hodnoty potřebuji. Musel jsem si nejprve položit otázku... Uvedu příklad:
Např. v hladině D=25–50 nmol/L je srdeční choroba cca 1/3 populace v dané kategorii (nemoci). Tato hladina se v nemocné populaci vyskytuje nejčastěji a tvoří jednu třetinu všech hladin.
Oproti tomu vitamín D v hladině ≥125 má jen 1/53 populace (nemocných!).
Co by se stalo, kdyby všichni měli stejné množství vitamínu D? Byl by 3x větší nárůst srdeční choroby?. Vezmu tedy tu hodnotu v procentech (tj. 33) a vynásobím ji podílem a dostanu: cca 100,46
Co by se stalo, kdyby všichni začali brát vysoké dávky vitamínu D? Byl by 53x větší nárůst srdeční choroby(?) 1,85* podíl 53,15 = 98,59
Když tedy porovnáme tyto dvě hladiny, tak nám vychází, že u vyšších hladin je menší pravděpodobnost.
Jenomže je otázka, jak velká je vlastně toleranční odchylka? Čím méně vzorků, tím větší nepřesnost. Musel jsem tedy počítat i s tou možností, že u malého množství vzorků odhadem pod 300 je větší tolerance chyby. Jak jsem již napsal, tak nejsem statistik, takže jsem si nakonec nastavil jako limit číslo 1000. To je dostatečný počet vzorků. Po výpočtu pokud je hodnota menší než 1, je počet vzorků dostatečný a výsledek je přesný. Čím větší tím horší...
VÝPOČET TOLERANČNÍ ODCHYLKY (1000 považuji za dostatečně široký vzorek) | ||||||||
<12.5 | 12.5–25 | 25–50 | 50–75 | 75–100 | 100–125 | ≥125 | ||
Charlsonův index komorbidity nula | 0,104877 | 0,034576 | 0,017538 | 0,02145 | 0,04774 | 0,15569 | 0,3058104 | |
Cukrovka | 0,722543 | 0,263922 | 0,168039 | 0,26688 | 0,692521 | 2,364066 | 4,85436893 | |
Srdeční choroba | 0,783085 | 0,276091 | 0,144051 | 0,184298 | 0,411353 | 1,201923 | 2,57731959 | |
Plicní nemoc | 0,631313 | 0,226552 | 0,126231 | 0,165153 | 0,361795 | 1,055966 | 2,0242915 | |
Neurologická onemocnění a demence | 0,754148 | 0,260349 | 0,130685 | 0,160823 | 0,346021 | 0,97371 | 1,97628458 | |
Nemoc jater | 2,03252 | 0,869565 | 0,600962 | 0,942507 | 2,164502 | 6,369427 | 10,6382979 | |
Nemoc ledvin | 4,608295 | 1,592357 | 0,962464 | 1,187648 | 2,777778 | 8,064516 | 15,8730159 | |
Vřed | 1,689189 | 0,697837 | 0,404204 | 0,586854 | 1,373626 | 3,802281 | 6,99300699 | |
Reumatologická nemoc | 5,524862 | 1,453488 | 0,635324 | 0,637349 | 1,149425 | 3,021148 | 5,23560209 | |
HIV | 200 | 40 | 26,31579 | 33,33333 | 66,66667 | 200 | #HODNOTA! |
Z toho mně vychází že pokud mě zajímá nejvyšší hladina, tak bohužel nemám dostatek dat, abych bylo možné s jistotou říct, že na ty čísla se dá spolehnout.
Plicní onemocnění a nemoc neurologická vyšly dost negativně. Čím vyšší hladina vitamínu D, tím horší. Nicméně příčina těchto výsledků není jasná. Zajimavější pro mě je diabetes. Tam jsou čísla spolehlivá až do hladiny 30-40 g/L, tj. v normě, ale pod optimální hladinou.
Nyní k samotným výsledkům. Jak jsem napsal výše, co by se stalo, pokud by všichni v dané kategorii/kohortě měli stejnou hladinu?
PŘEPOČET ODHAD* [%] | <12.5 | 12.5–25 | 25–50 | 50–75 | 75–100 | 100–125 | ≥125 |
Charlsonův index komorbidity nula | 97,75 | 98,92 | 99,93 | 101 | 100,9 | 98,81 | 100,6 |
Cukrovka | 144,6 | 132,1 | 106,3 | 82,73 | 70,91 | 66,34 | 64,62 |
Srdeční choroba | 108,1 | 102,3 | 100,5 | 97,04 | 96,7 | 105,7 | 98,59 |
Plicní nemoc | 116 | 107,9 | 99,18 | 93,68 | 95,11 | 104,1 | 108,6 |
Neurologická onemoc. a demence | 100,1 | 96,73 | 98,74 | 99,15 | 102,5 | 116,3 | 114,6 |
Nemoc jater | 171,5 | 133,7 | 99,16 | 78,13 | 75,67 | 82,13 | 98,36 |
Nemoc ledvin | 117,4 | 113,4 | 96,1 | 96,24 | 91,52 | 100,7 | 102,3 |
Vřed | 142,9 | 115,4 | 102,1 | 86,88 | 82,56 | 95,26 | 103,6 |
Reumatologická nemoc | 59,31 | 75,22 | 88,17 | 108,6 | 134 | 162,8 | 187,9 |
HIV | 75,03 | 125,2 | 97,49 | 95,11 | 105,8 | 112,6 | - |
* Odhad při změně poměru lidí v dané hladině na 1. |
PŘEPOČET ODHAD* [%] | <12.5 | 12.5–25 | 25–50 | 50–75 | 75–100 | 100–125 | ≥125 |
Charlsonův index komorbidity nula | 97,75 | 98,92 | 99,93 | 101 | 100,9 | 98,81 | 100,6 |
Cukrovka | 144,6 | 132,1 | 106,3 | 82,73 | 70,91 | 66,34 | 64,62 |
Srdeční choroba | 108,1 | 102,3 | 100,5 | 97,04 | 96,7 | 105,7 | 98,59 |
Plicní nemoc | 116 | 107,9 | 99,18 | 93,68 | 95,11 | 104,1 | 108,6 |
Neurologická onemoc. a demence | 100,1 | 96,73 | 98,74 | 99,15 | 102,5 | 116,3 | 114,6 |
Nemoc jater | 171,5 | 133,7 | 99,16 | 78,13 | 75,67 | 82,13 | 98,36 |
Nemoc ledvin | 117,4 | 113,4 | 96,1 | 96,24 | 91,52 | 100,7 | 102,3 |
Vřed | 142,9 | 115,4 | 102,1 | 86,88 | 82,56 | 95,26 | 103,6 |
Reumatologická nemoc | 59,31 | 75,22 | 88,17 | 108,6 | 134 | 162,8 | 187,9 |
HIV | 75,03 | 125,2 | 97,49 | 95,11 | 105,8 | 112,6 | - |
* Odhad při změně poměru lidí v dané hladině na 1. |
Čím nižší hodnota tím lepší.
Z toho jsem odvodil tabulku s rozdílem oproti 100:
μg/L | 5 | 5-10 | 10-20 | 20-30 | 30-40 | 40-50 | 50> |
ROZDÍL [%] | <12.5 | 12.5–25 | 25–50 | 50–75 | 75–100 | 100–125 | ≥125 |
Charlsonův index komorbidity nula | 2,251 | 1,076 | 0,071 | -0,962 | -0,902 | 1,186 | -0,626 |
Cukrovka | -44,64 | -32,12 | -6,323 | 17,27 | 29,09 | 33,66 | 35,38 |
Srdeční choroba | -8,105 | -2,302 | -0,466 | 2,964 | 3,301 | -5,697 | 1,405 |
Plicní nemoc | -16 | -7,852 | 0,818 | 6,323 | 4,887 | -4,077 | -8,596 |
Neurologická onemoc. a demence | -0,09 | 3,268 | 1,258 | 0,849 | -2,501 | -16,33 | -14,65 |
Nemoc jater | -71,5 | -33,75 | 0,839 | 21,87 | 24,33 | 17,87 | 1,643 |
Nemoc ledvin | -17,41 | -13,36 | 3,9 | 3,764 | 8,48 | -0,678 | -2,315 |
Vřed | -42,88 | -15,39 | -2,079 | 13,12 | 17,44 | 4,743 | -3,6 |
Reumatologická nemoc | 40,69 | 24,78 | 11,83 | -8,613 | -33,96 | -62,77 | -87,87 |
HIV | 24,97 | -25,17 | 2,511 | 4,893 | -5,771 | -12,6 | - |
Čím vyšší hodnota tím lepší.
Když se tedy podíváte na tu cukrovku, tak vydíte, že čím více vitamínu D, tím menší pravděpodobnost, že onemocníte. Což by mělo jít ruku v ruce s následující studií:
Pro přidání komentáře se musíš přihlásit nebo registrovat na signály.cz.