Misplaced Pages

MM algorithm

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.

The MM algorithm is an iterative optimization method which exploits the convexity of a function in order to find its maxima or minima. The MM stands for “Majorize-Minimization” or “Minorize-Maximization”, depending on whether the desired optimization is a minimization or a maximization. Despite the name, MM itself is not an algorithm, but a description of how to construct an optimization algorithm.

The expectation–maximization algorithm can be treated as a special case of the MM algorithm. However, in the EM algorithm conditional expectations are usually involved, while in the MM algorithm convexity and inequalities are the main focus, and it is easier to understand and apply in most cases.

History

The historical basis for the MM algorithm can be dated back to at least 1970, when Ortega and Rheinboldt were performing studies related to line search methods. The same concept continued to reappear in different areas in different forms. In 2000, Hunter and Lange put forth "MM" as a general framework. Recent studies have applied the method in a wide range of subject areas, such as mathematics, statistics, machine learning and engineering.

Algorithm

MM algorithm

The MM algorithm works by finding a surrogate function that minorizes or majorizes the objective function. Optimizing the surrogate function will either improve the value of the objective function or leave it unchanged.

Taking the minorize-maximization version, let f ( θ ) {\displaystyle f(\theta )} be the objective concave function to be maximized. At the m step of the algorithm, m = 0 , 1... {\displaystyle m=0,1...} , the constructed function g ( θ | θ m ) {\displaystyle g(\theta |\theta _{m})} will be called the minorized version of the objective function (the surrogate function) at θ m {\displaystyle \theta _{m}} if

g ( θ | θ m ) ≤ f ( θ )  for all  θ {\displaystyle g(\theta |\theta _{m})\leq f(\theta ){\text{ for all }}\theta }
g ( θ m | θ m ) = f ( θ m ) {\displaystyle g(\theta _{m}|\theta _{m})=f(\theta _{m})}

Then, maximize g ( θ | θ m ) {\displaystyle g(\theta |\theta _{m})} instead of f ( θ ) {\displaystyle f(\theta )} , and let

θ m + 1 = arg ⁡ max θ g ( θ | θ m ) {\displaystyle \theta _{m+1}=\arg \max _{\theta }g(\theta |\theta _{m})}

The above iterative method will guarantee that f ( θ m ) {\displaystyle f(\theta _{m})} will converge to a local optimum or a saddle point as m goes to infinity. By the above construction

f ( θ m + 1 ) ≥ g ( θ m + 1 | θ m ) ≥ g ( θ m | θ m ) = f ( θ m ) {\displaystyle f(\theta _{m+1})\geq g(\theta _{m+1}|\theta _{m})\geq g(\theta _{m}|\theta _{m})=f(\theta _{m})}

The marching of θ m {\displaystyle \theta _{m}} and the surrogate functions relative to the objective function is shown in the figure.

Majorize-Minimization is the same procedure but with a convex objective to be minimised.

Constructing the surrogate function

One can use any inequality to construct the desired majorized/minorized version of the objective function. Typical choices include

References

  1. Lange, Kenneth. "The MM Algorithm" (PDF).
  2. Lange, Kenneth (2016). MM Optimization Algorithms. SIAM. doi:10.1137/1.9781611974409. ISBN 978-1-61197-439-3.
  3. Lange, K.; Zhou, H. (2022). "A legacy of EM algorithms". International Statistical Review. 90 (Suppl 1): S52–S66. doi:10.1111/insr.12526. PMC 10191373. PMID 37204987.
  4. Ortega, J.M.; Rheinboldt, W.C. (1970). Iterative Solutions of Nonlinear Equations in Several Variables. New York: Academic. pp. 253–255. ISBN 9780898719468.
  5. Hunter, D.R.; Lange, K. (2000). "Quantile Regression via an MM Algorithm". Journal of Computational and Graphical Statistics. 9 (1): 60–77. CiteSeerX 10.1.1.206.1351. doi:10.2307/1390613. JSTOR 1390613.
  6. Wu, C. F. Jeff (1983). "On the Convergence Properties of the EM Algorithm". Annals of Statistics. 11 (1): 95–103. doi:10.1214/aos/1176346060. JSTOR 2240463.
Category: