Odborné textyRepetitorium - klinická genetikaRepetitorium - statistika

Chí-kvadrát test

JZAAB

c2-test patří mezi neparametrické testy hypotéz o rozložení pracující s daty na nominální nebo ordinální škále, v nejjednodušší podobě ve dvou třídách: pozitivní/negativní, úspěch/neúspěch a podobně. Obecně se jedná o testování shody teoretického rozložení v kontingenčních tabulkách. K výkladu je připojen příklad 006.

 

 

Použití c2-testu pro testování hypotéz v kontingenčních tabulkách.

Kontigenční tabulkou se rozumí tabulka o r řádcích a s sloupcích. Pomocí výpočtu hodnoty χ˛ se testují experimentální (o) a teoretické četnosti (e). Používá se tedy v situacích, kdy se porovnává více populací a více než jeden sledovaný kvalitativní znak. V nejjednodušší modifikaci, pro dvě populace a dvě třídy jednoho kvalitativního znaku jsou pozorované a teoretické četnosti následující:

 

 

 

Soubor I.

Soubor II.

Celkem

Pozitivní

o1,1 

resp. e1,1

o1,2

resp. e1,2

o1,. =

e1,.

Negativní

o2,1

resp. e2,1

o2,2

resp. e2,2

o2,. =

e2,.

Celkem

o.,1 =

e.,1

o.,2 =

e.,2

N

 

Symboly s tečkou označují příslušné řádkové, resp. sloupcové součty. Ty musejí být totožné pro experimentální i teoretické četnosti.

 

Hodnota c2 se počítá analogicky podle již uvedeného vzorce, ale pro všechny buňky tabulky:

 

c2 = (o1,1 – e1,1)2 / e1,1 + (o1,2 – e1,2)2 / e1,2 + (o2,1 – e2,1)2 / e2,1 + (o2,2 – e2,2)2 / e2,2

 

Pro výpočet teoretických četností e platí obecné vzorce, ve kterých se využívá příslušných řádkových a sloupcových součtů a celkového počtu pozorování:

e1,1 = o.,1 * o1,. / N

e1,2 = o.,2 * o1,. / N

e2,1 = o.,1 * o2,. / N

e2,2 = o.,2 * o2,. / N,

přitom o.,1=x11+x21 , o1,.=x11+x12  a podobně.

 

Počet stupňů volnosti se počítá jako součin (počet proměnných ve sloupcích – 1) * (počet proměnných v řádcích – 1), tedy v tomto případě (2 – 1)*(2 – 1) = 1.

 

Zcela analogicky se postupuje u kontingenční tabulky s r řádky a s sloupci, počet stupňů volnosti je (r-1)*(s-1).

 

 

 

Použití c2-testu pro posouzení rozložení experimentálních a teoretických četností

Podobný výpočet hodnoty c2 se používá při posuzování rozložení experimentálních a teoretických četností. Zjednodušeně lze říci, že se pracuje s kontingenční tabulkou o jednom řádku a k sloupcích. Postup se používá pro testování shody experimentálního a libovolného teoretického rozložení. Pro k-1 stupňů volnosti je výpočet následující:

 

c2 = (o1 – e1)2 / e1 + (o2 – e2)2 / e2 + .... + (ok – ek)2 / ek

 

Základní podmínkou je, že součet všech experimentálních četností a četností teoretických musí být totožný.

 

 

Použití c2-test pro testování hypotéz ve čtyřpolní tabulce

Čtyřpolní tabulka je speciálním případem kontingenční tabulky se dvěma řádky a dvěma sloupci. Pro výpočet hodnoty se používá zjednodušený vzorec, počet stupňů volnosti je 1. Musí být splněna podmínka, že četnost v každé buňce čtyřpolní tabulky je nejméně 5. Pro nižší počty v jednotlivých buňkách je nutné použít Yatesovu korekci nebo Fisherův přímý exaktní test.

 

Postup

Zjistíme počty pozitivních a negativních osob v souboru I. a II. a četnosti spořádáme do čtyřpolní tabulky.

 

 

 

 

Soubor I.

Soubor II.

Celkem

Pozitivní

a

b

a+b

Negativní

c

d

c+d

Celkem

a+c

b+d

N

 

Formulujeme nulovou hypotézu (soubory pocházejí ze stejné výchozí populace, frekvence pozitivních se neliší).

 

Vypočteme hodnotu c2 (chí-kvadrát):

 

c2 = [(ad – bc)2 * N] / [(a+c)*(b+d)*(a+b)*(c+d)]

 

 

Zjištěnou hodnotu c2 porovnáme s kritickou hodnotou (pro čtyřpolní tabulku vždy pro 1 stupeň volnosti):

 

 

p

0,05

0,025

0,01

0,005

0,001

χ2

3,84

5,02

6,64

7,88

10,82

 

 

Pokud je vypočítaná hodnota c2 vyšší než některá z kritických hodnot uvedených v tabulce, nulovou hypotézu zamítáme (pravděpodobnost omylu odečteme z tabulky nebo aproximujeme).

 

 

 

Yatesova korekce

Pro situace, kdy není splněna podmínka dostatečných četností v buňkách čtyřpolní tabulky, lze použít modifikovaný výpočet testovací charakteristiky, nazývaný Yatesova korekce. Je rovněž možné použít  Fisherův přímý exaktní test.

 

 

Literatura

Kubánková, V., Hendl, J.: Statistika pro zdravotníky. Zdravotnické aktuality 208, Praha, Avicenum 1986.

Geigy Scientific Tables. Volum 2, Ciba-Geigy, Basle, 1982.

 

 

Další informace

Fisherův přímý exaktní test

Yatesova korekce

Příklad 006: chí-kvadrát test

 

Antonín Jabor, recenzoval Zdenek Kubíček