A new smoothing-regularization approach for a maximum-likelihood estimation problem |
| |
Authors: | Alfredo Noel Iusem B. F. Svaiter |
| |
Affiliation: | (1) Instituto de Matemática Pura e Aplicada, Estrada Dona Castorina 110, Jardim Botânico, CEP 22.460 Rio de Janeiro, Brazil |
| |
Abstract: | We consider the problem min i=1m( ai,x –bilog ai, z ) subject tox 0 which occurs as a maximum-likelihood estimation problem in several areas, and particularly in positron emission tomography. After noticing that this problem is equivalent to mind(b, Ax) subject tox 0, whered is the Kullback-Leibler information divergence andA, b are the matrix and vector with rows and entriesai,bi, respectively, we suggest a regularized problem mind(b, Ax) + d(v, Sx), where is the regularization parameter,S is a smoothing matrix, andv is a fixed vector. We present a computationally attractive algorithm for the regularized problem, establish its convergence, and show that the regularized solutions, as goes to 0, converge to the solution of the original problem which minimizes a convex function related tod(v, Sx). We give convergence-rate results both for the regularized solutions and for their functional values.The research of A. N. Iusem was partially supported by CNPq Grant No. 301280/86-MA. |
| |
Keywords: | Convex programming EM algorithm Maximum likelihood Statistical estimation Regularization Positron emission tomography |
本文献已被 SpringerLink 等数据库收录! |
|