#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Analýza dat v neurologii
L. Vybrané komentáře k odhadům a interpretaci poměru šancí


Autoři: L. Dušek;  T. Pavlík;  J. Jarkovský ;  J. Koptíková
Působiště autorů: Masarykova univerzita, Brno ;  Institut bio­statistiky a analýz
Vyšlo v časopise: Cesk Slov Neurol N 2015; 78/111(2): 240-243
Kategorie: Okénko statistika

Tento jubilejní díl seriálu budeme věnovat dílčím komentářům zaměřeným na vybrané problémy spojené s odhady šance (O) a poměru šancí (OR) v asociačních studiích. Uzavíráme tak kapitoly, na které nebyl dostatečný prostor v předchozích dílech seriálu, a také odpovídáme na nejčastější otázky, které v této souvislosti zaznívají v literatuře či na přednáškách.

Pro srozumitelnost v úvodu opakujeme nejjednodušší tabulku četností 2 × 2 s typickým značením četností v jednotlivých polích a, b, c, d (viz též díl seriálu 49) (tab. 1).

Tab. 1. Výpočet relativního rizika (RR) a poměru šancí (OR) na základě tabulky četností 2 × 2.
Výpočet relativního rizika (RR) a poměru šancí (OR) na základě tabulky četností 2 × 2.

Jak postupovat v případě, kdy se v tabulce četností objeví políčko s nulou?

Pro výpočet poměru šancí je tato situace vážnou komplikací, což je patrné již z výše uvedených vztahů pod tabulkou četností. Pokud nulová četnost nastane v buňkách b, c, pak poměr šancí nelze spočítat, neboť nemůžeme dělit nulou. Jsou‑li nulové buňky a, d, pak OR sice spočítat můžeme, ale nevyhnutelně vyjde nulové, což nedává interpretační smysl. Do publikací v obou případech píšeme poznámku „OR not defined due to zero cells in frequency table“. Nastane‑li taková situace, je třeba prověřit, zda skutečně jde o kategorie vzácné až nemožné, kde je nulová četnost objektivní fakt, anebo jde o chybu vzorkování. S takovou chybou se setkáváme především u malých vzorků, které snadno podhodnotí výskyt ně­kte­ré vzácnější kategorie.

Jak se u asociačních prospektivních a retrospektivních studií liší cílový parametr, endpoint?

V zásadě se neliší, oba experimentální plány sledují stejný cíl. Jednoduše řečeno, jde o to posoudit vliv expozice (tedy experimentálního zásahu, léčby, vliv nějakého znaku, prediktoru, rizikového faktoru) na výskyt zkoumaného jevu (typicky nemoci). Výskyt nemoci je tedy cílovým parametrem, přičemž ho měříme jako relativní četnost ve skupině exponovaných a neexponovaných jedinců. Rozdíl mezi prospektivním a retrospektivním plánem je pouze ten, že se na daný problém dívají z opačné strany. Pouze prospektivní plán ovšem umožňuje reprezentativní odhad prevalence zkoumaného jevu, proto je možné u něj odhadovat relativní riziko (RR) jeho výskytu. Retrospektivní plány již při zařazování jedinců do studie rozhodují o poměru případů a kontrol, relativní výskyt případů tedy není nezávislý odhad, ale výsledek náběru probandů.

Existuje nějaké jednoduché pravidlo, kdy použít OR a kdy RR?

Ano, existuje. Jak bylo vysvětleno v dílech 48 a  49 našeho seriálu, nejlepší doporučení jsou tato:

  • RR nemá smysl a může být silně zavádějící u retrospektivních studií (studie případů a kontrol, retrospektivní kohortové studie); zde je naopak dobře interpretovatelný odhad OR.
  • U prospektivních observačních studií (kohortové, průřezové) lze sice odhadnout ORRR, z hlediska interpretace je ale lepší RR, neboť pracuje s reprezentativním odhadem rizika (pravděpodobnosti) výskytu jevu v exponované a neexponované skupině; to samé platí u prospektivních randomizovaných studií.
  • Rovněž následné metaanalýzy, je‑li to možné, srovnávají preferenčně odhady RR. U studií, kde to jde, je možné hodnotu RR aproximovat (přepočítat) z odhadu OR (např. při splnění předpokladu vzácných onemocnění).

Proč je v literatuře stále častěji používán OR, a to i u prospektivních kohortových studií?

Důvodů je více. V první řadě můžeme uvést snadný způsob výpočtu a jeho dobrou dostupnost v běžných softwarových nástrojích. Více sofistikovaný důvod je, že odhad OR je běžný výstup složitějších modelovacích technik, které jsou velmi efektivní pro hodnocení komplikovaných tabulek četností. Zejména jde o logistickou regresi, které se budeme věnovat v ně­kte­rém z dalších dílů našeho seriálu. Tato modelovací technika umožňuje odhady i vícerozměrně adjustovaného poměru šancí, a tedy odfiltrování současného vlivu více zavádějících faktorů. Právě proto bývá využívána i pro kohortové prospektivní studie.

Odhad OR má také řadu dalších numerických výhod, které ho často favorizují před RR. Například platí následující vztah:

  • ORriziko  = 1/ORprotekce nebo jiným označením ORneúspěch = 1/ ORúspěch

Jinými slovy jednoduchou inverzí můžeme přepínat odhad OR pro hodnocení neúspěchu nebo úspěchu, a to nejen vlastní bodový odhad OR, ale i hraniční hodnoty jeho intervalu spolehlivosti. Výhodou je, že můžeme velmi lehce měnit interpretaci stejného odhadu z pohledu rizikovosti nebo naopak protekce dané sledovaným faktorem. Při publikaci výsledků, kde zpracováváme více odhadů OR, pak můžeme snadno nastavit jejich „směr“ tak, aby vypovídaly konzistentně. Odhad RR takovou vlastnost nemá, a tudíž musíme již při primárním výpočtu velmi pozorně sledovat, zda hodnotíme „úspěch“ či „neúspěch“, nebo „rizikový vliv faktoru vs. protektivní vliv faktoru“. Ukázku takových recipročních výpočtů dokumentujeme na příkladu 1.

Příklad 1. Ukázka recipročního přepočtu poměru šancí a relativního rizika v asociační studii.
Příklad 1. Ukázka recipročního přepočtu poměru šancí a relativního rizika v asociační studii.

V příkladu 2 dále popisujeme situaci, kdy hodnotíme asociaci více potenciálních prediktorů vůči cílovému parametru „nemoc: ano/ ne“. Hodnotíme tedy více tabulek četností a výsledkem je více odhadů OR. Pro čtenáře je vhodné a přehledné, pokud zkoumané prediktivní faktory mají stejný „směr“, tudíž vyjadřují buď rizikový, nebo naopak protektivní účinek. Toto sladění provedeme u odhadů OR velmi snadno pomocí výše uvedeného vztahu mezi OR pro úspěch a OR pro neúspěch. Příklad 2 také dokumentuje, že stejně lze postupovat i u hranic intervalu spolehlivosti.

Příklad 2. Přepočet hodnot poměru šancí (OR) z „šance protekce“ na „šanci rizika“.
Příklad 2. Přepočet hodnot poměru šancí (OR) z „šance protekce“ na „šanci rizika“.

Je možné jednoduše srovnat hodnoty RROR? A za jakých situací?

Zjednodušeně shrnuto má takové přímé srovnání smysl v podstatě jen ve dvou modelových situacích:

  • U retrospektivních studií (typicky studie případů a kontrol) při platnosti předpokladu vzácného onemocnění (viz díl 48 a  49 seriálu), kdy lze hodnotu OR aproximovat na RR.
  • U prospektivních studií (typicky observační kohortové nebo průřezové studie, případně i randomizované klinické studie), kde můžeme současně odhadovat RROR; v takovém případě:
    • hodnoty ORRR se absolutně shodují pouze v případě, kdy žádná asociace mezi zkoumaným jevem a expozicí neexistuje, potom platí OR = RR = 1,
    • hodnoty OR jsou vždy extrémnější (tj. vzdálenější od hodnoty 1, a to oběma směry) než odpovídající hodnoty RR; tento rozdíl narůstá s rostoucí prevalencí zkoumaného znaku ve zdrojové tabulce četností –  u běžně se vyskytujících znaků (prevalence větší než 5– 10 %) nelze odhady RROR vzájemně zaměňovat.

V jakém vztahu je odhad rizika a odhad šance?

V literatuře je většinou pojednáváno o srovnání a vzájemné aproximaci ORRR (viz odstavec výše). Srovnání šance výskytu jevu A (OA) a rizika jevu A (PA) je ale stejně důležité. Připomeňme, že rizikem zde označujeme pravděpodobnost nastání nějakého jevu či události. Šance nastání jevu a jeho riziko ve smyslu pravděpodobnosti jsou vzájemně spojeny již definičním vztahem:

  • OA PA/ (1 –  PA), a tedy recipročně PA = OA/ (1 + OA).

Z těchto vztahů opět vidíme, že u vzácných jevů, kde je PA nízké, bude hodnota OA blízká až rovna hodnotě PA. U běžných jevů budou hodnoty číselně nesouměřitelné, neboť platí:

  • PA je pravděpodobnost, a tudíž může nabývat hodnot pouze v intervalu 0 až 1;
  • naopak OA nabývá hodnot od 0 až do nekonečna.

Výše uvedené omezení možných hodnot rizika (PA) se následně promítá i do omezení rozsahu hodnot poměru rizik, tedy relativního rizika RR. Relativní riziko je totiž poměrem rizik, tedy pravděpodobností, a ty jsou limitovány rozsahem možných hodnot 0 až 1. Tudíž i hodnoty RR jsou omezeny, a proto také hodnoty OR vycházejí „extrémněji“, tedy jsou číselně vzdálenější od referenční hodnoty 1 než hodnoty RR. Uvažujme hodnotu rizika v kontrolní skupině stejnou jako v příkladu 1, tedy 0,4. Potom je maximální možná vypočtená hodnota RR = 2,5, protože riziko v čitateli poměru RR nemůže překročit 1, a tedy 1/ 0,4 = 2,5. Hodnota rizika v kontrolní skupině tak determinuje maximální možnou hodnotu RR; s klesajícím rizikem v kontrolní skupině roste maximální možná hodnota RR. Hodnoty poměru šancí takové omezení nemají. Srovnání různých hodnot RROR uvádíme v tab. 2, ze které je patrné, že:

  • s rostoucí prevalencí sledovaného rizikového jevu narůstá rozdíl mezi OR RR; OR je přitom vždy vzdálenější od hodnoty 1;
  • i v případě, kdy je hodnota RR téměř rovna 1, může OR nabývat hodnoty > 1;
  • hodnota OR na rozdíl od RR souvisí s absolutním rozdílem rizika mezi kontrolní a experimentální skupinou, s rostoucí diferencí rizik roste i OR;
  • RR nesouvisí s absolutním rozdílem rizik; tedy i při stejné absolutní diferenci rizika mezi kontrolní a experimentální skupinou může RR nabývat různých hodnot.

Tab. 2. Srovnání různých odhadů relativního rizika (RR) a poměru šancí (OR) z asociačních tabulek četností 2 x 2.
Srovnání různých odhadů relativního rizika (RR) a poměru šancí (OR) z asociačních tabulek četností 2 x 2.
Modelové příklady odvozené z prospektivních kohortových studií, u kterých lze smysluplně hodnotit RR i OR.

Velký číselný rozsah možných hodnot OR (tj. od 0 až do nekonečna) představuje problém při vzájemném srovnávání výstupů různých studií. Pro tento účel existují standardizované postupy, díky kterým je interpretace hodnot OR snazší. Tomuto problému se budeme věnovat v příštím díle seriálu.

doc. RNDr. Ladislav Dušek, Ph.D.

Institut bio­statistiky a analýz

MU, Brno

e‑mail: dusek@iba.muni.cz


Štítky
Dětská neurologie Neurochirurgie Neurologie

Článek vyšel v časopise

Česká a slovenská neurologie a neurochirurgie

Číslo 2

2015 Číslo 2

Nejčtenější v tomto čísle
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#