XI. Miery závislosti dvoch znakov

11.1 Druhy dát

Doteraz sme uvažovali len kvantitatívne štatistické znaky, t.j. tie, ktoré majú číselne vyjadriteľné hodnoty. V praxi sa však často stretávame aj s inými druhmi dát. Z hľadiska ich povahy rozlišujeme

  1. kvalitatívne (kategoriálne) dáta: popisujú vlastnosti štatistickej jednotky slovne, definíciou alebo symbolmi. Nadobúdajú vždy diskrétne hodnoty a nemusia sa dať jednoznačne merať. Také sú
    1. nominálne dáta: jednotlivé kategórie vieme pomenovať, ale nie usporiadať. Nedajú sa robiť žiadne matematické operácie okrem určovania početností (napr. pohlavie, farba očí).
    2. ordinálne (poradové) dáta: jednotlivé kategórie vieme usporiadať. Ani tu nemožno robiť matematické operácie (napr. kategórie odpovedí: veľmi nesúhlasím - skôr nesúhlasím - skôr súhlasím - veľmi súhlasím).
  2. kvantitatívne (kardinálne) dáta: vyjadrujú vlastnosti štatistickej jednotky priamo číslami. Môžu byť diskrétneho i spojitého charakteru.
    1. intervalové dáta: jednotlivé hodnoty sú dobre usporiadané, ale dáta nemajú absolútnu nulu. Hodnoty je možné sčítať a odčítať, ale nie násobiť a deliť (napr. časové údaje).
    2. racionálne dáta: dobre usporiadané, spravidla spojité dáta s absolútnou nulou. Je možné robiť všetky matematické operácie (napr. hmotnosť, veľkosť obratu).
Z hľadiska množiny možných hodnôt rozlišujeme diskrétne (t.j. oddeliteľné) a spojité dáta. Oproti nášmu doterajšiemu chápaniu je kategória diskrétnych dát bohatšia o kvalitatívne dáta. U kvalitatívnych štatistických znakov sa podľa počtu možných kategórií rozlišujú
  1. dichotomické (alternatívne, binárne, dvojné) dáta, ktoré majú iba dve možné kategórie (napr. pohlavie) a
  2. multinomické dáta, ktoré majú viac možných kategórií (napr. farba očí).

11.2 Závislosť dvoch kvalitatívnych znakov

Základom pre skúmanie závislosti či nezávislosti dvoch kvalitatívnych znakov je podobne ako pri jednorozmerných štatistických súboroch zistenie početností. Skúmané znaky označíme symbolmi A a B a budeme ich nazývať riadkový a stĺpcový faktor. Nech znak A môže nadobúdať a znak B rôznych hodnôt. Predpokladáme, že máme k dispozícii štatistických jednotiek a na každej z nich bola určená hodnota oboch znakov. Jednotlivé štatistické jednotky zatriedime do kategórií zodpovedajúcich všetkým možným kombináciám hodnôt oboch faktorov. Zistené početnosti zhŕňa kontingenčná tabuľka:

Ak , hovoríme tiež o asociačnej tabuľke.

Z počtu pravdepodobnosti vyplýva, že v prípade nezávislosti faktorov A a B (t.j. príslušnosť do kategórie nemá žiadny vplyv na pravdepodobnosť príslušnosti do a naopak) by v ideálnom prípade mali pozorované početnosti spĺňať vzťah (V prípade nezávislosti totiž pre príslušné pravdepodobnosti platí a najlepšie odhady týchto pravdepodobností sú , a .) Tieto početnosti označíme a budeme ich nazývať očakávané početnosti.

Definícia 11.2.1: Pre kontingenčnú tabuľku s početnosťami a očakávanými početnosťami definujeme veličinu chí-kvadrát vzťahom

Veličina je základom pre test hypotézy o nezávislosti faktorov A a B. Jej malé hodnoty hovoria v prospech hypotézy, veľké hodnoty v neprospech hypotézy.

Tvrdenie 11.2.2: Ak faktory A a B kontingenčnej tabuľky sú stochasticky nezávislé, potom rozdelenie pravdepodobnosti veličiny sa pre blíži k rozdeleniu stupňami voľnosti.

V praxi sa zhoda s asymptotickým rozdelením považuje za dostatočnú, ak . Ak je , hypotézu o nezávislosti faktorov A, B zamietame (na hladine ).

Podobnou veličinou je . Má asymptoticky rovnaké rozdelenie ako . Niekedy sa nazýva -pomer vierohodností. V praxi sa veličiny a väčšinou líšia iba málo.

Aplikácia -testu naráža v praxi na viacero problémov. Okrem toho skutočná hladina testu sa od nominálnej (asymptotickej) môže dosť líšiť. Preto je na testovanie závislosti lepšie použiť USP-test, ktorý je založený na štatistike Tento test je pre vždy definovaný a presne dodržiava nominálnu hladinu . Jeho jedinou nevýhodou je, že je založený na náhodných permutáciách; preto sa odporúča zvoliť ich počet čo najvyšší.

Hodnota veličiny nám hovorí, či môžeme závislosť medzi faktormi A a B považovať za štatisticky významnú alebo nie. Nehovorí však o miere závislosti týchto znakov. Preto boli definované viaceré miery stupňa závislosti:

  1. koeficient Platí, že ; pritom hodnoty blízke 0 hovoria o vysokej miere nezávislosti. Nevýhodou je neexistencia horného ohraničenia.
  2. kontingenčný koeficient (Pearsonovo C) Platí . Hodnoty blízke 0 značia malú závislosť, hodnoty blízke 1 veľkú závislosť znakov. Jeho maximálna hodnota však závisí od rozmerov tabuľky.
  3. upravený kontingenčný koeficient (Sakodov koeficient) Je . Tento koeficient nezávisí od rozmerov tabuľky. Jeho interpretácia je podobná ako u CC.
  4. Cramerovo V Je . Ani tento koeficient nezávisí od rozmerov tabuľky. Pre asociačnú tabuľku platí . Interpretácia je rovnaká ako u .

Koeficient asociácie založený na USP štatistike je

  1. koeficient D Platí .

Ďalšie miery závislosti sú odvodené od iných princípov. Goodman a Kruskal definovali niekoľko koeficientov vyjadrujúcich, do akej miery nám znalosť hodnoty jedného znaku pomáha predpovedať hodnotu druhého znaku:

  1. koeficienty Je . Vyjadruje to relatívne zmenšenie chyby predpovede druhého znaku.
  2. koeficienty neurčitosti kde , a sú entropie (miery neurčitosti) príslušných znakov. Je . Koeficienty predstavujú zníženie neurčitosti určenia jedného znaku pri znalosti hodnoty druhého.

Závislosť ordinálnych znakov:

V tomto oddieli budeme predpokladať, že hodnoty znakov A a B sú v tabuľke usporiadané (rovnakým spôsobom). Nech označuje pozorovanie patriace do kategórie , t.j. pozorovali sme na ňom hodnoty a . Spoločným základom pre výpočet rôznych mier závislosti je zistenie nasledujúcich údajov:

     = počet dvojíc pozorovaní , pre ktoré platí  a alebo  a (konkordantné páry)
     = počet dvojíc pozorovaní , pre ktoré platí  a alebo  a (diskordantné páry)
     = počet dvojíc pozorovaní , pre ktoré platí
     = počet dvojíc pozorovaní , pre ktoré platí .
Používané miery závislosti sú:
  1. koeficient Platí . Je to vlastne rozdiel pravdepodobností konkordantného a diskordantného usporiadania náhodne vybranej dvojice pozorovaní za podmienky, že neexistujú zhodné poradia. Ak sú faktory A, B nezávislé, je blízke 0, ak sú silne závislé, sa blíži +1 alebo -1. Nevýhodou je, že zo vzťahu nevyplýva nezávislosť faktorov. Túto mieru navrhli Goodman a Kruskal.
  2. Kendallovo Aj tu je . Je to vlastne modifikácia Kendallovho koeficientu poradovej korelácie pre prípad zhodných poradí. Interpretácia je podobná ako u .
  3. Kendallovo Niekedy sa nazýva tiež Stuartovo . Tiež je . Berie do úvahy celkový počet pozorovaní a rozmery tabuľky. Interpretácia je podobná ako u .
  4. Somersovo d Sú to asymetrické koeficienty predpokladajúce jednosmernú závislosť (v prvom prípade je B závislý a A nezávislý faktor, v druhom naopak). Je . Je to vlastne rozdiel pravdepodobností konkordantného a diskordantného poradia pri náhodnom výbere dvoch pozorovaní za podmienky, že nezávislá premenná nemá zhodné poradia. Interpretácia je podobná ako u .

Na inom prístupe je založený

  1. Spearmanov koeficient poradovej korelácie
    Nech sú pôvodné pozorovania; označme ich poradia v usporiadaní podľa hodnôt faktoru A a ich poradia v usporiadaní podľa hodnôt faktoru B (t.j. vektory i   sú permutáciou čísel ). Pritom keď sú niektoré pozorovania zhodné (t.j. patria do tej istej kategórie resp. ), priradíme im všetkým rovnaké priemerné poradie.

    Veta 11.2.3: Majme kontingenčnú tabuľku. Ak všetky pozorovania usporiadame podľa hodnôt faktora A, potom priemerné poradie všetkých hodnôt patriacich do triedy je Ak všetky pozorovania usporiadame podľa hodnôt faktora B, potom priemerné poradie všetkých hodnôt patriacich do triedy je

    Dôkaz: Zrejme poradie prvej hodnoty v triede je a poradie poslednej hodnoty je . Všetky poradia v danej triede tvoria aritmetickú postupnosť. Ich priemer musí byť priemer čísel , ktorý je . Podobne to platí pre faktor B. QED.

    Môžeme si ešte všimnúť, že priemerné poradie pozorovaní z triedy sa dá tiež napísať v tvare , kde a sú kumulatívne a absolútne početnosti príslušných tried.

    Definícia 11.2.4: Spearmanov koeficient poradovej korelácie je

    Zrejme je . Je to vlastne Pearsonov korelačný koeficient poradí . Interpretácia je preto rovnaká ako u korelačného koeficientu.

11.3 Ďalšie miery závislosti

Ak obidva merané znaky sú kvantitatívne, používame ako mieru ich závislosti (Pearsonov) korelačný koeficient. Tento pojem bol vysvetlený v kapitole 6. Tu uvedieme vetu, ktorá nám umožňuje testovať hypotézu, že korelačný koeficient je nulový.

Tvrdenie 11.3.1: Nech je náhodný výber z normálneho rozdelenia s kladnými rozptylmi a nulovou koreláciou. Ak , má veličina Studentovo rozdelenie .

Ak teda je , kde je zvolená hladina testu, zamietame hypotézu o nulovosti korelačného koeficientu. V opačnom prípade ju nezamietame.

Ak jeden znak je kvalitatívny (nominálny) a druhý kvantitatívny, meriame ich závislosť koeficientom nelineárnej korelácie . Nech kvalitatívny znak je ; označíme hodnoty znaku príslušné -tej hodnote znaku , t.j. ležiace v kategórii . Ďalej označíme a , kde . Potom definujeme

Definícia 11.3.2: Koeficient nelineárnej korelácie je , kde

Je to vlastne percento celkového rozptylu závisle premennej, ktoré je vysvetlené rozptylom medzi jednotlivými kategóriami nezávisle premennej. Zrejme je . Koeficient je asymetrický. Ak sú kategórie kvalitatívnej premennej numericky kódované, je možné vyrátať aj symetrický koeficient ako priemer dvoch asymetrických, ale jeho interpretácia je nejasná. nezávisí od poradia kategórií, ak je kvalitatívna premenná ordinálna.

Niekedy je zadeľovanie do kategórií subjektívne. V takej situácii nás môže zaujímať miera zhody dvoch hodnotiteľov. Vtedy sú zrejme znaky A a B v podstate totožné, t.j. majú rovnaké kategórie a . Mierou zhody je Cohenovo :

Definícia 11.3.3: Pre kontingenčnú tabuľku s rovnakými kategóriami znakov A i B definujeme Cohenovo vzťahom

Platí . Hodnota 0 znamená, že zhoda hodnotiteľov je iba na náhodnej úrovni. Preto nás za normálnych okolností záporné hodnoty nezaujímajú. Hodnota 1 znamená úplnú zhodu hodnotiteľov.