Inconsistency of Template Estimation with the Fréchet mean in Quotient Space

Loïc Devilliers (1) , Xavier Pennec (1) , Stéphanie Allassonnière (2)

1 ASCLEPIOS - Analysis and Simulation of Biomedical Images
2 CRC - Centre de Recherche des Cordeliers

We tackle the problem of template estimation when data have been randomly transformed under an isometric group action in the presence of noise. In order to estimate the template, one often minimizes the variance when the influence of the transformations have been removed (computation of the Fréchet mean in quotient space). The consistency bias is defined as the distance (possibly zero) between the orbit of the template and the orbit of one element which minimizes the variance. In this article we establish an asymptotic behavior of the consistency bias with respect to the noise level. This behavior is linear with respect to the noise level. As a result the inconsistency is unavoidable as soon as the noise is large enough. In practice, the template estimation with a finite sample is often done with an algorithm called max-max. We show the convergence of this algorithm to an empirical Karcher mean. Finally, our numerical experiments show that the bias observed in practice cannot be attributed to the small sample size or to a convergence problem but is indeed due to the previously studied inconsistency.

Domaines: Statistiques [math.ST]

Loïc Devilliers, Xavier Pennec, Stéphanie Allassonnière. Inconsistency of Template Estimation with the Fréchet mean in Quotient Space. Information Processing in Medical Imaging. IPMI 2017., Martin Styner; Marc Niethammer; Dinggang Shen; Stephen Aylward; Ipek Oguz; Hongtu Zhu, Jun 2017, Boone, United States. pp.16-27, ⟨10.1007/978-3-319-59050-9_2⟩. ⟨hal-01481074v2⟩