| Von | → Nach | β | ×W (Moderator) | β_int |
|---|
| Variable | Verteilung | Mittelwert / p | SD / – |
|---|
Generiere die simulierten Daten, dann verändere die Modellformel und beobachte was passiert:
· Confounder weglassen (aus brms/glmmTMB-Code entfernen) → β(X) wird verzerrt (Confounding Bias)
· Collider aufnehmen (zur Formel hinzufügen) → öffnet Scheinassoziation zwischen X und Y (Berkson Bias)
· Mediator kontrollieren bei totalem Effekt → unterschätzt β(X), weil der indirekte Pfad blockiert wird
Der wahre β steht im Simulationscode — vergleiche ihn mit dem Schätzer im gefitteten Modell!
Eine Instrumentalvariable Z erlaubt kausale Identifikation auch dann, wenn ein ungemessener Confounder U existiert. Z muss drei Bedingungen erfüllen: (1) Relevanz: Z→X (Z beeinflusst die Exposition), (2) Exklusion: Z hat keinen direkten Effekt auf Y außer über X, (3) Unabhängigkeit: Z teilt keine gemeinsamen Ursachen mit Y (Z ist unabhängig von U).
Frequentistisch: Two-Stage Least Squares (2SLS) via ivreg(Y ~ X | Z). Bayesianisch: Joint-Modell beider Gleichungen mit korrelierten Residuen — rescor = TRUE schätzt ρ (Konfundierungsstärke durch U) direkt als Posterior.
Angrist & Pischke (2009). Mostly Harmless Econometrics. Princeton UP. — McElreath (2020). Statistical Rethinking, Kap. 15.
Wenn U X und Y konfundiert, aber ein Mediator M existiert mit (1) M auf allen Pfaden X→Y, (2) kein Backdoor X→M, (3) alle Backdoors M→Y durch X blockiert — dann ist der kausale Effekt ohne U zu messen identifizierbar.
Front Door Formel: P(Y|do(X)) = Σm P(M|X) · Σx' P(Y|M,X=x') · P(X=x'). Im linearen Fall: β_FD = β_XM × β_MY, wobei β_MY durch Kontrolle für X geschätzt wird (X blockt M←X←U→Y).
Bayesianisch: multivariates brms-Modell mit rescor = FALSE liefert die gemeinsame Posteriori von β_XM und β_MY — deren Produkt ist der kausale Gesamteffekt als volle Posterior-Verteilung.
Pearl (2009). Causality, Kap. 3.3. Cambridge UP. — Pearl, Glymour & Jewell (2016). Causal Inference in Statistics: A Primer. Wiley.