ISSN:
0170-4214
Keywords:
Mathematics and Statistics
;
Applied Mathematics
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mathematics
Notes:
The expectation maximization (EM) algorithm is an iterative procedure used to determine maximum likelihood estimators in situations of incomplete data. In the case of independent Poisson variables it converges to a solution of a problem of the form min ∑[〈ai,x〉 - bi log 〈ai, x〉] such that x ≥0. Thus, it can be used to solve systems of the form Ax = b, x≥0 (with A stochastic and b positive.) It converges to a feasible solution if it exists and to an approximate one otherwise (the one that minimizes d (b, Ax), where d is the Kullback-Leibler information divergence). We study the convergence of the multiplicatively relaxed version proposed by Tanaka for use in positron emission tomography. We prove global convergence in the underrelaxed and unrelaxed cases. In the overrelaxed case we present a local convergence theorem together with two partial global results: the sequence generated by the algorithm is bounded and, if it converges, its limit point is a solution of the aforementioned problem.
Additional Material:
1 Ill.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/mma.1670140805
