Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
16/2013

Computing medians and means in Hadamard spaces

Miroslav Bačák

Abstract

The geometric median as well as the Fréchet mean of points in an Hadamard space are important in both theory and applications. Surprisingly, no algorithms for their computation are hitherto known. To address this issue, we use a split version of the proximal point algorithm for minimizing a sum of convex functions and prove that this algorithm produces a sequence converging to a minimizer of the objective function, which extends a recent result of D. Bertsekas (2001) into Hadamard spaces. The method is quite robust and not only does it yield algorithms for the median and the mean, but it also applies to various other optimization problems. We moreover show that another algorithm for computing the Fréchet mean can be derived from the law of large numbers due to K.-T. Sturm (2002).

In applications, computing medians and means is probably most needed in tree space, which is an instance of an Hadamard space, invented by Billera, Holmes, and Vogtmann (2001) as a tool for averaging phylogenetic trees. It turns out, however, that it can be also used to model numerous other tree-like structures. Since there now exists a polynomial-time algorithm for computing geodesics in tree space due to M. Owen and S. Provan (2011), we obtain efficient algorithms for computing medians and means, which can be directly used in practice.

Received:
Feb 1, 2013
Published:
Feb 1, 2013
MSC Codes:
51F99, 60B99, 92B05, 49M27, 62E99
Keywords:
Hadamard space, tree space, Frechet mean, law of large numbers, proximal point algorithm, median and mean

Related publications

inJournal
2014 Repository Open Access
Miroslav Bačák

Computing medians and means in Hadamard spaces

In: SIAM journal on optimization, 24 (2014) 3, pp. 1542-1566