Prädiktionsmatrix

In der Statistik ist die Prädiktionsmatrix (Vorlage:EnS prediction matrix) eine symmetrische und idempotente Matrix und damit eine Projektionsmatrix. Die Prädiktionsmatrix wird gelegentlich Hut-Matrix oder Dach-Matrix genannt, da sie $y$ auf $\hat{y}$ abbildet. Dementsprechend wird sie entweder mit $𝐏$ oder $𝐇$ notiert. Der Begriff „Prädiktionsmatrix“ bzw. „Vorhersagematrix“ wurde von Hoaglin & Welsh (1978) sowie Chatterjee & Hadi (1986) geprägt und rührt daher, dass wenn man die Matrix auf die $y$-Werte anwendet sie die vorhergesagten Werte ($\hat{y}$-Werte) generiert. Eine weitere in der Statistik wichtige Matrix ist die Residualmatrix, die durch die Prädiktionsmatrix definiert wird und ebenfalls eine Projektionsmatrix ist.

Definition

Gegeben ein typisches multiples lineares Regressionsmodell $𝐲=𝐗\mathit{\beta }+\mathit{\epsilon }$, mit $\mathit{\beta }$ dem $p\times 1$ Vektor der unbekannten Regressionsparameter, der $n\times p$ Versuchsplanmatrix $𝐗$, dem $n\times 1$ Vektor der abhängigen Variablen $𝐲$ und dem $n\times 1$ Vektor der Störgrößen $\mathit{\epsilon }$. Dann ist die Prädiktionsmatrix definiert durch

$𝐏\equiv 𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}$ mit $𝐏\in {ℝ}^{n\times n}$.

Die Matrix ${𝐗}^{+}={\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}$ wird auch Moore-Penrose-Inverse von $𝐗$ genannt.

Die mithilfe der Methode der kleinsten Quadrate geschätzte Regressions(hyper)ebene ist dann gegeben durch die Stichproben-Regressionsfunktion $\widehat{𝐲}=\widehat{\mathrm{E}(𝐲)}=𝐗\hat{\mathit{\beta }}$, wobei $\hat{\mathit{\beta }}={\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}𝐲$ der Kleinste-Quadrate-Schätzvektor ist. Die Prädiktionsmatrix $𝐏$ ist die Matrix der Orthogonalprojektion auf den Spaltenraum von $𝐗$ und hat maximal den Rang $p$ ($p=k+1$ ist die Anzahl der Parameter des Regressionsmodells). Falls $𝐗$ eine $(n\times p)$ Matrix mit $\mathrm{Rang}(𝐗)=p$ ist, dann ist $\mathrm{Rang}(𝐏)=p$. Da $𝐏$ eine Projektionsmatrix ist, gilt $\mathrm{Rang}(𝐏)=\mathrm{Spur}(𝐏)=p$. Die Idempotenz- und die Symmetrieeigenschaft ($𝐏\cdot 𝐏=𝐏$ und ${𝐏}^{\mathrm{\top }}=𝐏$) implizieren, dass $𝐏$ ein orthogonaler Projektor auf den Spaltenraum $S(𝐗)=S(𝐏)$ ist. Die Projektionsrichtung ergibt sich aus der Matrix $(𝐈-𝐏)$, deren Spalten senkrecht auf $S(𝐗)$ stehen. Die Matrix $𝐏$ wird Prädiktionsmatrix genannt, da sich die Vorhersagewerte $\widehat{𝐲}$ durch die linksseitige Multiplikation des Vektors $𝐲$ mit dieser Matrix ergeben. Dies kann durch Einsetzen des KQ-Parameterschätzers wie folgt gezeigt werden:

$\widehat{𝐲}=𝐗\hat{\mathit{\beta }}=\underset{=𝐏}{\underbrace{𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}}}𝐲=𝐏𝐲$.

Die Vorhersagewerte von $y$ (die $\hat{y}$-Werte) können also als eine Funktion der beobachteten $y$-Werte verstanden werden. Zahlreiche statistische Resultate lassen sich auch mit der Prädiktionsmatrix darstellen. Beispielsweise lässt sich der Residualvektor mittels der Prädiktionsmatrix darstellen als: $\hat{\mathit{\epsilon }}=𝐲-\widehat{𝐲}=𝐲-𝐗\hat{\mathit{\beta }}=(𝐈-𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }})𝐲=(𝐈-𝐏)𝐲$. Die (nichttriviale) Kovarianzmatrix des Residualvektors lautet $\mathrm{Cov}(\hat{\mathit{\epsilon }})={\sigma }^2(𝐈-𝐏)$ und spielt für die Analyse von Hebelwerten eine Rolle.

Eigenschaften

Idempotenz

Die Prädiktionsmatrix ist idempotent. Dies kann so interpretiert werden, dass „zweimaliges Anwenden der Regression zum gleichen Ergebnis führt“. Die Idempotenzeigenschaft der Prädiktionsmatrix kann wie folgt gezeigt werden:

${𝐏}^2=𝐏\cdot 𝐏=𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}=𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}𝐈{𝐗}^{\mathrm{\top }}=𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}=𝐏$,

wobei $𝐈$ die Einheitsmatrix ist.

Symmetrie

Die Prädiktionsmatrix ist symmetrisch. Die Symmetrieeigenschaft der Prädiktionsmatrix kann wie folgt gezeigt werden

${𝐏}^{\mathrm{\top }}={\left(𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}\right)}^{\mathrm{\top }}={\left(\left(𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}\right)\left({𝐗}^{\mathrm{\top }}\right)\right)}^{\mathrm{\top }}={\left({𝐗}^{\mathrm{\top }}\right)}^{\mathrm{\top }}{\left(𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}\right)}^{\mathrm{\top }}=𝐗{\left({\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}\right)}^{\mathrm{\top }}{𝐗}^{\mathrm{\top }}=𝐗{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐗}^{\mathrm{\top }}=𝐏$

Hebelwerte

Die Diagonalelemente $p_{ii}$ der Prädiktionsmatrix $𝐏$ können als Hebelwerte interpretiert werden und spielen in der Regressionsdiagnostik eine große Rolle. Sie sind gegeben durch

$p_{ii}={𝐱}_i^{\mathrm{\top }}{\left({𝐗}^{\mathrm{\top }}𝐗\right)}^{-1}{𝐱}_i$.

Diese Hebelwerte werden bei der Berechnung des Cook-Abstands verwendet und können genutzt werden, um einflussreiche Beobachtungen zu identifizieren. Es gilt $\frac{1}{n}\le p_{ii}\le \frac{1}{r}$, wobei $r$ die Anzahl der Zeilen in der Versuchsplanmatrix $𝐗$ darstellt, die unterschiedlich sind. Wenn alle Zeilen unterschiedlich sind, dann gilt $\frac{1}{n}\le p_{ii}\le 1$.