Weyl Prior and Bayesian Statistics

Abstract

When using Bayesian inference, one needs to choose a prior distribution for parameters. The well-known Jeffreys prior is based on the Riemann metric tensor on a statistical manifold. Takeuchi and Amari defined the $\alpha$-parallel prior, which generalized the Jeffreys prior by exploiting a higher-order geometric object, known as a Chentsov–Amari tensor. In this paper, we propose a new prior based on the Weyl structure on a statistical manifold. It turns out that our prior is a special case of the $\alpha$-parallel prior with the parameter $\alpha$ equaling $-n$, where n is the dimension of the underlying statistical manifold and the minus sign is a result of conventions used in the definition of $\alpha$-connections. This makes the choice for the parameter $\alpha$ more canonical. We calculated the Weyl prior for univariate Gaussian and multivariate Gaussian distribution. The Weyl prior of the univariate Gaussian turns out to be the uniform prior.

1. Introduction

In Bayesian inference, a parameter is regarded as a random variable $\mathsf{\Theta }$. A density of $\mathsf{\Theta }$ is called, by abuse of terminology, a prior distribution $p\left(\theta \right)$. After collecting some data, one obtains a conditional density $p\left(x\right|\theta )$, referred to as the likelihood function. The Bayes’ theorem then computes the posterior distribution $p\left(\theta \right|x)$ using $p\left(\theta \right)$ and $p\left(x\right|\theta )$. This is interpreted as an update of the information about the unknown parameter $\mathsf{\Theta }$ in Bayesian inference. One such choice of the prior distribution $p\left(\theta \right)$ is Jeffreys prior ${}^J\omega$, which is the correct choice of uniform distribution. Here, the word “uniform” means uninformative, not favorable of any particular choice of the parameter.

Information geometry, in its narrowest sense, is an attempt to use differential geometry to study statistical inference. It has found applications in statistical inference, signal processing, and machine learning [1]. References [2,3] are two elementary introductions. In information geometry, geometric structures, for example metric tensors g and affine connections ∇, can be put on the set of prior distributions $\mathcal{P}\left(\theta \right)$. These geometric structures help to single out some particular prior distributions, for example, the Jeffreys prior ${}^J\omega ,$ which, by the fundamental theorem of Riemannian geometry, is the unique volume form parallel with respect to the Levi–Civita connection ${}^{LC}\nabla$. Since the Jeffreys prior ${}^J\omega$ is provided by geometry, it is automatically invariant under reparametrization, which reflects the opinion that information can be at best not lost during a transformation of parameters and this is encoded in the notion of sufficient statistics. Similarly, if one can find a unique prior distribution satisfying some specified geometric conditions, then that prior distribution is called canonically chosen. Matsuzoe, Takeuchi, and Amari used information geometry to define the $\alpha$-parallel prior ${}^{\alpha }\omega$ such that, when $\alpha =0$, it reduces to the Jeffreys prior [4].

Historically, Weyl proposed a generalization of general relativity to unify gravity and electromagnetism. Einstein soon pointed out that Weyl’s theory predicted substantial broadening of the characteristic length of atoms, which is contradictory to the well-observed thin atomic spectra. Even though Weyl geometry failed the unification of gravity and electromagnetism, which is still an open problem, Weyl geometry has found applications in possible generalization of general relativity [5] and the differential-geometric study of defects in continuum mechanics [6]. Weyl geometry is kept in mathematics. The relation of affine differential geometry, Weyl geometry, and Riemannian geometry are shown below. Let $\pi :E\to B$ be a fibre bundle with base B and let each fibre ${\pi }^{-1}\left(x\right),$ $x\in B$, be a Lie group G, called the structure group of the fibre bundle. For different G, we obtain different geometries as follows:

• $GL\left(n\right),$ the general linear group, affine differential geometry;
• $C\left(n\right):=\{kA | k\in {\mathbb{R}}^{+}, A\in O\left(n\right)\},$ the conformal group, Weyl geometry;
• $O\left(n\right),$ the orthogonal group, Riemannian geometry.

With the reduction of structure groups $O\left(n\right)\subseteq C\left(n\right)\subseteq GL\left(n\right)$, and the fact that smaller the structure group, the more geometric properties are expected. In our case, the reduction of group gives rise to a canonical choice for the parameter $\alpha$ of the $\alpha$-parallel prior. For more about bundle-theoretic differential geometry, see [7].

In this paper, we will use Weyl geometry to define a prior distribution for Bayesian inference, which we call the Weyl prior. We will elucidate the relation between the dimension of a statistical manifold and the parameter $\alpha$ in Takeuchi and Amari’s $\alpha$-parallel prior.

The organization of the paper is as follows: In Section 2, we review information geometry and the $\alpha$-parallel prior. We discuss Weyl geometry in Section 3. We define the Weyl prior, and elucidate the relation between the Weyl prior and the $\alpha$-parallel prior in Section 4. We calculate the Weyl prior for the univariate Gaussian distribution as an example in Section 5 and the multivariate Gaussian distribution in Section 6. All functions in the paper are real-valued and smooth, all connections are torsion-free, and the Einstein summation rule is used.

2. Information Geometry and $\mathbf{\alpha }$-Priors

In this section, we review some basics of information geometry. For more details, see [1].

Let us consider a statistical model $\mathcal{P},$ which is a set of parametric densities $\mathcal{P}=\left\{p\right(x\left|\theta \right)\}.$ $\mathcal{P}$ can be geometrized as follows: first, we introduce the Fisher metric tensor, which is a 2nd order tensor,

Definition 1.

The Fisher metric tensor is defined by

$$g_{ij}=E_{\theta }\left[{\partial }_{i}l{\partial }_{j}l\right]$$

(1)

where $E_{\theta }$ is the the transition kernel $X\times \mathsf{\Theta }\to [0,\infty )$, l is the log likelihood function, and ${\partial }_i$ is the partial derivative with respect to coordinate i.

Then, we introduce the Amari–Chentsov tensor, which is a 3rd order tensor.

Definition 2

(Amari–Chentsov Tensor). The Amari–Chentsov tensor C is defined by

$$C_{ijk}=E_{\theta }\left[{\partial }_{i}l{\partial }_{j}l{\partial }_{k}l\right].$$

(2)

Remark 1.

The Amari–Chentsov tensor C defined above satisfies

$$C=\nabla g.$$

In other words, C is the covariant derivative of the metric tensor $g.$

In Riemannian geometry, C vanishes everywhere, which is required if the length of a tangent vector is to be preserved under the parallel transport. In information geometry, this requirement is dropped and thus a duality theory arises. Let ∇ be an arbitrary torsion-free affine connection on a Riemannian manifold $(M,g).$ The dual connection ${\nabla }^{*}$ of ∇ plays an important role in information geometry.

Definition 3

(Dual Connection). The dual connection ${\nabla }^{*}$ on a Riemannian manifold $(M,g)$ with affine connection ∇ is defined as the unique affine connection satisfying the following equation:

$$X{g}_p(Y,Z)=g_p({\nabla }_{X}Y,Z)+g_p(Y,{\nabla }_X^{*}Z),$$

(3)

where $p\in M$ and $X,Y,Z\in T_{p}M.$

Remark 2.

The dual connection ${\nabla }^{*}$ preserves the metric tensor g together with $\nabla :$

$$g_p(X,Y)=g_q(\mathsf{\Pi }X,{\mathsf{\Pi }}^{*}Y),$$

where $X,Y\in T_{p}M,$ and Π and ${\mathsf{\Pi }}^{*}$ are parallel transports induced by ∇ and ${\nabla }^{*}$, respectively, along some curve from p to q. In general, $g_p(X,Y)\ne g_q(\mathsf{\Pi }X,\mathsf{\Pi }Y)$ and $g_p(X,Y)\ne g_q({\mathsf{\Pi }}^{*}X,{\mathsf{\Pi }}^{*}Y)$, unless $\nabla ={\nabla }^{*}{=}^{LC}\nabla$. See [1].

Now, we introduce $\alpha$-connections.

Definition 4.

The α-connections are defined in terms of Christoffel symbols by

$${}^{\alpha }{\mathsf{\Gamma }}_{jk}^i{=}^{LC}{\mathsf{\Gamma }}_{jk}^i-\frac{\alpha }{2}{g}^{il}{C}_{ljk},$$

(4)

where $\alpha \in \mathbb{R}$ and LC stands for Levi–Civita.

Remark 3.

The dual connection of ${\nabla }^{\alpha }$ is then given by

$${\nabla }^{\alpha *}={\nabla }^{-\alpha }.$$

(5)

Remark 4.

The α-parallel prior ${}^{\alpha }\omega$ is the volume form parallel with respect to ${\nabla }^{\alpha }.$ Unlike the Jeffreys prior, which always exists, the α-parallel prior do not necessarily exist. An α-parallel prior exists if and only if the Ricci curvature tensor is symmetric [4]. However, if ${}^{\alpha }\omega$ exists for one $\alpha \in \mathbb{R},$ then it exists for all α [8].

The following characterization will be used in Section 5 to obtain the relation between the $\alpha$-parallel prior and the Weyl prior defined therein.

Proposition 1.

[4] Let $(M,g{,}^{\alpha }\nabla )$ be a statistical manifold. If there exists an exact 1-form $T=d\mathsf{\Omega }$ for some function Ω determined by ∇ and g, then the α-parallel prior is ${}^{\alpha }\omega =exp\{-\frac{\alpha }{2}\mathsf{\Omega }\}\sqrt{detg}.$

Remark 5.

$d\mathsf{\Omega }$ is known as the Chebyshev 1-form. A differential form ϕ is called closed if the exterior derivative vanishes i.e., $d\varphi =0$, and is called exact if there exists a differential form φ such that $\varphi =d\phi$. By definition, every exact form is closed. By Poincare’s lemma, every closed form is locally exact. Because statistical manifolds are simply connected, closedness implies exactness.

3. Weyl Geometry

In this section, we review some concepts of Weyl geometry which are needed in the next section. For more details, see [9].

Two Riemannian metrics g and $g^{\prime }$ on a manifold M are said to be conformally equivalent if $g^{\prime }={\mathrm{e}}^{\lambda }g$ for some smooth function $\lambda$ on M.

A conformal structure $\mathcal{C}$ on M is an equivalent class of conformally equivalent Riemannian metrics, i.e., $\mathcal{C}:=\left\{g^{\prime }\right|g^{\prime }={\mathrm{e}}^{\lambda }g\}$.

A Weyl structure is a map $F:\mathcal{C}\to {\mathsf{\Lambda }}^1\left(M\right)$ from the conformal structure $\mathcal{C}$ to the set of 1-forms on $M,$ satisfying

$$F\left({\mathrm{e}}^{\lambda }g\right)=F\left(g\right)-d\lambda .$$

The image of g under F is called the Weyl 1-form $F\left(g\right):=\phi .$

A Weyl structure enables us to translate a scalar product ${( , )}_p$ at p to ${( , )}_q$ at q along a curve $c:[0,1]\to M:$

$${( , )}_q=exp\left[{\int }_0^1{c}^{*}\phi \right]g_q,$$

(6)

where $c^{*}\phi$ is the pullback of the Weyl $1-$form $\phi$ along curve $c.$ A Weyl manifold is a manifold with a Weyl structure.

Remark 6.

The meaning of this equation is: If we start with a scalar product ${(,)}_p$ at a point p arising from the conformal class $\mathcal{C}$, then there exists a metric tensor $g\in \mathcal{C}$ extending ${(,)}_p$, i.e., $g_p={(,)}_p$. The value of this particular choice of g at another point q is $g_q$. However, different choice of g gives rise to different $g_q$. The scalar product ${(,)}_q$ determined by Weyl translation is proven to be independent of g [9]. Hence, by Weyl translation, we can compare lengths of vectors at different points on a Weyl manifold, whereas, with only the conformal structure $\mathcal{C}$, we can only compare ratios of lengths.

An affine connection ∇ is said to be a Weyl connection if the parallel transport of a scalar product under ∇ coincides with the Weyl translation.

The Weyl connection is characterized by the following propositions.

Proposition 2

([9]). An affine connection ∇ is a Weyl connection if and only if $\nabla g+\phi \otimes g=0$ for all $g\in \mathcal{C}.$

Proposition 3

(Fundamental Theorem of Weyl Geometry [9]). There exists a unique torsion-free Weyl connection ${}^W\nabla$ on a Weyl manifold $M.$ The Christoffel symbols of ${}^W\nabla$ are given by

$${}^W{\mathsf{\Gamma }}_{jk}^i{=}^{LC}{\mathsf{\Gamma }}_{jk}^i+\frac{1}{2}\left({\delta }_j^i{\phi }_k+{\delta }_k^i{\phi }_j-g^{im}{g}_{jk}{\phi }_m\right),$$

(7)

where ${\delta }_j^i$ is the Kronecker delta.

4. Weyl Prior

In this section, we define the Weyl prior and show its relation to the $\alpha$-parallel prior.

First, we define the Weyl prior as follows.

Definition 5

(Weyl Prior). Let $(M,g)$ be an n-dimensional Riemannian manifold with the conformal structure $\mathcal{C}=\left[g\right]$ and the Weyl structure $F.$ Let ${}^W\nabla$ be the Weyl connection. The Weyl prior ${}^W\omega$ is defined as the unique volume form parallel with respect to ${}^W\nabla .$

Remark 7.

The uniqueness of the Weyl prior is the result of the uniqueness of the Weyl connection.

Now, we prove the main result of this paper.

Theorem 1.

Let $(M,g)$ be a Riemannian manifold. Let ${}^W\nabla$ and ${}^{-n}\nabla$ be the Weyl connection and the $-n$-connection, i.e., the α-connection with $\alpha =-n$, where n is the dimension of $M.$ Suppose that the $-n$-prior ${}^{-n}\omega$ exists, then

$${}^W\omega {=}^{-n}\omega .$$

Proof. 

Consider an arbitrary volume form $f\sqrt{detg},$ where f is a positive function on $M.$ For $f\sqrt{detg}$ to be parallel with respect to ${}^W\omega ,$ it is necessary and sufficient that

$$\begin{array}{c}\hfill {}^W\nabla \left(f\sqrt{detg}\right)=f^W\nabla \sqrt{detg}{+}^W\nabla f\sqrt{detg}=0.\end{array}$$

(8)

Componentwise, Equation (8) becomes

$$\begin{array}{c}\hfill f^W{\nabla }_j\sqrt{detg}{+}^W{\nabla }_{j}f\sqrt{detg}=f^W{\nabla }_j\sqrt{detg}+{\partial }_{j}f\sqrt{detg},\end{array}$$

(9)

since covariant derivative coincides with partial derivative for functions.

Since $\sqrt{detg}$ is a scalar density of weight $1,$ its covariant derivative is given by

$${}^W{\nabla }_j\sqrt{detg}={\partial }_j\sqrt{detg}{-}^W{\mathsf{\Gamma }}_j\sqrt{detg},$$

(10)

where ${\mathsf{\Gamma }}_j$ is obtained by the contraction of Equation (7) over i and $k:$

$$\begin{array}{cc}\hfill {}^W{\mathsf{\Gamma }}_j{=}^W{\mathsf{\Gamma }}_{ji}^i& {=}^{LC}{\mathsf{\Gamma }}_{ji}^i+\frac{1}{2}\left({\delta }_j^i{\phi }_i+{\delta }_i^i{\phi }_j-g^{im}{g}_{ji}{\phi }_m\right)\hfill \\ & ={\partial }_{j}ln\sqrt{detg}+\frac{1}{2}\left({\phi }_j+n{\phi }_j-{\delta }_j^m{\phi }_m\right)\hfill \\ & ={\partial }_{j}ln\sqrt{detg}+\frac{n}{2}{\phi }_j.\hfill \end{array}$$

Substituting Equation (10) into Equation (9), we obtain

$${\partial }_{j}f=\frac{n}{2}{\phi }_{j}f.$$

(11)

Since the covariant derivative coincides with exterior derivative for functions, collect indices in Equation (11)

$$\phi =\frac{2}{n}\mathrm{d}lnf.$$

(12)

Assume for now that the Weyl 1-form $\phi$ is exact, that is, $\phi =\mathrm{d}\mathsf{\Omega }$ for some function $\mathsf{\Omega }$ on $M.$ Then, from Equation (29), the Weyl prior is given by

$${}^W\omega =exp\left\{\frac{n}{2}\mathsf{\Omega }\right\}\sqrt{detg}.$$

(13)

By comparison of Equation (13) with Proposition 1, the theorem is proved under the assumption of the exactness of the Weyl 1-form.

Since we proved that the Weyl prior ${}^W\omega$ is the $\alpha -$prior with $\alpha =-n,$ and we required the existence of ${}^n\omega ,$ our assumption of the exactness of the Weyl $1-$form $\phi$ is indeed true by Remark 4. □

Remark 8.

The minus sign in $\alpha =-n$ is a result of the definition of α-connection. By Remark 3, the dual connection of ${}^{\alpha }\nabla$ is ${}^{-\alpha }\nabla ,$ we would have $\alpha =n$ here, had we defined the α-connection to be its dual connection in Definition 4. It would seem more natural to consider the dual prior of the Weyl prior.

5. Weyl Prior for Gaussian Family

In this section, we calculate the Weyl prior of the Gaussian family as an example.

Example 1

(Gaussian Family). Consider the Gaussian family

$$\mathcal{P}=\left\{p\left(x\right|\mu ,{\sigma }^2)=\frac{1}{\sqrt{2\pi }\sigma }exp\left[-\frac{1}{2{\sigma }^2}{(x-\mu )}^2\right] | (\mu ,\sigma )\in M\right\}.$$

Choosing $(\mu ,{\sigma }^2)$ as a coordinate system, we have

$${\partial }_{\mu }l=\frac{x-\mu }{{\sigma }^2},$$

$${\partial }_{{\sigma }^2}l=\frac{{(x-\mu )}^2}{2{\sigma }^4}-\frac{1}{2{\sigma }^2},$$

where l is the log likelihood function.

The first element of the Fisher metric tensor g in the $(\mu ,{\sigma }^2)$-coordinate is given by

$$\begin{array}{cc}\hfill g_{\mu \mu }& ={\mathrm{E}}_{\theta }\left[{\partial }_{\mu }l{\partial }_{\mu }l\right]\hfill \\ & ={\int }_{-\infty }^{\infty }\frac{{(x-\mu )}^2}{{\sigma }^4}\frac{1}{\sqrt{2\pi }\sigma }exp\left[-\frac{1}{2{\sigma }^2}{(x-\mu )}^2\right]dx\hfill \\ & =\frac{1}{{\sigma }^2},\hfill \end{array}$$

where $E_{\theta }$ is the conditional expectation of X given $\mathsf{\Theta }.$ The other elements of the Fisher metric tensor are

$$g_{\mu {\sigma }^2}=g_{{\sigma }^2\mu }=0,$$

and

$$g_{{\sigma }^2{\sigma }^2}=\frac{1}{2{\sigma }^4}.$$

Hence,

$$\sqrt{detg}=\frac{1}{\sqrt{2}{\sigma }^3}.$$

(14)

To calculate the Weyl 1-form, we first calculate the Amari–Chentsov tensor $C,$

$$\begin{array}{cc}\hfill C_{\mu \mu \mu }& ={\mathrm{E}}_{\theta }\left[{\partial }_{\mu }l{\partial }_{\mu }l{\partial }_{\mu }l\right]\hfill \\ & ={\int }_{-\infty }^{\infty }\frac{{(x-\mu )}^3}{{\sigma }^6}\frac{1}{\sqrt{2\pi }\sigma }exp\left[-\frac{1}{2{\sigma }^2}{(x-\mu )}^2\right]dx\hfill \\ & =0.\hfill \end{array}$$

Similarly,

$$C_{{\sigma }^2\mu \mu }=C_{\mu {\sigma }^2\mu }=C_{\mu \mu {\sigma }^2}=\frac{1}{{\sigma }^4},$$

$$C_{{\sigma }^2{\sigma }^2\mu }=C_{{\sigma }^2\mu {\sigma }^2}=C_{\mu {\sigma }^2{\sigma }^2}=0,$$

and

$$C_{{\sigma }^2{\sigma }^2{\sigma }^2}=\frac{1}{{\sigma }^6}.$$

Hence, the Weyl 1-form is given by

$$\begin{array}{cc}\hfill \phi & =\frac{1}{2}{C}_{ijk}{g}^{jk}d{\theta }^i\hfill \\ & =\frac{3}{2{\sigma }^2}d{\sigma }^2.\hfill \end{array}$$

(15)

Now, it is easy to check that $\phi =\mathrm{d}(\frac{3}{2}ln{\sigma }^2)$ is an exact form. Hence, for Gaussian family $\mathcal{P},$ a Weyl prior exists and is given by

$$\begin{array}{cc}\hfill {}^W\omega & =exp\left\{\frac{2}{2}\frac{3}{2}ln{\sigma }^2\right\}\frac{1}{\sqrt{2}{\sigma }^3}\hfill \\ & =\frac{1}{\sqrt{2}}.\hfill \end{array}$$

(16)

Remark 9.

Based on our calculation, we find that the Weyl prior for the univariate Gaussian distribution with unknown mean and unknown variance is just the uniform prior. This shows that the uniform prior is in fact a uninformative prior. This counter-intuitive result is related to the fact that every two-dimensional manifold is conformally-flat, which can be proved using the existence of isothermal coordinates in two dimensions [10].

6. Multivariate Gaussian

The above example can be extended to the multivariate case. Consider the multivariate Gaussian distribution

$$f\left(x\right|\mu ,\mathsf{\Sigma })=\frac{1}{{\left(2\pi \right)}^{n/2}\sqrt{det\mathsf{\Sigma }}}exp\{-\frac{1}{2}{(\mathit{x}-\mathbf{\mu })}^{⊺}{\mathsf{\Sigma }}^{-1}(\mathbf{x}-\mathbf{\mu })\},$$

(17)

where $\mathbf{\mu }$ is the mean vector and $\mathsf{\Sigma }$ is the covariance matrix.

Using matrix calculus, we have

$${\partial }_{\mathbf{\mu }}l={\mathsf{\Sigma }}^{-1}(\mathit{x}-\mathbf{\mu })$$

(18)

and

$$\begin{array}{cc}\hfill {\partial }_{\mathsf{\Sigma }}l& =-\frac{1}{2}{\mathsf{\Sigma }}^{-1}+\frac{1}{2}{\mathsf{\Sigma }}^{-1}(\mathit{x}-\mathbf{\mu }){(\mathit{x}-\mathbf{\mu })}^{⊺}{\mathsf{\Sigma }}^{-1}\hfill \\ & =-\frac{1}{2}{\mathsf{\Sigma }}^{-1}+\frac{1}{2}\left[{\mathsf{\Sigma }}^{-1}(\mathit{x}-\mathbf{\mu })\right]\otimes \left[{\mathsf{\Sigma }}^{-1}(\mathit{x}-\mathbf{\mu })\right]\hfill \\ & =-\frac{1}{2}{\mathsf{\Sigma }}^{-1}+\frac{1}{2}({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1})[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })]\hfill \end{array}$$

(19)

We can now compute the Fisher information matrix.

$$\begin{array}{cc}\hfill g_{\mathbf{\mu }\mathbf{\mu }}& ={\mathrm{E}}_{\theta }[{\partial }_{\mathbf{\mu }}\otimes {\partial }_{\mathbf{\mu }}]\hfill \\ & =({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}){\mathrm{E}}_{\theta }[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })]\hfill \\ & =({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1})\mathsf{\Sigma }\hfill \\ & ={\mathsf{\Sigma }}^{-1}\mathsf{\Sigma }{\mathsf{\Sigma }}^{-⊺}\hfill \\ & ={\mathsf{\Sigma }}^{-1}\hfill \end{array}$$

(20)

where the second last line is by the action of matrix tensor product and the last line is by the definition of covariance matrix.

Similarly,

$$g_{\mathbf{\mu }\mathsf{\Sigma }}=g_{\mathsf{\Sigma }\mathbf{\mu }}=\mathbf{0},$$

(21)

and

$$g_{\mathsf{\Sigma }\mathsf{\Sigma }}=\frac{1}{2}{\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}.$$

(22)

The Amari–Chentsov tensor can be computed in the same way:

$$\begin{array}{c}\hfill C_{\mathbf{\mu }\mathbf{\mu }\mathbf{\mu }}={\mathrm{E}}_{\theta }[{\partial }_{\mathbf{\mu }}l\otimes {\partial }_{\mathbf{\mu }}l\otimes {\partial }_{\mathbf{\mu }}l]=\mathbf{0}.\end{array}$$

(23)

$$C_{\mathsf{\Sigma }\mathbf{\mu }\mathbf{\mu }}=C_{\mathbf{\mu }\mathsf{\Sigma }\mathbf{\mu }}=C_{\mathbf{\mu }\mathbf{\mu }\mathsf{\Sigma }}={\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}.$$

(24)

$$C_{\mathbf{\mu }\mathsf{\Sigma }\mathsf{\Sigma }}=C_{\mathsf{\Sigma }\mathbf{\mu }\mathsf{\Sigma }}=C_{\mathsf{\Sigma }\mathsf{\Sigma }\mathbf{\mu }}=\mathbf{0}.$$

(25)

$$C_{\mathsf{\Sigma }\mathsf{\Sigma }\mathsf{\Sigma }}={\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}.$$

(26)

The detail computation of Equation (26) is as follows:

$$\begin{array}{cc}\hfill C_{\mathsf{\Sigma }\mathsf{\Sigma }\mathsf{\Sigma }}& ={\mathrm{E}}_{\theta }[{\partial }_{\mathsf{\Sigma }}l\otimes {\partial }_{\mathsf{\Sigma }}l\otimes {\partial }_{\mathsf{\Sigma }}l]\hfill \\ & ={\mathrm{E}}_{\theta }\{-\frac{1}{8}{\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}+\frac{1}{8}{\mathsf{\Sigma }}^{-1}\otimes ({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1})\left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]\otimes {\mathsf{\Sigma }}^{-1}\hfill \\ & +\frac{1}{8}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]\otimes ({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1})-\frac{1}{8}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\hfill \\ & (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\otimes {\mathsf{\Sigma }}^{-1}+\frac{1}{8}{\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes \left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\hfill \\ & \left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]-\frac{1}{8}{\mathsf{\Sigma }}^{-1}\otimes \left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]\otimes \left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\hfill \\ & \left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]-\frac{1}{8}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]\otimes {\mathsf{\Sigma }}^{-1}\otimes \left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\hfill \\ & \left[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\right]+\frac{1}{8}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)[(\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\hfill \\ & \otimes (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\otimes (\mathit{x}-\mathbf{\mu })\left]\right\}={\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\hfill \end{array}$$

(27)

The above expression can be evaluated by 4th and 6th moments of multivariate Gaussian.

The Weyl prior is then given by:

$$\begin{array}{cc}\hfill \phi & =\frac{1}{2}{C}_{ijk}{g}^{jk}d{\theta }^i\hfill \\ & =\frac{1}{2}{C}_{\mathsf{\Sigma }\mathbf{\mu }\mathbf{\mu }}{g}^{\mathbf{\mu }\mathbf{\mu }}d\mathsf{\Sigma }+\frac{1}{2}{C}_{\mathsf{\Sigma }\mathsf{\Sigma }\mathsf{\Sigma }}{g}^{\mathsf{\Sigma }\mathsf{\Sigma }}d\mathsf{\Sigma }\hfill \\ & =\frac{1}{2}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)\mathsf{\Sigma }d\mathsf{\Sigma }+\frac{1}{2}\left({\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right)2\left(\mathsf{\Sigma }\otimes \mathsf{\Sigma }\right)d\mathsf{\Sigma }\hfill \\ & =\frac{3}{2}{\mathsf{\Sigma }}^{-1}d\mathsf{\Sigma }\hfill \\ & =\mathrm{d}\left(\frac{3}{2}lndet\mathsf{\Sigma }\right).\hfill \end{array}$$

(28)

The Weyl prior is thus given by:

$$\begin{array}{cc}\hfill {}^W\omega & =exp\left\{\frac{n+(n+1)n/2}{2}\frac{3}{2}lndet\mathsf{\Sigma }\right\}\sqrt{det{\mathsf{\Sigma }}^{-1}det\left[\frac{1}{2}{\mathsf{\Sigma }}^{-1}\otimes {\mathsf{\Sigma }}^{-1}\right]}\hfill \\ & ={\left(det\mathsf{\Sigma }\right)}^{(3n^2+9n)/8}\frac{{\left(det\mathsf{\Sigma }\right)}^{(-2n-1)/2}}{2^{n/2}}\hfill \\ & =\frac{{\left(det\mathsf{\Sigma }\right)}^{(n-1)(3n+4)/8}}{2^{n/2}},\hfill \end{array}$$

(29)

where, in the first line, $n+(n+1)n/2$ is the dimension of the statistical manifold for the multivariate Gaussian distribution.

Remark 10.

Our calculation of the Weyl prior of the multivariate Gaussian distribution is generally not a uniform prior. However, Equation (29) shows that, when $n=1$, that is, the univariate case, the Weyl prior is indeed the uniform prior. This is in accordance with our direct calculation for the univariate case.

7. Discussion and Conclusions

We discussed Weyl geometry and Weyl prior in this paper. We also calculated Weyl prior for the Gaussian family as an example.

The underlying principle of Jeffreys prior, $\alpha$-parallel prior, and Weyl prior is the concept of invariance in statistics. Jeffreys prior is invariant under a change of the coordinate of parameters. Weyl prior and $\alpha$-parallel prior, as generalizations of Jeffreys prior, automatically satisfy this invariance. Moreover, Weyl prior, as a volume form defined on a Weyl manifold, is also invariant under a gauge transformation [11]. In addition, invariant under the gauge transformation is the generalized conjugate connection [11].

One possible use of the Weyl prior is using the uniform prior for distributions with two parameters. This is because any two-dimensional manifold is conformally-flat.