Enhet 5 / 11

Regressionsdiagnostik och kränkningar: kollinearitet, heteroskedasticitet, autokorrelation

Vinster:

  • Förmåga att upptäcka överträdelser som multikollinearitet, heteroskedasticitet och autokorrelation med diagnostiska tester och grafik med stöd för artificiell intelligens
  • Förmåga att särskilja om varje överträdelse förvränger koefficientuppskattningar eller standardfel och välja lämplig korrigering (robust standardfel, transformation, modellrevision)
  • Förmåga att verifiera att de korrigeringar som föreslås av AI löser problemet snarare än att dölja det och att diagnosen är baserad på en förståelse av processen som producerar data

Att läsa utdata från en regression är lätt; Det är svårt att lita på honom. Eftersom koefficienterna är opartiska, är standardfelen korrekta och p-värdena är giltiga, beroende på att de antaganden som vi introducerade i föregående enhet är uppfyllda. I denna enhet kommer vi att täcka de tre vanligaste kränkningarna: multikollinearitet, heteroskedasticitet och autokorrelation. Vi kommer att svara på tre frågor för varje: hur man diagnostiserar det, vad det går sönder (koefficient eller standardfel) och hur man fixar det. Artificiell intelligens (AI) genererar koden för diagnostisk testning och korrigering; men du bestämmer vilken överträdelse som verkligen är problemet och om åtgärden löser problemet.

Låt oss göra den mest kritiska skillnaden i förväg: Förvränger en överträdelse koefficientuppskattningen eller bara standardfelen? Denna distinktion avgör lösningen. Ett problem som förvränger koefficienten (t.ex. exogenitetsöverträdelse) kräver att modellen omprövas; Ett problem som bara förvränger standardfelet (heteroscedasticitet, autokorrelation) löses ofta genom standardfelkorrigering. Ett vanligt misstag för AI är att tillämpa en teknik som korrigerar standardfelet och sedan förklara problemet "löst"; medan den underliggande strukturen fortfarande kan finnas där.

Överträdelse 1: Multikollinearitet

Varför? Två eller flera förklaringsvariabler är starkt relaterade till varandra; till exempel att lägga in både "total års erfarenhet" och "ålder" i modellen. Eftersom dessa variabler bär nästan samma information har modellen svårt att skilja dem åt.

Vad går det sönder? Koefficientuppskattningar förblir opartiska men blir instabila: standardfel blir uppblåsta, konfidensintervall vidgas, individuella koefficienter kan visa sig vara insignifikanta medan modellen kan vara signifikant totalt sett (F-test). En liten dataändring flyttar koefficienten avsevärt.

Hur diagnostiseras det? VIF (Variance Inflation Factor): för varje variabel mäter den hur mycket den variabeln förklaras av de andra. Som en grov tröskel, VIF > 5 (ungefär 10) motiverar försiktighet. Dessutom undersöks korrelationsmatrisen mellan förklaringsvariabler.

Hur fixar man? Släpp en av de korrelerade variablerna, kombinera dem (gör ett index) eller (om teorin kräver) behålla båda och undvika att tolka individuella koefficienter. Vilken variabel som stannar är ett teoretiskt beslut, inte ett statistiskt – AI föreslår eliminering, du ger motiveringen.

Överträdelse 2: Heteroskedasticitet

Varför? Variansen av feltermen är inte konstant över observationer. Typiskt exempel: när inkomsten ökar, ökar också spridningen kring utgifterna; Ett "tratt"-mönster bildas, smalt vid låg inkomst och brett vid hög inkomst.

Vad går det sönder? Koefficientuppskattningar är återigen opartiska — detta är viktigt. Det som blir förvrängt är standardfelen: de beräknas felaktigt, så p-värden och konfidensintervall blir opålitliga. Så din koefficient är i rätt mitt, men svaret på frågan "hur säker är du" är fel.

Hur diagnostiseras det? Spridningsdiagram av residualer kontra förutsagda värden (sökt efter trattmönster) och formella tester: Breusch-Pagan-test, Vitt test. Litet p-värde säger "ingen konstant varians".

Hur fixar man? Den vanligaste och mest praktiska lösningen är att använda robusta standardfel (robust / heteroskedasticitetskonsistent, HC-standardfel): det ändrar inte koefficienten, det korrigerar standardfelet för överträdelse. Alternativ: transformering av den beroende variabeln (som log) eller viktad LCM. Robusta standardfel har blivit nästan standard i empiriskt arbete idag.

Överträdelse 3: Autokorrelation

Varför? Feltermerna är inte oberoende av varandra; Det är vanligast i tidsserier: felet i en period påverkar nästa (t.ex. återstoden förblir positiv i perioder efter en chock).

Vad går det sönder? Återigen, det förvränger i första hand standardfel (underskattar dem ofta, vilket skapar falsk betydelse); Koefficienterna förblir opartiska i LCC men förlorar sin effektivitet.

Hur diagnostiseras det? Durbin-Watson-test (för första ordningens autokorrelation), Breusch-Godfrey-test (mer allmänt) och plotter av residualer kontra fördröjda värden.

Hur fixar man? Newey-West (HAC — resistent mot autokorrelation och heteroskedasticitet) standardfel, lägga till fördröjda termer till modellen eller byta till lämplig tidsseriemodell (enhet 8).

Varning: Heteroskedasticitet och autokorrelation förvränger standardfelet, inte koefficienten. Innan du säger "koefficienten var signifikant", se därför till att standardfelet beräknas korrekt; Annars kan meningsfullhet vara en illusion.

jämförelsediagram

kränkning

Vad går sönder

Diagnos

Vanlig lösning

multilänk

Blåser upp standardfel, gör koefficienten instabil

VIF, korrelationsmatris

Subtrahera/kombinera variabel (genom teori)

heteroskedasticitet

Standardfel (koefficient opartisk)

Breusch-hednisk, vit, resttomt

Robust (HC) standardfel, konvertering

Autokorrelation

Standardfel (koefficient opartisk)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), försenad termin

Fyra kopierbara uppmaningar

1. Komplett diagnostikpaket:

Din roll: regressionsdiagnostikassistent. Jag vill ha R-kod, jag delar inte data. Modell: lm(utgifter ~ inkomst + ålder + region, data = data). Uppgift: (1) Beräkna VIF, (2) Testa heteroskedasticitet med Breusch-Pagan och residual-estimat plot, (3) Testa autokorrelation med Durbin-Watson. Förklara i kommentarsraden vilken överträdelse varje test mäter och hur man tolkar p-värdet. Rättelse ANSÖKAN; ställa diagnos.

2. Robust standardfel:

Skriv R-kod som återger koefficienttabellen med heteroskedasticity-robust (HC3) standardfel för samma modell (sandwich + lmtest-paket). Ta fram en tabell som visar de klassiska och robusta standardfelen sida vid sida så att jag kan se skillnaden. Betona i kommentarsraden att koefficienterna INTE ändras, det gör bara standardfelen.

3. Granskning av flera länkar:

Ge mig en TÄNK-lista över variabler med hög VIF: vilka variabler skulle teoretiskt kunna mäta samma sak, vilka som har ett starkt argument för att behålla dem. Kvalificera INTE automatiskt; Jag kommer att fatta beslutet baserat på teori. Förklara också varför det kan vara fel att bara titta på VIF och tilldela variabler.

4. Autokorrelationskorrigering:

I min tidsserieregression var Breusch-Godfrey p-värdet 0,001. Skriv R-kod som producerar koefficienttabellen med Newey-West (HAC) standardfel och förklara hur man väljer eftersläpning. Observera att denna korrigering inte löser ORSAKEN till autokorrelation, den korrigerar bara slutsatsen.

Svag prompt / Stark prompt

Svag uppmaning:

Finns det något problem med min regression? I så fall, fixa det.

Denna prompt sätter AI i ett blindt "fix"-läge: den kan hoppa över tester och tillämpa en direkt transformation eller variabel eliminering, utan att visa dig vilken överträdelse som faktiskt finns och vad den gick sönder.

Kraftfull uppmaning:

Din roll: diagnostisk assistent, inte autokorrigerare. Första testet för TRE överträdelser separat (multikolinjäritet, heteroskedasticitet, autokorrelation) och för varje: vilket test, hur man läser resultatet, bryter det COEFFICIENT eller STANDARD ERROR. Föreslå sedan bara en åtgärd för överträdelsen som FAKTISKT upptäcktes och förklara hur jag kan verifiera att korrigeringen löste problemet.

Skillnad: stark vilja tvingar fram diagnos före korrigering och kräver tydlig åtskillnad mellan "vad det går sönder".

tre minifodral

Fall 1 — Falsk betydelse (heteroscedasticitet). En analytiker fann att intäktskoefficienten i utgifts~inkomstmodellen var "signifikant" vid p=0,01. Men nu fanns det ett distinkt trattmönster på hans sjökort. När robusta standardfel tillämpades ökade p-värdet till 0,09; meningsfullheten går förlorad. Det första resultatet var en illusion skapad av ett felräknat standardfel. Lärdom: signifikans kan inte litas på utan att korrigera standardfelet.

Fall 2 — Blind variabel eliminering. En elev tappade variabeln "erfarenhet" från modellen eftersom hans VIF var hög; Koefficienterna ”rensades upp” men modellens teoretiska betydelse förvrängdes eftersom erfarenhet var den huvudsakliga variabeln av intresse. Rätt sätt: subtrahera ålder och behålla erfarenhet, eller kombinera de två. Lektion: VIF-tröskeln ensam är inte skäl för eliminering; teorin avgör.

Fall 3 — Osäkerhet dold av autokorrelation. I en tidsseriestudie ignorerades Durbin-Watson; koefficienten verkade mycket "precis" (snävt konfidensintervall). När Newey-West standardfel tillämpades fördubblades konfidensintervallet och policyrekommendationen blev mycket mer konservativ. Lärdom: autokorrelation kan underskatta osäkerhet.

Vanliga misstag

  • Att missförstå överträdelsen som förvränger standardfelet som ett koefficientproblem. Heteroskedasticitet och autokorrelation förvränger inte koefficienten; Få inte panik och bygg om modellen i onödan, korrigera standardfelet först.
  • Att titta på VIF och blint tilldela variabler. Att ta bort en teoretiskt nödvändig variabel bara för att VIF är hög gör modellen meningslös.
  • Verifierar inte korrigeringen. Efter att ha tillämpat ett robust standardfel, kontrollera att överträdelsen faktiskt korrigerar slutsatsen; Säg inte "Jag tillämpade det, det är klart".
  • Fixa utan diagnos. Att tillämpa transformation/eliminering utan att testa vilken överträdelse som finns kan "lösa" problemet där det inte finns och dölja den som gör det.
  • Att fixa varför. Newey-West löser inte orsaken till autokorrelation; det fixar bara slutsatsen. Om det finns ett strukturellt problem måste modellen ändras.
Tips: För varje överträdelse, fråga dig själv "förstör detta koefficienten eller mitt förtroende för den?" Om svaret är "förtroende" är lösningen nästan alltid standardfelkorrigering, inte att bygga om modellen.

Sammanfattningsvis

Regressionsdiagnostik är en förutsättning för att lita på resultatet. Av de tre vanliga överträdelserna gör multikollinearitet koefficienten instabil och blåser upp standardfelet; Heteroskedasticitet och autokorrelation, å andra sidan, lämnar koefficienten opartisk och förvränger bara standardfelet. AI:n genererar diagnostik- och fixkoden, men du bestämmer vilken överträdelse som är det verkliga problemet, vilken variant som förblir teoretisk och om fixen löser problemet. Varken panik eller blind korrigering är korrekt utan att förtydliga skillnaden mellan "vad som går sönder".

Applikationsuppgift

Ställ in en multivariat regression och be AI om en kod som bara ger diagnoser (inga korrigeringar): VIF, Breusch-Pagan/White och Durbin-Watson/Breusch-Godfrey. För varje test, tolka resultatet och ange om överträdelsen förvränger koefficienten eller standardfelet. För en verklig överträdelse du upptäcker (t.ex. heteroskedasticitet), tillämpa robusta standardfel och jämför klassiska och robusta resultat; Visa att koefficienten inte ändras men standardfelet ändras. Rapportera i ett stycke hur ditt signifikansresultat påverkades av korrigeringen.

checklista

  • [ ] Jag testade tre överträdelser (multikolinjäritet, heteroskedasticitet, autokorrelation) separat.
  • [ ] För varje överträdelse särskiljde jag om det förvränger koefficienten eller standardfelet.
  • [ ] Jag baserade variabeleliminering i multikollinearitet på teori, inte bara VIF.
  • [ ] Jag använde standardfel med robusthet mot heteroskedasticitet/autokorrelation (HC/HAC).
  • [ ] Efter korrigering kontrollerade och verifierade jag överträdelsens inverkan på min slutsats.
  • [ ] Jag rapporterar hur mitt signifikansresultat påverkades av standardfelkorrigering.