Iteratively reweighted least squares

Regression analysis
Part of a series on Statistics

Models
Linear regression Simple regression Ordinary least squares Polynomial regression General linear model
Generalized linear model Discrete choice Logistic regression Multinomial logit Mixed logit Probit Multinomial probit Ordered logit Ordered probit Poisson
Multilevel model Fixed effects Random effects Mixed model
Nonlinear regression Nonparametric Semiparametric Robust Quantile Isotonic Principal components Least angle Local Segmented
Errors-in-variables
Estimation
Least squares Ordinary least squares Linear (math) Partial Total Generalized Weighted Non-linear Non-negative Iteratively reweighted Ridge regression
Least absolute deviations Bayesian Bayesian multivariate
Background
Regression model validation Mean and predicted response Errors and residuals Goodness of fit Studentized residual Gauss–Markov theorem
Statistics portal

The method of iteratively reweighted least squares (IRLS) is used to solve certain optimization problems with objective functions of the form:

{\underset {{\boldsymbol \beta }}{\operatorname {arg\,min}}}\sum _{{i=1}}^{n}{\big |}y_{i}-f_{i}({\boldsymbol \beta }){\big |}^{p},

by an iterative method in which each step involves solving a weighted least squares problem of the form:^[1]

{\boldsymbol \beta }^{{(t+1)}}={\underset {{\boldsymbol \beta }}{\operatorname {arg\,min}}}\sum _{{i=1}}^{n}w_{i}({\boldsymbol \beta }^{{(t)}}){\big |}y_{i}-f_{i}({\boldsymbol \beta }){\big |}^{2}.

IRLS is used to find the maximum likelihood estimates of a generalized linear model, and in robust regression to find an M-estimator, as a way of mitigating the influence of outliers in an otherwise normally-distributed data set. For example, by minimizing the least absolute error rather than the least square error.

Although not a linear regression problem, Weiszfeld's algorithm for approximating the geometric median can also be viewed as a special case of iteratively reweighted least squares, in which the objective function is the sum of distances of the estimator from the samples.

One of the advantages of IRLS over linear programming and convex programming is that it can be used with Gauss–Newton and Levenberg–Marquardt numerical algorithms.

Examples

L₁ minimization for sparse recovery

IRLS can be used for ℓ₁ minimization and smoothed ℓ_p minimization, p < 1, in the compressed sensing problems. It has been proved that the algorithm has a linear rate of convergence for ℓ₁ norm and superlinear for ℓ_t with t < 1, under the restricted isometry property, which is generally a sufficient condition for sparse solutions.^[2]^[3] However in most practical situations, the restricted isometry property is not satisfied.

L^p norm linear regression

To find the parameters β = (β₁, …,β_k)^T which minimize the L^p norm for the linear regression problem,

{\underset {{\boldsymbol \beta }}{\operatorname {arg\,min}}}{\big \|}{\mathbf y}-X{\boldsymbol \beta }\|_{p}={\underset {{\boldsymbol \beta }}{\operatorname {arg\,min}}}\sum _{{i=1}}^{n}\left|y_{i}-X_{i}{\boldsymbol \beta }\right|^{p},

the IRLS algorithm at step t + 1 involves solving the weighted linear least squares problem:^[4]

\boldsymbol\beta^{(t+1)} = \underset{\boldsymbol\beta}{ \operatorname{arg\,min} } \sum_{i=1}^n w_i^{(t)} \left| y_i - X_i \boldsymbol\beta \right|^2 = (X^{\rm T} W^{(t)} X)^{-1} X^{\rm T} W^{(t)} \mathbf{y},

where W^(t) is the diagonal matrix of weights, usually with all elements set initially to:

w_{i}^{{(0)}}=1

and updated after each iteration to:

w_{i}^{{(t)}}={\big |}y_{i}-X_{i}{\boldsymbol \beta }^{{(t)}}{\big |}^{{p-2}}.

In the case p = 1, this corresponds to least absolute deviation regression (in this case, the problem would be better approached by use of linear programming methods,^[5] so the result would be exact) and the formula is:

w_{i}^{{(t)}}={\frac {1}{{\big |}y_{i}-X_{i}{\boldsymbol \beta }^{{(t)}}{\big |}}}.

To avoid dividing by zero, regularization must be done, so in practice the formula is:

w_{i}^{(t)}={\frac {1}{\max \left\{\delta ,\left|y_{i}-X_{i}{\boldsymbol {\beta }}^{(t)}\right|\right\}}}.

where $\delta$ is some small value, like 0.0001.^[5] Note the use of $\delta$ in the weighting function is equivalent to the Huber loss function in robust estimation.

Notes

↑ C. Sidney Burrus, Iterative Reweighted Least Squares
↑ Chartrand, R.; Yin, W. (March 31 – April 4, 2008). "Iteratively reweighted algorithms for compressive sensing". IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), 2008. pp. 3869–3872.
↑ Daubechies, I.; Devore, R.; Fornasier, M.; Güntürk, C. S. N. (2010). "Iteratively reweighted least squares minimization for sparse recovery". Communications on Pure and Applied Mathematics. 63: 1. doi:10.1002/cpa.20303.
↑ Gentle, James (2007). "6.8.1 Solutions that Minimize Other Norms of the Residuals". Matrix algebra. New York: Springer. ISBN 978-0-387-70872-0. doi:10.1007/978-0-387-70873-7.
1 2 William A. Pfeil, Statistical Teaching Aids, Bachelor of Science thesis, Worcester Polytechnic Institute, 2006

References

Numerical Methods for Least Squares Problems by Åke Björck (Chapter 4: Generalized Least Squares Problems.)
Practical Least-Squares for Computer Graphics. SIGGRAPH Course 11

External links

Solve under-determined linear systems iteratively

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.