Deming-Regression

In der Statistik wird mit der Deming-Regression eine Ausgleichsgerade für eine endliche Menge metrisch skalierter Datenpaare ($x_{i},y_{i}$) nach der Methode der kleinsten Quadrate bestimmt. Es handelt sich um eine Variante der linearen Regression. Bei der Deming-Regression werden die Residuen (Messfehler) sowohl für die $x$- als auch für die $y$-Werte in das Modell einbezogen.

Die Deming-Regression ist somit ein Spezialfall der Regressionsanalyse; sie beruht auf einer Maximum-Likelihood-Schätzung der Regressionsparameter, bei der die Residuen beider Variablen als unabhängig und normalverteilt angenommen werden und der Quotient $\delta$ ihrer Varianzen als bekannt unterstellt wird.

Die Deming-Regression geht auf eine Arbeit von C.H. Kummell (1879) zurück; 1937 wurde die Methode von T.C. Koopmans wieder aufgegriffen und in allgemeinerem Rahmen 1943 von W. E. Deming für technische und ökonomische Anwendungen bekannt gemacht.

Die orthogonale Regression ist ein wichtiger Spezialfall der Deming-Regression; sie behandelt den Fall $\delta =1$. Die Deming-Regression wiederum ist ein Spezialfall der York-Regression.

Rechenweg

Die gemessenen Werte $x_{i}$ und $y_{i}$ werden als Summen der „wahren Werte“ $x_{i}^{*}$ bzw. $y_{i}^{*}$ und der „Fehler“ $\eta _{i}$ bzw. $\varepsilon _{i}$ aufgefasst, d. h. $(x_{i},y_{i})=(x_{i}^{*}+\eta _{i},y_{i}^{*}+\varepsilon _{i})$ Die Datenpaare ($x_{i}^{*},y_{i}^{*}$) liegen auf der zu berechnenden Geraden. $\eta _{i}$ und $\varepsilon _{i}$ seien unabhängig mit bekanntem Quotienten der Fehlervarianzen $\delta ={\frac {\sigma _{\varepsilon }^{2}}{\sigma _{\eta }^{2}}}$.

Es wird eine Gerade

$y=\beta _{0}+\beta _{1}x$

gesucht, die die gewichtete Residuenquadratsumme minimiert:

$SQR=\sum _{i=1}^{n}{\bigg (}{\frac {\eta _{i}^{2}}{\sigma _{\eta }^{2}}}+{\frac {\varepsilon _{i}^{2}}{\sigma _{\varepsilon }^{2}}}{\bigg )}={\frac {1}{\sigma _{\varepsilon }^{2}}}\sum _{i=1}^{n}{\Big (}(y_{i}-\beta _{0}-\beta _{1}x_{i}^{*})^{2}+\delta (x_{i}-x_{i}^{*})^{2}{\Big )}\ \to \ \min _{\beta _{0},\beta _{1},x_{i}^{*}}SQR$

Für die weitere Rechnung werden die folgenden Hilfswerte benötigt:

${\overline {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}$     (arithmetisches Mittel der $x_{i}$)

${\overline {y}}={\frac {1}{n}}\sum _{i=1}^{n}y_{i}$     (arithmetisches Mittel der $y_{i}$)

$s_{x}^{2}={\tfrac {1}{n-1}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{2}$     (Stichprobenvarianz der $x_{i}$)

$s_{y}^{2}={\tfrac {1}{n-1}}\sum _{i=1}^{n}(y_{i}-{\overline {y}})^{2}$     (Stichprobenvarianz der $y_{i}$)

$s_{xy}={\tfrac {1}{n-1}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})(y_{i}-{\overline {y}})$     (Stichprobenkovarianz der $(x_{i},y_{i})$).

Damit ergeben sich die Parameter zur Lösung des Minimierungsproblems:

$\beta _{1}={\frac {s_{y}^{2}-\delta s_{x}^{2}+{\sqrt {(s_{y}^{2}-\delta s_{x}^{2})^{2}+4\delta s_{xy}^{2}}}}{2s_{xy}}}$

$\beta _{0}={\overline {y}}-\beta _{1}{\overline {x}}$.

Die $x_{i}^{*}$-Koordinaten berechnet man mit

$x_{i}^{*}=x_{i}+{\frac {\beta _{1}}{\beta _{1}^{2}+\delta }}(y_{i}-\beta _{0}-\beta _{1}x_{i})$.