The method of iteratively reweighted least squares (IRLS) is used to solve certain optimization problems with objective functions of the form of a p-norm,
a r g m i n β ∑ i = 1 n | y i − f i ( β ) | p , {\displaystyle \operatorname {arg\,min} _{\boldsymbol {\beta }}\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
β ( t + 1 ) = a r g m i n β ∑ i = 1 n w i ( β ( t ) ) | y i − f i ( β ) | 2 . {\displaystyle {\boldsymbol {\beta }}^{(t+1)}=\operatorname {arg\,min} _{\boldsymbol {\beta }}\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 errors rather than the least square errors. 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
L1 minimization for sparse recovery IRLS can be used for ℓ1 minimization and smoothed ℓp minimization, p < 1, in compressed sensing problems. It has been proved that the algorithm has a linear rate of convergence for ℓ1 norm and superlinear for ℓt with t < 1, under the restricted isometry property, which is generally a sufficient condition for sparse solutions.
Lp norm linear regression To find the parameters β = (β1, …,βk)T which minimize the Lp norm for the linear regression problem,
a r g m i n β ‖ y − X β ‖ p = a r g m i n β ∑ i = 1 n | y i − X i β | p , {\displaystyle {\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
… excerpt ends here. Continue reading the full article.
