Üksus 5 / 11

Regressioonidiagnostika ja rikkumised: kollineaarsus, heteroskedastilisus, autokorrelatsioon

Kasu:

  • Võimalus tuvastada rikkumisi, nagu multikollineaarsus, heteroskedastilisus ja autokorrelatsioon tehisintellekti toega diagnostiliste testide ja graafikaga
  • Võimalus eristada, kas iga rikkumine moonutab koefitsientide hinnanguid või standardvigu, ja valida sobiv parandus (tugev standardviga, teisendus, mudeli revisjon)
  • Võimalus kontrollida, kas tehisintellekti pakutud parandused lahendavad probleemi, mitte ei varja seda ja et diagnoos põhineb andmete loomise protsessi mõistmisel

Regressiooni väljundi lugemine on lihtne; Teda on raske usaldada. Kuna koefitsiendid on erapooletud, on standardvead õiged ja p-väärtused kehtivad, olenevalt sellest, kas eelmises ühikus kasutatud eeldused on täidetud. Selles üksuses käsitleme kolme kõige levinumat rikkumist: multikollineaarsus, heteroskedastilisus ja autokorrelatsioon. Vastame igaühe puhul kolmele küsimusele: kuidas seda diagnoosida, mida see rikub (koefitsient või standardviga) ja kuidas seda parandada. Tehisintellekt (AI) genereerib koodi diagnostiliseks testimiseks ja korrigeerimiseks; kuid otsustate, milline rikkumine tegelikult probleem on ja kas parandus lahendab probleemi.

Teeme kohe kõige kriitilisema vahe: kas rikkumine moonutab koefitsiendi hinnangut või ainult standardvigu? See eristus määrab lahenduse. Koefitsienti moonutav probleem (nt eksogeensuse rikkumine) nõuab mudeli ümbermõtestamist; Probleem, mis moonutab ainult standardviga (heteroskedastilisus, autokorrelatsioon), lahendatakse sageli standardvea paranduse abil. AI levinud viga on kasutada tehnikat, mis parandab standardvea ja kuulutab seejärel probleemi lahendatuks; alusstruktuur võib siiski alles olla.

Rikkumine 1: Multikollineaarsus

Miks? Kaks või enam selgitavat muutujat on üksteisega tugevalt seotud; näiteks panna mudelisse nii "aastate kogukogemus" kui ka "vanus". Kuna need muutujad kannavad peaaegu sama teavet, on mudelil raskusi nende eraldamisega.

Mis see katki läheb? Koefitsientide hinnangud jäävad erapooletuks, kuid muutuvad ebastabiilseks: standardvead suurenevad, usaldusvahemikud laienevad, üksikud koefitsiendid võivad osutuda ebaolulisteks, samas kui mudel võib üldiselt olla oluline (F-test). Väike andmete muutus muudab koefitsienti oluliselt.

Kuidas seda diagnoositakse? VIF (Variance Inflation Factor): iga muutuja puhul mõõdab see, kui palju see muutuja on teistega seletatav. Ligikaudne läviväärtus, VIF > 5 (umbes 10), nõuab ettevaatust. Lisaks uuritakse selgitavate muutujate vahelist korrelatsioonimaatriksit.

Kuidas parandada? Ühe korrelatsioonimuutuja mahajätmine, nende kombineerimine (indeksi tegemine) või (kui teooria nõuab) mõlema hoidmine ja üksikute koefitsientide tõlgendamise vältimine. Milline muutuja jääb, on teoreetiline, mitte statistiline otsus – AI soovitab kõrvaldada, teie esitate põhjenduse.

Rikkumine 2: Heteroskedastilisus

Miks? Vealiikme dispersioon ei ole vaatluste lõikes konstantne. Tüüpiline näide: sissetulekute kasvades suureneb hajuvus kulutuste ümber; Moodustub "lehtri" muster, väikese sissetuleku korral kitsas ja suure sissetuleku korral lai.

Mis see katki läheb? Koefitsientide hinnangud on jällegi erapooletud – see on oluline. Mis moonutatakse, on standardvead: need on valesti arvutatud, nii et p-väärtused ja usaldusvahemikud muutuvad ebausaldusväärseks. Seega on teie koefitsient õiges keskel, kuid vastus küsimusele "kui enesekindel sa oled" on vale.

Kuidas seda diagnoositakse? Jääkide hajuvusdiagramm võrreldes prognoositud väärtustega (otsitud lehtri muster) ja ametlikud testid: Breusch-Pagan test, White test. Väike p-väärtus ütleb "konstantset dispersiooni pole".

Kuidas parandada? Levinuim ja praktilisem lahendus on kasutada robustseid standardvigu (robust / heteroskedasticity-consistent, HC standard errors): see ei muuda koefitsienti, parandab standardvea rikkumise korral. Alternatiiv: sõltuva muutuja (nt log) või kaalutud LCM-i teisendamine. Tugevad standardvead on tänapäeval muutunud empiirilises töös peaaegu vaikimisi.

3. rikkumine: autokorrelatsioon

Miks? Veaterminid ei ole üksteisest sõltumatud; Kõige sagedamini esineb see aegridades: ühe perioodi viga mõjutab järgmist (nt jääk jääb šokile järgnevatel perioodidel positiivseks).

Mis see katki läheb? Jällegi, see moonutab peamiselt standardvigu (sageli alahinnab neid, luues vale tähenduse); Koefitsiendid jäävad LCC-s erapooletuks, kuid kaotavad oma tõhususe.

Kuidas seda diagnoositakse? Durbin-Watsoni test (esimest järku autokorrelatsiooni jaoks), Breusch-Godfrey test (üldisem) ja jääkide ja mahajäänud väärtuste graafikud.

Kuidas parandada? Newey-Westi (HAC — autokorrelatsioonile ja heteroskedastilisusele vastupidav) standardvead, hilinenud terminite lisamine mudelisse või üleminek sobivale aegrea mudelile (ühik 8).

Ettevaatust: Heteroskedastilisus ja autokorrelatsioon moonutavad standardviga, mitte koefitsienti. Seetõttu veenduge enne "koefitsient oli oluline" ütlemist, et standardviga on õigesti arvutatud; Vastasel juhul võib tähenduslikkus olla illusioon.

võrdlustabel

rikkumine

Mis puruneb

Diagnoos

Ühine lahendus

mitme lingiga

Paisutab standardvead, muudab koefitsiendi ebastabiilseks

VIF, korrelatsioonimaatriks

Muutuja lahutamine/ühendamine (teooria järgi)

heteroskedastilisus

Standardvead (koefitsient erapooletu)

Breusch-Pagan, valge, jääkkrunt

Tugev (HC) standardviga, teisendus

Autokorrelatsioon

Standardvead (koefitsient erapooletu)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), hilinenud tähtaeg

Neli kopeeritavat viipa

1. Täielik diagnostikapakett:

Sinu roll: regressioonidiagnostika assistent. Ma tahan R-koodi, ma ei jaga andmeid. Mudel: lm(kulu ~ tulu + vanus + piirkond, andmed = andmed). Ülesanne: (1) Arvutage VIF, (2) Testige heteroskedastilisust Breusch-Pagani ja jääkhinnangu diagrammiga, (3) Testige autokorrelatsiooni Durbin-Watsoniga. Selgitage kommentaarireal, millist rikkumist iga test mõõdab ja kuidas p-väärtust tõlgendada. Parandustaotlus; diagnoosi panna.

2. Tugev standardviga:

Kirjutage R-kood, mis reprodutseerib sama mudeli (sandwich + lmtest paketid) heteroskedastilisuse-robust (HC3) standardvigadega koefitsientide tabeli. Koostage tabel, mis näitab kõrvuti klassikalisi ja tugevaid standardvigu, et saaksin näha erinevust. Rõhutage kommentaarireas, et koefitsiendid EI muutu, muutuvad ainult standardvead.

3. Mitme lingi ülevaade:

Andke mulle THINK loend kõrge VIF-iga muutujatest: millised muutujad võiksid teoreetiliselt mõõta sama asja, milliste puhul on nende hoidmine põhjendatud. ÄRGE kvalifitseeruge automaatselt; Otsuse teen teooria põhjal. Samuti selgitage, miks võib olla vale vaadata lihtsalt VIF-i ja määrata muutujaid.

4. Autokorrelatsiooni parandus:

Minu aegrea regressioonis oli Breusch-Godfrey p-väärtus 0,001. Kirjutage R-kood, mis loob Newey-Westi (HAC) standardvigadega koefitsientide tabeli, ja selgitage, kuidas viivitust valida. Pange tähele, et see parandus ei lahenda autokorrelatsiooni PÕHJUST, vaid parandab lihtsalt järelduse.

Nõrk viip / Tugev viip

Nõrk viip:

Kas minu regressioonis on probleem? Kui jah, siis parandage see.

See viip paneb tehisintellekti pimedasse "parandusrežiimi": see võib testid vahele jätta ja rakendada otsest teisendust või muutuja kõrvaldamist, näitamata teile, milline rikkumine tegelikult eksisteerib ja mida see rikkus.

Võimas viip:

Teie roll: diagnostikaassistent, mitte automaatkorrektor. Esimene test KOLME rikkumise kohta eraldi (multikollineaarsus, heteroskedastilisus, autokorrelatsioon) ja igaühe kohta: milline test, kuidas tulemust lugeda, kas see rikub KOEFITSIENT või STANDARDVIGA. Seejärel soovitage ainult TEGELIKULT tuvastatud rikkumise parandamist ja selgitage, kuidas saan kontrollida, kas parandus probleemi lahendas.

Erinevus: tugev tahe sunnib diagnoosi enne korrigeerimist ja nõuab selget vahet "mida see rikub".

kolm minikarpi

Juhtum 1 – võlts tähendus (heteroskedastilisus). Üks analüütik leidis, et tulukoefitsient kulu-tulu mudelis on "oluline" p=0,01. Kuid nüüd oli tema diagrammil selge lehtermuster. Tugevate standardvigade rakendamisel tõusis p-väärtus 0,09-ni; tähenduslikkus on kadunud. Esimene tulemus oli illusioon, mille lõi valesti arvutatud standardvea. Õppetund: olulisust ei saa usaldada ilma standardviga parandamata.

Juhtum 2 – pimemuutuja kõrvaldamine. Üks õpilane jättis mudelist välja "kogemuse" muutuja, kuna tema VIF oli kõrge; Koefitsiendid "puhastati", kuid mudeli teoreetiline tähendus moonutati, kuna kogemus oli peamine huvipakkuv muutuja. Õige viis: lahutage vanus ja säilitage kogemus või ühendage need kaks. Õppetund: VIF-i lävi üksi ei anna põhjust kõrvaldamiseks; teooria määrab.

Juhtum 3 – autokorrelatsiooniga peidetud määramatus. Aegridade uuringus jäeti Durbin-Watson tähelepanuta; koefitsient tundus väga "täpne" (kitsas usaldusvahemik). Newey-Westi standardvigade rakendamisel usaldusvahemik kahekordistus ja poliitikasoovitus muutus palju konservatiivsemaks. Õppetund: autokorrelatsioon võib ebakindlust alahinnata.

Levinud vead

  • Standardviga moonutava rikkumise eksitamine koefitsiendi probleemiks. Heteroskedastilisus ja autokorrelatsioon ei moonuta koefitsienti; Ärge sattuge paanikasse ja ehitage mudel asjatult ümber, parandage esmalt standardviga.
  • VIF-i vaatamine ja muutujate pimesi määramine. Teoreetiliselt vajaliku muutuja eemaldamine lihtsalt seetõttu, et VIF on kõrge, muudab mudeli mõttetuks.
  • Parandust ei kontrollita. Pärast tugeva standardvea rakendamist kontrollige, kas rikkumine parandab järelduse; Ärge öelge: "Ma panin selle peale, see on tehtud".
  • Parandage ilma diagnoosita. Transformatsiooni/elimineerimise rakendamine ilma esineva rikkumise testimata võib "lahendada" probleemi seal, kus seda pole, ja peita probleemi, mis esineb.
  • Paranduse paika panemine, miks. Newey-West ei lahenda autokorrelatsiooni põhjust; see lihtsalt parandab järelduse. Kui on struktuurne probleem, tuleb mudelit muuta.
Näpunäide: iga rikkumise puhul küsige endalt: "Kas see hävitab koefitsiendi või minu usalduse selle vastu?" Kui vastus on "usaldus", on peaaegu alati lahendus standardvea parandamine, mitte mudeli taastamine.

Kokkuvõttes

Regressioondiagnostika on väljundi usaldamise eeldus. Kolmest tavalisest rikkumisest muudab multikollineaarsus koefitsiendi ebastabiilseks ja suurendab standardvea; Heteroskedastilisus ja autokorrelatsioon seevastu jätavad koefitsiendi erapooletuks ja moonutavad ainult standardviga. AI genereerib diagnostika- ja paranduskoodi, kuid teie otsustate, milline rikkumine on tegelik probleem, milline variant jääb teoreetiliseks ja kas parandus lahendab probleemi. Ei paanika ega pime korrigeerimine pole õiged, kui ei tee selgeks vahet "mis puruneb".

Rakenduse ülesanne

Seadistage mitme muutujaga regressioon ja küsige AI-lt koodi, mis toodab ainult diagnoose (parandusi pole): VIF, Breusch-Pagan/White ja Durbin-Watson/Breusch-Godfrey. Iga testi puhul tõlgendage tulemust ja märkige, kas rikkumine moonutab koefitsienti või standardviga. Tõelise tuvastatud rikkumise (nt heteroskedastilisus) korral rakendage tugevaid standardvigu ning kõrvutage klassikalisi ja tugevaid tulemusi; Näidake, et koefitsient ei muutu, kuid standardviga muutub. Ühes lõigus teatage, kuidas parandus teie olulisuse tulemust mõjutas.

kontrollnimekiri

  • [ ] Testisin kolme rikkumist (multikollineaarsus, heteroskedastilisus, autokorrelatsioon) eraldi.
  • [ ] Iga rikkumise puhul eristasin, kas see moonutab koefitsienti või standardviga.
  • [ ] Ma põhinesin multikollineaarsuse muutujate elimineerimisel teoorial, mitte ainult VIF-il.
  • [ ] Kasutasin standardviga, mis on vastupidav heteroskedastilisusele/autokorrelatsioonile (HC/HAC).
  • [ ] Pärast parandamist kontrollisin ja kontrollisin rikkumise mõju minu järeldusele.
  • [ ] Teatan, kuidas standardvea parandus mõjutas minu olulisuse tulemust.