entropy
Review

An Elementary Introduction to Information Geometry
Frank Nielsen
Sony Computer Science Laboratories, Tokyo 141-0022, Japan; Frank.Nielsen@acm.org
Received: 6 September 2020; Accepted: 25 September 2020; Published: 29 September 2020


Abstract: In this survey, we describe the fundamental differential-geometric structures of information
manifolds, state the fundamental theorem of information geometry, and illustrate some use cases ofthese information manifolds in information sciences. The exposition is self-contained by concisely
introducing the necessary concepts of differential geometry. Proofs are omitted for brevity.
Keywords: differential geometry; metric tensor; affine connection; metric compatibility; conjugate
connections; dual metric-compatible parallel transport; information manifold; statistical manifold;curvature and flatness; dually flat manifolds; Hessian manifolds; exponential family; mixture family;
statistical divergence; parameter divergence; separable divergence; Fisher–Rao distance; statisticalinvariance; Bayesian hypothesis testing; mixture clustering; α-embeddings; mixed parameterization;
gauge freedom

1. Introduction
1.1. Overview of Information Geometry
We present a concise and modern view of the basic structures lying at the heart of Information
Geometry (IG), and report some applications of those information-geometric manifolds (herein termed“information manifolds”) in statistics (Bayesian hypothesis testing) and machine learning (statisticalmixture clustering).
By analogy to Information Theory (IT) (pioneered by Claude Shannon in his celebrated 1948
paper [1]) which considers primarily the communication of messages over noisy transmission
channels, we may define Information Sciences (IS) as the fields that study “communication” between(noisy/imperfect) data and families of models (postulated as a priori knowledge). In short, information
sciences seek methods to distill information from data to models. Thus information sciences encompass
information theory but also include the fields of Probability and Statistics, Machine Learning (ML),
Artificial Intelligence (AI), Mathematical Programming, just to name a few.
We review some key milestones of information geometry and report some definitions of the fieldby its pioneers in Section 5.2. Professor Shun-ichi Amari, the founder of modern information geometry,defined information geometry in the preface of his latest textbook [2] as follows: “Information geometry
is a method of exploring the world of information by means of modern geometry.” In short, informationgeometry geometrically investigates information sciences. It is a mathematical endeavour to define andbound the term geometry itself as geometry is open-ended. Often, we start by studying the invarianceof a problem (e.g., invariance of distance between probability distributions) and get as a result a novel
geometric structure (e.g., a “statistical manifold”). However, a geometric structure is “pure” andthus may be applied to other application areas beyond the scope of the original problem (e.g., useof the dualistic structure of statistical manifolds in mathematical programming [3]): the method ofgeometry [4] thus yields a pattern of abduction [5,6].
A narrower definition of information geometry can be stated as the field that studies the geometryof decision making. This definition also includes model fitting (inference) which can be interpreted

Entropy 2020, 22, 1100; doi:10.3390/e22101100

www.mdpi.com/journal/entropyEntropy 2020, 22, 1100

2 of 61

as a decision problem as illustrated in Figure 1; namely, deciding which model parameter to choose
from a family of parametric models. This framework was advocated by Abraham Wald [7–9] who
considered all statistical problems as statistical decision problems. Dissimilarities (also loosely called
distances among others) play a crucial role not only for measuring the goodness-of-fit of data tomodel (say, likelihood in statistics, classifier loss functions in ML, objective functions in mathematicalprogramming or operations research, etc.) but also for measuring the discrepancy (or deviance)
between models.

mθ̂n (D)

mθ2

mθ1
M
Figure 1. The parameter inference θ̂ of a model from data D can also be interpreted as a decision
making problem: decide which parameter of a parametric family of models M = {mθ }θ ∈Θ suits the“best” the data. Information geometry provides a differential-geometric structure on manifold M whichuseful for designing and studying statistical decision rules.

One may ponder why adopting a geometric approach? Geometry allows one to study invarianceof “figures” in a coordinate-free framework. The geometric language (e.g., line, ball or projection) also
provides affordances that help us reason intuitively about problems. Note that although figures canbe visualized (i.e., plotted in coordinate charts), they should be thought of as purely abstract objects,
namely, geometric figures.
Geometry also allows one to study equivariance: For example, the centroid c( T ) of a triangle
is equivariant under any affine transformation A: c( A.T ) = A.c( T ). In Statistics, the Maximum
Likelihood Estimator (MLE) is equivariant under a monotonic transformation g of the model parameter
θ: (\g(θ )) = g(θ̂ ), where the MLE of θ is denoted by θ̂.
1.2. Rationale and Outline of the Survey
The goal of this survey on information geometry [2] is to describe the core dualistic structures
on manifolds without assuming any prior background on differential geometry [10], and explain
several important related principles and concepts like invariance, covariance, projections, flatness andcurvature, information monotonicity, etc. In doing so, we shall illustrate the basic underlying concepts
with selected examples and applications, and shall make clear of some potential sources of confusion(e.g., a geometric statistical structure can be used in non-statistical applications [3], untangle themeaning of α in the α-connections, the α-divergences, and the α-representations, etc.). In particular,
we shall name and state the fundamental theorem of information geometry in Section 3.5. We refer
the reader to the books [2,4,11–17] for an indepth treatment of the field with its applications in
information sciences.
This survey is organized as follows:
In the first part (Section 2), we start by concisely introducing the necessary background on
differential geometry in order to define a manifold structure ( M, g, ∇), i.e., a manifold M equippedwith a metric tensor field g and an affine connection ∇. We explain how this framework generalizes
the Riemannian manifolds ( M, g) by stating the fundamental theorem of Riemannian geometry thatdefines a unique torsion-free metric-compatible Levi–Civita connection which can be derived from themetric tensor.

Entropy 2020, 22, 1100

3 of 61

In the second part (Section 3), we explain the dualistic structures of information manifolds: we
present the conjugate connection manifolds ( M, g, ∇, ∇∗ ), the statistical manifolds ( M, g, C ) where C
denotes a cubic tensor, and show how to derive a family of information manifolds ( M, g, ∇−α , ∇α ) forα ∈ R provided any given pair (∇ = ∇−1 , ∇∗ = ∇1 ) of conjugate connections. We explain how to get
conjugate connections ∇ and ∇∗ coupled to the metric g from any smooth (potentially asymmetric)
distances (called divergences), present the dually flat manifolds obtained when considering Bregman
divergences, and define, when dealing with parametric family of probability models, the exponential
connection e ∇ and the mixture connection m ∇ that are dual connections coupled to the Fisher
information metric. We discuss the concept of statistical invariance for the metric tensor and the notionof information monotonicity for statistical divergences [2,18]. It follows that the Fisher information
metric is the unique invariant metric (up to a scaling factor), and that the f -divergences are the uniqueseparable invariant divergences.
In the third part (Section 4), we illustrate how to use these information-geometric structuresin simple applications: First, we described the natural gradient descent method in Section 4.1 andits relationships with the Riemannian gradient descent and the Bregman mirror descent. Second,
we consider two applications in dually flat spaces in Section 4.2: In the first application, we consider
the problem of Bayesian hypothesis testing and show how Chernoff information (which defines thebest error exponent) can be geometrically characterized on the dually flat structure of an exponential
family manifold. In the second application, we show how to cluster statistical mixtures sharing thesame component distributions on the dually flat mixture family manifold.
Finally, we conclude in Section 5 by summarizing the important concepts and structures of
information geometry, and by providing further references and textbooks [2,16] for further readings tomore advanced structures and applications of information geometry. We also mention recent studies
of generic classes of principled distances and divergences.
In Appendix A, we show how to estimate the statistical f -divergences between two probability
distributions in order to ensure that the estimates are non-negative in Appendix B, and report the
canonical decomposition of the multivariate Gaussian family, an example of exponential family which
admits a dually flat structure.
At the beginning of each part, we start by outlining its contents. A summary of the notations usedthroughout this survey is provided in Appendix C.
2. Prerequisite: Basics of Differential Geometry
In Section 2.1, we review the very basics of Differential Geometry (DG) for defining a manifold( M, g, ∇) equipped with both a metric tensor field g and an affine connection ∇. We explainthese two independent metric/connection structures in Sections 2.2 and 2.3, respectively. Froman affine connection ∇, we show how to derive the notion of covariant derivative in Section 2.3.1,
parallel transport in Section 2.3.2 and geodesics in Section 2.3.3. We further explain the intrinsic
curvature and torsion of manifolds induced by the connection in Section 2.3.4, and state thefundamental theorem of Riemannian geometry in Section 2.4: the existence of a unique torsion-freeLevi–Civita connection LC ∇ compatible with the metric (metric connection) that can be derived fromthe metric tensor g. Thus the Riemannian geometry ( M, g) is obtained as a special case of the moregeneral manifold structure ( M, g, LC ∇): ( M, g) ≡ ( M, g, LC ∇). Information geometry shall further
consider a dual structure ( M, g, ∇∗ ) associated to ( M, g, ∇), and the pair of dual structures shall form
an information manifold ( M, g, ∇, ∇∗ ).
2.1. Overview of Differential Geometry: Manifold ( M, g, ∇)Informally speaking, a smooth D-dimensional manifold M is a topological space that locallybehaves like the D-dimensional Euclidean space RD . Geometric objects (e.g., points, balls, and vector
fields) and entities (e.g., functions and differential operators) live on M, and are coordinate-free but canconveniently be expressed in any local coordinate system of an atlas A = {(Ui , xi )}i of charts(Ui , xi )’sEntropy 2020, 22, 1100

4 of 61(fully covering the manifold) for calculations. Historically, René Descartes (1596–1650) allegedly
invented the global Cartesian coordinate system while wondering how to locate a fly on the ceilingfrom his bed. In practice, we shall use the most expedient coordinate system to facilitate calculations.
In information geometry, we usually handle a single chart fully covering the manifold.
A C k manifold is obtained when the change of chart transformations are C k . The manifold is saidsmooth when it is C ∞ . At each point p ∈ M, a tangent plane Tp locally best linearizes the manifold.
On any smooth manifold M, we can define two independent structures:
1.
2.

A metric tensor g, andAn affine connection ∇.

The metric tensor g induces on each tangent plane Tp an inner product space that allows one tomeasure vector magnitudes (vector “lengths”) and angles/orthogonality between vectors. The affine
connection ∇ is a differential operator that allows one to define:
1.
2.
3.
4.

The covariant derivative operator which provides a way to calculate differentials of a vector fieldY with respect to another vector field X: namely, the covariant derivative ∇ X Y,
The parallel transport ∏∇c which defines a way to transport vectors between tangent planes along
any smooth curve c,
The notion of ∇-geodesics γ∇ which are defined as autoparallel curves, thus extending theordinary notion of Euclidean straightness,
The intrinsic curvature and torsion of the manifold.

2.2. Metric Tensor Fields g
The tangent bundle of M is defined as the “union” of all tangent spaces:
TM := t p Tp := {( p, v),

p ∈ M, v ∈ Tp }.

(1)

Thus the tangent bundle TM of a D-dimensional manifold M is of dimension 2D (the tangent
bundle is a particular example of a fiber bundle with base manifold M).
Informally speaking, a tangent vector v plays the role of a directional derivative, with v f
informally meaning the derivative of a smooth function f (belonging to the space of smooth functions
F( M )) along the direction v. Since the manifolds are abstract and not embedded in some Euclidean
space, we do not view a vector as an “arrow” anchored on the manifold. Rather, vectors can beunderstood in several ways in differential geometry like directional derivatives or equivalent class ofsmooth curves at a point. That is, tangent spaces shall be considered as the manifold abstract too.
A smooth vector field X is defined as a “cross-section” of the tangent bundle: X ∈ X( M ) = Γ( TM),
where X( M ) or Γ( TM) denote the space of smooth vector fields. A basis B = {b1 , . . . , bD } of afinite D-dimensional vector space is a maximal linearly independent set of vectors: A set of vectors
B = {b1 , . . . , bD } is linearly independent if and only if ∑iD=1 λi bi = 0 iff λi = 0 for all i ∈ [ D ]. That is,
in a linearly independent vector set, no vector of the set can be represented as a linear combinationof the remaining vectors. A vector set is linearly independent maximal when we cannot add anotherlinearly independent vector. Tangent spaces carry algebraic structures of vector spaces. Furthermore,
to any vector space V, we can associate a dual covector space V ∗ which is the vector space of real-valued
linear mappings. We do not enter into details here to preserve this gentle introduction to informationgeometry with as little intricacy as possible. Using local coordinates on a chart (U , x ), the vector fieldΣX can be expressed as X = ∑iD=1 X i ei = X i ei using Einstein summation convention on dummy indicesΣ(using notation =), where ( X ) B := ( X i ) denotes the contravariant vector components (manipulated as“column vectors” in algebra) in the natural basis B = {e1 = ∂1 , . . . , e D = ∂ D } with ∂i :=: ∂x∂ . A tangentiplane (vector space) equipped with an inner product h·, ·i yields an inner product space. We definea reciprocal basis B∗ = {e∗ i = ∂i }i of B = {ei = ∂i }i so that vectors can also be expressed usingEntropy 2020, 22, 1100

5 of 61

the covariant vector components in the natural reciprocal basis. The primal and reciprocal basis aremutually orthogonal by construction as illustrated in Figure 2.

x2
hei , ej i = δij
e2

e2

e1

x1

e1
Figure 2.
Primal basis (red) and reciprocal basis (blue) of an inner product h·, ·i space.
The primal/reciprocal basis are mutually orthogonal: e1 is orthogonal to e2 , and e1 is orthogonal
to e2 .

For any vector v, its contravariant components vi ’s (superscript notation) and its covariant
components vi ’s (subscript notation) can be retrieved from v using the inner product with the use of
the reciprocal and primal basis, respectively:vi= hv, e∗ i i,

(2)

vi

= hv, ei i.

(3)

The inner product defines a metric tensor g and a dual metric tensor g∗ :
gij
g

∗ ij

:=

h ei , e j i,

(4)

:=

∗i

(5)

∗j

h e , e i.

Technically speaking, the metric tensor g p : Tp M × Tp M → R is a 2-covariant tensor field:Σg = gij dxi ⊗ dx j ,(6)where ⊗ is the dyadic tensor product performed on pairwise covector basis {dxi }i (the covectors
corresponding to the reciprocal vector basis). We do not describe tensors in details for sake of brevity.
A tensor is a geometric entity of a tensor space that can also be interpreted as a multilinear map.
A contravariant vector lives in a vector space while a covariant vector lives in the dual covector space.
We recommend the textbook [19] for a concise and well-explained description of tensors.
Let G = [ gij ] and G ∗ = [ g∗ ij ] denote the D × D matrices. It follows by construction of thereciprocal basis that G ∗ = G −1 . The reciprocal basis vectors e∗ i ’s and primal basis vectors ei ’s can be
expressed using the dual metric g∗ and metric g on the primal basis vectors e j ’s and reciprocal basis
vectors e∗ j ’s, respectively:
e∗ i

=

Σ

g∗ ij e j ,

(7)

ei

=

Σ

gij e∗ j .

(8)

The metric tensor field g (“metric tensor” or “metric” for short) defines a smooth symmetricpositive-definite bilinear form on the tangent bundle so that for u, v ∈ Tp , g(u, v) ≥ 0 ∈ R. We canalso write equivalently g p (u, v):=:hu, vi p :=:hu, vi g( p) :=:hu, vi. Two vectors u and v are said orthogonal,
denoted by u ⊥ v, iff hu, vi = 0. The length of a vector is induced from the norm kuk p :=:kuk g( p) =qhu, ui g( p) . Using local coordinates of a chart (U , x ), we get the vector contravariant/covariantEntropy 2020, 22, 1100

6 of 61

components, and compute the metric tensor using matrix algebra (with column vectors by convention)
as follows:>−1g(u, v) = (u)>(9)
B × Gx ( p ) × ( v ) B = ( u ) B ∗ × Gx ( p ) × ( v ) B ∗ ,since it follows from the primal/reciprocal basis that G × G ∗ = I, the identity matrix. Thus on any
tangent plane Tp , we get a Mahalanobis distance:v
uD DuMG (u, v) := ku − vkG = t ∑ ∑ Gij (ui − vi )(u j − v j ).

(10)

i =1 j =1

The inner product of two vectors u and v is a scalar (a 0-tensor) that can be equivalently
calculated as:ΣΣhu, vi := g(u, v) = ui vi = ui vi .
(11)
A metric tensor g of manifold M is said conformal when h·, ·i p = κ ( p)h·, ·iEuclidean . That is,when the inner product is a scalar function κ (·) of the Euclidean dot product. More precisely, we define
the notion of a metric g0 conformal to another metric g when these metrics define the same angles
between vectors u and v of a tangent plane Tp :q

g0p (u, v)
g p (u, v)
q
q
= q
.
g0p (u, u) g0p (v, v)
g p (u, u) g p (v, v)(12)Usually g0 is chosen as the Euclidean metric. In conformal geometry, we can measure angles between
vectors in tangent planes as if we were in an Euclidean space, without any deformation. This is handy
for checking orthogonality in charts. For example, the Poincaré disk model of hyperbolic geometry is
conformal but Klein disk model is not conformal (except at the origin), see [20].
2.3. Affine Connections ∇An affine connection ∇ is a differential operator defined on a manifold that allows us to define(1) a covariant derivative of vector fields, (2) a parallel transport of vectors on tangent planes along asmooth curve, and (3) geodesics. Furthermore, an affine connection fully characterizes the curvatureand torsion of a manifold.
2.3.1. Covariant Derivatives ∇ X Y of Vector Fields
A connection defines a covariant derivative operator that tells us how to differentiate a vector
field Y according to another vector field X. The covariant derivative operator is denoted using the
traditional gradient symbol ∇. Thus a covariate derivative ∇ is a function:∇ : X( M ) × X( M ) → X( M ),(13)

that has its own special subscript notation ∇ X Y:=:∇( X, Y ) for indicating that it is differentiating avector field Y according to another vector field X.
By prescribing D3 smooth functions Γijk = Γijk ( p), called the Christoffel symbols of the second
kind, we define the unique affine connection ∇ that satisfies in local coordinates of chart (U , x ) the
following equations:∇∂i ∂ j = Γijk ∂k .
(14)

Entropy 2020, 22, 1100

7 of 61

The Christoffel symbols can also be written as Γijk := (∇∂i ∂ j )k , where (·)k denotes the k-th
coordinate. The k-th component (∇ X Y )k of the covariant derivative of vector field Y with respect tovector field X is given by:
k Σ

k Σ

(∇ X Y ) = X i (∇i Y ) = X i∂Y k+ Γijk Y j
∂xi

!
.

(15)

The Christoffel symbols are not tensors (fields) because the transformation rules induced by achange of basis do not obey the tensor contravariant/covariant rules.
2.3.2. Parallel Transport ∏∇c along a Smooth Curve c
Since the manifold is not embedded in a Euclidean space, we cannot add a vector v ∈ Tp toa vector v0 ∈ Tp0 as the tangent vector spaces are unrelated to each others without a connection(the Whitney embedding theorem [21] states that any D-dimensional Riemannian manifold can beembedded into R2D ; when embedded, we can implicitly use the ambient Euclidean connection Euc ∇on the manifold, see [22]). Thus a connection ∇ defines how to associate vectors between infinitesimallyclose tangent planes Tp and Tp+dp . Then the connection allows us to smoothly transport a vectorv ∈ Tp by sliding it (with infinitesimal moves) along a smooth curve c(t) (with c(0) = p and c(1) = q),
so that the vector v p ∈ Tp “corresponds” to a vector vq ∈ Tq : this is called the parallel transport. This
mathematical prescription is necessary in order to study dynamics on manifolds (e.g., study the motionof a particle on the manifold). We can express the parallel transport along the smooth curve c as:∇∀v ∈ Tp , ∀t ∈ [0, 1],

vc(t) =∏v ∈ Tc(t)(16)c(0)→c(t)The parallel transport is schematically illustrated in Figure 3.

vc(t)
vp

vq

q
c(t)
Q∇
v
=
q
c vp
p

M
Figure 3. Illustration of the parallel transport of vectors on tangent planes along a smooth curve. For asmooth curve c, with c(0) = p and c(1) = q, a vector v p ∈ Tp is parallel transported smoothly to avector vq ∈ Tq such that for any t ∈ [0, 1], we have vc(t) ∈ Tc(t) .

Elie Cartan introduced the notion of affine connections [23,24] in the 1920s motivated by theprinciple of inertia in mechanics: a point particle, without any force acting on it, shall move along a
straight line with constant velocity.
2.3.3. ∇-Geodesics γ∇ : Autoparallel Curves
A connection ∇ allows one to define ∇-geodesics as autoparallel curves, that are curves γ suchthat we have:∇γ̇ γ̇ = 0.
(17)
That is, the velocity vector γ̇ is moving along the curve parallel to itself (and all tangent
vectors on the geodesics are mutually parallel): In other words, ∇-geodesics generalize the notion ofEntropy 2020, 22, 1100

8 of 61“straight Euclidean” lines. In local coordinates (U , x ), γ(t) = (γk (t))k , the autoparallelism amounts tosolve the following second-order Ordinary Differential Equations (ODEs):γ̈(t) + Γijk γ̇(t)γ̇(t) = 0,

γ l ( t ) = x l ◦ γ ( t ),(18)

where Γijk are the Christoffel symbols of the second kind, with:Σ

Γijk = Γij,l glk ,ΣΓij,k = glk Γijl ,(19)

where Γij,l the Christoffel symbols of the first kind. Geodesics are 1D autoparallel submanifolds and∇-hyperplanes are defined similarly as autoparallel submanifolds of dimension D − 1. We may specify
in subscript the connection that yields the geodesic γ: γ∇ .
The geodesic equation ∇γ̇(t) γ̇(t) = 0 may be either solved as an Initial Value Problem (IVP) or asa Boundary Value Problem (BVP):••

Initial Value Problem (IVP): fix the conditions γ(0) = p and γ̇(0) = v for some vector v ∈ Tp .
Boundary Value Problem (BVP): fix the geodesic extremities γ(0) = p and γ(1) = q.

2.3.4. Curvature and Torsion of a Manifold
An affine connection ∇ defines a 4D curvature tensor R (expressed using components Rijkl of a(1, 3)-tensor). The coordinate-free equation of the curvature tensor is given by:
R( X, Y ) Z := ∇ X ∇Y X − ∇Y ∇ X Z − ∇[ X,Y ] Z,

(20)

where [ X, Y ]( f ) = X (Y ( f )) − Y ( X ( f )) (∀ f ∈ F( M )) is the Lie bracket of vector fields. When theconnection is the metric Levi–Civita, the curvature is called Riemann–Christoffel curvature tensor. In a
local coordinate system, we have:ΣR(∂ j , ∂k )∂i = Rljki ∂l .

(21)

Informally speaking, the curvature tensor as defined in Equation (20) quantifies the amount of
non-commutativity of the covariant derivative. It follows from symmetry constraints that the number
D 2 ( D 2 −1)of independent components of the Riemann tensor isin D dimensions.
12
A manifold M equipped with a connection ∇ is said flat (meaning ∇-flat) when R = 0. This holds
in particular when finding a particular coordinate system x of a chart (U , x ) such that Γijk = 0, i.e., when
all connection coefficients vanish. For example, the Christoffel symbols vanish in a rectangularcoordinate system of a plane but not in the polar coordinate system of it.
A manifold is torsion-free when the connection is symmetric. A symmetric connection satisfies
the following coordinate-free equation:∇ X Y − ∇Y X = [ X, Y ].

(22)

Using local chart coordinates, this amounts to check that Γijk = Γkji . The torsion tensor is a (1, 2)-tensor
defined by:
T ( X, Y ) := ∇ X Y − ∇Y X − [ X, Y ].
(23)
For a torsion-free connection, we have the first Bianchi identity:
R( X, Y ) Z + R( Z, X )Y + R(Y, Z ) X = 0,(24)

and the second Bianchi identity:(∇V R)( X, Y ) Z + (∇ X R)(Y, V ) Z + (∇Y R)(V, X ) Z = 0.

(25)

Entropy 2020, 22, 1100

9 of 61

In general, the parallel transport is path-dependent. The angle defect of a vector transported on
an infinitesimal closed loop (a smooth curve with coinciding extremities) is related to the curvature.
However for a flat connection, the parallel transport does not depend on the path, and yields absoluteparallelism geometry [25]. Figure 4 illustrates the parallel transport along a loop curve for a curved
manifold (the sphere manifold) and a flat manifold (the cylinder manifold).

Figure 4. Parallel transport with respect to the metric connection: the curvature effect can be visualized
as the angle defect along the parallel transport on smooth (infinitesimal) loops. For a sphere manifold,a vector parallel-transported along a loop does not coincide with itself (e.g., a sphere), while it always
conside with itself for a (flat) manifold (e.g., a cylinder).

Historically, the Gaussian curvature at of point of a manifold has been defined as the
product of the minimal and maximal sectional curvatures: κ G := κmin κmax . For a cylinder, sinceκmin = 0, it follows that the Gaussian curvature of a cylinder is 0. Gauss’s Theorema Egregium(meaning “remarkable theorem”) proved that the Gaussian curvature is intrinsic and does not dependon how the surface is embedded into the ambient Euclidean space.
An affine connection is a torsion-free linear connection. Figure 5 summarizes the various conceptsof differential geometry induced by an affine connection ∇ and a metric tensor g.
Affine connection∇
Curvature
∇ i
R jkl

ParallelQtransport∇c(t) v

Geodesic∇γ(∇γ̇ γ̇ = 0)

Volume form∇
ω = f dv(∇ ∇ ω = 0)

Covariantderivative∇X Y

Levi-civita connectiong∇
Ricci curvature
Metric tensor∇
Ric(Y, Z) = tr(X 7→ ∇ R(X, Y, Z))g

Scalar curvature
Scal = g ij Rij
Manifold (M, g, ∇)Divergence
div(X) = tr(Y 7→ g(∇Y X, Y ))Gradientg(gradf, X) = df (X)HessianHessf (x)[v] = ∇v gradf

Laplacian∆f = div(grad(f ))Figure 5. Differential-geometric concepts associated to an affine connection ∇ and a metric tensor g.

Entropy 2020, 22, 1100

10 of 61

Curvature is a fundamental concept inherent to geometry [26]: there are several notions ofcurvatures in differential geometry: scalar curvature, sectional curvature, Gaussian curvature ofsurfaces to Riemannian–Christoffel 4-tensor, Ricci symmetric 2-tensor, synthetic Ricci curvature in
Alexandrov geometry, etc.
For example, the real-valued Gaussian curvature sec Tp M on a 2D Riemannian manifold ( M, g)with Riemann curvature (1, 3)-tensor R at a point p (with local basis {∂1 , ∂2 } on its tangent plane Tp M)
is defined by:

sec Tp M



=
=

R2112
,
2
g11 g22 − g12∇ ∂2 ∇ ∂1 ∂ 1 − ∇ ∂1 ∇ ∂2 ∂ 1 , ∂ 2 pdet( g)

(26)
.

(27)

In general, the sectional curvatures are real values defined for 2-dimensional subspaces π p of thetangent plane Tp M (called tangent 2-planes) as:
sec p (π ) :=h R( X, Y )Y, X i p,
Q p ( X, Y )(28)

where X and Y are linearly independent vectors of Tp , andQ p ( X, Y ) := h X, X i p hY, Y i p − h X, Y i2p ,(29)denotes the squared area of the parallelogram spanned by vectors X and Y of Tp M. It can be shownthat sec p (π ) is independent of the chosen basis X and Y. In a local basis {∂i }i of D-dimensionaltangent plane Tp , we thus get the sectional curvatures at point p ∈ M as the following real values:κij := sec p span ∂i , ∂ j ,i 6= j.

(30)

A Riemannian manifold ( M, g) is said of constant curvature κ if and only if sec p (π ) = κ for allp ∈ M and π p ⊂ Tp M. In particular, the Riemannian manifold is said flat when it is of constantcurvature 0. Notice that the definition of sectional curvatures relies on the metric tensor g but the
Riemann–Christoffel curvature tensor is defined with respect to an affine connection (which can betaken as the default Levi–Civita metric connection induced by the metric g).
2.4. The Fundamental Theorem of Riemannian Geometry: The Levi–Civita Metric Connection
By definition, an affine connection ∇ is said metric compatible with g when it satisfies for any
triple ( X, Y, Z ) of vector fields the following equation:
X hY, Z i = h∇ X Y, Z i + hY, ∇ X Z i,(31)

which can be written equivalently as:
Xg(Y, Z ) = g(∇ X Y, Z ) + g(Y, ∇ X Z )(32)

Using local coordinates and natural basis {∂i } for vector fields, the metric-compatibility property
amounts to check that we have:∂k gij = h∇∂k ∂i , ∂ j i + h∂i , ∇∂k ∂ j i(33)

Entropy 2020, 22, 1100

11 of 61

A property of using a metric-compatible connection is that the parallel transport ∏∇ of vectors
preserve the metric:
*
+
∇

∇

∏

hu, vic(0) =∏u,c(0)→c(t)∀t.

v

c(0)→c(t)

(34)

c(t)

That is, the parallel transport preserves angles (and orthogonality) and lengths of vectors in tangent
planes when transported along a smooth curve.
The fundamental theorem of Riemannian geometry states the existence of a unique torsion-free
metric compatible connection:
Theorem 1 (Levi–Civita metric connection). There exists a unique torsion-free affine connection compatible
with the metric called the Levi–Civita connection LC ∇.
The Christoffel symbols of the Levi–Civita connection can be expressed from the metric tensor g
as follows:LC k Σ 1 kl
Γij = g ∂i gil + ∂ j gil − ∂l gij ,(35)
2
where gij denote the matrix elements of the inverse matrix g−1 .
The Levi–Civita connection can also be defined coordinate-free with the Koszul formula:2g(∇ X Y, Z ) = X ( g(Y, Z )) + Y ( g( X, Z )) − Z ( g( X, Y )) + g([ X, Y ], Z ) − g([ X, Z ], Y ) − g([Y, Z ], X ).

(36)

There exists metric-compatible connections with torsions studied in theoretical physics. See forexample the flat Weitzenböck connection [27].
The metric tensor g induces the torsion-free metric-compatible Levi–Civita connection that
determines the local structure of the manifold. However, the metric g does not fix the global topologicalstructure: For example, although a cone and a cylinder have locally the same flat Euclidean metric,
they exhibit different global structures.
2.5. Preview: Information Geometry versus Riemannian Geometry
In information geometry, we consider a pair of conjugate affine connections ∇ and ∇∗ (often but
not necessarily torsion-free) that are coupled to the metric g: the structure is conventionally written
as ( M, g, ∇, ∇∗ ). The key property is that those conjugate connections are metric compatible,
and therefore the induced dual parallel transport preserves the metric:*hu, vic(0) =

∇∗

∇

∏

c(0)→c(t)u,∏c(0)→c(t)

+
v

.

(37)

c(t)

Thus the Riemannian manifold ( M, g) can be interpreted as the self-dual information-geometricmanifold obtained for ∇ = ∇∗ = LC ∇ the unique torsion-free Levi–Civita metric connection: ( M, g) ≡∗( M, g, LC ∇, LC ∇ = LC ∇). However, let us point out that for a pair of self-dual Levi–Civita conjugate
connections, the information-geometric manifold does not induce a distance. This contrasts with the
Riemannian modeling ( M, g) which provides a Riemmanian metric distance Dρ ( p, q) defined by thelength of the geodesic γ connecting the two points p = γ(0) and q = γ(1):Dρ ( p, q)

:=

=

Z 1

Z 1q
kγ0 (t)kγ(t) dt =gγ(t) (γ̇(t), γ̇(t))dt,00
Z 1qγ̇(t)> gγ(t) γ̇(t)dt.
0

(38)
(39)

Entropy 2020, 22, 1100

12 of 61

This geodesic length distance Dρ ( p, q) can also be interpreted as the shortest path linking point p to
R1
point q: Dρ ( p, q) = infγ 0 kγ0 (t)kγ(t) dt (with p = γ(0) and q = γ(1)).
Usually, this Riemannian geodesic distance is not available in closed-form (and need to beapproximated or bounded) because the geodesics cannot be explicitly parameterized (see geodesic
shooting methods [28]).
We are now ready to introduce the key geometric structures of information geometry.
3. Information Manifolds
3.1. Overview
In this part, we explain the dualistic structures of manifolds in information geometry. In Section 3.2,
we first present the core Conjugate Connection Manifolds (CCMs) ( M, g, ∇, ∇∗ ), and show how tobuild Statistical Manifolds (SMs) ( M, g, C ) from a CCM in Section 3.3. From any statistical manifold,we can build a 1-parameter family ( M, g, ∇−α , ∇α ) of CCMs, the information α-manifolds. We state thefundamental theorem of information geometry in Section 3.5. These CCMs and SMs structures are notrelated to any distance a priori but require at first a pair (∇, ∇∗ ) of conjugate connections coupled to a
metric tensor g. We show two methods to build an initial pair of conjugate connections. A first method
consists of building a pair of conjugate connections ( D ∇, D ∇∗ ) from any divergence D in Section 3.6.
Thus we obtain self-conjugate connections when the divergence is symmetric: D (θ1 : θ2 ) = D (θ2 : θ1 ).
When the divergences are Bregman divergences (i.e., D = BF for a strictly convex and differentiable
Bregman generator), we obtain Dually Flat Manifolds (DFMs) ( M, ∇2 F, F ∇, F ∇∗ ) in Section 3.7.
DFMs nicely generalize the Euclidean geometry and exhibit Pythagorean theorems. We further
characterize when orthogonal F ∇-projections and dual F ∇∗ -projections of a point on submanifold
a is unique. In Euclidean geometry, the orthogonal projection of a point p onto an affine subspace S
is proved to be unique using the Pythagorean theorem. A second method to get a pair of conjugate
connections (e ∇, m ∇) consists of defining these connections from a regular parametric family of
probability distributions P = { pθ ( x )}θ . In that case, these ‘e’xponential connection e ∇ and ‘m’ixture
connection m ∇ are coupled to the Fisher information metric P g. A statistical manifold (P , P g, P C )can be recovered by considering the skewness Amari–Chentsov cubic tensor P C, and it follows
a 1-parameter family of CCMs, (P , P g, P ∇−α , P ∇+α ), the statistical expected α-manifolds. In thisparametric statistical context, these information manifolds are called expected information manifolds
because the various quantities are expressed from statistical expectations E· [·]. Notice that theseinformation manifolds can be used in information sciences in general, beyond the traditional fieldsof statistics. In statistics, we motivate the choice of the connections, metric tensors and divergencesby studying statistical invariance criteria, in Section 3.10. We explain how to recover the expected
α-connections from standard f -divergences that are the only separable divergences that satisfy theproperty of information monotonicity. Finally, in Section 3.11, the recall the Fisher–Rao expected
Riemannian manifolds that are Riemannian manifolds (P , P g) equipped with a geodesic metricdistance called the Fisher–Rao distance, or Rao distance for short.
3.2. Conjugate Connection Manifolds: ( M, g, ∇, ∇∗ )We begin with a definition:
Definition 1 (Conjugate connections). A connection ∇∗ is said to be conjugate to a connection ∇ withrespect to the metric tensor g if and only if we have for any triple ( X, Y, Z ) of smooth vector fields the following
identity satisfied:
X hY, Z i = h∇ X Y, Z i + hY, ∇∗X Z i, ∀ X, Y, Z ∈ X( M ).
(40)

Entropy 2020, 22, 1100

13 of 61

We can notationally rewrite Equation (40) as:
Xg(Y, Z ) = g(∇ X Y, Z ) + g(Y, ∇∗X Z ),(41)

and further explicit that for each point p ∈ M, we have:X p g p (Yp , Z p ) = g p ((∇ X Y ) p , Z p ) + g p (Yp , (∇∗X Z ) p ).

(42)

We check that the right-hand-side is a scalar and that the left-hand-side is a directional derivative of areal-valued function, that is also a scalar.
Conjugation is an involution: (∇∗ )∗ = ∇.
Definition 2 (Conjugate Connection Manifold). The structure of the Conjugate Connection Manifold(CCM) is denoted by ( M, g, ∇, ∇∗ ), where (∇, ∇∗ ) are conjugate connections with respect to the metric g.
A remarkable property is that the dual parallel transport of vectors preserves the metric. That is,for any smooth curve c(t), the inner product is conserved when we transport one of the vector u using∇∗
the primal parallel transport ∏∇c and the other vector v using the dual parallel transport ∏c .
*+∗∇hu, vic(0) =

∏

∇

u,

c(0)→c(t)∏c(0)→c(t)

v

.

(43)

c(t)

Property 1 (Dual parallel transport preserves the metric). A pair (∇, ∇∗ ) of conjugate connections
preserves the metric g if and only if:*∀t ∈ [0, 1],

∇∗

∇

∏

c(0)→c(t)u,∏c(0)→c(t)+= hu, vic(0) .

v

(44)

c(t)

Property 2. Given a connection ∇ on ( M, g) (i.e., a structure ( M, g, ∇)), there exists a unique conjugate
connection ∇∗ (i.e., a dual structure ( M, g, ∇∗ )).
We consider a manifold M equipped with a pair of conjugate connections ∇ and ∇∗ that arecoupled with the metric tensor g so that the dual parallel transport preserves the metric. We define the¯:mean connection ∇∗¯ = ∇+∇ ,(45)∇
2with corresponding Christoffel coefficients denoted by Γ̄. This mean connection coincides with theLevi–Civita metric connection:¯ = LC ∇.
∇(46)¯ is self-conjugate, and coincide with the Levi–Civita metric connection.
Property 3. The mean connection ∇3.3. Statistical Manifolds: ( M, g, C )Lauritzen introduced this corner structure [29] of information geometry in 1987. Beware that
although it bears the name “statistical manifold”, it is a purely geometric construction that may be
used outside of the field of Statistics. However, as we shall mention later, we can always find astatistical model P corresponding to a statistical manifold [30]. We shall see how we can convert a
conjugate connection manifold into such a statistical manifold, and how we can subsequently derive
an infinite family of CCMs from a statistical manifold. In other words, once we have a pair of conjugate
connections, we will be able to build a family of pairs of conjugate connections.

Entropy 2020, 22, 1100

14 of 61

We define a cubic (0, 3)-tensor (i.e., 3-covariant tensor) called the Amari–Chentsov tensor:
Cijk := Γijk − Γ∗ ijk ,(47)

C ( X, Y, Z ) := h∇ X Y − ∇∗X Y, Z i.

(48)

or in coordinate-free equation:

The cubic tensor is totally symmetric, meaning that Cijk = Cσ(i)σ( j)σ(k) for any permutation σ.
The metric tensor is totally symmetric.
Using the local basis, this cubic tensor can be expressed as:
Cijk = C (∂i , ∂ j , ∂k ) = h∇∂i ∂ j − ∇∗∂i ∂ j , ∂k i

(49)

Definition 3 (Statistical manifold [29]). A statistical manifold ( M, g, C ) is a manifold M equipped with ametric tensor g and a totally symmetric cubic tensor C.
3.4. A Family {( M, g, ∇−α , ∇α = (∇−α )∗ )}α∈R of Conjugate Connection Manifolds
For any pair (∇, ∇∗ ) of conjugate connections, we can define a 1-parameter family of connections{∇α }α∈R , called the α-connections such that (∇−α , ∇α ) are dually coupled to the metric, with ∇0 =¯ = LC ∇, ∇1 = ∇ and ∇−1 = ∇∗ . By observing that the scaled cubic tensor αC is also a totally∇symmetric cubic 3-covariant tensor, we can derive the α-connections from a statistical manifold( M, g, C ) as:
α
Γij,k

−α
Γij,k

α
= Γ0ij,k − Cij,k ,2α0= Γij,k + Cij,k ,
2

(50)
(51)

Σ

where Γ0ij,k are the Levi–Civita Christoffel symbols, and Γki,j = Γijl glk (by index juggling).
The α-connection ∇α can also be defined as follows:g(∇αX Y, Z ) = g(LC ∇ X Y, Z ) +αC ( X, Y, Z ), ∀ X, Y, Z ∈ X( M ).
2

(52)

Theorem 2 (Family of information α-manifolds). For any α ∈ R, ( M, g, ∇−α , ∇α = (∇−α )∗ ) is aconjugate connection manifold.
The α-connections ∇α can also be constructed directly from a pair (∇, ∇∗ ) of conjugate
connections by taking the following weighted combination:
α
Γij,k
=

1+α
1−α ∗
Γij,k +
Γij,k .
2
2

(53)

3.5. The Fundamental Theorem of Information Geometry: ∇ κ-Curved ⇔ ∇∗ κ-Curved
We now state the fundamental theorem of information geometry and its corollaries:
Theorem 3 (Dually constant curvature manifolds). If a torsion-free affine connection ∇ has constantcurvature κ then its conjugate torsion-free connection ∇∗ has necessarily the same constant curvature κ.
The proof is reported in [16] (Proposition 8.1.4, page 226).
A statistical manifold ( M, g, C ) is said α-flat if its induced α-connection is flat. It can be shownthat Rα = − R−α .
We get the following two corollaries:

Entropy 2020, 22, 1100

15 of 61

Corollary 1 (Dually α-flat manifolds). A manifold ( M, g, ∇−α , ∇α ) is ∇α -flat if and only if it is ∇−α -flat.
Corollary 2 (Dually flat manifolds (α = ±1)). A manifold ( M, g, ∇, ∇∗ ) is ∇-flat if and only if it is ∇∗ -flat.
Refer to Theorem 3.3 of [4] for a proof of this corollary.
Let us now define the notion of constant curvature of a statistical structure [31]:Definition 4 (Constant curvature κ). A statistical structure ( M, g, ∇) is said of constant curvature κ whenR∇ ( X, Y ) Z = κ { g(Y, Z ) X − g( X, Z )Y },∀ X, Y, Z ∈ Γ( TM),

where Γ( TM ) denote the space of smooth vector fields.
It can be proved that the Riemann–Christoffel (RC) 4-tensors of conjugate α-connections [16] arerelated as follows:g R(α) ( X, Y ) Z, W + g Z, R(−α) ( X, Y )W = 0.
(54)
 ∗


We have g R∇ ( X, Y ) Z, W = − g Z, R∇ ( X, Y )W .
Thus once we are given a pair of conjugate connections, we can always build a 1-parametricfamily of manifolds. Manifolds with constant curvature κ are interesting from the computational
viewpoint as dual geodesics have simple closed-form expressions.
∗3.6. Conjugate Connections from Divergences: ( M, D ) ≡ ( M, D g, D ∇, D ∇∗ = D ∇)Loosely speaking, a divergence D (· : ·) is a smooth distance [32], potentially asymmetric. In orderto define precisely a divergence, let us first introduce the following handy notations: ∂i,· f ( x, y) =22∂
f ( x, y), ∂·,j f ( x, y) = ∂y∂ j f ( x, y), ∂ij,k f ( x, y) = ∂x∂i ∂x j ∂y∂ k f ( x, y) and ∂i,jk f ( x, y) = ∂x∂ i ∂y∂j ∂yk f ( x, y), etc.
∂xi
Definition 5 (Divergence). A divergence D : M × M → [0, ∞) on a manifold M with respect to a local chartΘ ⊂ RD is a C3 -function satisfying the following properties:
1.
2.
3.

D (θ : θ 0 ) ≥ 0 for all θ, θ 0 ∈ Θ with equality holding iff θ = θ 0 (law of the indiscernibles),∂i,· D (θ : θ 0 )|θ =θ 0 = ∂·,j D (θ : θ 0 )|θ =θ 0 = 0 for all i, j ∈ [ D ],−∂·,i ∂·,j D (θ : θ 0 )|θ =θ 0 is positive-definite.
The dual divergence is defined by swapping the arguments:
D ∗ ( θ : θ 0 ) : = D ( θ 0 : θ ),(55)

and is also called the reverse divergence (reference duality in information geometry). Reference dualityof divergences is an involution: ( D ∗ )∗ = D.
The Euclidean distance is a metric distance but not a divergence. The squared Euclidean distanceis a non-metric symmetric divergence. The metric tensor g yields Riemannian metric distance Dρ but it
is never a divergence.
From any given divergence D, we can define a conjugate connection manifold following theconstruction of Eguchi [33,34] (1983):∗

Theorem 4 (Manifold from divergence). ( M, D g, D ∇, D ∇) is an information manifold with:

D
D∗

∗

g

:=

−∂i,j D (θ : θ 0 )|θ =θ 0 = D g,

(56)

Γijk

:=

−∂ij,k D (θ : θ 0 )|θ =θ 0 ,

(57)

D

Γijk

:=

0

−∂k,ij D (θ : θ )|θ =θ 0 .

(58)

Entropy 2020, 22, 1100

16 of 61

The associated statistical manifold is ( M, D g, D C ) with:D∗

Cijk = D Γijk − D Γijk .

(59)

Since α D C is a totally symmetric cubic tensor for any α ∈ R, we can derive a one-parameter familyof conjugate connection manifolds:nα( M, D g, D C ) ≡ ( M, D g, D ∇−α, (D ∇−α ∗α) = D∇ )

o
α ∈R

.

(60)

In the remainder, we use the shortcut ( M, D ) to denote the divergence-induced information∗
manifold ( M, D g, D ∇, D ∇ ). Notice that it follows from construction that:D∗∗∇ = D ∇.

(61)

3.7. Dually Flat Manifolds (Bregman Geometry): ( M, F ) ≡ ( M, BF g, BF ∇, BF ∇∗ = BF∗ ∇)We consider dually flat manifolds that satisfy asymmetric Pythagorean theorems. These flat
manifolds can be obtained from a canonical Bregman divergence.
Consider a strictly convex smooth function F (θ ) called a potential function, with θ ∈ Θ where Θ isan open convex domain. Notice that the function convexity does not change by an affine transformation.
We associate to the potential function F a corresponding Bregman divergence (parameter divergence):
B F ( θ : θ 0 ) : = F ( θ ) − F ( θ 0 ) − ( θ − θ 0 ) > ∇ F ( θ 0 ).

(62)

We write also the Bregman divergence between point P and point Q as D ( P : Q) := BF (θ ( P) :θ ( Q)), where θ ( P) denotes the coordinates of a point P.
The information-geometric structure induced by a Bregman generator is ( M, F g, F C ) :=( M, BF g, BF C ) with:
F
F
F

g
Γ

Cijk

:=

BF



g = − ∂i ∂ j BF (θ : θ 0 )|θ 0 =θ = ∇2 F (θ ),

(63)

:=

BF

Γij,k (θ ) = 0,

(64)

:=

BF

Cijk = ∂i ∂ j ∂k F (θ ).

(65)

Here, we define a Bregman generator as a proper, lower semi-continuous, and strictly convex andC3 differentiable real-valued function.
Since all coefficients of the Christoffel symbols vanish (Equation (64)), the information manifold
is F ∇-flat. The Levi–Civita connection LC ∇ is obtained from the metric tensor F g (usually not flat),
and we get the conjugate connection ( F ∇)∗ = F ∇1 from ( M, F g, F C ).
The Legendre–Fenchel transformation yields the convex conjugate F ∗ that is interpreted as thedual potential function:F ∗ (η ) := sup{θ > η − F (θ )}.
(66)
θ ∈Θ

A function f is lower semicontinous (lsc) at x0 iff f ( x0 ) ≤ limx→ x0 inf f ( x ). A function f is lsc if it
is lsc at x for all x in the function domain. The following theorem states that the conjugation of lower
semicontinuous and convex functions is an involution:
Theorem 5 (Fenchel–Moreau biconjugation [35]). If F is a lower semicontinuous and convex function,then its Legendre–Fenchel transformation is involutive: ( F ∗ )∗ = F (biconjugation).
In a dually flat manifold, there exists two global dual affine coordinate systems η = ∇ F (θ ) andθ = ∇ F ∗ (η ), and therefore the manifold can be covered by a single chart. Thus if a probability family

Entropy 2020, 22, 1100

17 of 61

belongs to an exponential family then its natural parameters cannot belong to, say, a spherical space(that requires at least two charts).
We have the Crouzeix [36] identity relating the Hessians of the potential functions:∇2 F (θ )∇2 F ∗ (η ) = I,

(67)

where I denote the D × D identity matrix. This Crouzeix identity reveals that B = {∂i }i and B∗ = {∂ j } jare the primal and reciprocal basis, respectively.
The Bregman divergence can be reinterpreted using Young–Fenchel (in)equality as the canonicaldivergence A F,F∗ [37]:BF (θ : θ 0 ) = A F,F∗ (θ : η 0 ) = F (θ ) + F ∗ (η 0 ) − θ > η 0 = A F∗ ,F (η 0 : θ ).

(68)

The dual Bregman divergence BF ∗ (θ : θ 0 ) := BF (θ 0 : θ ) = BF∗ (η : η 0 ) yields
F ij

g (η )F ∗ ijkΓ= ∂ i ∂ j F ∗ ( η ),(η ) = 0,F∂l :=:

∂
∂η l

(69)

Cijk = ∂i ∂ j ∂k F ∗ (η )(70)

Thus the information manifold is both F ∇-flat and F ∇∗ -flat: This structure is called a dually flat
manifold (DFM). In a DFM, we have two global affine coordinate systems θ (·) and η (·) related by the
Legendre–Fenchel transformation of a pair of potential functions F and F ∗ . That is, ( M, F ) ≡ ( M, F ∗ ),and the dual atlases are A = {( M, θ )} and A∗ = {( M, η )}.
In a dually flat manifold, any pair of points P and Q can either be linked using the ∇-geodesic(that is θ-straight) or the ∇∗ -geodesic (that is η-straight). In general, there are 23 = 8 types of geodesic
triangles in a dually flat manifold.
On a Bregman manifold, the primal parallel transport of a vector does not change the contravariant
vector components, and the dual parallel transport does not change the covariant vector components.
Because the dual connections are flat, the dual parallel transports are path-independent.
Moreover, the dual Pythagorean theorems [38] illustrated in Figure 6 holds. Let γ( P, Q) =γ∇ ( P, Q) denote the ∇-geodesic passing through points P and Q, and γ∗ ( P, Q) = γ∇∗ ( P, Q) denotethe ∇∗ -geodesic passing through points P and Q. Two curves γ1 and γ2 are orthogonal at pointp = γ1 (t1 ) = γ2 (t2 ) with respect to the metric tensor g when g(γ̇1 (t1 ), γ̇2 (t2 )) = 0.
γ ∗ (P, Q) ⊥F γ(Q, R)

γ(P, Q) ⊥F γ ∗ (Q, R)

P

P

R
Q
D(P : R) = D(P : Q) + D(Q : R)
BF (θ(P ) : θ(R)) = BF (θ(P ) : θ(Q)) + BF (θ(Q) : θ(R))

R
Q
D∗ (P : R) = D∗ (P : Q) + D∗ (Q : R)
BF ∗ (η(P ) : η(R)) = BF ∗ (η(P ) : η(Q)) + BF ∗ (η(Q) : η(R))

Figure 6. Dual Pythagorean theorems in a dually flat space.

Entropy 2020, 22, 1100

18 of 61

Theorem 6 (Dual Pythagorean identities).
Σγ∗ ( P, Q) ⊥ γ( Q, R)⇔ (η ( P) − η ( Q))> (θ ( Q) − θ ( R)) = (ηi ( P) − ηi ( Q))(θi ( Q) − θi ( R)) = 0,γ( P, Q) ⊥ γ∗ ( Q, R)⇔ (θ ( P) − θ ( Q))> (η ( Q) − η ( R)) = (θi ( P) − θi ( Q))> (ηi ( Q) − ηi ( R)) = 0.

Σ

We can define dual Bregman projections and characterize when these projections are unique:A submanifold S ⊂ M is said ∇-flat (∇∗ -flat) iff. it corresponds to an affine subspace in the θ-coordinate
system (in the η-coordinate system, respectively).
Theorem 7 (Uniqueness of projections). The ∇-projection PS of P on S is unique if S is ∇∗ -flat andminimizes the divergence D (θ ( P) : θ ( Q)):∇-projection: PS = arg min D (θ ( P) : θ ( Q)).

(71)

Q∈S

The dual ∇∗ -projection PS∗ is unique if M ⊆ S is ∇-flat and minimizes the divergence D (θ ( Q) : θ ( P)):∇∗ -projection: PS∗ = arg min D (θ ( Q) : θ ( P)).

(72)

Q∈S

Let S ⊂ M and S0 ⊂ M, then we define the divergence between S and S0 as
D (S : S0 ) :=min D (s : s0 ).

(73)

s∈S,s0 ∈S0When S is a ∇-flat submanifold and S0 ∇∗ -flat submanifold, the divergence D (S : S0 ) between
submanifold S and submanifold S0 can be calculated using the method of alternating projections [2].
Let us remark that Kurose [39] reported a Pythagorean theorem for dually constant curvature manifolds
that generalizes the Pythagorean theorems of dually flat spaces.
We shall concisely explain the space of Bregman spheres explained in details in [40]. Let D denotethe dimension of Θ. We define the lifting of primal coordinates θ to the primal potential function
F = {θ̂ = (θ, θ D+1 = F (θ )) : θ ∈ Θ} using an extra dimension θ D+1 . A Bregman ball Σ
BallF (C : r ) := { P such that F (θ ( P)) + F ∗ (η (C )) − hθ ( P), η (C )i ≤ r }(74)

can then be lifted to F : Σ̂ = {θ̂ ( P) : P ∈ σ }. The boundary Bregman sphere σ = ∂Σ is lifted to∂Σ̂ = σ̂, and the lifted points are all supported by a supporting ( D + 1)-dimensional hyperplane (ofdimension D):
Hσ̂ : θ D+1 = hθ − θ (C ), η (C )i + F (θ (C )) + r.
(75)
Let Hσ̂− denotes the halfspaces bounded by Hσ̂ and containing θ̂ (C ) = (θ (C ), F (θ (C ))). A point P
belongs to a Bregman ball Σ iff θ (ˆP) ∈ Hσ̂− , see [40]. Reciprocally, a ( D + 1)-dimensional hyperplaneH : θ D+1 = hθ, ηa i + b cutting the potential function F yields a Bregman sphere σH of center Cwith θ (C ) = ∇ F ∗ (ηa ) and radius r = h∇ F ∗ (ηa ), ηa i − F (θ a ) + b = F ∗ (ηa ) + b, where θ a = ∇ F ∗ (ηa ).
It follows that the intersection of k Bregman balls is a ( D − k)-dimensional Bregman ball, and thata Bregman sphere can be defined by D + 1 points in general position since an hyperplane in theaugmented space is defined by D + 1 points. We can test whether a point P belongs to a Bregman ballwith bounding Bregman sphere passing through D + 1 points P1 , . . . , PD+1 or not by checking the signof a ( D + 2) × ( D + 2) determinant:1InBregmanBallF ( P1 , . . . , Pd+1 ; P) := sign  θ ( P1 )F (θ ( P1 ))

...
...
...

1
θ ( PD+1 )F (θ ( PD+1 ))1θ ( P)
.
F (θ ( P))(76)Entropy 2020, 22, 1100

19 of 61

We have: = −1
InBregmanBallF ( P1 , . . . , Pd+1 ; P) :=0
 = +1⇔⇔⇔P ∈ InBregmanBall◦F ( P1 , . . . , PD+1 ; P)P ∈ ∂InBregmanBallF ( P1 , . . . , PD+1 ; P)P 6∈ InBregmanBallF ( P1 , . . . , PD+1 ; P)

(77)

Similarly, a dual-type Bregman ball Σ∗ can be defined byBall∗F (C : r ) := { P such that F (θ (C )) + F ∗ (η ( P)) − hθ (C ), η ( P)i ≤ r },(78)

and be lifted to the dual potential function F ∗ . Notice that Ball∗F (C : r ) = BallF∗ (C : r ).
Figure 7 displays five concentric pairs of dual Itakura–Saito circles obtained for the separable Burgnegentropy generator F ( x, y) = − log( x ) − log(y) (with corresponding Bregman divergence the
Itakura–Saito divergence).

Figure 7. Five concentric pairs of dual Itakura–Saito circles.

Using the space of spheres, it is easy to design algorithms for calculating the union or intersectionof Bregman spheres [40], or data-structures for proximity queries [41] (relying on the radical axis of
two Bregman spheres). The Bregman spheres are considered for building Bregman Voronoi diagramsin [40,42].
The smallest enclosing Bregman ball [43,44] (SEBB) of a set of points P1 , . . . , Pn (with respective
θ-coordinates θ1 , . . . , θn ) can also be modeled as a convex program; indeed, point Pi belongs to the lower
halfspace H − of equation θ D+1 ≤ hηa , θ i + b (parameterized by vector ηa ∈ RD and scalar b ∈ R) iff
hηa , θ i + b ≥ F (θi ). Thus we seek to minimize minηa ,b r = F ∗ (ηa ) + b such that hθi , ηa i + b − F (θi ) ≥ 0
for all i ∈ {1, . . . , n}. This is a convex program since F ∗ is the convex conjugate of a convex generator F.
When F (θ ) = 12 θ > θ (i.e., Euclidean geometry), we recover the fact that the smallest enclosing ball of apoint set in Euclidean geometry can be solved using quadratic programming [45]. Faster approximation
algorithms for the smallest enclosing Bregman ball can be built based on core-sets [43].
In general, we have the following quadrilateral relation for Bregman divergences:
Property 4 (Bregman 4-parameter property [46]). For any four points P1 , P2 , Q1 , Q2 , we have thefollowing identity:

Entropy 2020, 22, 1100

20 of 61

BF (θ ( P1 ) : θ ( Q1 )) + BF (θ ( P2 ) : θ ( Q2 ))− BF (θ ( P1 ) : θ ( Q2 )) − BF (θ ( P2 ) : θ ( Q1 ))−(θ ( P2 ) − θ ( P1 ))> (η ( Q1 ) − η ( Q2 )) = 0.

(79)

In summary, to define a dually flat space, we need a convex Bregman generator. When the
α-geometries are neither dually flat (e.g., Cauchy manifolds [47], we may still build a dually flat
structure on the manifold by considering some Bregman generator (e.g., Bregman- -Tsallis generator
for the dually flat Cauchy manifold [47]). The dually flat geometry can be investigated under the wider
scope of Hessian manifolds [48] which consider locally potential functions. In general, a dually flatspace can be built from any smooth strictly convex generator F. For example, a dually flat geometry canbe built on homogeneous cones with the characteristic function F of the cone [48]. Figure 8 illustrates
several common constructions of dually flat spaces.
Exponential family
F : cumulant function
strictly convex anddifferentiable C 3 function

Homogeneous convex cone
F : characteristic function
Mixture family
Shannon negentropy

Dually flat space(Hessian structure)(M, g, ∇, ∇∗ )Bregman divergence

Figure 8. Common dually flat spaces associated to smooth and strictly convex generators.

3.8. Hessian α-Geometry: ( M, F, α) ≡ ( M, F g, F ∇−α , F ∇α )The dually flat manifold is also called a manifold with a Hessian structure [48] induced by a convex∗
potential function F. Since we built two dual affine connections BF ∇ = F ∇ and BF ∇∗ = F ∇∗ = F ∇,we can build a family of α-geometry as follows:F

andgij (θ ) = ∂i ∂ j F (θ ),1−αF αΓijk (θ ) =∂ i ∂ j ∂ k F ( θ ),2

F ij

g ( η ) = ∂ i ∂ j F ( η ),

(80)

∗ α
1+α i j k ∗F α ∗Γijk (η ) = F Γijk(η ) =∂ ∂ ∂ F ( η ).
2

(81)
∗

−1
Thus when α = ±1, the Hessian α-geometry is dually flat since F Γ1ijk (θ ) = 0 and F Γijk(η ) = 0.
We now consider information manifolds induced by parametric statistical models.

3.9. Expected α-Manifolds of a Family of Parametric Probability Distributions: (P , P g, P ∇−α , P ∇α )Informally speaking, an expected manifold is an information manifold built on a regularparametric family of distributions. It is sometimes called “expected” manifold or “expected” geometry
in the literature [49] because the components of the metric tensor g and the Amari–Chentsov cubic
tensor C are expressed using statistical expectations.
Let P be a parametric family of probability distributions:

P := { pθ ( x )}θ ∈Θ ,(82)

with θ belonging to the open parameter space Θ. The order of the family is the dimension of its
parameter space. We define the likelihood function L(θ; x ) := pθ ( x ) as a function of θ, and its
corresponding log-likelihood function:l (θ; x ) := log L(θ; x ) = log pθ ( x ).

(83)

More precisely, the likelihood function is an equivalence class of functions defined modulo apositive scaling factor.

Entropy 2020, 22, 1100

21 of 61

The score vector:s θ = ∇ θ l = ( ∂i l )i ,(84)

indicates the sensitivity of the likelihood ∂i l:=: ∂θ∂ l (θ; x ).
i
The Fisher information matrix (FIM) of D × D for dim(Θ) = D is defined by:
P I ( θ ) : = Eθ∂i l∂ j l ij  0,(85)

where  denotes the Löwner order. That is, for two symmetric positive-definite matrices A and B,
A  B if and only if matrix A − B is positive semidefinite. For regular models [16], the FIM is positive
definite: P I (θ )  0, where A  B if and only if matrix A − B is positive-definite.
The FIM is invariant by reparameterization of the sample space X , and covariant by
reparameterization of the parameter space Θ, see [16]. That is, let p̄( x; η ) = p(θ (η ); x ). Then we have:"Ī (η ) =Matrix Jij =

h

∂θi
∂η j

#>
ij

"

∂θi
× I (θ (η )) ×
∂η j

#
.

(86)

ij

i

∂θi
∂η j ij is the Jacobian matrix.

Let us give illustrate the covariance of the Fisher information matrix with the following example:
Example 1. Consider the familyN =( x − µ )2)) : (µ, σ) ∈ R × R++p( x; µ, σ ) = √
exp(−
2σ2
2πσ
1


(87)

of univariate normal distributions. The 2D parameter vector is λ = (µ, σ ) with µ denoting the mean andσ the standard deviation. Another common parameterization of the normal family is λ0 = (µ, σ2 ). The λ0
parameterization extends naturally to d-variance normal distributions with λ0 = (µ, Σ), where Σ denotes thecovariance matrix (with Σ = σ2 when d = 1). For multivariate normal distributions, the λ-parameterization
can be interpreted as λ = (µ, L> ) where L> is the upper triangular matrix in the Cholesky decomposition(when d = 1, L> = σ). We have the following Fisher information matrices in the λ-parameterization and
λ0 -parameterization:
 "

#
1
1
0
2
0
λ
2
σ
=
(88)
Iλ (λ) =  2
0 λ22
0 σ22
2

and

"
Iλ0 λ

0



=

1
λ2

0

#

0

"

=

1
2λ22

1
σ2

0

#

0

(89)

1
2σ4

Since the FIM is covariant, we have the following the change of transformation:>0Iλ0 λ0 = Jλ,λ× Jλ,λ0 ,0 × Iλ λ λwith"
Jλ0 ,λ =

1
0

0
2σ

(90)

#
(91)

Thus we check that"Iλ (λ) =

1
0

0
2σ

#"

1
σ2

0

0
1
2σ4

#"

1
0

0
2σ

#

"

=

Notice that the infinitesimal length elements are invariant: dsλ = dsλ0 .

1
σ2

0

0
2
σ2

#
(92)

Entropy 2020, 22, 1100

22 of 61

As a corollary, notice that we can recognize the Euclidean metric in any other coordinate system> J 0 . For example, the Riemannian geometry induced by aif the metric tensor g can be written Jλ,λ0 λ,λdually flat space with a separable potential function is Euclidean [50].
In statistics, the FIM plays a role in the attainable precision of unbiased estimators. For any
unbiased estimator, the Cramér–Rao lower bound [51] on the variance of the estimator is:
Varθ [θ̂n ( X )] 

1 −1
I ( θ ).
nP

(93)

Figure 9 illustrates the Cramér–Rao lower bound (CRLB) for the univariate distributions:
At regular grid locations (µ, σ ) of the upper space of normal parameters, we repeat 200 runs (trials)
of estimating the normal parameters (\µ, σ ) using the MLE on 100 iid samples x1 , . . . , xn ∼ N (µ, σ ).
The sample mean and the sample covariance matrix are calculated for the number of trials anddisplayed as back ellipses. The Fisher information matrix is plotted as red ellipses at the grid locations:the red ellipses have semi-axes parallel to the coordinate system since the parameters µ and σ are
orthogonal (diagonal FIM). This is not true anymore for the sample covariance matrix of the MLE
estimator, and the centers of the sample covariance matrices deviate from the grid locations.

Figure 9. Visualizing the Cramér–Rao lower bound: the red ellipses display the Fisher information
matrix of normal distributions N (µ, σ2 ) at grid locations. The black ellipses are sample covariancematrices centered at the sample means calculated by repeating 200 runs of sampling 100 iid variates forthe normal parameters of the grid.

We report the expression of the FIM for two important generic parametric family of probability
distributions: (1) an exponential family (with its prominent multivariate normal family), and (2) amixture family.

Entropy 2020, 22, 1100

23 of 61

Example 2 (FIM of an exponential family E ). An exponential family [52] E is defined for a sufficient statistic
vector t( x ) = (t1 ( x ), . . . , t D ( x )), and an auxiliary carrier measure k ( x ) by the following canonical density:
(

E=

D

pθ ( x ) = exp∑ ti ( x ) θi − F ( θ ) + k ( x )!

)such that θ ∈ Θ

,

(94)

i =1

where F is the strictly convex cumulant function (also called log-normalizer, and log partition function or freeenergy in statistical mechanics). Exponential families include the Gaussian family, the Gamma and Beta families,
the probability simplex ∆, etc. The FIM of an exponential family is given by:
2
2 ∗
−1
 0.
E I ( θ ) = CovX ∼ pθ ( x ) [ t ( x )] = ∇ F ( θ ) = (∇ F ( η ))(95)

Indeed, under mild conditions [2], we have I (θ ) = − E pθ [∇2 log pθ ( x )]. Since −∇2 log pθ ( x ) = ∇2 F (θ ),it follows that E I (θ ) = ∇2 F (θ ). Natural parameters beyond vector types can also be used in the canonicaldecomposition of the density of an exponential family. For example, we may use a matrix type for defining
the zero-centered multivariate Gaussian family or the Wishart family, a complex numbers for defining thecomplex-valued Gaussian distribution family, etc. We then replace the term ∑iD=1 ti ( x )θi in Equation (94) byan inner product defined for the natural parameter type (e.g., dot product for vectors, matrix product trace
for matrices, etc.). Furthermore, natural parameters can be of compound types: For example, the multivariate
Gaussian distribution can be written using θ = (θv , θ M ) where θv is a vector part and θ M a matrix part, see [52].
Let Σ = [σij ] denote the covariance matrix and Σ−1 = [σij ] the precision matrix of a multivariate normal
distribution. The Fisher information matrix of the multivariate Gaussian [53,54] N (µ, Σ) is given byµ ijσ
I (µ, Σ) =0

Σ = [σij ]

0
µ
il
jk
ik
jl
σ σ +σ σ
Σ = [σkl ](96)

Notice that the lower right block matrix is a 4D tensor of dimension d × d × d × d. The zero subblock matrices
in the FIM indicate that the parameters µ and Σ are orthogonal to each other. In particular, when d = 1,
since σ11 = σ12 , we recover the Fisher information matrix of the univariate Gaussian:"I (µ, Σ) =

1
σ2

0

0

#
(97)

1
2σ4

We refer to [55] for the FIM of a Gaussian distribution using other canonical parameterizations(natural/expectation parameters of exponential family).
Example 3 (FIM of a mixture family M). A mixture family is defined for D + 1 functions F1 , . . . , FD and
C as:
(
)
D

M=

pθ ( x ) = ∑ θi Fi ( x ) + C ( x ) such that θ ∈ Θ

,

(98)

i =1

R
where the functions { Fi ( x )}i are linearly independent on the common support X and satisfying Fi ( x )dµ( x ) =R
0. Function C is such that C ( x )dµ( x ) = 1. Mixture families include statistical mixtures with prescribedcomponent distributions and the probability simplex ∆. The FIM of a mixture family is given by:
M I ( θ ) = EX ∼ p θ ( x ) Z
Fi ( x ) Fj ( x )Fi ( x ) Fj ( x )=dµ( x )  0.
2
pθ ( x )( pθ ( x ))X

(99)

The family of Gaussian mixture model (GMM) with prescribed component distributions (i.e., convex weightcombinations of D + 1 Gaussian densities) form a mixture family [56].

Entropy 2020, 22, 1100

24 of 61

Notice that the probability simplex of discrete distributions can be both modeled as an exponential
family or a mixture family [2].
The expected α-geometry is built from the expected dual ±α-connections. The Fisher“information metric” tensor is built from the FIM as follows:>
P g ( u, v ) : = ( u )θ P I ( θ ) (v )θ

(100)

The expected exponential connection and expected mixture connection are given by
e
P∇
m
P∇

:=
:=



Eθ (∂i ∂ j l )(∂k l ) ,Eθ (∂i ∂ j l + ∂i l∂ j l )(∂k l ) .

(101)
(102)

e
The dualistic structure is denoted by (P , P g, m
P ∇, P ∇) with Amari–Chentsov cubic tensor called
the skewness tensor:
Cijk := Eθ ∂i l∂ j l∂k l .
(103)

It follows that we can build a one-family of expected information α-manifolds:
(P , P g, P ∇−α , P ∇+α ) α∈R ,

(104)

with
α
P Γ ij,k ( θ )

:=

=


 1−α
Eθ ∂i ∂ j l∂k l +Cijk (θ ),
2



1−α
Eθ
∂i ∂ j l +∂i l∂ j l (∂k l ) .
2

(105)
(106)

The Levi–Civita metric connection is recovered as follows:
¯

P∇ =

P∇

−α +

2

P∇

α

LC
= LC
∇(P g)
P ∇ :=(107)

The α-Riemann–Christoffel curvature tensor is:ααrs
P Rijkl = ∂i Γ jk,l − ∂ j Γik,l + g




α
α
Γik,r
Γαjs,l − Γαjk,r Γis,l
,

(108)

α = − R−α . We check that the expected ± α-connections are coupled with the metric: ∂ g =with Rijkl
i jk
ijlk

α + Γ−α .
Γij,k
ik,j
In case of an exponential family E or a mixture family M equipped with the dual
exponential/mixture connection, we get dually flat manifolds (Bregman geometry).
Indeed, for the exponential/mixture families, it is easy to check that the Christoffel symbols of ∇eand ∇m vanish:
e
m
e
m
(109)
M Γ = M Γ = E Γ = E Γ = 0.

3.10. Criteria for Statistical Invariance
So far we have explained how to build an information manifold (or information α-manifold)from a pair of conjugate connections. Then we reported two ways to obtain such a pair ofconjugate connections: (1) from a parametric divergence, or (2) by using the predefined expected
exponential/mixture connections. We now ask the following question: which information manifold
makes sense in Statistics? We can refine the question as follows:•••

Which metric tensors g make sense in statistics?
Which affine connections ∇ make sense in statistics?
Which statistical divergences make sense in statistics (from which we can get the metric tensor
and dual connections)?

Entropy 2020, 22, 1100

25 of 61

By definition, an invariant metric tensor g shall preserve the inner product under importantstatistical mappings called Markov embeddings. Informally, we embed ∆ D into ∆ D0 with D 0 > D and
the induced metric should be preserved (see [2], page 62).
Theorem 8 (Uniqueness of Fisher information metric [57,58]). The Fisher information metric is the uniqueinvariant metric tensor under Markov embeddings up to a scaling constant.
A D-dimensional parameter (discrete) divergence satisfies the information monotonicity if andonly if:0
D (θĀ : θĀ) ≤ D (θ : θ 0 )(110)
for any coarse-grained partition A = {Ai }iE=1 of [ D ] = {1, . . . , D } (A-lumping [59]) withiE ≤ D, where θĀ= ∑ j∈Ai θ j for i ∈ [ E]. This concept of coarse-graining is illustrated in
Figure 10. This information monotonicity property could be renamed as the “distance coarse-binninginequality property.”

p1

p2 p3

p4 p5 p6 p7

p8pcoarse graining
p1 + p2 p3 + p4 + p5 p6 p7 + p8

pA0 ) ≤ D ( θ : θ 0 ).
Figure 10. A divergence satisfies the property of information monotonicity iff D (θĀ : θĀ
Here, parameter θ represents a discrete distribution.

A separable divergence D (θ1 : θ2 ) is a divergence that can be expressed as the sum of elementary
scalar divergences d( x : y):j
D (θ1 : θ2 ) := ∑ d(θ1i : θ2 ).
(111)
i

For example, the squared Euclidean distance D (θ1 : θ2 ) = ∑i (θ1i − θ2i )2 is a separable divergence for theqscalar Euclidean divergence d( x : y) = ( x − y)2 . The Euclidean distance DE (θ1 , θ2 ) = ∑i (θ1i − θ2i )2
is not separable because of the square root operation.
The only invariant and decomposable divergences when D > 1 are f -divergences [60] defined fora convex functional generator f : 0
θi
I f ( θ : θ ) : = ∑ θi f
≥ f (1),
θi
i =1
0

D

f (1) = 0.

(112)

The standard f -divergences are defined for f -generators satisfying f 0 (1) = 0 (choose f λ (u) :=f (u) + λ(u − 1) since I f λ = I f ), and f 00 (u) = 1 (scale fixed).
Statistical f -divergences are invariant [61] under one-to-one/sufficient statistic transformations
y = t( x ) of sample space: p( x; θ ) = q(y( x ); θ ):I f [ p( x; θ ) : p( x; θ 0 )]

=
=
=


p( x; θ 0 )dµ( x ),
p( x; θ )XZq(y; θ 0 )q(y; θ ) fdµ(y),
q(y; θ )YZp( x; θ ) f

I f [q(y; θ ) : q(y; θ 0 )].

The dual f -divergences for reference duality is
I f ∗ [ p( x; θ ) : p( x; θ 0 )] = I f [ p( x; θ 0 ) : p( x; θ )] = I f  [ p( x; θ ) : p( x; θ 0 )](113)

Entropy 2020, 22, 1100

26 of 61

for the standard conjugate f -generator (diamond f  generator) with:f  (u) := u f

 
1
.
u

(114)

One can check that f  is a standard f -generator when f is standard.
Let us report some common examples of f -divergences:•

The family of α-divergences:
Iα [ p : q] :=4
1 − α2


1−

Z


1− α
1+ α
p 2 ( x )q 2 ( x )dµ( x ) ,

(115)

1+ α

obtained for f (u) = 1−4α2 (1 − u 2 ). The α-divergences include:–

The Kullback–Leibler when α → 1:
KL[ p : q] =–Z

H [ p : q] :=•q( x ) log

q( x )dµ( x ) = KL[q : p],p( x )

(117)

Z q

p( x ) −q

2q( x )dµ( x ),

(118)√
for f (u) = ( u − 1)2 (corresponding to α = 0)
The Pearson and Neyman chi-squared divergences [62], etc.

The Jensen–Shannon divergence:JS[ p : q] :=

•

(116)

for f (u) = u log u.
The symmetric squared Hellinger divergence:2–p( x )dµ( x ),q( x )p( x ) log

for f (u) = − log u.
The reverse Kullback–Leibler α → −1:
KL∗ [ p : q] :=

–

Z

1
2

Z 

p( x ) log

2p( x )
2q( x )+ q( x ) logp( x ) + q( x )p( x ) + q( x )for f (u) = −(u + 1) log 1+2 u + u log u.
The Total Variation
Z
1
TV[ p : q] :=| p( x ) − q( x )| dµ( x ),2dµ( x ),

(119)

(120)

for f (u) = 12 |u − 1|. The total variation distance is the only metric f -divergence (up to ascaling factor).
The f -topology is the topology generated by open f -balls, open balls with respect to f -divergences.
A topology T is said to be stronger than a topology T 0 if T contains all the open sets of T 0 .
Csiszar’s theorem [63] states that when |α| < 1, the α-topology is equivalent to the topology inducedby the total variation metric distance. Otherwise, the α-topology is stronger than the TV topology.
Let us state an important feature of f divergences:
Theorem 9. The f -divergences are invariant by diffeomorphisms m( x ) of the sample space X : Let Y = m( X ),and Xi ∼ pi with Yi = m( Xi ) ∼ qi . Then we have I f [q1 : q2 ] = I f [ p1 : p2 ].

Entropy 2020, 22, 1100

27 of 61

Example 4. Consider the exponential distributions and the Rayleigh distributions which are related by:√

X ∼ Exponential(λ) ⇔ Y = m( X ) =1
X ∼ Rayleigh σ = √.
2λThe densities of the exponential distributions are defined bypλ ( x ) = λ exp(−λx ) with support X = [0, ∞),and the densities of the Rayleigh distributions are defined byx2xqσ ( x ) = 2 exp − 2 with support X = [0, ∞).
σ
2σWe have

λ22
λ21

DKL [qσ1 : qσ2 ] = log!

+

σ12 − σ22.
σ22

It follows that

DKL pλ1 : pλ2





=
=
=


1
1
DKL q √
: q√
2λ1
2λ1


1
2λ1
1
log
+ 2λ2
−
2λ2
2λ1
λ2
 
λ1
λ
log
+ 2 − 1.
λ2
λ1

A remarkable property is that invariant standard f -divergences yield the Fisher information
matrix and the α-connections. Indeed, the invariant standard f -divergences is related infinitesimallyto the Fisher metric as follows:Zp( x; θ + dθ )dµ( x )
(121)
I f [ p( x; θ ) : p( x; θ + dθ )] =p( x; θ ) fp( x; θ )
1
Σ
i
j
=
(122)
F g ( θ )dθ dθ .
2 ij
A statistical parameter divergence D on a parametric family of distributions P yields an equivalentparameter divergence P D:
0
0
(123)
P D ( θ : θ ) : = D [ p ( x; θ ) : p ( x; θ )].
Thus we can build the information manifold induced by this parameter divergence P D (· : ·). For P D (· :I( I )∗

I∗·) = I f [· : ·], the induced ±1-divergence connections Pf ∇ := P I f ∇ and P f ∇ := P f ∇ are precisely
the expected ±α-connections (derived from the exponential/mixture connections) with:
α = 2 f 000 (1) + 3.

(124)

Thus the invariant connections which coincide with the connections induced by the invariant
statistical divergences are the expected α-connections. Note that the curvature of an expected
α-connection depends both on α and on the considered statistical model [64].
3.11. Fisher–Rao Expected Riemannian Manifolds: (P , P g)
Historically, a first manifold modeling of a regular parametric family of distributions P ={ pθ ( x )}θ was to consider the Fisher Information Matrix (FIM) as the Riemannian metric tensor g(see [65,66]), with:

Entropy 2020, 22, 1100

28 of 61

P I ( θ ) : = E pθ∂i l∂ j l ,(125)where ∂i l:=: ∂θ∂ log p( x; θ ). Under some regularity conditions, we can rewrite the FIM:i

P I ( θ ) : = − E pθ∂i ∂ j l .

(126)

The Riemannian geodesic metric distance Dρ is commonly called the Fisher–Rao distance:Dρ ( p θ 1 , p θ 2 ) =Z 1q0γ̇(t)> gγ(t) γ̇(t)dt,(127)

where γ denotes the geodesic passing through γ(0) = θ1 and γ(1)q= θ2 . The Fisher–Rao distance can
R1also be defined as the shortest path length: Dρ ( pθ1 , pθ2 ) = infγ 0 γ̇(t)> gγ(t) γ̇(t)dt.
Definition 6 (Fisher–Rao distance). The Fisher–Rao distance is the geodesic metric distance of theFisher–Riemannian manifold (P , P g).
Let us give some examples of Fisher–Riemannian manifolds:•••

The Fisher–Riemannian manifold of the family of categorical distributions (also called finite
discrete distributions in [2]) amount to the spherical geometry [14] (spherical manifold).
The Fisher–Riemannian manifold of the family of bivariate location-scale families amount tohyperbolic geometry (hyperbolic manifold).
The Fisher–Riemannian manifold of the family of location families amount to Euclidean geometry(Euclidean manifold).
ΣThe first fundamental form of the Riemannian geometry is ds2 = hdx, dx i = gij dxi dx j where ds
denotes the line element. Let us give an example of Fisher–Rao geometry for location-scale families:
Example 5. Consider the location-scale family induced by a symmetric probability density f ( x ) with respect toRRR
0 such that X f ( x )dµ( x ) = 1, X x f ( x )dµ( x ) = 0 and X x2 f ( x )dµ( x ) = 1 (with support X = R):P=pθ ( x ) =

1
f
θ2



x − θ1θ2
θ = (θ1 , θ2 ) ∈ R × R++

,

.

(128)

The density f ( x ) is called the standard density of the location-scale family, and corresponds to the parameterθ0 = (0, 1): p(0,1) ( x ) = f ( x ). The parameter space Θ = R × R++ corresponds to the upper plane, and theFisher information matrix can be structurally calculated [67] as the following diagonal matrix:"
I (θ ) =

a2
0

0
b2

#
,

(129)

f ( x )dµ( x ),(130)f ( x )dµ( x ).

(131)

with scalars:
2

a :=
2

b :=

Z 

2
Z  0
f (x)
f (x)
f 0 (x)
+1
x
f (x)

2

Entropy 2020, 22, 1100

29 of 61

By rescaling θ = (θ1 , θ2 ) as θ 0 = (θ10 , θ20 ) with θ10 =a√θ and θ20 = θ2 , we get the FIM with respect tob 2 1

θ 0 expressed as:b2
I (θ ) = 0 2( θ2 )

"

0

1
0

0
1

#
.

(132)

We recognize that this metric is a constant time the metric of the Poincaré upper plane. Thus the Fisher–Rao
manifold of a location-scale family (with symmetric standard probability density f ) is isometric to the planar
hyperbolic space of negative curvature κ = − b12 . In practice, the Klein non-conformal model of hyperbolic
geometry is often used to implement computational geometric algorithms [20].
This Riemannian geometric structure applied to a family of parametric probability distributions
was first proposed by Harold Hotelling [65] (in a handwritten note of 1929, reprinted typeset in [68])and independently later by C. R. Rao [66] (1945, reprinted in [69]). In a similar vein, Jeffreys [70]proposed to use the volume element of a manifold as an invariant prior to 1946.
Notice that for a parametric family of probability distributions P , the Riemannian structure
I

I(P , P g) coincides with the self-dual conjugate connection manifold (P , P g, Pf ∇, Pf ∇∗ ) induced by asymmetric f -divergence like the squared Hellinger divergence.
The exponential map exp p at point p ∈ M provides a way to map back a vector v ∈ Tp to a point
exp p (v) ∈ M (when well-defined). The exponential map can be used to parameterize a geodesic γwith γ(0) = p and unit tangent vector γ̇(0) = v: t 7→ exp p (tv). For geodesically complete manifolds,
the exponential map is defined everywhere.
3.12. The Monotone α-Embeddings and the Gauge Freedom of the Metric
Another common mathematically equivalent expression of the FIM [16] is given by:
Iij (θ ) := 4

Z

q
∂i

p( x; θ )∂ jqp( x; θ )dµ( x ).

(133)

This form of the FIM is well-suited to prove that the FIM is always a positive semi-definite matrix [16](I (θ )  0). It turns out that we can define a family of equivalent representations of the FIM using the
α-embedding [71] of the parametric family.
First, we define the α-representation of densities l α ( x; θ ) := k α ( p( x; θ )) with:
(
1− α
2
2 ,
if α 6= 1,
1− α uk α (u) :=(134)
log u,
if α = 1.
The function l α ( x; θ ) is called the α-likelihood function. Then the α-representation of the FIM,
the α-FIM for short, is expressed as:
Iijα (θ ) :=Z∂i l α ( x; θ )∂ j l −α ( x; θ )dµ( x ).

We can rewrite compactly the α-FIM, as Iijα (θ ) =(Iijα (θ ) =

1
R1−α2

R

1− α

R

(135)

∂i l α ∂ j l −α dµ( x ). Expanding the α-FIM, we get:1+ α∂i p( x; θ ) 2 ∂ j p( x; θ ) 2 dµ( x )∂i log p( x; θ )∂ j p( x; θ )dµ( x )for α 6= ±1
for α ∈ {−1, 1}.

(136)

The 1-representation of the density is called the logarithmic representation (or e-representation),the −1-representation the mixture representation (or m-representation), and its 0-representation iscalled the square root representation. The set of α-scores vectors Bα := {∂i l α }i are interpreted as thetangent basis vectors of the α-base Bα . Thus the FIM is α-independent.

Entropy 2020, 22, 1100

30 of 61

Furthermore, the α-representation of the FIM can be rewritten under mild conditions [16] as:
Iijα (θ ) = −

2
1+α

Z

1+ α

p( x; θ ) 2 ∂i ∂ j l α ( x; θ )dµ( x ).

(137)


1−α
∂i l∂ j l ,∂i ∂ j l +
2

(138)

Since we have:α∂i ∂ j l ( x; θ ) = p
it follows that:
Iijα (θ ) = −

2
1+α

1− α
2





1−α
− Iij (θ ) +
I
2 ij



= Iij (θ ).

(139)

Notice that when α = 1, we recover the equivalent expression of the FIM (under mild conditions):
Iij1 (θ ) = − E[∇2 log p( x; θ )].

(140)

In particular, when the family is an exponential family [52] with cumulant function F (θ ) (satisfying themild conditions), we have:I ( θ ) = ∇2 F ( θ ).
(141)
Zhang [71,72] further discussed the representation/reference biduality which was confounded inthe α-geometry.
Gauge freedom of the Riemannian metric tensor has been investigated under the framework of(ρ, τ )-monotone embeddings [71–73] in information geometry: let ρ and τ be two strictly increasing
functions, and f a strictly convex function such that f 0 (ρ(u)) = τ (u) (with f ∗ denoting its convex
conjugate). Observe that the set of strictly increasing real-valued univariate functions has a groupstructure for the group operation chosen as the functional composition ◦. Let us write pθ ( x ) = p( x; θ ).
The (ρ, τ )-metric tensor ρ,τ g(θ ) = [ ρ,τ gij (θ )]ij can be derived from the (ρ, τ )-divergence:Dρ,τ ( p : q) =Z

( f (ρ( p( x ))) + f ∗ (τ (q( x ))) − ρ( p( x ))τ (q( x ))) dν( x ).

(142)

We have:ρ,τgij (θ )

=
=
=
=

Z


(∂i ρ( pθ ( x ))) ∂ j τ ( pθ ( x )) dν( x ),

(143)

Z


ρ0 ( pθ ( x ))τ 0 ( pθ ( x )) (∂i pθ ( x )) ∂ j pθ ( x ) dν( x ),

(144)

Z


f 00 (ρ( pθ ( x ))) (∂i ρ( pθ ( x ))) ∂ j ρ( pθ ( x )) dν( x ),

(145)

Z


( f ∗ )00 (τ ( pθ ( x ))) (∂i τ ( pθ ( x ))) ∂ j τ ( pθ ( x )) dν( x ).

(146)

3.13. Dually Flat Spaces and Canonical Bregman Divergences
We have described how to build a dually flat space from any strictly convex and smooth
generator F: A Hessian structure is built from F (θ ) with Riemannian Hessian metric ∇2 F (θ ), and theconvex conjugate F ∗ (η ) (obtained by the Legendre–Fenchel duality) yields the dual Hessian structurewith Riemannian Hessian metric ∇2 F ∗ (η ). The dual connections ∇ and ∇∗ are coupled with themetric. The connections are defined by their respective Christoffel symbols Γ(θ ) = 0 and Γ∗ (η ) = 0,
showing that they are flat connections.
Conversely, it can be proved [2] that given two dually flat connections ∇ and ∇∗ , we canreconstruct two dual canonical strictly convex potential functions F (θ ) and F ∗ (η ) such that η = ∇ F (θ )and θ = ∇ F ∗ (η ). The canonical divergence A F,F∗ yields the dual Bregman divergences BF and BF∗ .

Entropy 2020, 22, 1100

31 of 61

2 [40] with the
The only symmetric Bregman divergences are squared Mahalanobis distances MQ
Mahalanobis distance defined by:

MQ (θ, θ 0 ) =q( θ 0 − θ ) > Q ( θ 0 − θ ).

(147)

Let Q = LL> be the Cholesky decomposition of a positive-definite matrix Q  0. It is well-known
that the Mahalanobis distance MQ amounts to the Euclidean distance on affinely transformed points:
2
MQ(θ, θ 0 )= ∆θ > Q∆θ,

(148)

>

(149)

>

= ∆θ LL ∆θ,=

M2I ( L> θ, L> θ 0 ) = k L> θ − L> θ 0 k2 ,(150)

where ∆θ = θ 0 − θ.
2 does not satisfy the triangle inequality, but the
The squared Mahalanobis distance MQ
Mahalanobis distance MQ is a metric distance. We can convert a Mahalanobis distance MQ1 intoanother Mahalanobis distance MQ2 , and vice versa, as follows:
Proof. Let us write matrix Q = L> L  0 using the Cholesky decomposition. Then we have
MQ (θ1 , θ2 ) = M I ( L> θ1 , L> θ2 ) ⇔ M I (θ1 , θ2 ) = MQ (( L> )−1 θ1 , (( L> )−1 θ2 ).

(151)

Then we have for two symmetric positive-definite matrices Q1 = L1> L1  0 and Q2 = L2> L2  0:MQ1 (θ1 , θ2 ) = M I ( L1> θ1 , L1> θ2 ) = MQ2 (( L2> )−1 L1> θ1 , ( L2> )−1 L1> θ2 ).

(152)

It follows that we have:
MQ1 (θ1 , θ2 ) = MQ2 (( L2> )−1 L1> θ1 , ( L2> )−1 L1> θ2 ).

(153)

2 ( θ , θ ) = B ( θ , θ ) (Bregman divergence) with F ( θ ) = 1 θ > Qθ for a positive-definite
We have MQ
F 1 2
1 2
2

matrix Q  0. The convex conjugate F ∗ (η ) = 12 η > Q−1 η (with Q−1  0). We have η = Q−1 θ and
2 and M2 :
η = Qθ. We have the following identity between the dual Mahalanobis divergences MQ
Q −1
2
2
MQ
( θ1 , θ2 ) = MQ
− 1 ( η1 , η2 ) .

(154)

When the Bregman generator is based on an integral, i.e., the log-normalizer F (θ ) =R
log ( exp(ht( x ), θ idµ( x )) for exponential families E , or the negative Shannon entropy F (θ ) =Rmθ ( x ) log m(η )dµ( x ) for mixture families M, the associated Bregman divergences BF,E or BF,M canbe relaxed and interpreted as a statistical distance. We explain how to obtain the reconstruction below:•

Consider an exponential family E of order D with densities defined according to a dominating
measure µ:
E = { pθ ( x ) = exp(θ > t( x ) − F (θ )) : θ ∈ Θ},(155)
where the natural parameter θ and the sufficient statistic vector t( x ) belong to RD . We have theintegral-based Bregman generator:
F (θ ) = FE ( pθ ) = log

Z
exp(θ > t( x ))dµ( x ) ,(156)

Entropy 2020, 22, 1100

32 of 61

and the dual convex conjugateF ∗ (η ) = − h( pθ ) =where h( p) = −RZp( x ) log p( x )dµ( x ),(157)

p( x ) log p( x )dµ( x ) denotes Shannon’s entropy.

Let λ(i ) denotes the i-th coordinates of vector λ, and let us calculate the inner product θ1> η2 =∑i θ1 (i )η2 (i ) of the Legendre–Fenchel divergence. We have η2 (i ) = E pθ2 [ti ( x )]. Using the linear
property of the expectation E[·], we find that ∑i θ1 (i )η2 (i ) = E pθ [∑i θ1 (i )ti ( x )]. Moreover, we
2have ∑i θ1 (i )ti ( x ) = log pθ1 ( x ) + F (θ1 ). Thus we have:θ1> η2 = E pθ


2




log pθ1 + F (θ1 ) = F (θ1 ) + E pθ log pθ1 .
2

(158)

It follows that we get
BF,E [ pθ1 : pθ2 ]=F (θ1 ) + F ∗ (η2 ) − θ1> η2 ,

(159)

=

F (θ1 ) − h( pθ2 ) − E pθ [log pθ1 ] − F (θ1 ),
2


p θ2
E pθ log=: DKL∗ [ pθ1 : pθ2 ].
2
p θ1

(160)

=

(161)

By relaxing the exponential family densities pθ1 and pθ2 to be arbitrary densities p1 and p2 ,we obtain the reverse KL divergence between p1 and p2 from the dually flat structure induced by
the integral-based log-normalizer of an exponential family:

DKL∗[ p1 : p2 ] =
=

 Z
p (x)
p2
E p2 log= p2 ( x ) log 2dµ( x ),
p1
p1 ( x )DKL [ p2 : p1 ].

(162)
(163)

Thus we have recovered the reverse Kullback–Leibler divergence DKL∗ from BF,E .
The dual divergence D ∗ [ p1 : p2 ] := D [ p2 : p1 ] is obtained by swapping the distribution parameterorders. We have:


p1
∗
∗
DKL
[
p
:
p
]
:
=
D
[
p
:
p
]
=
E
log
=: DKL [ p1 : p2 ],
(164)
∗
p1
2
2
1
1
KL
p2
∗ [ p : p ] = D [ p : p ].
and DKL∗ [ p1 : p2 ] = DKL
∗
2
1
KL 2
1

To summarize, the canonical Legendre–Fenchel divergence associated with the log-normalizer ofan exponential family amounts to the statistical reverse Kullback–Leibler divergence between pθ1
and pθ1 (or the KL divergence between the swapped corresponding densities): DKL [ pθ1 : pθ2 ] =BF (θ2 : θ1 ) = A F,F∗ (θ2 : η1 ). Notice that it is easy to check that DKL [ pθ1 : pθ2 ] = BF (θ2 : θ1 ) [74,75].
Here, we took the opposite direction by constructing DKL from BF .
We may consider an auxiliary carrier term k( x ) so that the densities write pθ ( x ) = exp(θ > t( x ) −F (θ ) + k( x )). Then the dual convex conjugate writes [76] as F ∗ (η ) = −h( pθ ) + E pθ [k( x )].

•Notice that since the Bregman generator is defined up to an affine term, we may consider theequivalent generator F (θ ) = − log pθ (ω ) instead of the integral-based generator. This approachyields ways to build formula bypassing the explicit use of the log-normalizer for calculating
various statistical distances [77].
In this second example, we consider a mixture family
(

M=

D

D

i =1

i =1

m θ = ∑ θ i p i ( x ) + (1 − ∑ θ i ) p0 ( x )

)
,

(165)

Entropy 2020, 22, 1100

33 of 61

where p0 , . . . , p D are D + 1 linearly independent probability densities. The integral-based Bregman
generator F is chosen as Shannon negentropy:
F (θ ) = FM (mθ ) = −h(mθ ) =We haveZZmθ ( x ) log mθ ( x )dµ( x ).

(166)

( pi ( x ) − p0 ( x )) log mθ ( x )dµ( x ),(167)

p0 ( x ) log mθ ( x )dµ( x ) = h× ( p0 : mθ ),(168)

ηi = [∇ F (θ )]i =and the dual convex potential function isF ∗ (η ) = −Zi.e., the cross-entropy between the density p0 and the mixture mθ . Let us calculate the innerproduct θ1> η2 of the Legendre–Fenchel divergence as follows:∑ θ1 ( i )ZZ( pi ( x ) − p0 ( x )) log mθ2 ( x )dµ( x ) =i∑ θ1 (i) pi (x) log mθ2 (x)dµ(x)i− ∑ θ1 (i ) p0 ( x ) log mθ2 ( x )dµ( x ).

(169)

i

That isθ1> η2 =Z∑ θ1 (i) pi log mθ2 dµ − ∑ θ1 (i) p0 log mθ2 dµ.
i

(170)

i

Thus it follows that we have the following statistical distance:
BF,M [mθ1 : mθ2 ]:==F (θ1 ) + F ∗ (η2 ) − θ1> η2 ,− h ( m θ1 ) −

Z

(171)

p0 ( x ) log mθ2 ( x )dµ( x ) −Z∑ θ1 (i) pi (x) log mθ (x)dµ(x)2i+ ∑ θ1 (i ) p0 ( x ) log mθ2 ( x )dµ( x ),

(172)

i

=

− h ( m θ1 ) −Z((1 − ∑ θ1 (i )) p0 ( x ) + ∑ θ1 (i ) pi ( x )) log mθ2 ( x )dµ( x ),(173)mθ1 ( x ) log mθ2 ( x )dµ( x ),(174)

m θ1 ( x )dµ( x ),
m θ2 ( x )

(175)

i

=
=
=

− h ( m θ1 ) −ZZmθ1 ( x ) log

i

(176)

DKL [mθ1 : mθ2 ].

Thus we have DKL [mθ1 : mθ2 ] = BF (θ1 : θ2 ). By relaxing the mixture densities mθ1 and mθ2 toarbitrary densities m1 and m2 , we find that the dually flat geometry induced by the negentropy ofdensities of a mixture family induces a statistical distance which corresponds to the (forward) KL
divergence. That is, we have recovered the statistical distance DKL from BF,M . Note that in general
the entropy of a mixture is not available in closed-form (because of the log sum term), except when
the component distributions have pairwise disjoint supports. This latter case includes the case of
Dirac distributions whose mixtures represent the categorical distributions.
Let us consider the dually flat spaces induced by the family of discrete Poisson distributions andthe family of continuous gamma distributions:
Example 6. Consider the family P of Poisson distributions with rate parameter λ > 0:P=

λ x e−λpλ ( x ) =,x!

λ ∈ (0, ∞) .

(177)

Entropy 2020, 22, 1100

34 of 61

This family is a univariate discrete exponential family of order one (i.e., d = 1 and D = 1) with the following
canonical decomposition of its probability mass function pλ ( x ):••••µ( x )

Base measure: ν( x ) = x! = ek( x) µ( x ) where µ is the counting measure and k ( x ) = − log( x!)represents an auxiliary measure carrier term for defining the base measure ν,
Sufficient statistics: t( x ) = x,Natural parameter: θ = θ (λ) = log(λ) ∈ Θ = R,
Log-normalizer: F (θ ) = exp(θ ) since F (θ (λ)) = λ.

Thus we can rewrite the Poisson family as the following Discrete Exponential Family (DEF):

P = { pθ ( x ) = exp(θt( x ) − F (θ ))dν( x ) : θ ∈ Θ}.

(178)

The expectation is E pθ [ x ] = F 0 (θ ) = exp(θ ), or equivalently E pλ [ x ] = λ. The variance Var pθ [ x ] =00F (θ ) = exp(θ ), or or equivalently Var pλ [ x ] = λ. The Kullback–Leibler divergence between two Poisson
distributions pλ1 and pλ2 is:
DKL pλ1 : pλ2==

BF (θ (λ2 ) : θ (λ1 )),λλ1 log 1 + λ2 − λ1 .
λ2

(179)
(180)

We recognize the expression of the univariate Kullback–Leibler divergence extended to the positive scalars.
We have η = F 0 (θ ) = λ and Iη (η ) = ( F ∗ )00 (η ) = η1 where F ∗ (η ) = η log η − η is the convex conjugate
of F (θ ). Since η = λ, we deduce that the Fisher information is Iλ (λ) = λ1 . Notice that Iθ (θ ) = exp(θ ) = λ.
Thus we check the Crouzeix identity: F 00 (θ )( F ∗ )00 (η ) = λ × λ1 = 1. Beware that although Iθ (θ ) = λ, this is
not the FIM Iλ . Using the covariance equation of the FIM of Equation (86), we have:Iλ (λ)==dθdθI (θ (λ)) ,
dλ θ
dλ
1
1
1
exp(log(λ)) = .
λ
λ
λ

(181)
(182)

The Fisher–Rao distance [78] between two Poisson distributions pλ1 and pλ2 is:Dρ ( λ 1 , λ 2 ) = 2

p

λ1 −

p

λ2 .

(183)

In general, it is easy to get the Fisher–Rao distance of uniorder families because both the length elements and thegeodesics are available in closed forms.
The following example demonstrates the computational intractability of the Fisher–Rao distance.
Example 7. Consider the parametric family of Gamma distributions [79] with probability density:
pα,β ( x ) =βα x α−1 exp(− βx )
,
Γ(α)

(184)

for shape parameter α > 0, rate parameter β > 0 and support x ∈ X = (0, ∞). Function Γ(z) =R ∞ z −1 − x
e dx is the Gamma function defined for z > 0, and satisfying Γ(n) = (n − 1)! for integers n.
0 x
The Gamma distributions { pα,β , α, β > 0} form an univariate exponential family of order 2 (i.e., d = 1 and
D = 2) with the following canonical decomposition:•••

Natural parameters: θ = (θ1 , θ2 ) with θ (λ) = (− β, α − 1) with source parameter λ = (α, β),Sufficient statistics: t( x ) = ( x, log( x )),Log-normalizer: F (θ ) = − (θ2 + 1) log (−θ1 ) + log Γ (θ2 + 1),

Entropy 2020, 22, 1100•

35 of 61

Dual parameterization: η = (η1 , η2 ) = E pθ [t( x )] = ∇ F (θ ) =Γ0 ( x )θ2 +1−θ1 , − log (− θ1 ) + ψ ( θ2 + 1),dln(Γ( x )) = Γ( x) denotes the digamma function.
where ψ( x ) = dx

It follows that the Kullback–Leibler divergence between two Gamma distributions pλ1 and pλ2 withrespective source parameters λ1 = (α1 , β 1 ) and λ2 = (α2 , β 2 ) is:
DKL pα1 ,β1 : pα2 ,β2=

BF (θ (λ2 ) : θ (λ1 )),(185)= (α1 − α2 ) ψ (α1 ) − log Γ (α1 ) + log Γ (α2 )β − β1.
+α2 (log β 1 − log β 2 ) + α1 2
β1

(186)
(187)

The Fisher information matrix is Iθ (θ ) = ∇2 F (θ ). It can be expressed using the λ-parameterization [80] as:"
Iλ (α, β) =ψ1 (α)

− β1

− β1

α
β2

#
,

(188)

where ψ1 ( x ) is the trigamma function defined for x > 0 by:ψ1 ( x ) :=d2
log Γ( x ).
dx2

(189)

Because the Fisher information matrix is not diagonal using the λ-parameterization, we deduce that theparameters α and β are correlated (non-orthogonal). In general, by mixing the natural parameters θ with the
expectation parameters η of an exponential family, we obtain a block-diagonal Fisher information matrix [2,80].
That is, let δ = (θ1 , . . . , θl , η1 , . . . , ηD−l ) be a mixed primal/dual coordinate system for l ∈ {1, . . . , D − 1}
where D is the order of the family (or the dimension of the dually flat space). Then the Fisher information
matrix Iδ (δ) for the mixed parameterization is block diagonal. Thus we can always diagonalize the Fisher
information matrix of an exponential family of order 2. For example, for the gamma manifold, let us choose thereparameterization δ1 = αβ and δ2 = α. Then we have the Fisher information matrix that rewrites as:" δ

2Iδ (δ) =

#

0

δ12

ψ1 (δ2 ) − δ12

0

.

(190)

The parameters δ1 and δ2 are not correlated and orthogonal since the FIM is diagonal.
The numerical evaluation of the Fisher–Rao distance between two gamma distributions has been studied
in [81]. Let ω = log αβ . The length element is shown to be in this (α, ω ) parameterization:ds2 =
ψ1 (α) −1α(dα)2 + α(dω )2 .

(191)

However, no closed-form expression is known for the Fisher–Rao distance between two gamma distributions
because of the intractability of the geodesic equations on the gamma Fisher–Rao manifold [81]. This example
highlights the fact that computing the Fisher–Rao distance for simple family of distributions can be challenging.
In fact, we do not know the Fisher–Rao distance between any two multivariate Gaussian distibutions [82](except in a few cases including the univariate case).
In general, dually flat spaces can be built from any strictly convex C3 generator F. Vinberg and
Koszul [48] showed how to obtain such a convex generator for homogeneous cones. A cone C in avector space V yields a dual cone of positive linear functionals in the dual vector space V ∗ :C ∗ : = { ω ∈ V ∗ : ∀ v ∈ C , ω ( v ) ≥ 0} .

(192)

Entropy 2020, 22, 1100

36 of 61

The characteristic function of the cone is defined byχC (θ ) :=ZC∗

exp(−ω (θ ))dω ≥ 0,

(193)

and the function log χC (θ ) defines a Bregman generator which induces a Hessian structure and adually flat space.
Figure 11 displays the main types of information manifolds encountered in information geometry
with their relationships.
4. Some Applications of Information Geometry
Information geometry [2] found broad applications in information sciences. For example,
we can mention:
•
•
•
•
•

Statistics: Asymptotic inference, Expectation-Maximization (EM and the novel information-geometricem), time series (AutoRegressive Moving Average model, ARMA) models,Pattern recognition [83] and machine learning: Restricted Boltzmann machines [2] (RBMs),
neuromanifolds [84] and natural gradient [85],Signal processing: Principal Component Analysis (PCA), Independent Component Analysis(ICA), Non-negative Matrix Factorization (NMF),
Mathematical programming: Barrier function of interior point methods,
Game theory: Score functions.
Next, we shall describe a few applications, starting with the celebrated natural gradient descent.

4.1. Natural Gradient in Riemannian Space
The Natural Gradient [86] (NG) is an extension of the ordinary (Cartesian) gradient of Euclidean
geometry to the gradient in a Riemannian space analyzed in an arbitrary coordinate system. We explain
the natural gradient4.1.1. The Vanilla Gradient Descent Method
Given a real-valued function Lθ (θ ) parameterized by a a D-dimensional vector θ on parameter
space θ ∈ Θ ⊂ RD , we wish to minimize Lθ , i.e., solve minθ ∈Θ Lθ (θ ). The gradient descent (GD)
method, also called the steepest descent method, is a first-order local optimization procedure whichstarts by initializing the parameter to an arbitrary value (say, θ0 ∈ Θ), and then iteratively updates atstage t the current location of θt to θt+1 as follows:
GD :θ t +1 = θ t − α t ∇ θ L θ ( θ t ).

(194)

The scalar αt > 0 is called the step size or learning rate in machine learning. The ordinary gradient(OG) ∇θ Fθ (θ ) (vector of partial derivatives) represents the steepest vector at θ of the function graphLθ = {(θ, Lθ (θ )) : θ ∈ Θ}. The GD method was pioneered by Cauchy [87] (1847) and its convergence
proof to a stationary point was first reported in Curry [88] (1944).
If we reparameterize the function Lθ using a one-to-one and onto differentiable mapping η = η (θ )(with reciprocal inverse mapping θ = θ (η )), the GD update rule transforms as:
η t +1 = η t − α t ∇ η L η ( η t ),(195)

Lη (η ) := Lθ (θ (η )).

(196)

where

Entropy 2020, 22, 1100

37 of 61

Smooth Manifolds
Conjugate Connection Manifolds(M, g, ∇, ∇∗ )(M, g, C = Γ∗ − Γ)LC∗

∇ = ∇+∇
2

Riemannian Manifolds(M, g) = (M, g, LC ∇)Self-dual Manifold∇α = 1+α∇ + 1−α∇∗22Γ±α = Γ̄ ∓ αC
2(M, g, ∇−α , ∇α )(M, g, αC)
Divergence Manifold∗(M, Dg , D ∇, D ∇∗ = D ∇)D
∇ − flat ⇔ D ∇∗ − flat

Fisher-Riemannian
Manifold
Multinomial
familyI[pθ : p 0 ] = D(θ : θ 0 )θ

Parametric KL∗ on exponential families
KL on mixture families
families
Conformal divergences on deformed families
f -divergences Etc.
Bregman divergence
Expected Manifold(M, Fisher g, ∇−α , ∇α )α-geometry

Cubic skewness tensor
Cijk = E[∂i l∂j l∂k l]
αC = ∇αFisher gcanonical
divergence
Dually flat Manifolds(M, F, F ∗ )(Hessian Manifolds)
Dual Legendre potentials
Bregman Pythagorean theorem

Distance = Non-metric divergence

Location-scalefamilyLocationfamilySpherical Manifold

g = Fisher g
Fishergij = E[∂i l∂j l]Hyperbolic Manifold

Euclidean Manifold

Distance = Metric geodesic length

Frank Nielsen

Figure 11. Overview of the main types of information manifolds with their relationships in information geometry.

Entropy 2020, 22, 1100

38 of 61

Thus in general, the two gradient descent location sequences {θt }t and {ηt }t (initialized at
θ0 = θ (η0 ) and η0 = η (θ0 )) are different (because usually η (θ ) 6= θ), and the two GDs may
potentially reach different stationary points. In other words, the GD local optimization depends onthe choice of the parameterization of the function L (i.e., Lθ or Lη ). For example, minimizing with thegradient descent a temperature function Lθ (θ ) with respect to Celsius degrees θ may yield a different
result than minimizing the same temperature function Lη (η ) = Lθ (θ (η )) expressed with respect toFahrenheit degrees η. That is, the GD optimization is extrinsic since it depends on the choice of theparameterization of the function, and does not take into account the underlying geometry of theparameter space Θ.
The natural gradient precisely addresses this problem and solves it by choosing intrinsically the
steepest direction with respect to a Riemannian metric tensor field on the parameter manifold. We shall
explain the natural gradient descent method and highlight its connections with the Riemannian
gradient descent, the mirror descent and even the ordinary gradient descent when the parameter spaceis dually flat.
4.1.2. Natural Gradient and Its Connection with the Riemannian Gradient
Let ( M, g) be a D-dimensional Riemannian space [10] equipped with a metric tensor g, and L ∈C ∞ ( M ) a smooth function to minimize on the manifold M. The Riemannian gradient [89] uses the
Riemannian exponential map exp p : Tp → M to update the sequence of points pt ’s on the manifold
as follows:
RG : pt+1 = exp pt (−αt ∇ M L( pt )),(197)
where the Riemannian gradient ∇ M is defined according to a directional derivative ∇v by: 
∇ M L( p) := ∇v L exp p (v)with∇v L( p) := lim

h →0

v =0

,

L( p + hv) − L( p)
.
h

(198)

(199)

However, the Riemannian exponential mapping exp p (·) is often computationally intractable since
it requires to solve a system of second-order differential equations [10,22]. Thus instead of using exp p ,we shall rather use a computable Euclidean retraction R : Tp → RD of the exponential map expressedin a local θ-coordinate system as:
RetG :

θt+1 = Rθt (−αt ∇θ Lθ (θt )) .

(200)

Using the retraction [22] R p (v) = p + v which corresponds to a first-order Taylor approximationof the exponential map, we recover the natural gradient descent [86]:NG : θt+1 = θt − αt gθ−1 (θt )∇θ Lθ (θt ).

(201)

The natural gradient [86] (NG)NG

∇ Lθ (θ ) := gθ−1 (θ )∇θ Lθ (θ )(202)

encodes the Riemannian steepest descent vector, and the natural gradient descent method yields thefollowing update rule
NG : θt+1 = θt − αt NG ∇ Lθ (θt ).
(203)
Notice that the natural gradient is a contravariant vector while the ordinary gradient is a covariant
vector. Recall that a covariant vector [vi ] is transformed into a contravariant vector [vi ] by vi = ∑ j gij vi ,that is by using the dual Riemannian metric gη∗ (η ) = gθ (θ )−1 . The natural gradient is invariant underEntropy 2020, 22, 1100

39 of 61

an invertible smooth change of parameterization. However, the natural gradient descent does notguarantee that the locations θt ’s always stay on the manifold: Indeed, it may happen that for some t,θt 6∈ Θ when Θ 6= RD .
Property 5 ([89]). The natural gradient descent approximates the intrinsic Riemannian gradient descent using
a contravariant gradient vector induced by the Riemannian metric tensor g. The natural gradient is invariant tocoordinate transformations.
Next, we shall explain how the natural gradient descent is related to the mirror descent and theordinary gradient when the Riemannian space Θ is dually flat.
4.1.3. Natural Gradient in Dually Flat Spaces: Connections to Bregman Mirror Descent andOrdinary Gradient
Recall that a dually flat space ( M, g, ∇, ∇∗ ) is a manifold M equipped with a pair (∇, ∇∗ ) of dual
torsion-free flat connections which are coupled to the Riemannian metric tensor g [2,90] in the sense∗that ∇+∇= LC ∇, where LC ∇ denotes the unique metric torsion-free Levi–Civita connection.
2
On a dually flat space, there exists a pair of dual global Hessian structures [48] with dual canonicalBregman divergences [2,91]. The dual Riemannian metrics can be expressed as the Hessians ofdual convex potential functions F and F ∗ . Examples of Hessian manifolds are the manifolds ofexponential families or the manifolds of mixture families [92]. On a dually flat space induced
by a strictly convex and C3 function F (Bregman generator), we have two dual global coordinate
system: θ (η ) = ∇ F ∗ (η ) and η (θ ) = ∇ F (θ ), where F ∗ denotes the Legendre–Fenchel convex conjugate
function [51,93]. The Hessian metric expressed in the primal θ-coordinate system is gθ (θ ) = ∇2 F (θ ),and the dual Hessian metric expressed in the dual coordinate system is gη∗ (η ) = ∇2 F ∗ (η ). Crouzeix’s
identity [36] shows that gθ (θ ) gη (η ) = I, where I denotes the D × D matrix identity.
The ordinary gradient descent method can be extended using a proximity function Φ(·, ·)as follows:
1
PGD : θt+1 = arg min hθ, ∇ Lθ (θt )i + Φ(θ, θt ) .
(204)
αt
θ ∈Θ
When Φ(θ, θt ) = 12 kθ − θt k2 , the PGD update rule becomes the ordinary GD update rule.
Consider a Bregman divergence [91] BF for the proximity function Φ: Φ( p, q) = BF ( p : q).
Then the PGD yields the following mirror descent (MD):


MD :


1
θt+1 = arg min hθ, ∇ L(θt )i + BF (θ : θt ) .
αt
θ ∈Θ

(205)

This mirror descent can be interpreted as a natural gradient descent as follows:
Property 6 ([94]). Bregman mirror descent on the Hessian manifold ( M, g = ∇2 F (θ )) is equivalent to natural
gradient descent on the dual Hessian manifold ( M, g∗ = ∇2 F (η )), where F is a Bregman generator, η = ∇ F (θ )and θ = ∇ F ∗ (η ).
Indeed, the mirror descent rule yields the following natural gradient update rule:NG∗ : ηt+1

= ηt − αt ( gη∗ )−1 (ηt )∇η Lθ (θ (ηt )),(206)

ηt − αt ( gη∗ )−1 (ηt )∇η Lη (ηt ),

(207)

=

where gη∗ (η ) = ∇2 F ∗ (η ) = (∇2θ F (θ ))−1 and θ (η ) = ∇ F ∗ (θ ).
The method is called mirror descent [95] because it performs that gradient step in the dual
space (i.e., mirror space) H = {η = ∇ F (θ ) : θ ∈ Θ}, and thus solves the inconsistency

Entropy 2020, 22, 1100

40 of 61

contravariant/covariant type problem of subtracting a covariant vector from a contravariant vector of
the ordinary GD (Equation (194)).
Let us prove now the following property of the natural gradient in a dually flat space or Bregman
manifold [90]:Property 7 ([96]). In a dually flat space induced by potential convex function F, the natural gradient amounts
to the ordinary gradient on the dually parameterized function: NG ∇ Lθ (θ ) = ∇η Lη (η ) where η = ∇θ F (θ )and Lη (η ) = Lθ (θ (η )).
Proof. Let ( M, g, ∇, ∇∗ ) be a dually flat space. We have gθ (θ ) = ∇2 F (θ ) = ∇θ ∇θ F (θ ) = ∇θ η sinceη = ∇θ F (θ ). The function to minimize can be written either as Lθ (θ ) = Lθ (θ (η )) or as Lη (η ) =Lη (η (θ )). Recall the chain rule in the calculus of differentiation:∇θ Lθ (θ ) = ∇θ ( Lη (η (θ ))) = (∇θ η )(∇η Lη (η )).

(208)

Thus we have:NG∇ Lθ (θ ) :=gθ−1 (θ )∇θ Lθ (θ ),

(209)

−1

(210)

=

(∇θ η )(∇θ η )∇η Lη (η ),=∇ η L η ( η ).

(211)

It follows that the natural gradient descent on a loss function Lθ (θ ) amounts to an ordinary
gradient descent on the dually parameterized loss function Lη (η ) := Lθ (θ (η )). In short,NG ∇ L = ∇ L .
η ηθ θ
4.1.4. An Application of the Natural Gradient: Natural Evolution Strategies (NESs)
A nice family of applications of the natural gradient is the Natural Evolution Strategies (NESs)
for black-box minimization [97]. Let f ( x ) for x ∈ X ⊂ Rd be a real-valued function to minimize.
Berny [98] proposed to relax the optimization problem minx∈X f ( x ) by considering a parametric search
distribution pλ , and minimize instead:
min E pλ [ f ( x )],
λ∈Λ

(212)

where λ ∈ Λ ⊂ RD denotes the parameter space of the search distributions. Let J (λ) =E pλ [ f ( x )]. Minimizing J (λ) instead of f ( x ) is particularly useful when X is a discrete space: Indeed,
the combinatorial optimization [98] minx∈X f ( x ) is replaced by a continuous optimization minλ∈Λ J (λ)
when Λ is a continuous parameter, and the ordinary or natural GD methods can be used. The gradient∇ J (λ) is called the search gradient, and it can be approximated stochastically using the log-likelihood
trick [99] as
∼
1 n
(213)
∇ J (λ) := ∑ f ( xi )∇ log pλ ( xi ) ≈ ∇ J (λ),
n i =1
where x1 , . . . , xn ∼ pλ . Similarly, the Fisher information matrix (FIM) may be approximated by thefollowing empirical FIM:Ĩ (λ) =1 n∇λ lλ ( xi )(∇λ lλ ( xi ))> ≈ I (λ),
n i∑
=1

(214)

where lλ ( x ) := log pλ ( x ) denote the log-likelihood function. Notice that the approximated FIM
may potentially be degenerated and may not respect the structure of the true FIM. For example,

Entropy 2020, 22, 1100

41 of 61x −µ
( x − µ )2and∇
− 2σ1 2 . The non-diagonal of the approximate FIM Ĩ (λ)
2 =
2
σ
σ
2σ4


are close to but usually non-zero although the expected FIM is diagonal I (µ, σ2 ) = diag σ12 , 2σ1 4 .

we have ∇µ l ( x; µ, σ2 ) =Thus we may estimate the FIM until the non-diagonal elements have absolute values less than aprescribed e > 0. For multivariate normals, we have ∇µ l ( x; µ, Σ) = Σ−1 ( x − µ) and ∇Σ l ( x; µ, Σ) =1>
−1
2 (∇µ l ( x; µ, Σ )∇µ l ( x; µ, Σ ) − Σ ).
4.2. Some Illustrating Applications of Dually Flat Manifolds
In this part, we describe how to use the dually flat structures for handling an exponential family
E (in a hypothesis testing problem detailed in Section 4.3) and the mixture family M (clustering
statistical mixtures Section 4.4). Note that for a general divergence, neither (E , D ) nor (M, D ) is dually
flat. However, when D = KL, the Kullback–Leibler divergence, we get dually flat spaces that arecomputationally attractive since the primal/dual geodesics are straight lines in the corresponding
global affine coordinate system.
4.3. Hypothesis Testing in the Dually Flat Exponential Family Manifold (E , KL∗ )Given two probability distributions P0 ∼ p0 ( x ) and P1 ∼ p1 ( x ), we ask to classify a set of iid.
observations X1:n = { x1 , . . . , xn } as either sampled from P0 or from P1 ? This is a statistical decision
problem [100]. For example, P0 can represent the signal distribution and P1 the noise distribution.
Figure 12 displays the probability distributions and the unavoidable error that is made by any statisticaldecision rule (on observations x1 and x2 ).

p1 (x)

p0 (x)

x1

x2

xFigure 12. Statistical Bayesian hypothesis testing: the best Maximum A Posteriori (MAP) rule chooses
to classify an observation from the class that yields the maximum likelihood.

Assume that both distributions P0 ∼ Pθ0 and P1 ∼ Pθ1 belong to the same exponential family
E = { Pθ : θ ∈ Θ}, and consider the exponential family manifold with the dually flat structure(E , E g, E ∇e , E ∇m ). That is, the manifold equipped with the Fisher information metric tensor field andthe expected exponential connection and conjugate expected mixture connection. More generally,
the expected α-geometry of an exponential family E with cumulant function F is given by:gij (θ )αΓij,k= ∂ i ∂ j F ( θ ),1−α=∂ i ∂ j ∂ k F ( θ ).
2

(215)
(216)

α = 0 and ∇1 is flat, and so is ∇−1 by using the fundamental theorem of information
When α = 1, Γij,kgeometry.
The ±1-structure can also be derived from a divergence manifold structure by choosing thereverse Kullback–Leibler divergence KL∗ :(E , E g, E ∇e , E ∇m ) ≡ (E , KL∗ ).

(217)

Therefore, the Kullback–Leibler divergence KL[ Pθ : Pθ 0 ] amounts to a Bregman divergence (forthe cumulant function of the exponential family):KL∗ [ Pθ 0 : Pθ ] = KL[ Pθ : Pθ 0 ] = BF (θ 0 : θ ).

(218)

Entropy 2020, 22, 1100

42 of 61

The best exponent error α∗ of the best Maximum A Priori (MAP) decision rule is found by
minimizing the Bhattacharyya distance to get the Chernoff information [101]:
C [ P1 , P2 ] = − log minZα∈(0,1) x ∈X

p1α ( x ) p12−α ( x )dµ( x ) ≥ 0.

(219)

On the exponential family manifold E , the Bhattacharyya distance:
Bα [ p1 : p2 ] = − logZx ∈X

p1α ( x ) p12−α ( x )dµ( x ),(220)

amounts to a skew Jensen parameter divergence [102] (also called Burbea-Rao divergence):JFα (θ1 : θ2 ) = αF (θ1 ) + (1 − α) F (θ2 ) − F (θ1 + (1 − α)θ2 ).

(221)

It can be shown that the Chernoff information [100,103,104] (that minimizes α) is equivalent
to a Bregman divergence: Namely, the Bregman divergence for exponential families at the optimal
exponent value α∗ .
Theorem 10 (Chernoff information [100]). The Chernoff information between two distributions belonging to
the same exponential family amount to a Bregman divergence:
∗

∗

α
α
C [ Pθ1 : Pθ2 ] = B(θ1 : θ12) = B(θ2 : θ12
),

(222)

α = (1 − α ) θ + αθ , and α∗ denote the best exponent error.
where θ122
1∗

∗ : = θ α denote the best exponent error. The geometry [100] of the best error exponent can
Let θ1212
be explained on the dually flat exponential family manifold as follows:∗ = Ge ( P1 , P2 ) ∩ Bim ( P1 , P2 ),P∗ = Pθ12(223)

where Ge denotes the exponential geodesic γ∇e and Bim the m-bisector:
Bim ( P1 , P2 ) = { P : F (θ1 ) − F (θ2 ) + η ( P)> (θ2 − θ1 ) = 0}.

(224)

Figure 13 illustrates how to retrieve the best error exponent from an exponential arc (θ-geodesic)
intersecting the m-bisector.
m-bisector
Bim (Pθ1 , Pθ2 )η-coordinate systeme-geodesic Ge (Pθ1 , Pθ2 )
∗
pθ12
pθ2
Pθ ∗
12

pθ1
∗
C(θ1 : θ2 ) = B(θ1 : θ12)

Figure 13. Exact geometric characterization (not necessarily i closed-form) of the best exponent error
rate α∗ .

Furthermore, instead of considering two distributions for this statistical binary decision problem,
we may consider a set of n distributions of P1 , . . . , Pn ∈ E . The geometry of the error exponent in thismultiple hypothesis testing setting has been investigated in [105]. On the dually flat exponential family

Entropy 2020, 22, 1100

43 of 61

manifold, it corresponds to check the exponential arcs between natural neighbors (sharing Voronoi
subfaces) of a Bregman Voronoi diagram [40]. See Figure 14 for an illustration.

η-coordinate system

Chernoff distribution between
natural neighbours
Figure 14. Geometric characterization of the best exponent error rate in the multiple hypothesis
testing case.

4.4. Clustering Mixtures in the Dually Flat Mixture Family Manifold (M, KL)
Given a set of k prescribed statistical distributions p0 ( x ), . . . , pk−1 ( x ), all sharing the same support
X (say, R), a mixture family M of order D = k − 1 consists of all strictly convex combinations of thesecomponent distributions [56]:
(
!
)
k −1

M :=

m( x; θ ) = ∑ θi pi ( x ) +
i =1

k −1

k −1

1 − ∑ θip0 ( x ) such that θi > 0, ∑ θi < 1 .

i =1

(225)

i =1

Figure 15 displays two mixtures obtained as convex combinations of prescribed Laplacian,Gaussian and Cauchy component distributions (D = 2). When considering a set of prescribed
Gaussian component distributions, we obtain a w-Gaussian Mixture Model, or w-GMM for short.
0.5
M1
M2
Gaussian(-2,1)
Cauchy(2,1)
Laplace(0,1)

0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
-4

-2

0

2

4

6

8

Figure 15. Example of a mixture family of order D = 2 (3 components: Laplacian, Gaussian and
Cauchy prefixed distributions).

Entropy 2020, 22, 1100

44 of 61

We consider the expected information manifold (M, M g, M ∇m , M ∇e ) which is dually flatand equivalent to ( MΘ , KL). That is, the KL between two mixtures with prescribed components(w-mixtures, for short) is equivalent to a Bregman divergence for F (θ ) = −h(mθ ), where h( p) =R
p( x ) log p( x )dµ( x ) is the differential Shannon information (negative entropy) [56]:KL[mθ1 : mθ2 ] = BF (θ1 : θ2 ).

(226)

Consider a set {mθ1 , . . . , mθn } of n w-mixtures [56]. Because F (θ ) = −h(m( x; θ )) is the negative
differential entropy of a mixture (not available in closed form [106]), we approximate the untractable F
by another close tractable generator F̃. We use Monte Carlo stochastic sampling to get Monte-Carloconvex F̃S for an independent and identically distributed sample S .
Thus we can build a nested sequence (M, F̃S1 ), . . . , (M, F̃Sm ) of tractable dually flat manifolds
for nested sample sets S1 ⊂ . . . ⊂ Sm converging to the ideal mixture manifold (M, F ):limm→∞ (M, F̃Sm ) = (M, F ) (where convergence is defined with respect to the induced canonicalBregman divergence). A key advantage of this approach is that for a given sample S , all computations
carried inside the dually flat manifold (M, F̃S ) are consistent, see [56].
For example, we can apply Bregman k-means [107] on these Monte Carlo dually flat spaces [108]of w-GMMs (Gaussian Mixture Models) to cluster a set of w-GMMs. Figure 16 displays the result ofsuch a clustering.
0.25

0.2

0.15

0.1

0.05

0
-4

-2

0

2

4

Figure 16. Example of w-GMM clustering into k = 2 clusters.

We have briefly described two applications using dually flat manifolds: (1) the dually flat
exponential manifold induced by the statistical reverse Kullback–Leibler divergence on an exponential
family (structure (E , KL∗ )), and (2) the dually flat mixture manifold induced by the statistical
Kullback–Leibler divergence on a mixture family (structure (M, KL)). There are many other duallyflat structures that can be met in a statistical context. For example, two other dually flat structures forthe D-dimensional probability simplex ∆ D are reported in Amari’s textbook [2]: (1) the conformallydeforming of the α-geometry (page 88, Equation 4.95 of [2]), and (2) the χ-escort geometry (page 91,
Equation 4.114 of [2]).
5. Conclusions: Summary, Historical Background, and Perspectives
5.1. Summary
We explained the dualistic nature of information manifolds ( M, g, ∇, ∇∗ ) in information geometry.
The dualistic structure is defined by a pair of conjugate connections coupled with the metric tensor thatprovides a dual parallel transport that preserves the metric. We showed how to extend this structure
to a 1-parameter family of structures. From a pair of conjugate connections, the pipeline to build this
1-parameter family of structures can be informally summarized as:Entropy 2020, 22, 1100

45 of 61

( M, g, ∇, ∇∗ ) ⇒ ( M, g, C ) ⇒ ( M, g, αC ) ⇒ ( M, g, ∇−α , ∇α ),∀α ∈ R.

(227)

We stated the fundamental theorem of information geometry on dual constant-curvature
manifolds, including the special but important case of dually flat manifolds on which there exists
two potential functions and global affine coordinate systems related by the Legendre–Fencheltransformation. Although, information geometry historically started with the Riemannian modeling(P , P g) of a parametric family of probability distributions P by letting the metric tensor be the Fisher
information matrix, we have emphasized the dualistic view of information geometry which considers
non-Riemannian manifolds that can be derived from any divergence, and not necessarily tied toa statistical context (e.g., information manifold can be used in mathematical programming [109]).
Let us notice that for any symmetric divergence (e.g. any symmetrized f -divergence like the squared
Hellinger divergence), the induced conjugate connections coincide with the Levi–Civita connection
but the Fisher–Rao metric distance does not coincide with the squared Hellinger divergence.
On one hand, a Riemannian metric distance Dρ is never a divergence because the rooted distancefunctions fail to be smooth at the extremities but a squared Riemmanian metric distance is always
a divergence. On the other hand, taking the power δ of a divergence D (i.e., D δ ) for some δ > 0
may yield a metric distance (e.g., the square root of the Jensen–Shannon divergence [110]), but thismay not always be the case: the powered Jeffreys divergence J δ is never a metric distance (see [111],
page 889). Recently, the Optimal Transport (OT) theory [112] gained interest in statistics and machinelearning. However, the optimal transport between two members of a same elliptically-contoured
family has the same optimal transport formula distance (see [113] Eq. 16 and Eq. 17, although they
have different Kullback–Leibler divergences). Another essential difference is that the Fisher–Rao
manifold of location-scale families is hyperbolic but the Wasserstein manifold of location-scale families
has positive curvature [113,114].
Notice that we may convert back and forth a similarity S( p, q) ∈ (0, 1] to a dissimilarity D ( p, q) ∈[0, ∞) as follows:S( p, q)= exp (− D ( p, q)) ∈ (0, 1],(228)

D ( p, q)= − log (S( p, q)) ∈ [0, ∞).

(229)

When the dissimilarity satisfies the (additive) triangle inequality (i.e., D ( p, q) + D (q, r ) ≥ D ( p, r ) for
any triple ( p, q, r )) then the corresponding similarity satisfies the multiplicative triangle inequality:
S( p, q) × S(q, r ) ≤ S( p, r ). A metric transform on a metric distance D is a transformation T such that
T ( D ( p, q)) is a metric. The transformation T (u) = 1+1 u is a metric transform which bounds potentiallyD ( p,q)unbounded metric distances; that is, if D is an unbounded metric, then T ( D ( p, q)) = 1+ D( p,q) is abounded metric distance. The transformation S(u) = u2 is not a metric transform since the squared of
the Euclidean metric distance is not a metric distance.
5.2. A Brief Historical Review of Information Geometry
The field of Information Geometry (IG) was historically motivated by providing somedifferential-geometric structures to statistical models in order to reason geometrically about statisticalproblems with the endeavor goal of geometrizing mathematical statistics [4,12–14,115–118]: Professor
Harold Hotelling [65] first considered in the late 1920s the Fisher Information Matrix (FIM) I as a
Riemannian metric tensor g (i.e., the Fisher Information metric, FIm), and interpreted a parametricfamily of probability distributions M as a Riemannian manifold ( M, g). Historically speaking,Hotelling attended the American Mathematical Society’s Annual Meeting in Bethlehem (Bethlehem,
PA, USA) on 26–29 December 1929, but left before his scheduled talk on December 27. His handwritten
notes on the “Spaces of Statistical Parameters” was read by a colleague and are fully typeset in [68].
We warmly thank Professor Stigler for sending us the scanned handwritten notes and for discussing
by emails some historical aspects of the birth of information geometry. In this pioneering work,

Entropy 2020, 22, 1100

46 of 61

Hotelling mentioned that location-scale probability families yield Riemannian manifolds of constantnon-positive curvatures. This Riemannian modeling of parametric family of densities was further
independently studied by Calyampudi Radhakrishna Rao (C.R. Rao) in his celebrated paper [66](1945) that also includes the Cramér–Rao lower bound [51] and the Rao–Blackwellization techniqueused in statistics. Nowadays the induced Riemannian metric distance is often called the Fisher–Rao
distance [119] or Rao distance [81]. Yet another use of Riemannian geometry in statistics was pioneeredby Harold Jeffreys [70] that proposed to use as an invariant prior the normalized volume elementof the Fisher–Riemannian manifold. In those seminal papers, there was no theoretical justificationof using the Fisher information matrix as a metric tensor (besides the fact that it is a well-defined
positive-definite matrix for regular identifiable models). Nowadays, this Riemmanian metric tensor
is called the information metric for short. Modern information geometry considers a generalizationof this approach using a non-Riemannian dualistic modeling ( M, g, ∇, ∇∗ ) which coincides with the
Riemannian manifold when ∇ = ∇∗ = LC ∇, the Levi–Civita connection (the unique torsion-freeaffine connection compatible with the metric tensor). The Fisher–Rao geometry has also been exploredin thermodynamics yielding the Ruppeiner geometry [120], and the geometry of thermodynamics iscalled nowadays called geometrothermodynamics [121].
In the 1960s, Nikolai Chentsov (also commonly written Čencov) studied the algebraic categoryof all statistical decision rules with its induced geometric structures: Namely, the α-geometries(“equivalent differential geometry”) and the dually flat manifolds (“Nonsymmetric Pythagorean
geometry” of the exponential families with respect to the Kullback–Leibler divergence). In the prefaceof the english translation of his 1972 russian monograph [115], the field of investigation is definedas “geometrical statistics.” However in the original Russian monograph, Chentsov used the russian
term geometrostatistics. According to Professor Alexander Holevo, the geometrostatistics term wascoined by Andrey Kolmogorov to define the field of differential geometry of statistical models. In themonograph of Chentsov [115], the Fisher information metric is shown to be the unique metric tensor(up to a scaling factor) yielding statistical invariance under Markov morphisms (see [57] for a simpler
proof that generalizes to positive measures).
The dual nature of the information geometry was thoroughly investigated by Professor Shun-ichi
Amari [122]. In the preface of his 1985 monograph [116], Professor Amari coined the term informationgeometry as follows: “The differential-geometrical method developed in statistics is also applicable toother fields of sciences such as information theory and systems theory... They together will open a new
field, which I would like to call information geometry.” Professor Amari mentioned in [116] that heconsidered the Gaussian Riemannian manifold as a hyperbolic manifold in 1959, and was strongly
influenced by Efron’s paper on statistical curvature [123] (1975) to study the family of α-connections in
the 1980s [122,124]. Professor Amari prepared his PhD under the supervision of Professor Kondo [125],
an expert of differential geometry in touch with Professor Kawaguchi [126]. The role of differential
geometry in statistics has been discussed in [127].
Note that the dual affine connections of information geometry have also been investigatedindependently in affine differential geometry [128] which considers invariance undervolume-preserving affine transformations by defining a volume form instead of a metric form
for Riemannian geometry. The notion of dual parallel transport compatible with the metric isdue to Aleksandr Norden [129] and Rabindra Nath Sen [130–132] (see the Senian geometry in
http://insaindia.res.in/detail/N54-0728).

Entropy 2020, 22, 1100

47 of 61

We summarize the main fundamental structures of information manifolds below:( M, g)(P , P g)( M, g, ∇)(P , P g, P e ∇α )( M, g, ∇, ∇∗ )(P , P g, P ∇−α , P ∇α )( M, g, C )∗( M, D g, D ∇, D ∇)( M, F g, F C )Riemannian manifold
Fisher–Riemannian (expected) Riemannian manifold
Riemannian manifold ( M, g) with affine connection ∇Chentsov’s manifold with affine exponential α-connection
Amari’s dualistic information manifold
Amari’s (expected) information α-manifold, α-geometry
Lauritzen’s statistical manifold [29]Eguchi’s conjugate connection manifold induced by divergence D
Chentsov/Amari’s dually flat manifold induced by convex potential F

We use the ≡ symbol to denote the equivalence of geometric structures. For example, we have∗
( M, g) ≡ ( M, g, LC ∇, LC ∇ = LC ∇).
5.3. Perspectives
We recommend the two recent textbooks [2,16] for an indepth covering of (parametric) informationgeometry, and the book [133] for a thorough description of some infinite-dimensional statistical models.
(Japanese readers may refer to [134,135]) We did not report the various coefficients of the metric tensors,
Christoffel symbols and skewness tensors for the expected α-geometry of common parametric models
like the multivariate Gaussian distributions, the Gamma/Beta distributions, etc. They can be foundin [15,16] and in various articles dealing with less common family of distributions [15,64,136–140].
Although we have focused on the finite parametric setting, information geometry is also considering
non-parametric families of distributions [141], and quantum information geometry [142].
We have shown that we can always create an information manifold ( M, D ) from any divergence
function D. It is therefore important to consider generic classes of divergences in applications,that are ideally axiomatized and shown to have exhaustive characteristics. The α-skewed Jensen
divergences [102] are defined by for a real-valued strictly convex function F (θ ) by:JFα (θ1 : θ2 ) := (1 − α) F (θ1 ) + αF (θ2 ) − F ((1 − α)θ1 + αθ2 ) > 0,(230)

where both θ1 and θ2 belong to the parameter space Θ. Clearly, we have JFα (θ1 : θ2 ) = JF1−α (θ2 : θ1 ).
We have the following asymptotic properties of skewed Jensen divergences [102]:lim

1

α →0+ α ( 1 − α )lim1

α →1− α ( 1 − α )JFα (θ1 : θ2 )=B F ( θ1 : θ2 ),(231)

JFα (θ1 : θ2 )=B F ( θ2 : θ1 ),(232)

where BF (θ1 : θ2 ) is the Bregman divergence [91] induced by a strictly convex and differentiable
function F (θ ):B F ( θ1 : θ2 ) : = F ( θ1 ) − F ( θ2 ) − ( θ1 − θ2 ) F 0 ( θ2 ).
(233)
Appendix C further reports how to interpret geometrically these Jensen/Bregman divergences fromthe chordal slope theorem. Beyond the three main Bregman/Csiszár/Jensen classes (theses classes
overlap [143]), we may also mention the class of conformal divergences [73,144,145], the class ofprojective divergences [146,147], etc. Figure 17 illustrates the relationships between the principal
classes of distances.
There are many perspectives on information geometry as attested by the new Springer journal(see online at https://www.springer.com/mathematics/geometry/journal/41884), and the biannualinternational conference “Geometric Sciences of Information” (GSI) [148–150] with its collective
post-conference edited books [151,152]. We also mention the edited book [153] on the Occasionof Shun-ichi Amari’s 80th birthday.

Entropy 2020, 22, 1100

48 of 61

Additional materials are available online at https://FrankNielsen.github.io/SurveyIG/.
Dissimilarity measure
C3

Divergencev-Divergence Dv

Projective divergence
double sided

scaled conformal divergence CD,g (· : ·; ·)one-sided

γ-divergenceD(λp : λ0 p0 ) = D(p : p0 )Hÿvarinen SM/RM
D(λp : p0 ) = D(p : p0 )scaled Bregman divergence BF (· : ·; ·)conformal divergence CD,g (· : ·)total Bregman divergence tB(· : ·)Csiszár f -divergence If (· : ·)Bregman divergence BF (· : ·)tBF (P : Q) = √ BF (P :Q)

Dv (P : Q) = D(v(P ) : v(Q))


R
q(x)
If (P : Q) = p(x)f ( p(x)dν(x)

1+k∇F (Q)k2

BF (P : Q) = F (P ) − F (Q) − hP − Q, ∇F (Q)i

CD,g (P : Q) = g(Q)D(P : Q)QP
BF,g (P : Q; W ) = W BF Q
:W

Figure 17. Principled classes of distances/divergences.
Funding: This research received no external funding.
Acknowledgments: FN would like to thank the organizers of the Geometry In Machine Learning workshop in2018 (GiMLi, http://gimli.cc/2018/) for their kind keynote talk invitation, and specially Professor Søren Hauberg(Technical University of Denmark, DTU). This survey is based on the talk given at GiMLi. I am very thankful toProfessor Stigler (University of Chicago, USA) and Professor Holevo (Steklov Mathematical Institute, Russia) for
providing me feedback on some historical aspects of the field of information geometry. Finally, I would like toexpress my sincere thanks to Gaëtan Hadjeres (Sony Computer Science Laboratories Inc., Paris) for his careful
proofreading and feedback.
Conflicts of Interest: The author declare no conflict of interest.

Appendix A. Monte Carlo Estimations of f -Divergences
Let ( X, F, µ) be a probability space [154] with X denoting the sample space, F the σ-algebra, and µa reference positive measure. The f -divergence [62,63] between two probability measures P and Q
both absolutely continuous with respect to a positive measure µ for a convex generator f : (0, ∞) → R
strictly convex at 1 and satisfying f (1) = 0 isI f ( P : Q) = I f ( p : q) =Zp( x ) fq( x )p( x )dµ( x ),

(A1)

where P = pdµ and Q = qdµ (i.e., p and q denote the Radon-Nikodym derivatives with respect to µ).
We use the following conventions:
 
0
= 0,
0f
0

f (0) = lim f (u),u →0+

∀ a > 0, 0 fa
0= lim u f
u →0+

a
u

= a lim

u→∞

f (u)
.
u

(A2)

Entropy 2020, 22, 1100

49 of 61

When f (u) = − log u, we retrieve the Kullback–Leibler divergence (KLD):DKL ( p : q) =Z

p( x ) logp( x )dµ( x ).
q( x )(A3)

The KLD is usually difficult to calculate in closed-form, say, for example, between statisticalmixture models [155]. A common technique is to estimate the KLD using Monte Carlo sampling using
a proposal distribution r:nc n ( p : q) = 1 ∑ p( xi ) log p( xi ) ,KL(A4)n i =1 r ( x i )q ( xi )where x1 , . . . , xn ∼iid r. When r is chosen as p, the KLD can be estimated as:nc n ( p : q) = 1 ∑ log p( xi ) .
KLn i =1q ( xi )(A5)

c n ( p : q) = KL( p : q).
Monte Carlo estimators are consistent under mild conditions: limn→∞ KL
In practice, one problem when implementing Equation (A5), is that we may end up potentiallyc n ( p : q) < 0. This may have disastrous consequences as algorithms implemented by programswith KL
consider non-negative divergences to execute a correct workflow. The potential negative value problemof Equation (A5) comes from the fact that ∑i p( xi ) 6= 1 and ∑i q( xi ) 6= 1. One way to circumvent thisproblem is to consider the extended f -divergences:
Definition A1 (Extended f -divergence). The extended f -divergence for a convex generator f , strictly convex
at 1 and satisfying f (1) = 0 is defined by
I ef ( p : q) =

Z




 
q( x )q( x )0− f (1)−1dµ( x ).
p( x ) fp( x )p( x )(A6)

Indeed, for a strictly convex generator f , let us consider the scalar Bregman divergence [91]:B f ( a : b) = f ( a) − f (b) − ( a − b) f 0 (b) ≥ 0.

(A7)

q( x )

Setting a = p( x) and b = 1 in Equation (A7), and using the fact that f (1) = 0, we getfq( x )p( x )−q( x )−1p( x )

f 0 (1) ≥ 0.

(A8)

q( x ): 1 dµ( x ) ≥ 0.
p( x )(A9)

Therefore we define the extended f -divergences as
I ef ( p : q) =Zp( x ) B f

That is, the formula for the extended f -divergences is
I ef ( p : q) =Z

p( x ) q( x )q( x )f− f 0 (1)−1dµ( x ) ≥ 0.
p( x )p( x )(A10)

Then we estimate the extended f -divergence using importance sampling of the integral withrespect to distribution r, using n variates x1 , . . . , xn ∼iid p as:Î f ,n ( p : q) =

1 n
f
n i∑
=1



q ( xi )p ( xi )− f 0 (1)q ( xi )−1
p ( xi )



≥ 0.

(A11)

Entropy 2020, 22, 1100

50 of 61For example, for the KLD, we obtain the following Monte Carlo estimator:p ( xi )1 n
q ( xi )c
KLn ( p : q) = ∑ log+− 1 ≥ 0,n i =1q ( xi )p ( xi )(A12)

since the extended KLD is
DKLe ( p : q) =Z p( x ) logp( x )+ q( x ) − p( x ) dµ( x ).
q( x )(A13)

Equation (A12) can be interpreted as a sum of scalar Itakura–Saito divergences since thec n ( p : q) = 1 ∑n DIS ( p( xi ) : q( xi )) with the scalar
Itakura–Saito divergence is scale-invariant: KL
n i =1
Itakura–Saito divergence
DIS ( a : b) = DIS

a


a
a
: 1 = − log − 1 ≥ 0,
b
b
b

(A14)

a Bregman divergence obtained for the generator f (u) = − log u.
Notice that the extended f -divergence is a f -divergence for the generator
f e (u) = f (u) − f 0 (1)(u − 1).

(A15)

We check that the generator f e satisfies both f (1) = 0 and f 0 (1) = 0, and we have I ef ( p : q) = I f e ( p : q).
ee ( p : q ) with fThus DKLe ( p : q) = I fKL
KL ( u ) = − log u + u − 1.
Let us remark that we only need to have the scalar function strictly convex at 1 to ensure that
B f ba : 1 ≥ 0. Indeed, we may use the definition of Bregman divergences extended to strictly convex
functions but not necessarily smooth functions [156,157]:
B f ( x : y) =max

g(y)∈∂ f (y){ f ( x ) − f (y) − ( x − y) g(y)},(A16)

where ∂ f (y) denotes the subderivative of f at y.
Furthermore, noticing that Iλ f ( p : q) = λI f ( p : q) for λ > 0, we may enforce that f 00 (1) = 1,
and obtain a standard f -divergence [2] which enjoys the property that I f ( pθ ( x ) : pθ +dθ ( x )) =dθ > I (θ )dθ, where I (θ ) denotes the Fisher information matrix of the parameteric family { pθ }θ
of densities.
Appendix B. The Multivariate Gaussian Family: An Exponential Family
We report the canonical decomposition of the multivariate Gaussian [158] family{ N (µ, Σ) such that µ ∈ Rd , Σ  0} following [159]. The multivariate Gaussian family is also called theMultiVariate Normal family, or MVN family for short.
Let λ := (λv , λ M ) = (µ, Σ) denote the composite (vector,matrix) parameter of an MVN.
The d-dimensional MVN density is given by


1
1
> −1
pλ ( x; λ) :=exp − ( x − λv ) λ M ( x − λv ) ,
(A17)
dp
2
(2π ) 2 |λ M |where | · | denotes the matrix determinant. The natural parameters θ are also expressed using both avector parameter θv and a matrix parameter θ M in a compound object θ = (θv , θ M ). By defining thefollowing compound inner product on a composite (vector,matrix) object>hθ, θ 0 i := θv> θv0 + tr θ 0M θ M ,(A18)

Entropy 2020, 22, 1100

51 of 61

where tr(·) denotes the matrix trace, we rewrite the MVN density of Equation (A17) in the canonicalform of an exponential family [52]:pθ ( x; θ )exp (ht( x ), θ i − Fθ (θ )) = pλ ( x; λ(θ )),

:=

where


θ = (θv , θ M ) =

Σ

−1

1
µ, − Σ−12= θ (λ) =

(A19)

1 −1
1
λ−
M λv , − 2 λ M


,

(A20)

is the compound natural parameter andt( x ) = ( x, − xx > )(A21)

is the compound sufficient statistic. The function Fθ is the strictly convex and continuously
differentiable log-normalizer defined by:
1Fθ (θ ) =2


1 > −1d log π − log |θ M | + θv θ M θv ,
2

(A22)

The log-normalizer can be expressed using the ordinary parameters, λ = (µ, Σ), as:Fλ (λ)


1  > −1λv λ M λv + log |λ M | + d log 2π ,2
1  > −1µ Σ µ + log |Σ| + d log 2π .
2

=
=

(A23)
(A24)

The moment/expectation parameters [2] areη = (ηv , η M ) = E[t( x )] = ∇ F (θ ).

(A25)

We report the conversion formula between the three types of coordinate systems (namely theordinary parameter λ, the natural parameter θ and the moment parameter η) as follows:
(

(

1
−1
θv (λ) = λ−M λv = Σ µ
1 −1
1 −1
θ M (λ) = 2 λ M = 2 Σ1
ηv (θ ) = 12 θ −
M θv
1
−1
1 −1
>
η M (θ ) = − 2 θ M − 41 (θ −M θv )( θ M θv )(λ v ( η ) = ηv = µ
λ M (η ) = −η M − ηv ηv> = Σ

(

⇔
(

⇔
(

⇔

1
λv (θ ) = 12 θ −M θv = µ1
λ M (θ ) = 12 θ −
M =Σ

(A26)

θv (η ) = −(η M + ηv ηv> )−1 ηv
θ M (η ) = − 12 (η M + ηv ηv> )−1(A27)

ηv ( λ ) = λ v = µ> (A28)
η M (λ) = −λ M − λv λ>v = − Σ − µµ

The dual Legendre convex conjugate [2] isFη∗ (η ) = −1
−1log(1 + ηv> η Mηv ) + log | − η M | + d(1 + log 2π ) ,
2

(A29)

and θ = ∇η Fη∗ (η ).
We check the Fenchel-Young equality when η = ∇ F (θ ) and θ = ∇ F ∗ (η ):Fθ (θ ) + Fη∗ (η ) − hθ, η i = 0.

(A30)

The Kullback-Leibler divergence between two d-dimensional Gaussians distributions p(µ1 ,Σ1 ) andp(µ2 ,Σ2 ) (with ∆µ = µ2 − µ1 ) isKL( p(µ1 ,Σ1 ) : p(µ2 ,Σ2 ) )

=

1
2



| Σ2 |−1tr(Σ2−1 Σ1 ) + ∆>−d
µ Σ2 ∆µ + log| Σ1 |

= KL( pλ1 : pλ2 ). (A31)Entropy 2020, 22, 1100

52 of 61

We check that KL( p(µ,Σ) : p(µ,Σ) ) = 0 since ∆µ = 0 and tr(Σ−1 Σ) = tr( I ) = d. Notice that when
Σ1 = Σ2 = Σ, we haveKL( p(µ1 ,Σ) : p(µ2 ,Σ) ) =1 > −11∆ Σ ∆µ = DΣ2 −1 (µ1 , µ2 ),
2 µ
2

(A32)

that is half the squared Mahalanobis distance for the precision matrix Σ−1 (a positive-definite matrix:
Σ−1  0), where the Mahalanobis distance is defined for any positive matrix Q  0 as follows:
DQ ( p1 : p2 ) =q( p1 − p2 ) > Q ( p1 − p2 ).

(A33)

The Kullback–Leibler divergence between two probability densities of the same exponential
families amount to a Bregman divergence [2]:KL( p(µ1 ,Σ1 ) : p(µ2 ,Σ2 ) ) = KL( pλ1 : pλ2 ) = BF (θ2 : θ1 ) = BF∗ (η1 : η2 ),(A34)

where the Bregman divergence is defined by
BF (θ : θ 0 ) := F (θ ) − F (θ 0 ) − hθ − θ 0 , ∇ F (θ 0 )i,(A35)with η 0 = ∇ F (θ 0 ). Define the canonical divergence [2]A F ( θ 1 : η2 ) = F ( θ 1 ) + F ∗ ( η2 ) − h θ 1 , η2 i = A F ∗ ( η2 : θ 1 ) ,(A36)

since F ∗ ∗ = F. We have BF (θ1 : θ2 ) = A F (θ1 : η2 ).
Appendix C. Skew Jensen Divergences and Bregman Divergences
First, let us interpret geometrically the skewed Jensen divergences and the Bregman divergencesusing the chordal slope lemma of convexity theory [160] for univariate generators F:
Lemma A1 (Chordal slope lemma). Let F be a strictly convex function on interval I = ( a, b) for a < b.
For any θ1 < θ < θ2 in I, we have:F ( θ ) − F ( θ1 )F ( θ2 ) − F ( θ1 )F ( θ2 ) − F ( θ )<<.
θ − θ1θ2 − θ1θ2 − θ

(A37)

That is, define the points P1 = (θ1 , F (θ1 )), P = (θ, F (θ )), and P2 = (θ2 , F (θ2 )). Then the chordal
slope lemma states that the slope of the chord [ P1 P] is less than the slope of the chord [ P1 P2 ], which isless than the slope of [ PP2 ] (Figure A1).
slope(TP1 ) ≤ slope([P1 P ]) < slope([P1 P2 ]) < slope([P P2 ]) ≤ slope(TP2 )0

F (θ1 )

(θ1 )
(θ1 )
(θ)
F 0 (θ1 ) ≤ F (θ)−F< F (θθ22)−F< F (θθ22)−F≤ F 0 (θ2 )
θ−θ1
−θ1
−θ

F 0 (θ2 )F
P2 = (θ2 , F (θ2 ))P1 = (θ1 , F (θ1 ))

P = (θ, F (θ))

θ1

θ

θ2

Figure A1. Illustration of the chordal slope lemma.

Entropy 2020, 22, 1100

53 of 61Let θ = (1 − α)θ1 + αθ2 for α ∈ (0, 1). Then the inequalities of the chordal slope lemma can be
rewritten as:
F ( θ ) − F ( θ1 )α ( θ2 − θ1 )F ( θ2 ) − F ( θ1 )( θ2 − θ1 )F ( θ2 ) − F ( θ1 ),( θ2 − θ1 )F ( θ2 ) − F ( θ ),(1 − α)(θ2 − θ1 )

<
<

(A38)
(A39)

since θ − θ1 = α(θ2 − θ1 ) and θ2 − θ = (1 − α)(θ2 − θ1 ). Thus we have:F ( θ ) − F ( θ1 )< α( F (θ2 ) − F (θ1 )),(A40)

F ( θ2 ) − F ( θ ),(A41)

(1 − α)( F (θ2 ) − F (θ1 )) <or equivalently the following inequalities:α( F (θ2 ) − F (θ1 )) − F (θ ) + F (θ1 )
0

> 0,
<

(A42)

F (θ2 ) − F (θ ) + (1 − α)( F (θ1 ) − F (θ2 )).

(A43)

That is, we recover for θ1 < (1 − α)θ1 + αθ2 < θ2 (i.e., α ∈ (0, 1) and θ1 6= θ2 ) the α-skewed Jensen
divergence [102]:JFα (θ1 : θ2 ) := (1 − α) F (θ1 ) + αF (θ2 ) − F ((1 − α)θ1 + αθ2 ) > 0.

(A44)

Now, a consequence of the chordal slope lemma is that for a strictly convex and differentiable
function F, we have:F ( θ2 ) − F ( θ1 )F 0 ( θ1 ) ≤≤ F 0 ( θ2 ).
(A45)
θ2 − θ1
This can be geometrically visualized in Figure A1.
That is,F ( θ2 ) − F ( θ1 ) − ( θ2 − θ1 ) F 0 ( θ1 )0

F ( θ2 ) − F ( θ1 ) − ( θ2 − θ1 ) F ( θ2 )

≥ 0,

(A46)

≤ 0.

(A47)

We recognize the expressions of the Bregman divergences [91]:B F ( θ1 : θ2 ) : = F ( θ1 ) − F ( θ2 ) − ( θ1 − θ2 ) F 0 ( θ2 ),(A48)

and get:BF ( θ2 : θ1 )

≥ 0,

(A49)

BF ( θ1 : θ2 )

≥ 0.

(A50)

For multivariate strictly convex functions F, we observe that we can build a multivariate Bregman
divergence from a family of 1D Bregman divergences [161] induced by the 1D strictly convex functions
Fθ1 ,θ2 ((1 − λ)θ1 + θ2 ):Theorem A1. A multivariate Bregman divergence between two parameters θ1 and θ2 can be expressed as a
univariate Bregman divergence for the generator Fθ1 ,θ2 induced by the parameters:
BF (θ1 : θ2 ) = BFθ ,θ (0 : 1),
1 2∀ θ1 6 = θ2 ,

(A51)

where
Fθ1 ,θ2 (u) := F (θ1 + u(θ2 − θ1 )).

(A52)Entropy 2020, 22, 1100

54 of 61

Proof. The functions { Fθ1 ,θ2 }(θ1 ,θ2 )∈Θ2 are 1D Bregman generators. Consider the directional derivative:∇θ2 −θ1 Fθ1 ,θ2 (u) :==F (θ1 + (e + u)(θ2 − θ1 )) − F (θ1 + u(θ2 − θ1 )),ee →0(θ2 − θ1 )> ∇ F (θ1 + u(θ2 − θ1 )),
lim

(A53)
(A54)

Since Fθ1 ,θ2 (0) = F (θ1 ), Fθ1 ,θ2 (1) = F (θ2 ), and Fθ01 ,θ2 (u) = ∇θ2 −θ1 Fθ1 ,θ2 (u), it follows that
BFθ ,θ (0 : 1)
1 2=

Fθ1 ,θ2 (0) − Fθ1 ,θ2 (1) − (0 − 1)∇θ2 −θ1 Fθ1 ,θ2 (1),>=F ( θ1 ) − F ( θ2 ) + ( θ2 − θ1 ) ∇ F ( θ2 ),

(A56)

=

B F ( θ1 : θ2 ).

(A57)

Notations
Below is a list of notations we used in this document:[D]h·, ·i[ D ] := {1, . . . , D }inner productMQ (u, v) = ku − vkQ

Mahalanobis distance MQ (u, v) =D (θ : θ 0 )parameter divergencestatistical divergence
Divergence and dual (reverse)
divergenceD [ p( x ) : p0 ( x )]D, D ∗
Csiszár divergence I f
Bregman divergence BF
Canonical divergence A F,F∗
Bhattacharyya distance
Jensen/Burbea-Rao divergence
Chernoff information
F, F ∗Dρ ( p, q)
B, B∗
B = { e1 = ∂ 1 , . . . , e D = ∂ D }{dxi }i( v ) B : = ( vi )( v ) B∗ : = ( vi )u⊥vp
kvk = hv, vi
M, S
Tp
TM
F( M )X( M ) = Γ( TM)vf
X, Y, Z ∈ X( M)Σg = gij dxi ⊗ dx j
(U , x )(A55)

q∑i,j (ui − vi )(u j − v j ) Qij , Q  0

θi0
θi with f (1) = 0
0
BF ( θ : θ ) : = F ( θ ) − F ( θ 0 ) − ( θ − θ 0 )> ∇ F ( θ 0 )
A F,F∗ (θ : η 0 ) = F (θ ) + F ∗ (η 0 ) − θ > η 0RBα [ p1 : p2 ] = − log x∈X p1α ( x ) p12−α ( x )dµ( x )(α)
JF (θ1 : θ2 ) = αF (θ1 ) + (1 − α) F (θ2 ) − F (θ1 + (1 − α)θ2 )R
C [ P1 , P2 ] = − log minα∈(0,1) x∈X p1α ( x ) p12−α ( x )dµ( x )I f (θ : θ 0 ) := ∑iD=1 θi f

Potential functions related by Legendre–Fenchel transformation
R1
Riemannian distance Dρ ( p, q) := 0 kγ0 (t)kγ(t) dt
basis, reciprocal basis
natural basis
covector basis (one-forms)
contravariant components of vector v
covariant components of vector v
vector u is perpendicular to vector v (hu, vi = 0)
induced norm, length of a vector v
Manifold, submanifold
tangent plane at p
Tangent bundle TM = ∪ p Tp = {( p, v), p ∈ M, v ∈ Tp }space of smooth functions on M
space of smooth vector fields on M
direction derivative of f with respect to vector v
Vector fields
metric tensor (field)
local coordinates x in a chat U

Entropy 2020, 22, 1100

55 of 61∂i :=: ∂x∂i∂i :=: ∂x∂ i
∇
∇X Y
∏∇
c
∏∇
c v
γ, γ∇
Γij,l
Γijk

natural basis vector
natural reciprocal basis vector
affine connection
covariant derivative
parallel transport of vectors along a smooth curve cParallel transport of v ∈ Tc(0) along a smooth curve c
geodesic, geodesic with respect to connection ∇Christoffel symbols of the first kind (functions)
Christoffel symbols of the second kind (functions)R[ X, Y ]∇-projection∇∗ -projectionC
P = { pθ ( x )}θi nΘ
E , M, ∆ D
P I (θ )P I (θ )Pg

Riemann–Christoffel curvature tensor
Lie bracket [ X, Y ]( f ) = X (Y ( f )) − Y ( X ( f )), ∀ f ∈ F( M)
PS = arg minQ∈S D (θ ( P) : θ ( Q))PS∗ = arg minQ∈S D (θ ( Q) : θ ( P))Amari–Chentsov totally symmetric cubic 3-covariant tensor
parametric family of probability distributions
exponential family, mixture family, probability simplex
Fisher information matrix of family P
Fisher Information Matrix (FIM) for a parametric family P
Fisher information
metrichi tensor field

exponential connection eP ∇mixture connection mP∇e ∇ : = E ( ∂ ∂ l )( ∂ l )
θ
k
P
h i j
i
m ∇ : = E ( ∂ ∂ l + ∂ l∂ l )( ∂ l )
θ
i j
i j
k
P

expected skewness tensor Cijk

hi
Cijk := Eθ ∂i l∂ j l∂k lexpected α-connections

1+ ααk
P Γ ij : = − 2 Cijk = Eθ≡equivalence of geometric structureshi∂i ∂ j l + 1−2 α ∂i l∂ j l (∂k l )References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.

Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 623–656. [CrossRef]
Amari, S. Information Geometry and Its Applications; Applied Mathematical Sciences; Springer: Tokyo,
Japan, 2016.
Kakihara, S.; Ohara, A.; Tsuchiya, T. Information Geometry and Interior-Point Algorithms in Semidefinite
Programs and Symmetric Cone Programs. J. Optim. Theory Appl. 2013, 157, 749–780. [CrossRef]
Amari, S.; Nagaoka, H. Methods of Information Geometry; American Mathematical Society: Providence, RL,
USA 2007.
Peirce, C.S. Chance, Love, and Logic: Philosophical Essays; U of Nebraska Press: Lincoln, NE, USA, 1998.
Schurz, G. Patterns of abduction. Synthese 2008, 164, 201–234. [CrossRef]
Wald, A. Statistical decision functions. Ann. Math. Stat. 1949, 165–205. [CrossRef]
Wald, A. Statistical Decision Functions; Wiley: Chichester, UK, 1950.
Dabak, A.G. A Geometry for Detection Theory. Ph.D. Thesis, Rice University, Houston, TX, USA, 1993.
Do Carmo, M.P. Differential Geometry of Curves and Surfaces; Courier Dover Publications: New York, NY,
USA, 2016.
Amari, S.; Barndorff-Nielsen, O.E.; Kass, R.E.; Lauritzen, S.L.; Rao, C.R. Differential Geometry in Statistical
Inference; Institute of Mathematical Statistics: Hayward, CA, USA, 1987.
Dodson, C.T.J. (Ed.) Geometrization of Statistical Theory; ULDM Publications; University of Lancaster,
Department of Mathematics: Bailrigg, UK, 1987.
Murray, M.; Rice, J. Differential Geometry and Statistics; Number 48 in Monographs on Statistics and Applied
Probability; Chapman and Hall: London, UK, 1993.
Kass, R.E.; Vos, P.W. Geometrical Foundations of Asymptotic Inference; Wiley-Interscience: New York, NY,
USA, 1997.
Arwini, K.A.; Dodson, C.T.J. Information Geometry: Near Randomness and Near Independance; Springer: Berlin,
Germany, 2008.
Calin, O.; Udriste, C. Geometric Modeling in Probability and Statistics; Mathematics and Statistics, Springer
International Publishing: Cham, Switzerland, 2014.

Entropy 2020, 22, 1100

17.
18.
19.
20.

21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.

42.

43.
44.

56 of 61

Ay, N.; Jost, J.; Vân Lê, H.; Schwachhöfer, L. Information Geometry; Springer: Berlin, Germany, 2017;
Volume 64.
Corcuera, J.; Giummolè, F. A characterization of monotone and regular divergences. Ann. Inst. Stat. Math.
1998, 50, 433–450. [CrossRef]
Mühlich, U. Fundamentals of Tensor Calculus for Engineers with a Primer on Smooth Manifolds; Springer: Berlin,
Germany, 2017; Volume 230.
Nielsen, F.; Nock, R. Hyperbolic Voronoi diagrams made easy. In Proceedings of the IEEE International
Conference on Computational Science and Its Applications (ICCSA), Fukuoka, Japan, 23–26 March 2010;
pp. 74–80.
Whitney, H.; Eells, J.; Toledo, D. Collected Papers of Hassler Whitney; Nelson Thornes: London, UK, 1992.
Absil, P.A.; Mahony, R.; Sepulchre, R. Optimization Algorithms on Matrix Manifolds; Princeton University
Press: Princeton, NJ, USA, 2009.
Cartan, E.J. On Manifolds with an Affine Connection and the Theory of General Relativity; Bibliopolis; Humanities
Pr: London, UK, 1986.
Akivis, M.A.; Rosenfeld, B.A. Élie Cartan (1869–1951); American Mathematical Society: Cambridge, MA,
USA, 2011; Volume 123.
Wanas, M. Absolute parallelism geometry: Developments, applications and problems. arXiv 2002, arXiv:grqc/0209050.
Bourguignon, J.P. Ricci curvature and measures. Jpn. J. Math. 2009, 4, 27–45. [CrossRef]
Baez, J.C.; Wise, D.K. Teleparallel gravity as a higher gauge theory. Commun. Math. Phys. 2015, 333, 153–186.
[CrossRef]
Ashburner, J.; Friston, K.J. Diffeomorphic registration using geodesic shooting and Gauss-Newton
optimisation. NeuroImage 2011, 55, 954–967. [CrossRef]
Lauritzen, S.L. Statistical manifolds. Differ. Geom. Stat. Inference 1987, 10, 163–216.
Vân Lê, H. Statistical manifolds are statistical models. J. Geom. 2006, 84, 83–93.
Furuhata, H. Hypersurfaces in statistical manifolds. Differ. Geom. Its Appl. 2009, 27, 420–429. [CrossRef]
Zhang, J. Divergence functions and geometric structures they induce on a manifold. In Geometric Theory ofInformation; Nielsen, F., Ed.; Springer: Berlin, Germany, 2014; pp. 1–30.
Eguchi, S. Second order efficiency of minimum contrast estimators in a curved exponential family. Ann. Stat.
1983, 11, 793–803. [CrossRef]
Eguchi, S. A differential geometric approach to statistical inference on the basis of contrast functionals.
Hiroshima Math. J. 1985, 15, 341–391. [CrossRef]
Hiriart-Urruty, J.B.; Lemaréchal, C. Fundamentals of Convex Analysis; Springer Science & Business Media:
Berlin, Germany, 2012.
Crouzeix, J.P. A relationship between the second derivatives of a convex function and of its conjugate.
Math. Program. 1977, 13, 364–365. [CrossRef]
Ay, N.; Amari, S. A novel approach to canonical divergences within information geometry. Entropy 2015,
17, 8111–8129. [CrossRef]
Nielsen, F. What is ... an information projection? Not. AMS 2018, 65, 321–324. [CrossRef]
Kurose, T. On the divergences of 1-conformally flat statistical manifolds. Tohoku Math. J. Second Ser. 1994,
46, 427–433. [CrossRef]
Boissonnat, J.D.; Nielsen, F.; Nock, R. Bregman Voronoi diagrams. Discret. Comput. Geom. 2010, 44, 281–307.
[CrossRef]
Nielsen, F.; Piro, P.; Barlaud, M. Bregman vantage point trees for efficient nearest neighbor queries.
In Proceedings of the 2009 IEEE International Conference on Multimedia and Expo, New York, NY, USA,
28 June–3 July 2009; pp. 878–881.
Nielsen, F.; Boissonnat, J.D.; Nock, R. Visualizing Bregman Voronoi diagrams. In Proceedings of
the Twenty-Third Annual Symposium on Computational Geometry, Gyeongju, Korea, 6–8 June 2007;
pp. 121–122.
Nock, R.; Nielsen, F. Fitting the smallest enclosing Bregman ball. In Proceedings of the European Conference
on Machine Learning, Porto, Portugal, 3–7 October 2005; Springer: Berlin, Germany, 2005; pp. 649–656.
Nielsen, F.; Nock, R. On the smallest enclosing information disk. Inf. Process. Lett. 2008, 105, 93–97.
[CrossRef]

Entropy 2020, 22, 1100

45.

46.
47.
48.
49.
50.

51.
52.
53.
54.
55.

56.

57.
58.
59.

60.
61.
62.
63.
64.
65.
66.
67.
68.
69.

57 of 61

Fischer, K.; Gärtner, B.; Kutz, M. Fast smallest-enclosing-ball computation in high dimensions. In Proceedingsof the European Symposium on Algorithms, Budapest, Hungary, 16–19 September 2003; Springer: Berlin,
Germany, 2003; pp. 630–641.
Della Pietra, S.; Della Pietra, V.; Lafferty, J. Inducing features of random fields. IEEE Trans. Pattern Anal.
Mach. Intell. 1997, 19, 380–393. [CrossRef]
Nielsen, F. On Voronoi Diagrams on the Information-Geometric Cauchy Manifolds. Entropy 2020, 22, 713.
[CrossRef]
Shima, H. The Geometry of Hessian Structures; World Scientific: New Jersey, NJ, USA, 2007.
Zhang, J. Reference duality and representation duality in information geometry. AIP Conf. Proc. 2015,
1641, 130–146.
Gomes-Gonçalves, E.; Gzyl, H.; Nielsen, F. Geometry and Fixed-Rate Quantization in Riemannian MetricSpaces Induced by Separable Bregman Divergences. In Proceedings of the 4th International Conference onGeometric Science of Information (GSI), Toulouse, France, 27–29 August 2019; Nielsen, F., Barbaresco, F.,
Eds.; Lecture Notes in Computer Science; Springer: Berlin, Germany, 2019; Volume 11712, pp. 351–358.
[CrossRef]
Nielsen, F. Cramér-Rao lower bound and information geometry. In Connected at Infinity II; Springer: Berlin,
Germany, 2013; pp. 18–37.
Nielsen, F.; Garcia, V. Statistical exponential families: A digest with flash cards. arXiv 2009, arXiv:0911.4863.
Sato, Y.; Sugawa, K.; Kawaguchi, M. The geometrical structure of the parameter space of the two-dimensional
normal distribution. Rep. Math. Phys. 1979, 16, 111–119. [CrossRef]
Skovgaard, L.T. A Riemannian geometry of the multivariate normal model. Scand. J. Stat. 1984, 11, 211–223.
Malagò, L.; Pistone, G. Information geometry of the Gaussian distribution in view of stochastic optimization.
In Proceedings of the 2015 ACM Conference on Foundations of Genetic Algorithms XIII, Aberystwyth, UK,
17–20 January 2015; pp. 150–162.
Nielsen, F.; Nock, R. On the geometry of mixtures of prescribed distributions. In Proceedings of the 2018
IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, AB, Canada,
15–20 April 2018; pp. 2861–2865.
Campbell, L.L. An extended Čencov characterization of the information metric. Proc. Am. Math. Soc. 1986,
98, 135–141.
Vân Lê, H. The uniqueness of the Fisher metric as information metric. Ann. Inst. Stat. Math. 2017, 69,
879–896.
Csiszár, I.; Shields, P.C. Information Theory and Statistics: A Tutorial; Foundations and Trends R in
Communications and Information Theory; Now Publishers Inc.: Hanover, MA, USA, 2004; Volume 1,
pp. 417–528.
Jiao, J.; Courtade, T.A.; No, A.; Venkat, K.; Weissman, T. Information measures: The curious case of thebinary alphabet. IEEE Trans. Inf. Theory 2014, 60, 7616–7626. [CrossRef]
Qiao, Y.; Minematsu, N. A Study on Invariance of f -Divergence and Its Application to Speech Recognition.
IEEE Trans. Signal Process. 2010, 58, 3884–3890. [CrossRef]
Nielsen, F.; Nock, R. On the chi square and higher-order chi distances for approximating f -divergences.
IEEE Signal Process. Lett. 2013, 21, 10–13. [CrossRef]
Csiszár, I. Information-type measures of difference of probability distributions and indirect observation.
Stud. Sci. Math. Hung. 1967, 2, 229–318.
Mitchell, A.F.S. Statistical manifolds of univariate elliptic distributions. Int. Stat. Rev. 1988, 56, 1–16.
[CrossRef]
Hotelling, H. Spaces of statistical parameters. Bull. Am. Math. Soc. (AMS) 1930, 36, 191.
Rao, R.C. Information and the accuracy attainable in the estimation of statistical parameters. Bull. Calcutta
Math. Soc. 1945, 37, 81–91.
Komaki, F. Bayesian prediction based on a class of shrinkage priors for location-scale models. Ann. Inst.
Stat. Math. 2007, 59, 135–146. [CrossRef]
Stigler, S.M. The epic story of maximum likelihood. Stat. Sci. 2007, 22, 598–620. [CrossRef]
Rao, C.R. Information and the accuracy attainable in the estimation of statistical parameters. In Breakthroughsin Statistics; Springer: Berlin, Germany, 1992; pp. 235–247.

Entropy 2020, 22, 1100

70.
71.
72.
73.
74.
75.
76.
77.
78.
79.

80.
81.
82.
83.

84.
85.

86.
87.
88.
89.
90.
91.
92.
93.
94.
95.
96.

58 of 61

Jeffreys, H. An invariant form for the prior probability in estimation problems. Proc. R. Soc. Lond. A 1946,
186, 453–461.
Zhang, J. On monotone embedding in information geometry. Entropy 2015, 17, 4485–4499. [CrossRef]
Naudts, J.; Zhang, J. Rho–tau embedding and gauge freedom in information geometry. Inf. Geom. 2018.
[CrossRef]
Nock, R.; Nielsen, F.; Amari, S. On Conformal Divergences and Their Population Minimizers. IEEE TIT
2016, 62, 527–538. [CrossRef]
Azoury, K.S.; Warmuth, M.K. Relative loss bounds for on-line density estimation with the exponential familyof distributions. Mach. Learn. 2001, 43, 211–246. [CrossRef]
Banerjee, A.; Merugu, S.; Dhillon, I.S.; Ghosh, J. Clustering with Bregman divergences. J. Mach. Learn. Res.
2005, 6, 1705–1749.
Nielsen, F.; Nock, R. Entropies and cross-entropies of exponential families. In Proceedings of the 2010 IEEE
International Conference on Image Processing, Hong Kong, China, 26–29 September 2010; pp. 3621–3624.
Nielsen, F.; Nock, R. Cumulant-free closed-form formulas for some common (dis)similarities between
densities of an exponential family. Tech. Rep. 2020. [CrossRef]
Amari, S. Differential geometry of a parametric family of invertible linear systems: Riemannian metric, dual
affine connections, and divergence. Math. Syst. Theory 1987, 20, 53–82. [CrossRef]
Schwander, O.; Nielsen, F. Fast learning of Gamma mixture models with k-MLE. In Proceedings of theInternational Workshop on Similarity-Based Pattern Recognition, York, UK, 3–5 July 2013; Springer: Berlin,
Germany, 2013; pp. 235–249.
Miura, K. An introduction to maximum likelihood estimation and information geometry. Interdiscip. Inf. Sci.
2011, 17, 155–174. [CrossRef]
Reverter, F.; Oller, J.M. Computing the Rao distance for Gamma distributions. J. Comput. Appl. Math. 2003,
157, 155–167. [CrossRef]
Pinele, J.; Strapasson, J.E.; Costa, S.I. The Fisher-Rao Distance between Multivariate Normal Distributions:
Special Cases, Bounds and Applications. Entropy 2020, 22, 404. [CrossRef]
Nielsen, F. Pattern learning and recognition on statistical manifolds: An information-geometric review.
In Proceedings of the International Workshop on Similarity-Based Pattern Recognition, York, UK, 3–5 July
2013; Springer: Berlin, Germany, 2013; pp. 1–25.
Sun, K.; Nielsen, F. Lightlike Neuromanifolds, Occam’s Razor and Deep Learning.
arXiv 2019,arXiv:1905.11027.
Sun, K.; Nielsen, F. Relative Fisher Information and Natural Gradient for Learning Large Modular Models.
In Proceedings of the 34th International Conference on Machine Learning (ICML), Sydney, Australia,
6–11 August 2017; pp. 3289–3298.
Amari, S. Natural gradient works efficiently in learning. Neural Comput. 1998, 10, 251–276. [CrossRef]
Cauchy, A. Methode générale pour la résolution des systèmes d’équations simultanées. C. R. l’Académie Sci.
1847, 25, 536–538.
Curry, H.B. The method of steepest descent for non-linear minimization problems. Q. Appl. Math. 1944,
2, 258–261. [CrossRef]
Bonnabel, S. Stochastic gradient descent on Riemannian manifolds. IEEE Trans. Autom. Control 2013,
58, 2217–2229. [CrossRef]
Nielsen, F. On geodesic triangles with right angles in a dually flat space. arXiv 2019, arXiv:1910.03935.
Bregman, L.M. The relaxation method of finding the common point of convex sets and its application to thesolution of problems in convex programming. USSR Comput. Math. Math. Phys. 1967, 7, 200–217. [CrossRef]
Nielsen, F.; Hadjeres, G. Monte Carlo information-geometric structures. In Geometric Structures of Information;
Springer, Berlin, Germany, 2019; pp. 69–103.
Nielsen, F. Legendre Transformation and Information Geometry; Springer, Berlin/Heidelberg, Germany, 2010.
Raskutti, G.; Mukherjee, S. The information geometry of mirror descent. IEEE Trans. Inf. Theory 2015,
61, 1451–1457. [CrossRef]
Bubeck, S. Convex Optimization: Algorithms and Complexity; Foundations and Trends R in Machine Learning:
Hanover, MA, USA, 2015; Volume 8, pp. 231–357.
Zhang, G.; Sun, S.; Duvenaud, D.; Grosse, R. Noisy natural gradient as variational inference. In Proceedingsof the International Conference on Machine Learning, Stockholm, Sweden, 10–15 July 2018; pp. 5852–5861.

Entropy 2020, 22, 1100

59 of 61

Beyer, H.G.; Schwefel, H.P. Evolution strategies–A comprehensive introduction. Nat. Comput. 2002, 1, 3–52.
[CrossRef]
98. Berny, A. Selection and reinforcement learning for combinatorial optimization. In Proceedings of theInternational Conference on Parallel Problem Solving from Nature, Paris, France, 18–20 September 2000;
Springer: Berlin, Germany, 2000; pp. 601–610.
99. Wierstra, D.; Schaul, T.; Glasmachers, T.; Sun, Y.; Peters, J.; Schmidhuber, J. Natural evolution strategies.
J. Mach. Learn. Res. 2014, 15, 949–980.
100. Nielsen, F. An Information-Geometric Characterization of Chernoff Information. IEEE Sig. Proc. Lett. 2013,20, 269–272. [CrossRef]
101. Pham, G.; Boyer, R.; Nielsen, F. Computational Information Geometry for Binary Classification ofHigh-Dimensional Random Tensors. Entropy 2018, 20, 203. [CrossRef]
102. Nielsen, F.; Boltz, S. The Burbea-Rao and Bhattacharyya centroids. IEEE Trans. Inf. Theory 2011, 57, 5455–5466.
[CrossRef]
103. Nielsen, F. Chernoff Information of Exponential Families. arXiv 2011 arXiv:1102.2684.
104. Nielsen, F. Generalized Bhattacharyya and Chernoff upper bounds on Bayes error using quasi-arithmetic
means. Pattern Recognit. Lett. 2014, 42, 25–34. [CrossRef]
105. Nielsen, F. Hypothesis Testing, Information Divergence and Computational Geometry. In Proceedings ofthe International Conference on Geometric Science of Information Geometric Science of Information (GSI),
Paris, France, 28–30 August 2013; pp. 241–248.
106. Nielsen, F.; Sun, K. Guaranteed Bounds on Information-Theoretic Measures of Univariate Mixtures Using
Piecewise Log-Sum-Exp Inequalities. Entropy 2016, 18, 442. [CrossRef]
107. Nielsen, F.; Nock, R. Sided and symmetrized Bregman centroids. IEEE Trans. Inf. Theory 2009, 55, 2882–2904.
[CrossRef]
arXiv 2018,108. Nielsen, F.; Hadjeres, G. Monte Carlo Information Geometry: The dually flat case.
arXiv:1803.07225.
109. Ohara, A.; Tsuchiya, T. An Information Geometric Approach to Polynomial-Time Interior-Point Algorithms:
Complexity Bound via Curvature Integral; Research Memorandum; The Institute of Statistical Mathematics:
Tokyo, Japan, 2007; Volume 1055.
110. Fuglede, B.; Topsøe, F. Jensen-Shannon divergence and Hilbert space embedding. In Proceedings of theIEEE International Symposium on Information Theory (ISIT), Chicago, IL, USA, 27 June–2 July 2004; p. 31.
111. Vajda, I. On metric divergences of probability measures. Kybernetika 2009, 45, 885–900.
112. Villani, C. Optimal Transport: Old and New; Springer Science & Business Media: Berlin, Germany, 2008;
Volume 338.
113. Dowson, D.C.; Landau, B.V. The Fréchet distance between multivariate normal distributions. J. Multivar. Anal.
1982, 12, 450–455. [CrossRef]
114. Takatsu, A. Wasserstein geometry of Gaussian measures. Osaka J. Math. 2011, 48, 1005–1026.
115. Chentsov, N.N. Statistical Decision Rules and Optimal Inference; Monographs; American Mathematical Society:
Providence, RI, USA, 1982.
116. Amari, S. Differential-Geometrical Methods in Statistics; Lecture Notes on Statistics; Second Edition in 1990;
Springer: New York, NY, USA, 1985; Volume 28.
117. Amari, S.; Nagaoka, H. Methods of Information Geometry; Jouhou kika no Houhou; Iwanami Shoten: Tokyo,
Japan, 1993. (In Japanese)
118. Gibilisco, P.; Riccomagno, E.; Rogantin, M.P.; Wynn, H.P. (Eds.) Algebraic and Geometric Methods in Statistics;
Cambridge University Press: Cambridge, UK, 2009.
119. Srivastava, A.; Wu, W.; Kurtek, S.; Klassen, E.; Marron, J.S. Registration of Functional Data Using Fisher-Rao
Metric. arXiv 2011, arXiv:1103.3817.
120. Wei, S.W.; Liu, Y.X.; Mann, R.B. Ruppeiner geometry, phase transitions, and the microstructure of charged
AdS black holes. Phys. Rev. D 2019, 100, 124033. [CrossRef]
121. Quevedo, H. Geometrothermodynamics. J. Math. Phys. 2007, 48, 013506. [CrossRef]
122. Amari, S. Theory of Information Spaces: A Differential Geometrical Foundation of Statistics; Post RAAG Reports;
Tokyo, Japan, 1980. Available online: Available online: https://bsi-ni.brain.riken.jp/database/item/92(accessed on 8 May 2012). (accessed on 29 September 2020).

97.

Entropy 2020, 22, 1100

60 of 61

123. Efron, B. Defining the curvature of a statistical problem (with applications to second order efficiency).
Ann. Stat. 1975, 3, 1189–1242. [CrossRef]
124. Nagaoka, H.; Amari, S. Differential Geometry of Smooth Families of Probability Distributions; Technical Report;
METR 82-7; University of Tokyo: Tokyo, Japan, 1982.
125. Croll, G.J. The Natural Philosophy of Kazuo Kondo. arXiv 2007, arXiv:0712.0641.
126. Kawaguchi, M. An introduction to the theory of higher order spaces I. The theory of Kawaguchi spaces.
RAAG Memoirs 1960, 3, 718–734.
127. Barndorff-Nielsen, O.E.; Cox, D.R.; Reid, N. The role of differential geometry in statistical theory. Int. Stat. Rev.
1986, 54, 83–96. [CrossRef]
128. Nomizu, K.; Katsumi, N.; Sasaki, T. Affine Differential Geometry: Geometry of Affine Immersions; Cambridge
University Press: Cambridge, UK, 1994.
129. Norden, A.P. On Pairs of Conjugate Parallel Displacements in Multidimensional Spaces. In Doklady Akademii
nauk SSSR; Kazan State University, Comptes rendus de l’Académie des Sciences de l’URSS: Kazan, Russia,
1945; Volume 49, pp. 1345–1347.
130. Sen, R.N. On parallelism in Riemannian space I. Bull. Calcutta Math. Soc. 1944, 36, 102–107.
131. Sen, R.N. On parallelism in Riemannian space II. Bull. Calcutta Math. Soc. 1944, 37, 153–159.
132. Sen, R.N. On parallelism in Riemannian space III. Bull. Calcutta Math. Soc. 1946, 38, 161–167.
133. Giné, E.; Nickl, R. Mathematical Foundations of Infinite-Dimensional Statistical Models; Cambridge University
Press: Cambridge, UK, 2015; Volume 40.
134. Amari, S. New Developments of Information Geometry; Jouhou Kikagaku no Shintenkai; Saiensu’sha: Tokyo,
Japan, 2014. (In Japanese)
135. Fujiwara, A. Foundations of Information Geometry; Jouhou Kikagaku no Kisou; Makino Shoten: Tokyo, Japan,
2015; p. 223. (In Japanese)
136. Mitchell, A.F.S. The information matrix, skewness tensor and α-connections for the general multivariate
elliptic distribution. Ann. Inst. Stat. Math. 1989, 41, 289–304. [CrossRef]
137. Zhang, Z.; Sun, H.; Zhong, F. Information geometry of the power inverse Gaussian distribution. Appl. Sci.
2007, 9, 194–203.
138. Peng, T.L.L.; Sun, H. The geometric structure of the inverse gamma distribution. Contrib. Algebra Geom.
2008, 49, 217–225.
139. Zhong, F.; Sun, H.; Zhang, Z. The geometry of the Dirichlet manifold. J. Korean Math. Soc. 2008, 45, 859–870.
[CrossRef]
140. Peng, L.; Sun, H.; Jiu, L. The geometric structure of the Pareto distribution. Bol. Asoc. Mat. Venez. 2007,
14, 5–13.
141. Pistone, G. Nonparametric information geometry. In Geometric Science of Information; Springer: Berlin,
Germany, 2013; pp. 5–36.
142. Hayashi, M. Quantum Information; Springer: Berlin, Germany, 2006.
143. Pardo, M.d.C.; Vajda, I. About distances of discrete distributions satisfying the data processing theorem of
information theory. IEEE Trans. Inf. Theory 1997, 43, 1288–1293. [CrossRef]
144. Nielsen, F.; Nock, R. Total Jensen divergences: Definition, properties and k-means++ clustering. arXiv 2013,arXiv:1309.7109.
145. Nielsen, F.; Nock, R. Total Jensen divergences: Definition, properties and clustering. In Proceedings ofthe 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brisbane,
Australia, 19–24 April 2015; pp. 2016–2020.
146. Nielsen, F.; Nock, R. Patch Matching with Polynomial Exponential Families and Projective Divergences.
In Proceedings of the International Conference on Similarity Search and Applications (SISAP), Tokyo, Japan,
24–26 October 2016, pp. 109–116.
147. Nielsen, F.; Sun, K.; Marchand-Maillet, S. On Hölder Projective Divergences. Entropy 2017, 19, 122. [CrossRef]
148. Nielsen, F.; Barbaresco, F. (Eds.) Geometric Science of Information; Lecture Notes in Computer Science;
Springer: Berlin, Germany, 2013; Volume 8085. [CrossRef]
149. Nielsen, F.; Barbaresco, F. (Eds.) Geometric Science of Information; Lecture Notes in Computer Science;
Springer: Berlin, Germany, 2015; Volume 9389. [CrossRef]
150. Nielsen, F.; Barbaresco, F. (Eds.) Geometric Science of Information; Lecture Notes in Computer Science; Springer:
Berline, Germany, 2017; Volume 10589. [CrossRef]

Entropy 2020, 22, 1100

61 of 61

151. Nielsen, F. Geometric Structures of Information; Springer: Berlin, Germany, 2018.
152. Nielsen, F. Geometric Theory of Information; Springer: Berlin, Germany, 2014.
153. Ay, N.; Gibilisco, P.; Matús, F. Information Geometry and its Applications: On the Occasion of Shun-ichi Amari’s
80th Birthday, IGAIA IV Liblice, Czech Republic, 12–17 June 2016; Springer Proceedings in Mathematics &Statistics; Springer: Berlin, Germany, 2018; Volume 252.
154. Keener, R.W. Theoretical Statistics: Topics for a Core Course; Springer: Berlin, Germany, 2011.
155. Nielsen, F.; Sun, K. Guaranteed bounds on the Kullback–Leibler divergence of univariate mixtures.
IEEE Signal Process. Lett. 2016, 23, 1543–1546. [CrossRef]
156. Gordon, G.J. Approximate Solutions to Markov Decision Processes. Ph.D. Thesis, Carnegie Mellon
University, Pittsburgh, PA, USA, 1999.
157. Telgarsky, M.; Dasgupta, S. Agglomerative Bregman clustering. In Proceedings of the 29th International
Conference on Machine Learning, Edinburgh, UK, 26 June–1 July 2012; Omnipress: Madison, WI, USA,pp. 1011–1018.
158. Yoshizawa, S.; Tanabe, K. Dual differential geometry associated with Kullback–Leibler information on theGaussian distributions and its 2-parameter deformations. SUT J. Math. 1999, 35, 113–137.
159. Nielsen, F. On the Jensen–Shannon symmetrization of distances relying on abstract means. Entropy 2019,21, 485. [CrossRef]
160. Niculescu, C.; Persson, L.E. Convex Functions and Their Applications; Springer: Berlin, Germany, 2006.
161. Nielsen, F.; Nock, R. The Bregman chord divergence. In Proceedings of the International Conference
on Geometric Science of Information, Toulouse, France, 27–29 August 2019; Springer: Berlin, Germany,pp. 299–308.
c 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open accessarticle distributed under the terms and conditions of the Creative Commons Attribution(CC BY) license (http://creativecommons.org/licenses/by/4.0/).