Pseudo-Bestimmtheitsmaß

Für allgemeine Regressionsmodelle deren Parameter durch Maximum-Likelihood-Schätzung gefunden werden, lassen sich verschiedene Pseudo-Bestimmtheitsmaße (notiert als ${\text{Pseudo-R}}^2$) definieren, die auf der Likelihoodfunktion basieren.

Motivation

Im Falle einer linearen Regression beschreibt das Bestimmtheitsmaß den erklärten Anteil der Variabilität (Varianz) einer abhängigen Variablen $Y$ durch ein statistisches Modell. Bei einem nominalen oder ordinalen Skalenniveau von $Y$ (z. B. bei Klassifikationen) existiert jedoch kein Äquivalent, da man die Varianz und damit ein ${\mathrm{R}}^{\mathrm{2}}$ nicht berechnen kann.

Das Pseudo-Bestimmtheitsmaß

Pseudo-Bestimmtheitsmaße sind so konstruiert, dass sie den verschiedenen Interpretationen (z. B. erklärte Varianz, Verbesserung gegenüber dem Nullmodell oder als Quadrat der Korrelation) des Bestimmtheitsmaßes genügen. Sie sind dem ${\mathrm{R}}^{\mathrm{2}}$ in der Hinsicht ähnlich, dass dessen Werte ebenfalls im Intervall von 0 und 1 liegen und ein höherer Wert einer besseren Anpassung des Modells an die Daten entspricht.

Likelihood-basierte Maße

Maddalas / Cox & Snells Pseudo-R²

$R_{\mathrm{M}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{a}\mathrm{l}\mathrm{a}}^{\mathrm{2}}=1-{\left(\frac{L_0}{L_1}\right)}^{2/n}$,

mit

$L_0$: Nullmodell,

$L_1$: Modell mit erklärenden Variablen

$R_{\mathrm{M}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{a}\mathrm{l}\mathrm{a}}^{\mathrm{2}}\in [0,1)$

Vergleicht das Verhältnis der Werte $L_0$ der Wert der Likelihood-Funktion, in dem die völlige Unabhängigkeit aller Variablen angenommen wird (Nullmodell oder Leermodell) und $L_1$ der Likelihood-Funktionen unter Kenntnis des Zusammenhanges zwischen $Y$ und $X_i$ (volles Regressionsmodell). Je geringer dieses Verhältnis, desto größer die Verbesserung des ganzen Modells gegenüber dem Nullmodell. Maddalas ${\text{Pseudo-R}}^2$ kann auch bei perfekter Vorhersage nie den Wert 1 erreichen.

Nagelkerkes / Cragg & Uhlers Pseudo-R²

$R_{\mathrm{N}\mathrm{a}\mathrm{g}\mathrm{e}\mathrm{l}\mathrm{k}\mathrm{e}\mathrm{r}\mathrm{k}\mathrm{e}}^{\mathrm{2}}=\frac{1-{\left(\frac{L_0}{L_1}\right)}^{2/n}}{1-L_0^{2/n}}$,

mit

$L_0$: Nullmodell,

$L_1$: Modell mit erklärenden Variablen

$R_{\mathrm{N}\mathrm{a}\mathrm{g}\mathrm{e}\mathrm{l}\mathrm{k}\mathrm{e}\mathrm{r}\mathrm{k}\mathrm{e}}^{\mathrm{2}}\in [0,1]$

Nagelkerkes Pseudo-R² erweitert Maddalas Pseudo-R², sodass durch eine Reskalierung ein möglicher Wert von 1 erreicht werden kann, wenn das vollständige Modell eine perfekte Vorhersage mit einer Wahrscheinlichkeit von 1 trifft.

Nagelkerke gab auch allgemeine Bedingungen für ein Pseudo-Bestimmtheitsmaß an:

1. Ein Pseudo-Bestimmtheitsmaß sollte mit dem Bestimmtheitsmaß ${\mathrm{R}}^{\mathrm{2}}$ übereinstimmen, wenn beide berechnet werden können.
2. Es sollte ebenfalls mit der Maximum-Likelihood-Schätzung des Modells maximiert werden.
3. Es sollte, zumindest asymptotisch, unabhängig vom Stichprobenumfang sein.
4. Die Interpretation sollte die durch das Modell erklärte Variabilität von $Y$ sein.
5. Es sollte zwischen Null und Eins liegen. Bei einem Wert von Null sollte es keine Aussage über die Variabilität von $Y$ machen; bei einem Wert von Eins, sollte es die Variabilität von $Y$ vollständig erklären.
6. Es sollte keine Maßeinheit besitzen.

Log-Likelihood-basierte Maße

McFadden R²

$R_{\mathrm{M}\mathrm{c}\mathrm{F}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{e}\mathrm{n}}^{\mathrm{2}}=1-\frac{\ln L_1}{\ln L_0}=1-\frac{D_1}{D_0}$,

mit

$L_0$: Nullmodell,

$L_1$: Modell mit erklärenden Variablen

$D_i$ die Devianz (Statistik) der entsprechenden Modelle

$R_{\mathrm{M}\mathrm{c}\mathrm{F}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{e}\mathrm{n}}^{\mathrm{2}}\in [0,1)$

Das Verhältnis der logarithmierten der Werte $L_1$ und $L_0$ der Likelihood-Funktion (Wahrscheinlichkeiten) spiegelt den Grad der Verbesserung des vollständigen Modells mit Prädiktoren gegenüber dem Nullmodell wider. Ein Modell mit einem größeren McFaddens hat eine bessere Anpassung gegenüber einem anderen Modell mit einem geringeren Wert.

Daumenregel: Bereits $0,2<R_{\mathrm{M}\mathrm{c}\mathrm{F}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{e}\mathrm{n}}^{\mathrm{2}}<0,4$ stellt ein besonders gute Anpassung des Modells dar.

McFaddens korrigiertes R²

$R_{\mathrm{M}\mathrm{c}\mathrm{F}\mathrm{a}\mathrm{d}\mathrm{d}\mathrm{e}{\mathrm{n}}_{\mathrm{k}\mathrm{o}\mathrm{r}\mathrm{r}\mathrm{.}}}^{\mathrm{2}}=1-\frac{\ln L_1-K}{\ln L_0}$

Das korrigierte McFaddens bewertet die Anzahl der Prädiktoren $K$ für die Anpassungsgüte eines Modells. Ähnlich dem korrigierten Bestimmtheitsmaß ${\overline{\mathrm{R}}}^2$ verringern zu viele Prädiktoren, die dem Modell nicht genügend beitragen, die Effektivität eines Modells und schlagen sich negativ im korrigierten McFaddens ${\text{Pseudo-R}}^2$ nieder. Somit sind Werte kleiner 0 möglich.

Aldrich / Nelsons R²

$R_{\mathrm{A}\mathrm{l}\mathrm{d}\mathrm{r}\mathrm{i}\mathrm{c}\mathrm{h}\mathrm{N}\mathrm{e}\mathrm{l}\mathrm{s}\mathrm{o}\mathrm{n}}^{\mathrm{2}}=1-\frac{2\cdot (\ln L_1-L_0)}{2\cdot (\ln L_1-L_0)+c\cdot n}$, c = 1 (Probit-Modell), 3,29 (Logit-Modell)

$R_{\mathrm{A}\mathrm{l}\mathrm{d}\mathrm{r}\mathrm{i}\mathrm{c}\mathrm{h}\mathrm{N}\mathrm{e}\mathrm{l}\mathrm{s}\mathrm{o}\mathrm{n}}^{\mathrm{2}}\in [0,1)$

Aldrich / Nelsons setzt den Likelihood-Quotienten ins Verhältnis, der die Rate von Nullmodell und Alternativmodell bei eingetretenem Ereignis angibt. Es hat eine obere Grenze von weit unter 1.

Korrelations-basierte Maße

Lave / Efrons R²

$R_{\mathrm{L}\mathrm{a}\mathrm{v}\mathrm{e}}^{\mathrm{2}}=1-\frac{\frac{1}{n}\sum _{i=1}^n(Y_i-{\hat{P}}_i{)}^2}{\frac{1}{n}\sum _{i=1}^n(Y_i-\overline{Y}{)}^2}$

Lave / Efrons kann man ähnlich dem normalen Bestimmtheitsmaß als Quadrat der Korrelation und als erklärte Variabilität interpretieren. Es werden die quadrierten Residuen aufsummiert, wobei ${\hat{P}}_i$ eine vom Modell vorhergesagte Wahrscheinlichkeit für $Y_i=1$ ist, welche die diskrete abhängige Variable in eine stetige überführt (Hinweis: $Y_i$ kann nur die Werte 0 und 1 annehmen).

Auf der erklärten Variation basierend

McKelvey & Zavoinas R²

$R_{\mathrm{M}\mathrm{c}\mathrm{K}\mathrm{e}\mathrm{l}\mathrm{v}\mathrm{e}\mathrm{y}}^{\mathrm{2}}=\frac{\widehat{\mathrm{Var}}({\hat{y}}^{*})}{\widehat{\mathrm{Var}}({\hat{y}}^{*})+\mathrm{Var}(e)}$

$R_{\mathrm{M}\mathrm{c}\mathrm{K}\mathrm{e}\mathrm{l}\mathrm{v}\mathrm{e}\mathrm{y}}^{\mathrm{2}}\in [0,1]$

McKelvey & Zavoinas ist strukturell dem normalen Bestimmtheitsmaß nachempfunden. Es wird die geschätzte erklärte Quadratsumme der Regression mit der geschätzten erklärten und unerklärten Quadratsumme von Regression und Fehler ins Verhältnis gesetzt.

Vergleichbarkeit

Die Werte der verschiedenen Pseudo-Bestimmtheitsmaße können innerhalb eines Modells stark variieren. Somit können unterschiedliche Maße zwischen verschiedenen Datensätzen nicht miteinander verglichen und unabhängig interpretiert werden. Als beste Approximation hat sich McKelvey & Zavoinas erwiesen; Laves, McFaddens, Nagelkerkes unterschätzen das "wahre" ${\mathrm{R}}^{\mathrm{2}}$ einer Kleinste-Quadrate-Schätzung für ein Modell mit latenten Variablen stark.

Beispiel

Ein Wäscheklammerproduzent möchte seine neuartigen Wäscheklammern auf den Markt bringen und deswegen vorab die Wahrscheinlichkeit eines Kaufes berechnen. Er berät sich mit seinem Geschäftspartner, der ein Statistikprogramm besitzt. Dieser nimmt an, dass der Kauf nur von einem Attribut abhängt, dem Preis $X_{\text{Preis}}$. Der aggregierte Einfluss auf die Kaufentscheidung soll eine lineare Beziehung haben, $Z=b_0+b_1\cdot X_{\text{Preis}}$, auch Logit genannt. Der Wäscheklammerproduzent hingegen glaubt eher, dass die Kaufabsicht vom Preis, der Farbe und der Größe abhängt:$Z=b_0+b_1\cdot X_{\text{Preis}}+b_2\cdot Y_{\text{Farbe}}+b_3\cdot W_{\text{Größe}}$. Mittels Marktforschungsdaten sind die Regressionsparameter $b_0$, $b_1$, $b_2$ und $b_3$ nach der Maximum-Likelihood-Schätzung vom Computer iterativ ermittelt worden. Allerdings fragt sich nun der Wäscheklammerproduzent, welche Modellhypothese die Realität besser wiedergibt und auf welche man weitere Überlegungen stützen sollte. Zur Einschätzung der Anpassungsgüte der angenommenen Modelle an die vorhandenen Daten sollen verschiedene Pseudo-Bestimmtheitsmaße benutzt werden. Diese lassen sich die beiden Geschäftspartner vom Statistikprogramm ausgeben.

${\text{Pseudo-R}}^2$
Maße der Anpassungsgüte	Modell 1 ($b_0,b_1$)	Modell 2 ($b_0,b_1,b_2,b_3$)
McFadden R2	0,307	0,445
McFadden Adj R2	0,273	0,389
Cragg-Uhler (Nagelkerke) R2	0,436	0,578
McKelvey & Zavoina R2	0,519	0,643
Efron / Lave R2	0,330	0,472

Da die Pseudo-Bestimmtheitsmaße für Modell 2 durchweg höher sind, d. h., dass dieses Modell die Marktforschungsdaten besser abbildet, entscheidet man sich für dieses und schätzt damit die Kaufwahrscheinlichkeit bzw. den möglichen Marktanteil.

Siehe auch

• Devianz (Statistik)

Weblinks

• FAQ: What are pseudo R-squareds?

Literatur

• Cragg, J.G., Uhler, R. (1970), "The Demand for Automobiles", Canadian Journal of Economics 3, S. 386–406, Vorlage:JSTOR.
• Hagle, T. M., Mitchell II, G. E. (1992), "Goodness-of-Fit Measures for Probit and Logit", American Journal of Political Science 36, S. 762–784, Vorlage:JSTOR.
• McFadden, D. (1973), "Conditional Logit Analysis of Qualitative Choice Behavior" (PDF 1,77 MB), in: P. Zarembka (ed.) Frontiers in Econometrics, Academic Press: New York, ISBN 0-12-776150-0, S. 105–142.
• McKelvey, R., Zavoina, W. (1975), "A Statistical Model for the Analysis of Ordinal Level Dependent Variables", Journal of Mathematical Sociology 4, S. 103–120, Vorlage:DOI.
• Nagelkerke, N. J. D. (1991), "A Note on a General Definition of the Coefficient of Determination", Biometrika 78, Nr. 3, S. 691–692, Vorlage:DOI.
• Veall, M. R., Zimmermann, K. F. (1996), "Pseudo-R² Measures for Some Common Limited Dependent Variable Models", Sonderforschungsbereich 386, Paper 18, Vorlage:DOI.