Jednotka 5 / 11

Regresná diagnostika a porušenia: Kolinearita, Heteroscedasticita, Autokorelácia

zisky:

  • Schopnosť odhaliť porušenia, ako je multikolinearita, heteroskedasticita a autokorelácia pomocou diagnostických testov a grafiky s podporou umelej inteligencie
  • Schopnosť rozlíšiť, či každé porušenie skresľuje odhady koeficientov alebo štandardné chyby a zvoliť vhodnú korekciu (robustná štandardná chyba, transformácia, revízia modelu)
  • Schopnosť overiť, že opravy navrhnuté AI riešia problém a nie ho skrývajú a že diagnóza je založená na pochopení procesu, ktorý vytvára údaje

Čítanie výstupu regresie je jednoduché; Je ťažké mu dôverovať. Pretože koeficienty sú nezaujaté, štandardné chyby sú správne a p-hodnoty sú platné v závislosti od toho, či sú splnené predpoklady, ktoré sme zaviedli v predchádzajúcej jednotke. V tejto časti sa budeme zaoberať tromi najbežnejšími porušeniami: multikolinearitou, heteroskedasticitou a autokoreláciou. Ku každému odpovieme na tri otázky: ako ho diagnostikovať, čo pokazí (koeficient alebo štandardná chyba) a ako to opraviť. Umelá inteligencia (AI) generuje kód na diagnostické testovanie a korekciu; ale vy rozhodnete, ktoré porušenie je skutočne problémom a či oprava problém vyrieši.

Urobme najkritickejšie rozlíšenie dopredu: Skresľuje porušenie odhad koeficientu alebo len štandardné chyby? Toto rozlíšenie určuje riešenie. Problém, ktorý skresľuje koeficient (napr. porušenie exogenity), si vyžaduje prehodnotenie modelu; Problém, ktorý len skresľuje štandardnú chybu (heteroscedasticita, autokorelácia), sa často rieši štandardnou korekciou chýb. Bežnou chybou AI je použitie techniky, ktorá opraví štandardnú chybu a potom vyhlási problém za „vyriešený“; zatiaľ čo základná štruktúra tam stále môže byť.

Porušenie 1: Multikolinearita

prečo? Dve alebo viac vysvetľujúcich premenných spolu úzko súvisia; napríklad vloženie „celkového počtu rokov skúseností“ a „vek“ do modelu. Keďže tieto premenné nesú takmer rovnaké informácie, model má problém ich oddeliť.

Čo to láme? Odhady koeficientov zostávajú nezaujaté, ale stávajú sa nestabilnými: štandardné chyby sa zväčšujú, intervaly spoľahlivosti sa rozširujú, jednotlivé koeficienty sa môžu ukázať ako nevýznamné, zatiaľ čo model môže byť celkovo významný (F-test). Malá zmena údajov výrazne posunie koeficient.

Ako sa diagnostikuje? VIF (faktor variačnej inflácie): pre každú premennú meria, do akej miery je táto premenná vysvetlená ostatnými. Približná hranica VIF > 5 (asi 10) vyžaduje opatrnosť. Okrem toho sa skúma korelačná matica medzi vysvetľujúcimi premennými.

Ako opraviť? Vypustenie jednej z korelovaných premenných, ich kombinovanie (vytvorenie indexu) alebo (ak to teória vyžaduje) ponechanie oboch a vyhýbanie sa interpretácii jednotlivých koeficientov. Ktorá premenná zostane, je teoretické rozhodnutie, nie štatistické – AI navrhuje elimináciu, vy poskytnete odôvodnenie.

Porušenie 2: Heteroscedasticita

prečo? Rozptyl chybového termínu nie je konštantný naprieč pozorovaniami. Typický príklad: so zvyšujúcim sa príjmom sa zvyšuje aj rozptyl okolo výdavkov; Vytvára sa „lievik“ vzor, ​​úzky pri nízkom príjme a široký pri vysokom príjme.

Čo to láme? Odhady koeficientov sú opäť nezaujaté – to je dôležité. Čo sa skresľuje, sú štandardné chyby: sú vypočítané nesprávne, takže p-hodnoty a intervaly spoľahlivosti sa stávajú nespoľahlivými. Váš koeficient je teda v správnom strede, ale odpoveď na otázku „nakoľko ste si istý“ je nesprávna.

Ako sa diagnostikuje? Bodový graf rezíduí oproti predpokladaným hodnotám (hľadaný vzor lievika) a formálne testy: Breusch-Pagan test, White test. Malá p-hodnota hovorí, že „žiadny konštantný rozptyl“.

Ako opraviť? Najbežnejším a praktickým riešením je použitie robustných štandardných chýb (robustné / heteroskedasticity-konzistentné, HC štandardné chyby): nemení koeficient, opravuje štandardnú chybu za porušenie. Alternatíva: transformácia závislej premennej (ako log) alebo váženého LCM. Robustné štandardné chyby sa dnes v empirickej práci stali takmer štandardom.

Porušenie 3: Autokorelácia

prečo? Chybové pojmy nie sú navzájom nezávislé; Najčastejšie sa vyskytuje v časových radoch: chyba jedného obdobia ovplyvňuje ďalšie obdobie (napr. zvyšok zostáva kladný v obdobiach nasledujúcich po šoku).

Čo to láme? Opäť primárne skresľuje štandardné chyby (často ich podhodnocuje, čím vytvára falošný význam); Koeficienty zostávajú v LCC nezaujaté, ale strácajú svoju účinnosť.

Ako sa diagnostikuje? Durbin-Watsonov test (pre autokoreláciu prvého rádu), Breusch-Godfreyov test (všeobecnejší) a grafy rezíduí oproti oneskoreným hodnotám.

Ako opraviť? Štandardné chyby Newey-West (HAC – odolné voči autokorelácii a heteroskedasticite), pridanie oneskorených pojmov do modelu alebo prepnutie na vhodný model časového radu (jednotka 8).

Pozor: Heteroscedasticita a autokorelácia skresľujú štandardnú chybu, nie koeficient. Preto sa predtým, ako poviete „koeficient bol významný“, uistite, že štandardná chyba je vypočítaná správne; Inak môže byť zmysluplnosť ilúziou.

porovnávacia tabuľka

porušenie

Čo sa zlomí

Diagnóza

Spoločné riešenie

multi-link

Nafukuje štandardné chyby, robí koeficient nestabilným

VIF, korelačná matica

Odčítajte/kombinujte premennú (podľa teórie)

heteroskedasticita

Štandardné chyby (koeficient nezaujatý)

Breusch-Pagan, White, zvyškový pozemok

Robustná (HC) štandardná chyba, konverzia

Autokorelácia

Štandardné chyby (koeficient nezaujatý)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), oneskorený termín

Štyri kopírovateľné výzvy

1. Kompletný diagnostický balík:

Vaša úloha: asistent regresnej diagnostiky. Chcem kód R, údaje nezdieľam. Model: lm(výdavky ~ príjem + vek + región, údaje = údaje). Úloha: (1) Vypočítajte VIF, (2) Testujte heteroskedasticitu pomocou Breusch-Pagan a graf reziduálneho odhadu, (3) Testujte autokoreláciu s Durbin-Watsonom. Vysvetlite v riadku komentára, aké porušenie meria každý test a ako interpretovať hodnotu p. Oprava APLIKÁCIA; stanoviť diagnózu.

2. Robustná štandardná chyba:

Napíšte R kód, ktorý reprodukuje tabuľku koeficientov so štandardnými chybami odolnými voči heteroskedasticite (HC3) pre rovnaký model (balíky sendvič + lmtest). Vytvorte tabuľku zobrazujúcu klasické a robustné štandardné chyby vedľa seba, aby som videl rozdiel. V komentári zdôraznite, že koeficienty sa NEMENIA, menia sa iba štandardné chyby.

3. Kontrola viacerých odkazov:

Dajte mi zoznam premenných THINK s vysokým VIF: ktoré premenné by teoreticky mohli merať to isté, ktoré z nich majú silný dôvod na ich zachovanie. NEKvalifikujte sa automaticky; Rozhodnem sa na základe teórie. Vysvetlite tiež, prečo by mohlo byť nesprávne pozrieť sa na VIF a priradiť premenné.

4. Korekcia autokorelácie:

V mojej regresii časových radov bola Breusch-Godfreyova p-hodnota 0,001. Napíšte R kód, ktorý vytvorí tabuľku koeficientov so štandardnými chybami Newey-West (HAC) a vysvetlite, ako zvoliť oneskorenie. Všimnite si, že táto oprava nerieši PRÍČINU autokorelácie, iba opravuje záver.

Slabá výzva / Silná výzva

Slabá výzva:

Je problém s mojou regresiou? Ak áno, opravte to.

Táto výzva uvedie AI ​​do slepého režimu „opravy“: môže preskočiť testy a použiť priamu transformáciu alebo premennú elimináciu bez toho, aby vám ukázala, ktoré porušenie skutočne existuje a čo porušilo.

Výkonná výzva:

Vaša úloha: diagnostický asistent, nie autokorektor. Najprv otestujte TRI porušenia samostatne (multikolinearita, heteroskedasticita, autokorelácia) a pre každý: ktorý test, ako čítať výsledok, porušuje KOEFICIENT alebo ŠTANDARDNÚ CHYBU. Potom len navrhnite opravu porušenia, ktoré bolo SKUTOČNE zistené, a vysvetlite, ako môžem overiť, či oprava problém vyriešila.

Rozdiel: silná vôľa si vynucuje diagnostiku pred korekciou a vyžaduje jasné rozlíšenie medzi tým, „čo rozbije“.

tri mini prípady

Prípad 1 – Falošný význam (heteroscedasticita). Jeden analytik zistil, že koeficient príjmov v modeli výdavkov a príjmov je „významný“ pri p=0,01. Ale teraz bol na jeho grafe zreteľný vzor lievika. Keď sa použili robustné štandardné chyby, p-hodnota sa zvýšila na 0,09; zmysluplnosť sa stráca. Prvým výsledkom bola ilúzia vytvorená nesprávne vypočítanou štandardnou chybou. Ponaučenie: Významu nemožno dôverovať bez opravy štandardnej chyby.

Prípad 2 – eliminácia slepej premennej. Jeden študent vypustil z modelu premennú „skúsenosť“, pretože jeho VIF bolo vysoké; Koeficienty boli „vyčistené“, ale teoretický význam modelu bol skreslený, pretože hlavnou premennou záujmu bola skúsenosť. Správny spôsob: odpočítajte vek a ponechajte si skúsenosti, alebo oboje skombinujte. Poučenie: Samotný prah VIF nie je dôvodom na vylúčenie; teória určuje.

Prípad 3 – Neistota skrytá autokoreláciou. V štúdii časových radov bol Durbin-Watson ignorovaný; koeficient sa zdal veľmi "presný" (úzky interval spoľahlivosti). Keď sa použili štandardné chyby Newey-West, interval spoľahlivosti sa zdvojnásobil a politické odporúčanie sa stalo oveľa konzervatívnejším. Poučenie: autokorelácia môže podceňovať neistotu.

Časté chyby

  • Pomýlené porušenie, ktoré skresľuje štandardnú chybu, ako problém koeficientu. Heteroscedasticita a autokorelácia neskresľujú koeficient; Nepodliehajte panike a zbytočne prestavujte model, najskôr opravte štandardnú chybu.
  • Pohľad na VIF a slepé priraďovanie premenných. Odstránenie teoreticky potrebnej premennej len preto, že VIF je vysoké, robí model bezvýznamným.
  • Neoveruje sa oprava. Po použití robustnej štandardnej chyby skontrolujte, či porušenie skutočne opravuje záver; Nehovorte „aplikoval som to, hotovo“.
  • Oprava bez diagnostiky. Aplikácia transformácie/eliminácie bez testovania, ktoré porušenie existuje, môže „vyriešiť“ problém tam, kde neexistuje, a skryť ten, ktorý existuje.
  • Uvedenie opravy na miesto prečo. Newey-West nerieši príčinu autokorelácie; len opravuje záver. Ak sa vyskytne štrukturálny problém, model sa musí zmeniť.
Tip: Pri každom porušení si položte otázku: „ničí to koeficient alebo moju dôveru k nemu?“ Ak je odpoveď „dôvera“, riešením je takmer vždy štandardná oprava chýb, nie prestavba modelu.

V súhrne

Regresná diagnostika je predpokladom dôveryhodnosti výstupu. Z troch bežných porušení multikolinearita spôsobuje, že koeficient je nestabilný a zvyšuje štandardnú chybu; Na druhej strane heteroskedasticita a autokorelácia ponechávajú koeficient nezaujatý a iba skresľujú smerodajnú chybu. Umelá inteligencia generuje diagnostický a opravný kód, ale vy rozhodujete, ktoré porušenie je skutočným problémom, ktorý variant zostáva teoretický a či oprava problém vyrieši. Ani panika, ani slepá korekcia nie sú správne bez objasnenia rozdielu medzi tým, „čo sa láme“.

Aplikačná úloha

Nastavte viacrozmernú regresiu a požiadajte AI ​​o kód, ktorý vytvára iba diagnózy (bez opráv): VIF, Breusch-Pagan/White a Durbin-Watson/Breusch-Godfrey. Pre každý test interpretujte výsledok a uveďte, či porušenie skresľuje koeficient alebo štandardnú chybu. Pre skutočné porušenie, ktoré zistíte (napr. heteroskedasticita), použite robustné štandardné chyby a postavte klasické a robustné výsledky vedľa seba; Ukážte, že koeficient sa nemení, ale mení sa smerodajná chyba. V jednom odseku uveďte, ako oprava ovplyvnila váš výsledok významnosti.

kontrolný zoznam

  • [ ] Tri porušenia (multikolinearita, heteroskedasticita, autokorelácia) som testoval samostatne.
  • [ ] Pri každom porušení som rozlišoval, či skresľuje koeficient alebo smerodajnú chybu.
  • [ ] Elimináciu premenných v multikolinearite som založil na teórii, nielen na VIF.
  • [ ] Použil som štandardnú chybu s robustnosťou voči heteroskedasticite/autokorelácii (HC/HAC).
  • [ ] Po oprave som skontroloval a overil vplyv porušenia na môj záver.
  • [ ] Uvádzam, ako môj výsledok významnosti ovplyvnila štandardná oprava chýb.