Jedinica 5 / 11

Regresijska dijagnostika i poremećaji: kolinearnost, heteroskedastičnost, autokorelacija

Dobici:

  • Sposobnost otkrivanja kršenja kao što su multikolinearnost, heteroskedastičnost i autokorelacija s dijagnostičkim testovima i grafikom s podrškom za umjetnu inteligenciju
  • Sposobnost razlikovanja da li svako kršenje iskrivljuje procjene koeficijenata ili standardne pogreške i odabir odgovarajuće korekcije (robustna standardna pogreška, transformacija, revizija modela)
  • Sposobnost provjere da popravci koje predlaže AI rješavaju problem umjesto da ga skrivaju i da se dijagnoza temelji na razumijevanju procesa koji proizvodi podatke

Čitanje rezultata regresije je jednostavno; Teško mu je vjerovati. Budući da su koeficijenti nepristrani, standardne pogreške su točne, a p-vrijednosti valjane, ovisno o ispunjenju pretpostavki koje smo uveli u prethodnoj jedinici. U ovoj jedinici pokrit ćemo tri najčešća kršenja: multikolinearnost, heteroskedastičnost i autokorelaciju. Za svaki ćemo odgovoriti na tri pitanja: kako ga dijagnosticirati, što se kvari (koeficijent ili standardna pogreška) i kako ga popraviti. Umjetna inteligencija (AI) generira kod za dijagnostičko testiranje i ispravak; ali vi odlučujete koje je kršenje doista problem i rješava li popravak problem.

Napravimo najkritičniju razliku unaprijed: Iskrivljuje li kršenje procjenu koeficijenta ili samo standardne pogreške? Ova razlika određuje rješenje. Problem koji iskrivljuje koeficijent (npr. povreda egzogenosti) zahtijeva ponovno promišljanje modela; Problem koji iskrivljuje samo standardnu ​​pogrešku (heteroskedastičnost, autokorelacija) često se rješava standardnom korekcijom pogreške. Uobičajena pogreška umjetne inteligencije je primijeniti tehniku ​​koja ispravlja standardnu ​​pogrešku i zatim proglasiti problem "riješenim"; dok temeljna struktura može još uvijek biti tu.

Kršenje 1: Multikolinearnost

zašto Dvije ili više eksplanatornih varijabli snažno su povezane jedna s drugom; na primjer, stavljanje i "ukupnih godina iskustva" i "dobi" u model. Budući da te varijable nose gotovo iste informacije, model ih teško razdvaja.

Što kvari? Procjene koeficijenata ostaju nepristrane, ali postaju nestabilne: standardne pogreške postaju prenapuhane, intervali pouzdanosti se šire, pojedinačni koeficijenti mogu se pokazati beznačajnima, dok model može biti značajan u cjelini (F-test). Mala promjena podataka značajno pomiče koeficijent.

Kako se dijagnosticira? VIF (faktor inflacije varijance): za svaku varijablu mjeri koliko je ta varijabla objašnjena ostalima. Kao grubi prag, VIF > 5 (oko 10) zahtijeva oprez. Dodatno, ispituje se korelacijska matrica između eksplanatornih varijabli.

Kako popraviti? Izbacivanje jedne od koreliranih varijabli, njihovo kombiniranje (izrada indeksa) ili (ako teorija zahtijeva) zadržavanje obje i izbjegavanje tumačenja pojedinačnih koeficijenata. Koja varijabla ostaje teoretska je odluka, a ne statistička — AI predlaže eliminaciju, vi dajete obrazloženje.

Kršenje 2: Heteroskedastičnost

zašto Varijanca izraza pogreške nije konstantna u promatranjima. Tipičan primjer: kako se prihod povećava, tako raste i disperzija oko potrošnje; Formira se obrazac "lijevka", uzak kod niskih prihoda i širok kod visokih prihoda.

Što kvari? Procjene koeficijenata opet su nepristrane — ovo je važno. Ono što se iskrivljuje su standardne pogreške: one se netočno izračunavaju, pa p-vrijednosti i intervali pouzdanosti postaju nepouzdani. Dakle, vaš koeficijent je u desnoj sredini, ali odgovor na pitanje "koliko ste sigurni" je pogrešan.

Kako se dijagnosticira? Raspršeni dijagram reziduala u odnosu na predviđene vrijednosti (traži se uzorak lijevka) i formalni testovi: Breusch-Pagan test, White test. Mala p-vrijednost govori "nema konstantne varijance".

Kako popraviti? Najčešće i najpraktičnije rješenje je korištenje robusnih standardnih pogrešaka (robustne / heteroskedastičnosti konzistentne, HC standardne pogreške): ne mijenja koeficijent, ispravlja standardnu ​​pogrešku za kršenje. Alternativa: transformacija zavisne varijable (kao što je log) ili ponderirani LCM. Robusne standardne pogreške danas su postale gotovo standardne u empirijskom radu.

Kršenje 3: Autokorelacija

zašto Izrazi pogreške nisu neovisni jedni o drugima; Najčešći je u vremenskim serijama: pogreška jednog razdoblja utječe na sljedeće (npr. rezidual ostaje pozitivan u razdobljima nakon šoka).

Što kvari? Opet, prvenstveno iskrivljuje standardne pogreške (često ih podcjenjuje, stvarajući lažni značaj); Koeficijenti ostaju nepristrani u LCC-u, ali gube svoju učinkovitost.

Kako se dijagnosticira? Durbin-Watsonov test (za autokorelaciju prvog reda), Breusch-Godfreyjev test (općenitiji) i grafikoni reziduala u odnosu na zaostale vrijednosti.

Kako popraviti? Newey-West (HAC — otporan na autokorelaciju i heteroskedastičnost) standardne pogreške, dodavanje lagiranih članova u model ili prebacivanje na odgovarajući model vremenske serije (jedinica 8).

Oprez: heteroskedastičnost i autokorelacija iskrivljuju standardnu ​​pogrešku, a ne koeficijent. Stoga, prije nego kažete "koeficijent je bio značajan", provjerite je li standardna pogreška ispravno izračunata; Inače bi smislenost mogla biti iluzija.

usporedna tablica

kršenje

Što se lomi

Dijagnoza

Uobičajeno rješenje

multi-link

Povećava standardne pogreške, čini koeficijent nestabilnim

VIF, korelacijska matrica

Oduzmi/kombiniraj varijablu (po teoriji)

heteroskedastičnost

Standardne pogreške (nepristrani koeficijent)

Breusch-Pagan, Bijela, zaostala parcela

Robusna (HC) standardna pogreška, konverzija

Autokorelacija

Standardne pogreške (nepristrani koeficijent)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), odgođeni termin

Četiri upita za kopiranje

1. Kompletan dijagnostički paket:

Vaša uloga: pomoćnik regresijske dijagnostike. Želim R kod, ne dijelim podatke. Model: lm(rashod ~ prihod + dob + regija, podaci = podaci). Zadatak: (1) Izračunajte VIF, (2) Testirajte heteroskedastičnost s Breusch-Paganom i prikazom rezidualne procjene, (3) Testirajte autokorelaciju s Durbin-Watsonom. U retku komentara objasnite koje kršenje mjeri svaki test i kako protumačiti p-vrijednost. Ispravak PRIJAVA; postaviti dijagnozu.

2. Robusna standardna pogreška:

Napišite R kod koji reproducira tablicu koeficijenata sa standardnim pogreškama otpornim na heteroskedastičnost (HC3) za isti model (paketi sendvič + lmtest). Napravite tablicu koja prikazuje klasične i robusne standardne pogreške jednu pored druge kako bih mogao vidjeti razliku. Naglasite u retku komentara da se koeficijenti NE mijenjaju, već samo standardne pogreške.

3. Višestruki pregled:

Daj mi RAZMISLI popis varijabli s visokim VIF-om: koje varijable bi teoretski mogle mjeriti istu stvar, koje imaju jake razloge za zadržavanje. NEMOJTE se automatski kvalificirati; Odluku ću donijeti na temelju teorije. Također objasnite zašto bi moglo biti pogrešno samo pogledati VIF i dodijeliti varijable.

4. Ispravak autokorelacije:

U mojoj regresiji vremenske serije, Breusch-Godfreyjeva p-vrijednost bila je 0,001. Napišite R kod koji daje tablicu koeficijenata sa Newey-West (HAC) standardnim pogreškama i objasnite kako odabrati kašnjenje. Imajte na umu da ova korekcija ne rješava UZROK autokorelacije, ona samo ispravlja zaključak.

Slab upit / Jak upit

Slab upit:

Postoji li problem s mojom regresijom? Ako je tako, popravite to.

Ovaj upit stavlja AI u način rada slijepog "popravljanja": može preskočiti testove i primijeniti izravnu transformaciju ili eliminaciju varijable, bez da vam pokaže koje kršenje zapravo postoji i što je pokvarilo.

Snažan upit:

Vaša uloga: pomoćnik u dijagnostici, a ne autokorektor. Prvi test za TRI prekršaja zasebno (multikolinearnost, heteroskedastičnost, autokorelacija) i za svaki: koji test, kako očitati rezultat, krši li KOEFICIJENT ili STANDARDNU POGREŠKU. Zatim samo predložite popravak za kršenje koje je STVARNO otkriveno i objasnite kako mogu potvrditi da je popravak riješio problem.

Razlika: jaka volja forsira dijagnozu prije korekcije i zahtijeva jasnu razliku između "onog što kvari".

tri mini kućišta

Slučaj 1 — lažni značaj (heteroskedastičnost). Jedan je analitičar otkrio da je koeficijent prihoda u modelu potrošnje ~ prihoda "značajan" na p=0,01. Ali sada je na njegovom grafikonu bio jasan uzorak lijevka. Kada su primijenjene robusne standardne pogreške, p-vrijednost je porasla na 0,09; gubi se smisao. Prvi rezultat bila je iluzija koju je stvorila pogrešno izračunata standardna pogreška. Lekcija: značaju se ne može vjerovati bez ispravljanja standardne pogreške.

Slučaj 2 — Slijepa eliminacija varijable. Jedan je učenik izostavio varijablu "iskustvo" iz modela jer je njegov VIF bio visok; Koeficijenti su "očišćeni", ali je teorijsko značenje modela iskrivljeno jer je iskustvo bilo glavna varijabla od interesa. Ispravan način: oduzmite godine i zadržite iskustvo ili kombinirajte to dvoje. Lekcija: VIF prag sam po sebi nije osnova za eliminaciju; teorija određuje.

Slučaj 3 — Nesigurnost skrivena autokorelacijom. U studiji vremenskih serija Durbin-Watson je zanemaren; koeficijent se činio vrlo "precizan" (uzak interval pouzdanosti). Kada su primijenjene Newey-West standardne pogreške, interval pouzdanosti se udvostručio, a preporuka politike postala je mnogo konzervativnija. Lekcija: autokorelacija može podcijeniti neizvjesnost.

Uobičajene greške

  • Zamijeniti kršenje koje iskrivljuje standardnu pogrešku kao problem koeficijenta. Heteroskedastičnost i autokorelacija ne iskrivljuju koeficijent; Nemojte paničariti i nepotrebno obnavljati model, prvo ispravite standardnu ​​pogrešku.
  • Gledanje VIF-a i slijepo dodjeljivanje varijabli. Uklanjanje teoretski potrebne varijable samo zato što je VIF visok čini model besmislenim.
  • Ne provjerava se ispravak. Nakon primjene robusne standardne pogreške, provjerite ispravlja li kršenje zapravo zaključak; Nemojte reći "Primijenio sam, gotovo je".
  • Popravak bez dijagnoze. Primjena transformacije/eliminacije bez testiranja koje kršenje postoji može "riješiti" problem tamo gdje ne postoji i sakriti onaj koji postoji.
  • Postavljanje popravka na mjesto zašto. Newey-West ne rješava uzrok autokorelacije; samo popravlja zaključak. Ako postoji strukturni problem, model se mora promijeniti.
Savjet: Za svaki prekršaj zapitajte se "rušava li to koeficijent ili moje povjerenje u njega?" Ako je odgovor "povjerenje", rješenje je gotovo uvijek standardno ispravljanje pogrešaka, a ne ponovna izgradnja modela.

Ukratko

Regresijska dijagnostika je preduvjet za vjerovanje izlazu. Od tri uobičajene povrede, multikolinearnost čini koeficijent nestabilnim i povećava standardnu ​​pogrešku; Heteroskedastičnost i autokorelacija, s druge strane, ostavljaju koeficijent nepristranim i samo iskrivljuju standardnu ​​pogrešku. AI generira dijagnostički kod i kod za popravak, ali vi odlučujete koje je kršenje pravi problem, koja varijanta ostaje teoretska i rješava li popravak problem. Ni panika ni slijepa korekcija nisu ispravni bez pojašnjenja razlike između "onoga što se kvari".

Zadatak aplikacije

Postavite multivarijatnu regresiju i zatražite od umjetne inteligencije kod koji proizvodi samo dijagnoze (bez ispravaka): VIF, Breusch-Pagan/White i Durbin-Watson/Breusch-Godfrey. Za svaki test protumačite rezultat i naznačite iskrivljuje li kršenje koeficijent ili standardnu ​​pogrešku. Za stvarno kršenje koje otkrijete (npr. heteroskedastičnost), primijenite robusne standardne pogreške i usporedite klasične i robusne rezultate; Pokažite da se koeficijent ne mijenja, ali se mijenja standardna pogreška. U jednom odlomku navedite kako je ispravak utjecao na vaš rezultat značajnosti.

popis za provjeru

  • [ ] Zasebno sam testirao tri kršenja (multikolinearnost, heteroskedastičnost, autokorelacija).
  • [ ] Za svako kršenje razlikovao sam iskrivljuje li koeficijent ili standardnu ​​grešku.
  • [ ] Eliminaciju varijable u multikolinearnosti temeljio sam na teoriji, ne samo na VIF-u.
  • [ ] Koristio sam standardnu ​​pogrešku s robusnošću prema heteroskedastičnosti/autokorelaciji (HC/HAC).
  • [ ] Nakon ispravka, provjerio sam i potvrdio utjecaj kršenja na moj zaključak.
  • [ ] Izvještavam kako je na moj rezultat značajnosti utjecalo standardno ispravljanje pogrešaka.