Enhet 5 / 11

Regresjonsdiagnostikk og brudd: Kollinearitet, Heteroskedastisitet, Autokorrelasjon

Gevinster:

  • Evne til å oppdage brudd som multikollinearitet, heteroskedastisitet og autokorrelasjon med diagnostiske tester og grafikk med støtte for kunstig intelligens
  • Evne til å skille om hvert brudd forvrenger koeffisientestimater eller standardfeil og velge riktig korreksjon (robust standardfeil, transformasjon, modellrevisjon)
  • Evne til å verifisere at rettelsene foreslått av AI løser problemet i stedet for å skjule det, og at diagnosen er basert på en forståelse av prosessen som produserer dataene

Det er enkelt å lese resultatet av en regresjon; Det er vanskelig å stole på ham. Fordi koeffisientene er objektive, er standardfeilene korrekte, og p-verdiene er gyldige, avhengig av at forutsetningene vi introduserte i forrige enhet er oppfylt. I denne enheten vil vi dekke de tre vanligste bruddene: multikollinearitet, heteroskedastisitet og autokorrelasjon. Vi vil svare på tre spørsmål for hver: hvordan diagnostiseres det, hva det bryter (koeffisient eller standardfeil), og hvordan fikser det. Kunstig intelligens (AI) genererer koden for diagnostisk testing og korreksjon; men du bestemmer hvilket brudd som egentlig er problemet og om løsningen løser problemet.

La oss gjøre det mest kritiske skillet på forhånd: Forvrenger et brudd koeffisientestimatet eller bare standardfeilene? Denne forskjellen bestemmer løsningen. Et problem som forvrenger koeffisienten (f.eks. eksogenitetsbrudd) krever å tenke nytt om modellen; Et problem som bare forvrenger standardfeilen (heteroscedastisitet, autokorrelasjon) løses ofte ved standard feilretting. En vanlig feil ved AI er å bruke en teknikk som korrigerer standardfeilen og deretter erklære problemet "løst"; mens den underliggende strukturen fortsatt kan være der.

Brudd 1: Multikollinearitet

Hvorfor? To eller flere forklaringsvariabler er sterkt relatert til hverandre; for eksempel å legge både «totalt års erfaring» og «alder» inn i modellen. Siden disse variablene har nesten samme informasjon, har modellen vanskeligheter med å skille dem.

Hva bryter det? Koeffisientestimater forblir objektive, men blir ustabile: standardfeil blir oppblåst, konfidensintervaller utvides, individuelle koeffisienter kan vise seg å være ubetydelige mens modellen kan være signifikant totalt sett (F-test). En liten dataendring flytter koeffisienten betydelig.

Hvordan blir det diagnostisert? VIF (Variance Inflation Factor): for hver variabel måler den hvor mye den variabelen forklares av de andre. Som en grov terskel tilsier VIF > 5 (omtrent 10) forsiktighet. I tillegg undersøkes korrelasjonsmatrisen mellom forklaringsvariabler.

Hvordan fikse? Slippe en av de korrelerte variablene, kombinere dem (lage en indeks), eller (hvis teorien krever) beholde begge og unngå å tolke individuelle koeffisienter. Hvilken variabel som forblir er en teoretisk beslutning, ikke en statistisk – AI foreslår eliminering, du gir begrunnelsen.

Brudd 2: Heteroskedastisitet

Hvorfor? Variansen til feilleddet er ikke konstant på tvers av observasjoner. Typisk eksempel: ettersom inntekten øker, øker også spredningen rundt utgifter; Det dannes et "trakt"-mønster, smalt ved lav inntekt og bredt ved høy inntekt.

Hva bryter det? Koeffisientestimater er igjen objektive - dette er viktig. Det som blir forvrengt er standardfeilene: de beregnes feil, så p-verdier og konfidensintervaller blir upålitelige. Så koeffisienten din er i riktig senter, men svaret på spørsmålet "hvor sikker er du" er feil.

Hvordan er det diagnostisert? Spredningsplott av residualer kontra predikerte verdier (søkt etter traktmønster) og formelle tester: Breusch-hedensk test, hvit test. Liten p-verdi sier "ingen konstant varians".

Hvordan fikse? Den vanligste og mest praktiske løsningen er å bruke robuste standardfeil (robust / heteroskedastisitetskonsistent, HC standardfeil): det endrer ikke koeffisienten, det korrigerer standardfeilen for brudd. Alternativ: transformering av den avhengige variabelen (som log) eller vektet LCM. Robuste standardfeil har blitt nesten standard i empiri i dag.

Overtredelse 3: Autokorrelasjon

Hvorfor? Feilbegrepene er ikke uavhengige av hverandre; Det er mest vanlig i tidsserier: feilen i en periode påvirker den neste (f.eks. forblir resten positiv i perioder etter et sjokk).

Hva bryter det? Igjen, det forvrenger først og fremst standardfeil (undervurderer dem ofte, noe som skaper falsk betydning); Koeffisientene forblir objektive i LCC, men mister effektiviteten.

Hvordan er det diagnostisert? Durbin-Watson-test (for første-ordens autokorrelasjon), Breusch-Godfrey-test (mer generell), og plott av residualer versus etterskrevne verdier.

Hvordan fikse? Newey-West (HAC — motstandsdyktig mot autokorrelasjon og heteroskedastisitet) standardfeil, legger til etterslepende termer til modellen, eller bytter til riktig tidsseriemodell (enhet 8).

Forsiktig: Heteroskedastisitet og autokorrelasjon forvrenger standardfeilen, ikke koeffisienten. Derfor, før du sier "koeffisienten var signifikant", sørg for at standardfeilen er beregnet riktig; Ellers kan meningsfullhet være en illusjon.

sammenligningsdiagram

brudd

Hva går i stykker

Diagnose

Felles løsning

multi-link

Blåser opp standardfeil, gjør koeffisienten ustabil

VIF, korrelasjonsmatrise

Trekk fra/kombiner variabel (etter teori)

heteroskedastisitet

Standardfeil (koeffisient objektiv)

Breusch-hedensk, hvit, gjenværende plot

Robust (HC) standardfeil, konvertering

Autokorrelasjon

Standardfeil (koeffisient objektiv)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), forsinket termin

Fire kopierbare spørsmål

1. Komplett diagnosepakke:

Din rolle: regresjonsdiagnoseassistent. Jeg vil ha R-kode, jeg deler ikke dataene. Modell: lm(utgifter ~ inntekt + alder + region, data = data). Oppgave: (1) Beregn VIF, (2) Test heteroskedastisitet med Breusch-Pagan og residual-estimatplott, (3) Test autokorrelasjon med Durbin-Watson. Forklar i kommentarfeltet hvilket brudd hver test måler og hvordan p-verdien skal tolkes. Rettelse SØKNAD; stille diagnose.

2. Robust standardfeil:

Skriv R-kode som reproduserer koeffisienttabellen med heteroskedasticity-robust (HC3) standardfeil for samme modell (sandwich + lmtest-pakker). Lag en tabell som viser de klassiske og robuste standardfeilene side om side, slik at jeg kan se forskjellen. Understrek i kommentarfeltet at koeffisientene IKKE endres, kun standardfeilene gjør det.

3. Gjennomgang av flere lenker:

Gi meg en TENK-liste over variabler med høy VIF: hvilke variabler kan teoretisk måle det samme, hvilke som har en sterk sak for å beholde dem. IKKE autokvalifiser; Jeg tar avgjørelsen basert på teori. Forklar også hvorfor det kan være feil å bare se på VIF og tilordne variabler.

4. Autokorrelasjonskorreksjon:

I min tidsserieregresjon var Breusch-Godfrey p-verdien 0,001. Skriv R-kode som produserer koeffisienttabellen med Newey-West (HAC) standardfeil og forklar hvordan du velger etterslep. Merk at denne korreksjonen ikke løser ÅRSAKEN til autokorrelasjon, den korrigerer bare slutningen.

Svak forespørsel / Sterk forespørsel

Svak melding:

Er det et problem med regresjonen min? I så fall, fiks det.

Denne ledeteksten setter AI i en blind "fix"-modus: den kan hoppe over tester og bruke en direkte transformasjon eller variabel eliminering, uten å vise deg hvilket brudd som faktisk eksisterer og hva det brøt.

Kraftig ledetekst:

Din rolle: diagnoseassistent, ikke autokorrigerer. Første test for TRE brudd separat (multikolinearitet, heteroskedastisitet, autokorrelasjon) og for hver: hvilken test, hvordan man leser resultatet, bryter den KOEFFISK eller STANDARDFEIL. Foreslå deretter bare en løsning for bruddet som FAKTISK ble oppdaget, og forklar hvordan jeg kan bekrefte at løsningen løste problemet.

Forskjell: sterk vilje tvinger fram diagnose før korrigering og krever tydelig skille mellom "hva det går i stykker".

tre minisaker

Tilfelle 1 — falsk betydning (heteroscedastisitet). En analytiker fant at inntektskoeffisienten i forbruk~inntektsmodellen var "signifikant" ved p=0,01. Men nå var det et tydelig traktmønster på diagrammet hans. Når robuste standardfeil ble brukt, økte p-verdien til 0,09; meningsfullheten går tapt. Det første resultatet var en illusjon skapt av en feilberegnet standardfeil. Leksjon: Man kan ikke stole på betydning uten å korrigere standardfeilen.

Case 2 — Blind variabel eliminering. En elev droppet "erfaringsvariabelen" fra modellen fordi hans VIF var høy; Koeffisientene ble "ryddet opp", men den teoretiske betydningen av modellen ble forvrengt fordi erfaring var hovedvariabelen av interesse. Den riktige måten: trekk fra alder og behold erfaring, eller kombiner de to. Leksjon: VIF-terskelen alene er ikke grunnlag for eliminering; teori bestemmer.

Tilfelle 3 — Usikkerhet skjult av autokorrelasjon. I en tidsseriestudie ble Durbin-Watson ignorert; koeffisienten virket veldig "nøyaktig" (smalt konfidensintervall). Da Newey-West standardfeil ble brukt, doblet konfidensintervallet seg og policyanbefalingen ble mye mer konservativ. Leksjon: autokorrelasjon kan undervurdere usikkerhet.

Vanlige feil

  • Å ta feil av bruddet som forvrenger standardfeilen som et koeffisientproblem. Heteroskedastisitet og autokorrelasjon forvrenger ikke koeffisienten; Ikke få panikk og bygg om modellen unødvendig, rett standardfeilen først.
  • Ser på VIF og blindt tilordner variabler. Å fjerne en teoretisk nødvendig variabel bare fordi VIF er høy, gjør modellen meningsløs.
  • Verifiserer ikke rettelsen. Etter å ha brukt robust standardfeil, kontroller at bruddet faktisk korrigerer slutningen; Ikke si "Jeg brukte det, det er gjort".
  • Fiks uten diagnose. Å bruke transformasjon/eliminering uten å teste hvilken overtredelse som eksisterer kan "løse" problemet der det ikke eksisterer og skjule den som gjør det.
  • Setter fiksen på plass hvorfor. Newey-West løser ikke årsaken til autokorrelasjon; det fikser bare konklusjonen. Hvis det er et strukturelt problem, må modellen endres.
Tips: For hvert brudd, spør deg selv "ødelegger dette koeffisienten eller min tillit til den?" Hvis svaret er "tillit", er løsningen nesten alltid standard feilretting, ikke gjenoppbygging av modellen.

Oppsummert

Regresjonsdiagnostikk er en forutsetning for å stole på utgangen. Av de tre vanlige bruddene gjør multikollinearitet koeffisienten ustabil og blåser opp standardfeilen; Heteroskedastisitet og autokorrelasjon, derimot, etterlater koeffisienten objektiv og forvrenger bare standardfeilen. AI genererer diagnose- og fikskoden, men du bestemmer hvilket brudd som er det virkelige problemet, hvilken variant som forblir teoretisk, og om løsningen løser problemet. Verken panikk eller blind korreksjon er riktig uten å tydeliggjøre skillet mellom "hva som går i stykker".

Søknadsoppgave

Sett opp en multivariat regresjon og be AI om en kode som bare produserer diagnoser (ingen korreksjoner): VIF, Breusch-Pagan/White og Durbin-Watson/Breusch-Godfrey. For hver test skal du tolke resultatet og angi om bruddet forvrenger koeffisienten eller standardfeilen. For et reelt brudd du oppdager (f.eks. heteroskedastisitet), bruk robuste standardfeil og sett klassiske og robuste resultater sammen; Vis at koeffisienten ikke endres, men standardfeilen endres. Rapporter i ett avsnitt hvordan signifikansresultatet ditt ble påvirket av korrigeringen.

sjekkliste

  • [ ] Jeg testet tre brudd (multikolinearitet, heteroskedastisitet, autokorrelasjon) hver for seg.
  • [ ] For hvert brudd differensierte jeg om det forvrenger koeffisienten eller standardfeilen.
  • [ ] Jeg baserte variabel eliminering i multikollinearitet på teori, ikke bare VIF.
  • [ ] Jeg brukte standardfeil med robusthet til heteroskedastisitet/autokorrelasjon (HC/HAC).
  • [ ] Etter korrigering sjekket og bekreftet jeg virkningen av bruddet på slutningen min.
  • [ ] Jeg rapporterer hvordan mitt signifikansresultat ble påvirket av standard feilretting.