Triangulating Statistical Manifolds

Auteurs :
Publication MaxEnt 2014


Triangulating Statistical Manifolds


44.03 Ko


Creative Commons Aucune (Tous droits réservés)


Sponsors scientifique


Sponsors logistique


Sponsors financier

<resource  xmlns:xsi=""
        <identifier identifierType="DOI">10.23723/9603/11316</identifier><creators><creator><creatorName>Jean-Daniel Boissonnat</creatorName></creator><creator><creatorName>David Cohen-Steiner</creatorName></creator><creator><creatorName>Pooran Memari</creatorName></creator><creator><creatorName>Arnaud Poinas</creatorName></creator></creators><titles>
            <title>Triangulating Statistical Manifolds</title></titles>
        <resourceType resourceTypeGeneral="Text">Text</resourceType><dates>
	    <date dateType="Created">Sat 30 Aug 2014</date>
	    <date dateType="Updated">Mon 2 Oct 2017</date>
            <date dateType="Submitted">Thu 22 Mar 2018</date>
	    <alternateIdentifier alternateIdentifierType="bitstream">074d3057eefcaa959fa94016dda5dccda2d1f8c2</alternateIdentifier>
            <description descriptionType="Abstract"></description>

Triangulating Statistical Manifolds Jean-Daniel Boissonnat 1 David Cohen-Steiner 1 Pooran Memari 2 Arnaud Poinas 1 1 INRIA Sophia Antipolis - M´editerran´ee, 2 CNRS - LTCI, T´el´ecom ParisTech In this paper, we describe an algorithm to construct an intrinsic Delaunay triangulation of a smooth closed statistical manifold, i.e. a manifold of probability density functions (PDFs). In differ- ent applications (namely in data analysis, information geometry, biology, shape analysis, compression, scientific visualization), statistical data are often seen as points in a metric space, which is most fre- quently high dimensional. The generally accepted hypothesis is that, although it is embedded in spaces of high dimensions, data lives usually close to a much smaller structure with a low intrinsic dimension. Making use of this intrinsic geometric structure, this paper will develop a certified algo- rithm that can reconstruct statistical manifolds with a complexity that depends only linearly on the ambient dimension. We consider the Fisher information metric defined as a particular Riemannian metric on smooth statistical manifolds to calculate the informational difference between measurements. Interestingly, this metric can be understood as an infinitesimal form of the relative entropy of distributions; more specifically, it is the Hessian of the Kullback-Leibler divergence which is a particular instance of well-studied Bregman divergences for which efficient reconstruction algorithms have been proposed [1]. Our algorithm builds over existing algorithms for Bregman Delaunay triangulations [1] and mani- fold reconstruction using the tangential Delaunay complex proposed in [2]. The central idea is to define Bregman Delaunay triangulations locally and to glue these local triangulations together by removing inconsistencies between them. We view the inconsistencies as arising from instability in the Bregman Delaunay triangulations, and exploit and adapt the results presented in [3] for the Euclidean case. In particular, our technique heavily uses duality. Our main result is an algorithm whose complexity depends only linearly on the ambient dimension, and produces a Bregman Delaunay complex which is guaranteed to be a triangulation of the manifold under appropriate sampling conditions. We also demonstrate that the resulting triangulation coincides with the intrinsic Delaunay complex for the Fisher metric. References [1] Jean-Daniel Boissonnat, Frank Nielsen, Richard Nock, Bregman voronoi diagrams, Discrete Com- putational Geometry, 44, 2, 281-307, (2010). [2] Jean-Daniel Boissonnat, Arijit Ghosh, Manifold reconstruction using tangential Delaunay com- plexes, Discrete Computational Geometry, 51, 1, 221-267, (2014). [3] Jean-Daniel Boissonnat, Ramsay Dyer, and Arijit Ghosh, Constructing Intrinsic Delaunay Tri- angulations of Submanifolds, INRIA Research Report RR-8273, 2013. 1