Eenheid 5 / 11

Regressiediagnostiek en overtredingen: collineariteit, heteroscedasticiteit, autocorrelatie

Winst:

  • Mogelijkheid om overtredingen zoals multicollineariteit, heteroscedasticiteit en autocorrelatie te detecteren met diagnostische tests en grafische afbeeldingen met ondersteuning voor kunstmatige intelligentie
  • Mogelijkheid om te onderscheiden of elke overtreding coëfficiëntschattingen of standaardfouten vertekent en de juiste correctie te kiezen (robuuste standaardfout, transformatie, modelrevisie)
  • Mogelijkheid om te verifiëren dat de door de AI voorgestelde oplossingen het probleem oplossen in plaats van verbergen, en dat de diagnose gebaseerd is op inzicht in het proces dat de gegevens produceert

Het lezen van de uitkomst van een regressie is eenvoudig; Het is moeilijk om hem te vertrouwen. Omdat de coëfficiënten onbevooroordeeld zijn, zijn de standaardfouten correct en zijn de p-waarden geldig, afhankelijk van de aannames die we in de vorige unit hebben geïntroduceerd. In dit onderdeel behandelen we de drie meest voorkomende schendingen: multicollineariteit, heteroscedasticiteit en autocorrelatie. We zullen voor elk ervan drie vragen beantwoorden: hoe we het probleem kunnen diagnosticeren, wat het kapot maakt (coëfficiënt of standaardfout) en hoe we het kunnen repareren. Kunstmatige intelligentie (AI) genereert de code voor diagnostische tests en correctie; maar u beslist welke overtreding werkelijk het probleem is en of de oplossing het probleem oplost.

Laten we vooraf het meest kritische onderscheid maken: vertekent een overtreding de schatting van de coëfficiënten of alleen de standaardfouten? Dit onderscheid bepaalt de oplossing. Een probleem dat de coëfficiënt vertekent (bijvoorbeeld een schending van de exogeniteit) vereist een heroverweging van het model; Een probleem dat alleen de standaardfout vertekent (heteroscedasticiteit, autocorrelatie) wordt vaak opgelost door standaardfoutcorrectie. Een veel voorkomende fout bij AI is het toepassen van een techniek die de standaardfout corrigeert en vervolgens het probleem ‘opgelost’ verklaart; terwijl de onderliggende structuur er misschien nog steeds is.

Overtreding 1: Multicollineariteit

Waarom? Twee of meer verklarende variabelen zijn sterk met elkaar gerelateerd; door bijvoorbeeld zowel ‘totaal aantal ervaringsjaren’ als ‘leeftijd’ in het model te plaatsen. Omdat deze variabelen vrijwel dezelfde informatie bevatten, kan het model ze moeilijk scheiden.

Wat breekt het? Coëfficiëntschattingen blijven onbevooroordeeld, maar worden instabiel: standaardfouten worden groter, betrouwbaarheidsintervallen worden groter, individuele coëfficiënten kunnen onbeduidend blijken te zijn, terwijl het model over het geheel genomen significant kan zijn (F-test). Een kleine gegevensverandering verplaatst de coëfficiënt aanzienlijk.

Hoe wordt de diagnose gesteld? VIF (Variance Inflation Factor): voor elke variabele meet het hoeveel die variabele door de anderen wordt verklaard. Als ruwe drempel verdient VIF > 5 (ongeveer 10) voorzichtigheid. Daarnaast wordt de correlatiematrix tussen verklarende variabelen onderzocht.

Hoe op te lossen? Eén van de gecorreleerde variabelen laten vallen, ze combineren (een index maken), of (als de theorie vereist) beide behouden en de interpretatie van individuele coëfficiënten vermijden. Welke variabele blijft, is een theoretische beslissing, geen statistische beslissing. AI suggereert eliminatie, jij geeft de rechtvaardiging.

Overtreding 2: Heteroscedasticiteit

Waarom? De variantie van de foutterm is niet constant voor alle waarnemingen. Typisch voorbeeld: naarmate het inkomen toeneemt, neemt ook de spreiding rond de uitgaven toe; Er ontstaat een ‘trechterpatroon’: smal bij lage inkomens en breed bij hoge inkomens.

Wat breekt het? Coëfficiëntschattingen zijn wederom onbevooroordeeld – dit is belangrijk. Wat vervormd raakt zijn de standaardfouten: deze worden verkeerd berekend, waardoor p-waarden en betrouwbaarheidsintervallen onbetrouwbaar worden. Je coëfficiënt staat dus in het juiste midden, maar het antwoord op de vraag "hoe zeker ben je" is verkeerd.

Hoe wordt de diagnose gesteld? Spreidingsdiagram van residuen versus voorspelde waarden (gezocht trechterpatroon) en formele tests: Breusch-Pagan-test, White-test. Een kleine p-waarde zegt "geen constante variantie".

Hoe op te lossen? De meest voorkomende en praktische oplossing is het gebruik van robuuste standaardfouten (robuuste / heteroskedasticiteit-consistent, HC-standaardfouten): het verandert de coëfficiënt niet, het corrigeert de standaardfout voor overtreding. Alternatief: het transformeren van de afhankelijke variabele (zoals log) of gewogen LCM. Robuuste standaardfouten zijn tegenwoordig bijna de standaard geworden in empirisch werk.

Overtreding 3: Autocorrelatie

Waarom? De fouttermen zijn niet onafhankelijk van elkaar; Het komt het meest voor in tijdreeksen: de fout van de ene periode beïnvloedt de volgende (het residu blijft bijvoorbeeld positief in perioden na een schok).

Wat breekt het? Nogmaals, het verdraait voornamelijk standaardfouten (vaak onderschat ze, waardoor een valse betekenis ontstaat); De coëfficiënten blijven onbevooroordeeld in de LCC, maar verliezen hun effectiviteit.

Hoe wordt de diagnose gesteld? Durbin-Watson-test (voor autocorrelatie van de eerste orde), Breusch-Godfrey-test (meer algemeen) en plots van residuen versus vertraagde waarden.

Hoe op te lossen? Newey-West (HAC – resistent tegen autocorrelatie en heteroskedasticiteit) standaardfouten, het toevoegen van vertraagde termen aan het model, of het overschakelen naar het juiste tijdreeksmodel (eenheid 8).

Let op: heteroscedasticiteit en autocorrelatie vervormen de standaardfout, niet de coëfficiënt. Voordat u zegt: "de coëfficiënt was significant", moet u er daarom voor zorgen dat de standaardfout correct is berekend; Anders kan zinvolheid een illusie zijn.

vergelijkingstabel

overtreding

Wat breekt

Diagnose

Gemeenschappelijke oplossing

multi-link

Verhoogt standaardfouten, maakt de coëfficiënt onstabiel

VIF, correlatiematrix

Variabele aftrekken/combineren (volgens theorie)

heteroscedasticiteit

Standaardfouten (coëfficiënt onbevooroordeeld)

Breusch-Pagan, White, resterend plot

Robuuste (HC) standaardfout, conversie

Autocorrelatie

Standaardfouten (coëfficiënt onbevooroordeeld)

Durbin-Watson, Breusch-Godfrey

Newey-West (HAC), uitgestelde termijn

Vier kopieerbare aanwijzingen

1. Compleet diagnosepakket:

Jouw rol: regressiediagnostisch assistent. Ik wil R-code, ik deel de gegevens niet. Model: lm(uitgaven ~ inkomen + leeftijd + regio, data = data). Taak: (1) Bereken VIF, (2) Test heteroskedasticiteit met Breusch-Pagan en residuele schattingsplot, (3) Test autocorrelatie met Durbin-Watson. Leg in de commentaarregel uit welke overtreding elke test meet en hoe de p-waarde moet worden geïnterpreteerd. Correctie TOEPASSING; diagnose produceren.

2. Robuuste standaardfout:

Schrijf R-code die de coëfficiëntentabel reproduceert met heteroskedasticiteit-robuuste (HC3) standaardfouten voor hetzelfde model (sandwich + lmtest-pakketten). Maak een tabel waarin de klassieke en robuuste standaardfouten naast elkaar worden weergegeven, zodat ik het verschil kan zien. Benadruk in de commentaarregel dat de coëfficiënten NIET veranderen, alleen de standaardfouten.

3. Multi-link-beoordeling:

Geef me een THINK-lijst met variabelen met een hoge VIF: welke variabelen theoretisch hetzelfde zouden kunnen meten, en welke variabelen er goede argumenten voor hebben om ze te behouden. NIET automatisch kwalificeren; Ik zal de beslissing nemen op basis van de theorie. Leg ook uit waarom het verkeerd zou kunnen zijn om alleen naar VIF te kijken en variabelen toe te wijzen.

4. Autocorrelatiecorrectie:

In mijn tijdreeksregressie was de Breusch-Godfrey p-waarde 0,001. Schrijf R-code die de coëfficiëntentabel produceert met standaardfouten van Newey-West (HAC) en leg uit hoe u de vertraging selecteert. Merk op dat deze correctie de OORZAAK van autocorrelatie niet oplost, maar alleen de gevolgtrekking corrigeert.

Zwakke prompt/sterke prompt

Zwakke prompt:

Is er een probleem met mijn regressie? Zo ja, repareer het dan.

Deze prompt zet de AI in een blinde ‘fix’-modus: hij kan tests overslaan en een directe transformatie of eliminatie van variabelen toepassen, zonder je te laten zien welke overtreding daadwerkelijk bestaat en wat deze heeft verbroken.

Krachtige prompt:

Jouw rol: diagnostisch assistent, geen autocorrector. Test eerst op DRIE overtredingen afzonderlijk (multicollineariteit, heteroscedasticiteit, autocorrelatie) en voor elke: welke test, hoe het resultaat te lezen, doorbreekt deze COEFFICIENT of STANDARD ERROR. Stel dan alleen een oplossing voor voor de overtreding die WERKELIJK is gedetecteerd en leg uit hoe ik kan verifiëren dat de oplossing het probleem heeft opgelost.

Verschil: een sterke wil dwingt een diagnose af vóór correctie en vereist een duidelijk onderscheid tussen "wat het kapot maakt".

drie minikoffers

Geval 1 — Onechte betekenis (heteroscedasticiteit). Eén analist vond de omzetcoëfficiënt in het bestedingsinkomstenmodel “significant” met p=0,01. Maar nu was er een duidelijk trechterpatroon op zijn kaart te zien. Wanneer robuuste standaardfouten werden toegepast, nam de p-waarde toe tot 0,09; zingeving gaat verloren. Het eerste resultaat was een illusie gecreëerd door een verkeerd berekende standaardfout. Les: significantie kan niet worden vertrouwd zonder de standaardfout te corrigeren.

Geval 2 — Blinde eliminatie van variabelen. Eén student liet de variabele 'ervaring' uit het model vallen omdat zijn VIF hoog was; De coëfficiënten werden “opgeschoond”, maar de theoretische betekenis van het model werd vertekend omdat ervaring de belangrijkste variabele van belang was. De juiste manier: trek de leeftijd af en behoud de ervaring, of combineer de twee. Les: de VIF-drempel alleen is geen reden voor eliminatie; theorie bepaalt.

Geval 3 — Onzekerheid verborgen door autocorrelatie. In een tijdreeksstudie werd Durbin-Watson genegeerd; de coëfficiënt leek erg "precies" (nauw betrouwbaarheidsinterval). Toen Newey-West-standaardfouten werden toegepast, verdubbelde het betrouwbaarheidsinterval en werd de beleidsaanbeveling veel conservatiever. Les: autocorrelatie kan onzekerheid onderschatten.

Veel voorkomende fouten

  • De overtreding die de standaardfout vertekent, aanzien als een coëfficiëntenprobleem. Heteroscedasticiteit en autocorrelatie verstoren de coëfficiënt niet; Raak niet in paniek en bouw het model onnodig opnieuw op, corrigeer eerst de standaardfout.
  • Kijken naar VIF en blindelings variabelen toewijzen. Het verwijderen van een theoretisch noodzakelijke variabele alleen maar omdat de VIF hoog is, maakt het model zinloos.
  • De oplossing wordt niet geverifieerd. Controleer na het toepassen van een robuuste standaardfout of de overtreding de gevolgtrekking daadwerkelijk corrigeert; Zeg niet: "Ik heb het toegepast, het is klaar".
  • Fix zonder diagnose. Het toepassen van transformatie/eliminatie zonder te testen welke overtreding er bestaat, kan het probleem 'oplossen' waar het niet bestaat en het probleem dat dat wel doet, verbergen.
  • De oplossing aanbrengen waarom. Newey-West lost de oorzaak van autocorrelatie niet op; het corrigeert alleen de gevolgtrekking. Als er een structureel probleem is, moet het model worden gewijzigd.
Tip: Stel uzelf bij elke overtreding de vraag: 'Vernietigt dit de coëfficiënt of mijn vertrouwen daarin?' Als het antwoord ‘vertrouwen’ is, is de oplossing bijna altijd standaardfoutcorrectie, en niet het opnieuw opbouwen van het model.

Samengevat

Regressiediagnostiek is een voorwaarde voor het vertrouwen van de uitvoer. Van de drie veelvoorkomende overtredingen maakt multicollineariteit de coëfficiënt onstabiel en vergroot de standaardfout; Heteroskedasticiteit en autocorrelatie laten de coëfficiënt daarentegen onbevooroordeeld en vertekenen alleen de standaardfout. De AI genereert de diagnose- en fixcode, maar jij bepaalt welke overtreding het echte probleem is, welke variant theoretisch blijft en of de fix het probleem oplost. Noch paniek, noch blinde correctie zijn correct zonder het onderscheid tussen ‘wat kapot is’ te verduidelijken.

Applicatie taak

Zet een multivariate regressie op en vraag de AI om een ​​code die alleen diagnoses oplevert (geen correcties): VIF, Breusch-Pagan/White en Durbin-Watson/Breusch-Godfrey. Interpreteer voor elke test het resultaat en geef aan of de overtreding de coëfficiënt of de standaardfout vertekent. Voor een echte overtreding die u detecteert (bijvoorbeeld heteroscedasticiteit), past u robuuste standaardfouten toe en plaatst u klassieke en robuuste resultaten naast elkaar; Laat zien dat de coëfficiënt niet verandert, maar dat de standaardfout verandert. Vermeld in één alinea hoe uw significantieresultaat door de correctie is beïnvloed.

controlelijst

  • [ ] Ik heb drie overtredingen (multicollineariteit, heteroskedasticiteit, autocorrelatie) afzonderlijk getest.
  • [ ] Voor elke overtreding heb ik onderscheid gemaakt tussen de coëfficiënt of de standaardfout.
  • [ ] Ik heb de eliminatie van variabelen bij multicollineariteit gebaseerd op theorie, niet alleen op VIF.
  • [ ] Ik heb standaardfouten gebruikt met robuustheid voor heteroscedasticiteit/autocorrelatie (HC/HAC).
  • [ ] Na correctie heb ik de impact van de overtreding op mijn gevolgtrekking gecontroleerd en geverifieerd.
  • [ ] Ik rapporteer hoe mijn significantieresultaat werd beïnvloed door standaardfoutcorrectie.