Next Article in Journal
Application of the Self-Organization Phenomenon in the Development of Wear Resistant Materials—A Review
Next Article in Special Issue
Geometric Theory of Heat from Souriau Lie Groups Thermodynamics and Koszul Hessian Geometry: Applications in Information Geometry for Exponential Families
Previous Article in Journal
Bounds on Rényi and Shannon Entropies for Finite Mixtures of Multivariate Skew-Normal Distributions: Application to Swordfish (Xiphias gladius Linnaeus)
Previous Article in Special Issue
Non-Asymptotic Confidence Sets for Circular Means
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Explicit Formula of Koszul–Vinberg Characteristic Functions for a Wide Class of Regular Convex Cones

Graduate School of Mathematics, Nagoya University, Nagoya 464-8602, Japan
Entropy 2016, 18(11), 383; https://doi.org/10.3390/e18110383
Submission received: 16 September 2016 / Revised: 19 October 2016 / Accepted: 20 October 2016 / Published: 26 October 2016
(This article belongs to the Special Issue Differential Geometrical Theory of Statistics)

Abstract

:
The Koszul–Vinberg characteristic function plays a fundamental role in the theory of convex cones. We give an explicit description of the function and related integral formulas for a new class of convex cones, including homogeneous cones and cones associated with chordal (decomposable) graphs appearing in statistics. Furthermore, we discuss an application to maximum likelihood estimation for a certain exponential family over a cone of this class.

1. Introduction

Let Ω be an open convex cone in a vector space Z . The cone Ω is said to be regular if Ω contains no straight line, which is equivalent to the condition Ω ¯ ( Ω ¯ ) = { 0 } . In this paper, we always assume that a convex cone is open and regular. The dual cone Ω * with respect to an inner product ( · | · ) on Z is defined by:
Ω * : = ξ Z ; ( x | ξ ) > 0 ( x Ω ¯ \ { 0 } ) .
Then, Ω * is again a regular open convex cone, and we have ( Ω * ) * = Ω . The Koszul–Vinberg characteristic function φ Ω : Ω R > 0 defined by:
φ Ω ( x ) : = Ω * e ( x | ξ ) d ξ ( x Ω )
plays a fundamental role in the theory of regular convex cones [1,2,3,4].
In particular, φ Ω is an important function in the theory of convex programming [5], and it has also been studied recently in connection with thermodynamics [6,7]. There are several (not many) classes of cones for which an explicit formula of the Koszul–Vinberg characteristic function is known. Among them, the class of homogeneous cones [8,9,10] and the class of cones associated with chordal graphs [11] are particularly fruitful research objects. In this paper, we present a wide class of cones, including both of them, and give an explicit expression of the Koszul–Vinberg characteristic function (Section 3). Moreover, we get integral formulas involving the characteristic functions and the so-called generalized power functions, which are expressed as some product of powers of principal minors of real symmetric matrices (Section 4). After investigating the multiplicative Legendre transform of generalized power functions in Section 5, we study a maximum likelihood estimator for a Wishart-type natural exponential family constructed from the integral formula (Section 6).
A regular open convex cone Ω Z is said to be homogeneous if the linear automorphism group G L ( Ω ) : = α G L ( Z ) ; α Ω = Ω acts on Ω transitively. The cone P n of positive definite n × n real symmetric matrices is a typical example of homogeneous cones. It is known [12,13,14,15,16] that every homogeneous cone is linearly isomorphic to a cone P n Z with an appropriate subspace Z of the vector space Sym ( n , R ) of all n × n real symmetric matrices, where Z admits a specific block decomposition. Based on such results, our matrix realization method [15,17,18] has been developed for the purpose of the efficient study of homogeneous cones. In this paper, we present a generalization of matrix realization dealing with a wide class of convex cones, which turns out to include cones associated with chordal graphs. Actually, it was an enigma for the author that some formulas in [11,19] for the chordal graph resemble the formulas in [8,17] for homogeneous cones so much, and the mystery is now solved by the unified method in this paper to get the formulas. Furthermore, the techniques and ideas in the theory of homogeneous cones, such as Riesz distributions [8,20,21] and homogeneous Hessian metrics [4,18,22], will be applied to various cones to obtain new results in our future research.
Here, we fix some notation used in this paper. We denote by Mat ( p , q , R ) the vector space of p × q real matrices. For a matrix A, we write t A for the transpose of A. The identity matrix of size p is denoted by I p .

2. New Cones P V and P V *

2.1. Setting

We fix a partition n = n 1 + n 2 + + n r of a positive integer n. Let V = { V l k } 1 k < l r be a system of vector spaces V l k Mat ( n l , n k , R ) satisfying
(V1)
A V l k A t A R I n l ( 1 k < l r ) ,
(V2)
A V l j , B V k j A t B V l k ( 1 j < k < l r ) .
The integer r is called the rank of the system V . We denote by n l k the dimension of V l k . Note that some n l k can be zero. Let Z V be the space of real symmetric matrices x Sym ( n , R ) of the form:
x = X 11 t X 21 t X r 1 X 21 X 22 t X r 2 X r 1 X r 2 X r r X k k = x k k I n k , x k k R , k = 1 , , r X l k V l k , 1 k < l r ,
and P V the subset of Z V consisting of positive definite matrices. Then, P V is a regular open convex cone in Z V .
Example 1.
Let r = 3 , and set V 21 : = a 0 ; a R , V 31 : = 0 a ; a R , and V 32 : = R . Then, Z V is the space of symmetric matrices x of the form:
x = x 1 0 x 4 0 0 x 1 0 x 5 x 4 0 x 2 x 6 0 x 5 x 6 x 3 .
We shall see later that the cone P V = Z V P 4 is not homogeneous in this case, but admits various integral formulas, as well as explicit expression of the Koszul–Vinberg characteristic function.

2.2. Inductive Description of P V

If the system V = { V l k } 1 k < l r satisfies (V1) and (V2), any subsystem V I : = { V l k } k , l I with I { 1 , , r } also satisfies the same conditions. In particular, the cone corresponding to the subsystem { V l k } 2 k < l r will play an important role in this paper. Let us define V : = { V l k } 1 k < l r 1 by V l k : = V l + 1 , k + 1 . Then, V is a system of rank r 1 . Any x Z V is written as:
x = x 11 I n 1 t U U x ( x 11 R , U W , x Z V ) ,
where:
W : = U = X 21 X r 1 ; X l 1 V l 1 ( 1 < l r ) .
If x 11 0 , then we have:
x 11 I n 1 t U U x = I n 1 x 11 1 U I n n 1 x 11 I n 1 x x 11 1 U t U I n 1 x 11 1 t U I n n 1 .
Note that U t U belongs to Z V thanks to (V1) and (V2). Thus, we deduce the following lemma immediately from (5).
Lemma 1.
(i) Let x Z V as in (3). Then, x P V if and only if x 11 > 0 and x x 11 1 U t U P V .
(ii) For x P V , there exist unique U ˜ W and x ˜ P V for which:
x = I n 1 U ˜ I n n 1 x 11 I n 1 x ˜ I n 1 t U ˜ I n n 1 = x 11 I n 1 x 11 t U ˜ x 11 U ˜ x ˜ + x 11 U ˜ t U ˜ .
(iii) The closure P V ¯ of the cone P V is described as:
P V ¯ : = x 11 I n 1 x 11 t U ˜ x 11 U ˜ x ˜ + x 11 U ˜ t U ˜ ; x 11 0 , U ˜ W , x ˜ P V ¯ .

2.3. The Dual Cone P V *

We define an inner product on the space V l k by ( A | B ) V l k : = n l 1 tr A t B for A , B V l k . Then, we see from (V1) that:
A t B + B t A = 2 ( A | B ) V l k I n l .
Gathering these inner products ( · | · ) V l k , we introduce the standard inner product on the space Z V defined by:
( x | x ) : = k = 1 r x k k x k k + 2 1 k < l r ( X l k | X l k ) V l k
for x , x Z V of the form (1). When n 1 = n 2 = = n r = 1 (and only in this case), the standard inner product above equals the trace inner product tr ( x x ) .
Let W ˜ k ( k = 1 , , r ) be the vector space of W Mat ( n , n k , R ) of the form:
W = 0 n 1 + + n k 1 , n k X k k X k + 1 , k X r k ( X k k = x k k I n k , x k k R , X l k V l k , l > k ) .
Clearly, the space W ˜ k is isomorphic to R l > k V l k , which implies dim W ˜ k = 1 + q k with q k : = l > k n l k . Gathering orthogonal bases of V l k ’s, we take a basis of W ˜ k , so that we have an isomorphism W ˜ k W w = vect ( W ) R 1 + q k , where the first component w 1 of w is assumed to be x k k . Let us introduce a linear map ϕ k : Z V Sym ( 1 + q k , R ) defined in such a way that:
( W t W | ξ ) = t w ϕ k ( ξ ) w ( ξ Z V , W W ˜ k , w = vect ( W ) R 1 + q k ) .
It is easy to see that ϕ r ( ξ ) = ξ r r for ξ Z V .
Theorem 1.
The dual cone P V * Z V of P V with respect to the standard inner product is described as:
P V * = ξ Z V ; ϕ k ( ξ ) i s   p o s i t i v e   d e f i n i t e   f o r   a l l   k = 1 , , r = ξ Z V ; det ϕ k ( ξ ) > 0 f o r   a l l   k = 1 , , r .
Proof. 
We shall prove the statement by induction on the rank r. When r = 1 , we have ϕ 1 ( ξ ) = ξ 11 and ξ = ξ 11 I n 1 . Thus, (9) holds in this case.
Let us assume that (9) holds when the rank is smaller than r. In particular, the statement holds for P V * Z V , that is,
P V * = ξ Z V ; ϕ k ( ξ ) is   positive   definite   for   all k = 1 , , r 1 = ξ Z V ; det ϕ k ( ξ ) > 0 for   all k = 1 , , r 1 ,
where ϕ k is defined similarly to (8) for Z V . On the other hand, if:
ξ = ξ 11 I n 1 t V V ξ ( ξ 11 R , V W , ξ Z V ) ,
we observe that:
ϕ k ( ξ ) = ϕ k 1 ( ξ ) ( k = 2 , , r ) .
Therefore, in order to prove (9) for P V * of rank r, it suffices to show that:
P V * = ξ Z V ; ξ P V *  and  ϕ 1 ( ξ ) is   positive   definite = ξ Z V ; ξ P V *  and  det ϕ 1 ( ξ ) > 0 .
If q 1 = 0 , then any element ξ Z V is of the form:
ξ = ξ 11 I n 1 ξ ,
which belongs to P V if and only if ξ P V and ϕ 1 ( ξ ) = ξ 11 > 0 , so that (11) holds.
Assume q 1 > 0 . Keeping in mind that W ˜ 1 R W and W R q 1 by (4), we have for ξ Z V as in (10),
ϕ 1 ( ξ ) = ξ 11 t v v ψ ( ξ ) Sym ( 1 + q 1 , R ) ,
where v = vect ( V ) R q i and ψ : Z V Sym ( q i , R ) is defined in such a way that:
( U t U | ξ ) = t u ψ ( ξ ) u ( ξ Z V , U W , u = vect ( U ) R q 1 ) .
On the other hand, for x Z V as in (6), we have:
( x | ξ ) = x 11 ξ 11 + 2 x 11 t u ˜ v + x 11 t u ˜ ψ ( ξ ) u ˜ + ( x ˜ | ξ ) .
Owing to Lemma 1 (iii), the element ξ Z V belongs to P V * if and only if the right-hand side is strictly positive for all x 11 0 , U ˜ W and x ˜ P V ¯ with ( x 11 , x ˜ ) ( 0 , 0 ) . Assume ξ P V * . Considering the case x 11 = 0 , we have ( x ˜ | ξ ) > 0 for all x ˜ P V ¯ \ { 0 } , which means that ξ P V * . Then, the quantity in (13) is strictly positive for non-zero U because U t U belongs to P V ¯ \ { 0 } . Thus, ψ ( ξ ) is positive definite, and (14) is rewritten as:
( x | ξ ) = x 11 ( ξ 11 t v ψ ( ξ ) 1 v ) + x 11 t ( u ˜ + ψ ( ξ ) 1 v ) ψ ( ξ ) ( u ˜ + ψ ( ξ ) 1 v ) + ( x ˜ | ξ ) .
Therefore, we obtain:
P V * = ξ Z V ; ξ P V *  and  ξ 11 t v ψ ( ξ ) 1 v > 0 .
On the other hand, we see from (12) that:
ϕ 1 ( ξ ) = 1 t v ψ ( ξ ) 1 I q 1 ξ 11 t v ψ ( ξ ) 1 v ψ ( ξ ) 1 ψ ( ξ ) 1 v I q 1 .
Hence, we deduce (11) from (16) and (17). ☐
We note that, if q 1 > 0 , the ( 1 , 1 ) -component of the inverse matrix ϕ 1 ( ξ ) 1 is given by:
( ϕ 1 ( ξ ) 1 ) 11 = ( ξ 11 t v ψ ( ξ ) 1 v ) 1
thanks to (17).

3. Koszul–Vinberg Characteristic Function of P V *

We denote by φ V the Koszul–Vinberg characteristic function of P V * . In this section, we give an explicit formula of φ V .
Recall that the linear map ψ : Z V Sym ( q 1 , R ) plays an important role in the proof of Theorem 1. We shall introduce similar linear maps ψ k : Z V Sym ( q k , R ) for k such that q k > 0 . Let W k be the subspace of W ˜ k consisting of W W ˜ k for which w 1 = x k k = 0 . Then, clearly, W k l > k V l k and dim W k = q k . If q k > 0 , using the same orthogonal basis of V l k as in the previous section, we have the isomorphism W k W w = vect ( W ) R q k . Similarly to (8), we define ψ k by:
( W t W | ξ ) = t w ψ k ( ξ ) w ( ξ Z V , W W k , w = vect ( W ) R q k ) .
Then, we have:
ϕ k ( ξ ) = ξ k k t v k v k ψ k ( ξ ) ( ξ Z V ) ,
where v k R q k is a vector corresponding to the W k -component of ξ. If ξ P V * , we see from (19) that ψ k ( ξ ) is positive definite. In this case, we have:
ϕ k ( ξ ) = 1 t v k ψ k ( ξ ) 1 I q k ξ k k t v k ψ k ( ξ ) 1 v k ψ k ( ξ ) 1 ψ k ( ξ ) 1 v k I q k ,
so that we get a generalization of (18), that is,
( ϕ k ( ξ ) 1 ) 11 = ( ξ k k t v k ψ k ( ξ ) 1 v k ) 1 .
On the other hand, if q k = 0 , then ϕ k ( ξ ) 1 = ξ k k 1 .
We remark that ψ 1 ( ξ ) = ψ ( ξ ) , and that some part of the argument above is parallel to the proof of Theorem 1.
Theorem 2.
The Koszul–Vinberg characteristic function φ V of P V * is given by the following formula:
φ V ( ξ ) = C V k = 1 r ϕ k ( ξ ) 1 11 1 + q k / 2 q k > 0 ( det ψ k ( ξ ) ) 1 / 2 ( ξ P V * ) ,
where C V : = ( 2 π ) ( N r ) / 2 k = 1 r Γ ( 1 + q k 2 ) and N : = dim Z V .
Proof. 
We shall show the statement by induction on the rank as in the proof of Theorem 1. Then, it suffices to show that:
φ V ( ξ ) = ( 2 π ) q 1 / 2 Γ ( 1 + q 1 2 ) ( ϕ 1 ( ξ ) 1 ) 11 1 + q 1 / 2 ( det ψ 1 ( ξ ) ) sgn ( q 1 ) / 2 φ V ( ξ )
for ξ P V * as in (10), where ( det ψ 1 ( ξ ) ) sgn ( q 1 ) / 2 is interpreted as:
( det ψ 1 ( ξ ) ) sgn ( q 1 ) / 2 : = 1 ( q 1 = 0 ) , ( det ψ 1 ( ξ ) ) 1 / 2 ( q 1 > 0 ) .
When q 1 = 0 , we have:
φ V ( ξ ) = 0 P V e x 11 ξ 11 e ( x | ξ ) d x 11 d x = ξ 11 1 φ V ( ξ ) ,
which means (24).
When q 1 > 0 , the Euclidean measure d x equals 2 q 1 / 2 x 11 q 1 d x 11 d u ˜ d x ˜ by the change of variables in (6). Indeed, the coefficient 2 q 1 / 2 comes from the normalization of the inner product on W R q 1 regarded as a subspace of Z V . Then, we have by (15):
φ V ( ξ ) = 0 R q 1 P V e x 11 ( ξ 11 t v ψ ( ξ ) 1 v ) e x 11 t ( u ˜ + ψ ( ξ ) 1 v ) ψ ( ξ ) ( u ˜ + ψ ( ξ ) 1 v ) e ( x ˜ | ξ ) × 2 q 1 / 2 x 11 q 1 d x 11 d u ˜ d x ˜ .
By the Gaussian integral formula, we have:
R q 1 e x 11 t ( u ˜ + ψ ( ξ ) 1 v ) ψ ( ξ ) ( u ˜ + ψ ( ξ ) 1 v ) d u ˜ = π q 1 / 2 x 11 q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 .
Therefore, we get:
φ V ( ξ ) = ( 2 π ) q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 0 e x 11 ( ξ 11 t v ψ ( ξ ) 1 v ) x 11 q 1 / 2 d x 11 P V e ( x ˜ | ξ ) d x ˜ = ( 2 π ) q 1 / 2 ( det ψ 1 ( ξ ) ) 1 / 2 Γ ( 1 + q 1 2 ) ( ξ 11 t v ψ ( ξ ) 1 v ) 1 q k / 2 φ V ( ξ ) ,
which together with (18) leads us to (24). ☐
Example 2.
Let V = { V l k } 1 k < l 3 be as in Example 1. For:
ξ = ξ 1 0 ξ 4 0 0 ξ 1 0 ξ 5 ξ 4 0 ξ 2 ξ 6 0 ξ 5 ξ 6 ξ 3 Z V ,
we have:
ϕ 1 ( ξ ) = ξ 1 ξ 4 ξ 5 ξ 4 ξ 2 0 ξ 5 0 ξ 3 , ϕ 2 ( ξ ) = ξ 2 ξ 6 ξ 6 ξ 3 , ϕ 3 ( ξ ) = ξ 3 , ψ 1 ( ξ ) = ξ 2 0 0 ξ 3 , ψ 2 ( ξ ) = ξ 3 .
The cone P V * is described as:
P V * = ξ Z V ; ξ 1 ξ 4 ξ 5 ξ 4 ξ 2 0 ξ 5 0 ξ 3 > 0 , ξ 2 ξ 6 ξ 6 ξ 3 > 0 , ξ 3 > 0 ,
and its Koszul–Vinberg characteristic function φ V is expressed as:
φ V ( ξ ) = C V ξ 1 ξ 4 ξ 5 ξ 4 ξ 2 0 ξ 5 0 ξ 3 / ( ξ 2 ξ 3 ) 2 ξ 2 ξ 6 ξ 6 ξ 3 / ξ 3 3 / 2 ξ 3 1 · ( ξ 2 ξ 3 ) 1 / 2 ( ξ 3 ) 1 / 2 = C V ξ 1 ξ 4 ξ 5 ξ 4 ξ 2 0 ξ 5 0 ξ 3 2 ξ 2 ξ 6 ξ 6 ξ 3 3 / 2 ξ 2 3 / 2 ξ 3 3 / 2 ,
where C V = ( 2 π ) 3 / 2 Γ ( 2 ) Γ ( 3 / 2 ) Γ ( 1 ) = 2 π 2 .
Suppose that the cone P V is homogeneous. Then, P V * , as well as P V , is a homogeneous cone of rank 3, so that the Koszul–Vinberg characteristic function of P V * has at most three irreducible factors (see [8]). However, we have seen that there are four irreducible factors in the function φ V . Therefore, we conclude that neither P V , nor P V * is homogeneous.

4. Γ-Type Integral Formulas

For an n × n matrix A = ( A i j ) and 1 m n , we denote by A [ m ] the upper-left m × m submatrix ( A i j ) i , j m of A. Put M k : = i = 1 k n k ( k = 1 , , r ) . For s ̲ = ( s 1 , , s r ) C r , we define functions Δ s ̲ V on P V and δ s ̲ V on P V * respectively by:
Δ s ̲ V ( x ) : = ( det x [ M 1 ] ) s 1 / n 1 k = 2 r ( det x [ M k ] det x [ M k 1 ] ) s k / n k = ( det x ) s r / n r k = 1 r 1 ( det x [ M k ] ) s k / n k s k 1 / n k 1 ( x P V ) ,
δ s ̲ V ( ξ ) : = k = 1 r ( ϕ k ( ξ ) 1 ) 11 s k = q k = 0 ξ k k s k q k > 0 ( ξ k k t v k ψ k ( ξ ) 1 v k ) s k ( ξ P V * ) .
Recall (22) for the second equality of (27).
For a ̲ = ( a 1 , , a r ) R > 0 r , let D a ̲ denote the diagonal matrix defined by:
D a ̲ : = a 1 I n 1 a 2 I n 2 a r I n r G L ( n , R ) .
Then, the linear map Z V x D a x D a Z V preserves both P V and P V * , and we have:
Δ s ̲ V ( D a x D a ) = ( k = 1 r a k 2 s k ) Δ s ̲ V ( x ) ( x P V ) ,
δ s ̲ V ( D a ξ D a ) = ( k = 1 r a k 2 s k ) δ s ̲ V ( ξ ) ( ξ P V ) .
Assume q 1 > 0 . For B W , we denote by τ B the linear transform on Z V given by:
τ B x : = I n 1 B I n n 1 x 11 I n 1 t U U x I n 1 t B I n n 1 = x 11 I n 1 t U + x 11 t B U + x 11 B x + U t B + B t U + x 11 B t B ,
where x Z V is as in (3). Indeed, since:
U t B + B t U = ( U + B ) t ( U + B ) U t U B t B Z V ,
the matrix τ B x belongs to Z V . Clearly, τ B preserves P V , and we have:
Δ s ̲ V ( τ B x ) = Δ s ̲ V ( x ) ( x P V ) .
The formula (5) is rewritten as:
τ x 11 1 U ( x ) = x 11 I n 1 x x 11 1 U t U ,
which together with (30) tells us that:
Δ s ̲ V ( x ) = x 11 s 1 Δ s ̲ V ( x x 11 1 U t U ) ,
where s ̲ : = ( s 2 , , s r ) C r 1 .
Let us consider the adjoint map τ B * : Z V Z V of τ B with respect to the standard inner product. Let b R q 1 be the vector corresponding to B W . For x Z V and ξ Z V as in (3) and (10), respectively, we observe that:
( τ B x | ξ ) = x 11 ξ 11 + 2 t ( u + x 11 b ) v + ( x + U t B + B t U + x 11 B t B | ξ ) = x 11 ( ξ 11 + 2 t b v + t b ψ ( ξ ) b ) + 2 t u ( v + ψ ( ξ ) b ) + ( x | ξ ) .
Thus, if we write:
ι ( ξ 11 , v , ξ ) : = ξ 11 I n 1 t V V ξ ,
we have:
τ B * ι ( ξ 11 , v , ξ ) = ι ( ξ 11 + 2 t b v + t b ψ ( ξ ) b , v + ψ ( ξ ) b , ξ ) .
Furthermore, we see from (12) that ϕ 1 ( τ B * ι ( ξ 11 , v , ξ ) ) equals:
ξ 11 + 2 t b v + t b ψ ( ξ ) b t v + t b ψ ( ξ ) v + ψ ( ξ ) b ψ ( ξ ) = 1 t b I q 1 ξ 11 t v v ψ ( ξ ) 1 b I q 1 ,
so that we get for ξ = ι ( ξ 11 , v , ξ ) :
ϕ 1 ( τ B * ξ ) = 1 t b I q 1 ϕ 1 ( ξ ) 1 b I q 1 .
Therefore:
( ϕ 1 ( τ B * ξ ) 1 ) 11 = ( ϕ 1 ( ξ ) 1 ) 11 .
On the other hand, we have for ξ = ι ( ξ 11 , v , ξ ) P V * :
δ s ̲ V ( ξ ) = ( ϕ 1 ( ξ ) 1 ) 11 s 1 δ s ̲ V ( ξ ) .
Thus, we conclude that:
δ s ̲ V ( τ B * ξ ) = δ s ̲ V ( ξ ) .
Theorem 3.
When s k > 1 q k / 2 for k = 1 , , r , one has:
P V e ( x | ξ ) Δ s ̲ V ( x ) d x = C V 1 γ V ( s ̲ ) δ s ̲ V ( ξ ) φ V ( ξ ) ,
where γ V ( s ̲ ) : = ( 2 π ) ( N r ) / 2 k = 1 r Γ ( s k + 1 + q k 2 ) .
Proof. 
Recalling Theorem 2, we rewrite the right-hand side of (35) as:
( 2 π ) ( N r ) / 2 k = 1 r Γ ( s k + 1 + q k 2 ) k = 1 r ϕ k ( ξ ) 1 11 s k + 1 + q k / 2 q k > 0 ( det ψ k ( ξ ) ) 1 / 2 ,
which is similar to the right-hand side of (23). Thus, the proof is parallel to Theorem 2. Namely, by induction on the rank, it suffices to show that:
P V e ( x | ξ ) Δ s ̲ V ( x ) d x = ( 2 π ) q 1 / 2 Γ ( s 1 + 1 + q 1 2 ) ( ϕ 1 ( ξ ) 1 ) 11 s 1 + 1 + q 1 / 2 ( det ψ 1 ( ξ ) ) sgn ( q 1 ) / 2 × P V e ( x | ξ ) Δ s ̲ V ( x ) d x
thanks to (33).
When q 1 = 0 , we have ( x | ξ ) = x 11 ξ 11 + ( x | ξ ) and Δ s ̲ V ( x ) = x 11 s 1 Δ s ̲ V ( x ) . Thus:
P V e ( x | ξ ) Δ s ̲ V ( x ) d x = 0 e x 11 ξ 11 x 11 s 1 d x 11 × P V e ( x | ξ ) Δ s ̲ V ( x ) d x .
Since 0 e x 11 ξ 11 x 11 s 1 d x 11 = Γ ( s 1 + 1 ) ξ 11 s 1 1 , we get (36).
When q 1 > 1 , we use the change of variable (6). Since x ˜ = x x 11 1 U t U , we have Δ s ̲ V ( x ) = x 11 s 1 Δ s ̲ V ( x ˜ ) by (31). Therefore, by the same Gaussian integral formula as in the proof of Theorem 2, the integral P V e ( x | ξ ) Δ s ̲ V ( x ) d x equals:
0 W P V e x 11 ( ξ 11 t v ψ ( ξ ) 1 v ) e x 11 t ( u ˜ + ψ ( ξ ) 1 v ) ψ ( ξ ) ( u ˜ + ψ ( ξ ) 1 v ) e ( x ˜ | ξ ) x 11 s 1 Δ s ̲ V ( x ˜ ) × 2 q 1 / 2 x 11 q 1 d x 11 d u ˜ d x ˜ = ( 2 π ) q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 0 e x 11 ( ξ 11 t v ψ ( ξ ) 1 v ) x 11 s 1 + q 1 / 2 d x 11 × P V e ( x ˜ | ξ ) Δ s ̲ V ( x ˜ ) d x ˜ = ( 2 π ) q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 Γ ( s k + 1 + q 1 2 ) ( ξ 11 t v ψ ( ξ ) 1 v ) s k 1 q k / 2 × P V e ( x ˜ | ξ ) Δ s ̲ V ( x ˜ ) d x ˜ .
Hence, we get (36) by (18). ☐
We shall obtain an integral formula over P V * as follows.
Theorem 4.
When s k > q k / 2 for k = 1 , , r , one has:
P V * e ( x | ξ ) δ s ̲ V ( ξ ) φ V ( ξ ) d ξ = C V Γ V ( s ̲ ) Δ s ̲ V ( x ) ( x P V ) ,
where Γ V ( s ̲ ) : = ( 2 π ) ( N r ) / 2 k = 1 r Γ ( s k q k / 2 ) .
Proof. 
Using (24), (31) and (33), we rewrite (37) as:
P V * e ( x | ξ ) ( ϕ 1 ( ξ ) 1 ) 11 s 1 + 1 + q 1 / 2 ( det ψ 1 ( ξ ) ) sgn ( q 1 ) / 2 δ s ̲ V ( ξ ) φ V ( ξ ) d ξ = C V ( 2 π ) q 1 / 2 Γ ( s 1 q 1 / 2 ) Γ V ( s ̲ ) x 11 s 1 Δ s ̲ V ( x ˜ ) ,
where:
x ˜ : = x ( q 1 = 0 ) , x x 11 1 U t U ( q 1 > 0 ) .
Therefore, by induction on the rank, it suffices to show that the left-hand side of (38) equals:
( 2 π ) q 1 / 2 Γ ( s 1 q 1 / 2 ) x 11 s 1 P V * e ( x ˜ | ξ ) δ s ̲ V ( ξ ) φ V ( ξ ) d ξ .
When q 1 = 0 , since d ξ = d ξ 11 d ξ , the left-hand side of (38) equals:
0 e x 11 ξ 11 ξ 11 s 1 1 d ξ 11 P V * e ( x | ξ ) δ s ̲ V ( ξ ) φ V ( ξ ) d ξ ,
which coincides with (39) in this case.
Assume q 1 > 0 . Keeping (16) and (18) in mind, we put ξ ˜ 11 : = ξ 11 t v ψ ( ξ ) 1 v = ( ϕ 1 ( ξ ) 1 ) 11 1 > 0 . By the change of variables ξ = ι ( ξ ˜ 11 + t v ψ ( ξ ) 1 v , v , ξ ) , we have d ξ = 2 q 1 / 2 d ξ ˜ 11 d v d ξ . On the other hand, we observe:
( x | ξ ) = x 11 ( ξ ˜ 11 + t v ψ ( ξ ) 1 v ) + 2 t u v + ( x | ξ ) = x 11 ξ ˜ 11 + x 11 t ( v + x 11 1 ψ ( ξ ) u ) ψ ( ξ ) 1 ( v + x 11 1 ψ ( ξ ) u ) + ( x x 11 1 U t U | ξ ) .
Thus, the left-hand side of (39) equals:
0 R q 1 P V * e x 11 ξ ˜ 11 e x 11 t ( v + x 11 1 ψ ( ξ ) u ) ψ ( ξ ) 1 ( v + x 11 1 ψ ( ξ ) u ) e ( x x 11 1 U t U | ξ ) × ξ ˜ 11 s 1 1 q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 δ s ̲ V ( ξ ) φ V ( ξ ) 2 q 1 / 2 d ξ ˜ 11 d v d ξ .
By the Gaussian integral formula, we have:
R q 1 e x 11 t ( v + x 11 1 ψ ( ξ ) u ) ψ ( ξ ) 1 ( v + x 11 1 ψ ( ξ ) u ) d v = π q 1 / 2 x 11 q 1 / 2 ( det ψ ( ξ ) ) 1 / 2 ,
so that (40) equals:
( 2 π ) q 1 / 2 x 11 q 1 / 2 0 e x 11 ξ ˜ 11 ξ ˜ 11 s 1 1 q 1 / 2 d ξ ˜ 11 P V * e ( x x 11 1 U t U | ξ ) δ s ̲ V ( ξ ) φ V ( ξ ) d ξ ,
which coincides with (39) because: 0 e x 11 ξ ˜ 11 ξ ˜ 11 s 1 1 q 1 / 2 d ξ ˜ 11 = Γ ( s 1 q 1 / 2 ) x 11 s 1 + q 1 / 2 . ☐
Example 3.
Let Z V be as in Example 1, and let x P V and ξ P V * be as in (2) and (25), respectively. Then, we have for s ̲ = ( s 1 , s 2 , s 3 ) C 3 ,
Δ s ̲ V ( x ) = ( x 11 2 ) s 1 / 2 s 2 x 1 0 x 4 0 x 1 0 x 4 0 x 2 s 2 s 3 ( det x ) s 3 = x 11 s 1 s 2 s 3 x 1 x 4 x 4 x 2 s 2 s 3 ( det x ) s 3 ,
and:
δ s ̲ V ( ξ ) = ( ξ 1 ξ 4 2 ξ 2 ξ 5 2 ξ 3 ) s 1 ( ξ 2 ξ 6 2 ξ 3 ) s 2 ξ 3 s 3 = ξ 1 ξ 4 ξ 5 ξ 4 ξ 2 0 ξ 5 0 ξ 3 s 1 ξ 2 ξ 6 ξ 6 ξ 3 s 2 ξ 2 s 1 ξ 3 s 3 s 1 s 2 .
When s 1 > 2 , s 2 > 3 / 2 and s 3 > 1 , the integral formula (35) holds with:
γ V ( s ̲ ) = ( 2 π ) 3 / 2 Γ ( s 1 + 2 ) Γ ( s 2 + 3 / 2 ) Γ ( s 3 + 1 ) .
Furthermore, when s 1 > 1 , s 2 > 1 / 2 and s 3 > 0 , the integral formula (37) holds with:
Γ V ( s ̲ ) = ( 2 π ) 3 / 2 Γ ( s 1 1 ) Γ ( s 2 1 / 2 ) Γ ( s 3 ) .

5. Multiplicative Legendre Transform of Generalized Power Functions

For s ̲ R > 0 r , we see that log Δ s ̲ is a strictly convex function on the cone P V . In fact, Δ s ̲ is defined naturally on P n as a product of powers of principal minors, and it is well known that such log Δ s ̲ is strictly convex on the whole P n . In this section, we shall show that log Δ s ̲ V and log δ s ̲ V are related by the Fenchel–Legendre transform.
For x P V , we denote by I s ̲ V ( x ) the minus gradient log Δ s ̲ ( x ) at x with respect to the inner product. Namely, I s ̲ V ( x ) is an element of Z V for which:
( I s ̲ V ( x ) | y ) = d d t t = 0 log Δ s ̲ ( x + t y ) ( y Z V ) .
Similarly, J s ̲ V ( ξ ) : = log δ s ̲ ( ξ ) is defined for ξ P V * . If q 1 > 0 , then for any B W , we have:
I s ̲ V τ B = τ B * I s ̲ V ,
J s ̲ V τ B * = τ B J s ̲ V
owing to (30) and (34), respectively.
Theorem 5.
For any s ̲ R > 0 r , the map I s ̲ V : P V Z V gives a diffeomorphism from P V onto P V * , and J s ̲ V gives the inverse map.
Proof. 
We shall prove the statement by induction on the rank. When r = 1 , we have I s ̲ V ( x 11 I n 1 ) = s 1 x 11 I n 1 = J s ̲ V ( x 11 I n 1 ) for x 11 > 0 . Thus, the statement is true in this case.
When r > 1 , assume that the statement holds for the system of rank r 1 . Let Z V 0 be the subspace of Z V defined by:
Z V 0 : = x 11 I n 1 0 0 x ; x 11 R , x Z V .
By direct computation with (31) and (33), we have:
I s ̲ V x 11 I n 1 0 0 x = s 1 x 11 I n 1 0 0 I s ̲ V ( x ) ,
J s ̲ V ξ 11 I n 1 0 0 ξ = s 1 ξ 11 I n 1 0 0 J s ̲ V ( ξ )
for x 11 , ξ 11 > 0 , x P V and ξ P V * . By the induction hypothesis, we see that I s ̲ V : P V Z V 0 P V * Z V 0 is bijective with the inverse map J s ̲ V : P V * Z V 0 P V Z V 0 .
If q 1 = 0 , the statement holds because Z V = Z V 0 . Assume q 1 > 0 . Lemma 1 (ii) tells us that, for x P V , there exist unique x 0 Z V 0 P V and B W for which x = τ B x 0 . Similarly, we see from (32) that, for ξ P V * , there exist unique ξ 0 Z V 0 P V * and C W for which ξ = τ C * ξ 0 . Therefore, we deduce from (41) and (42) that I s ̲ V : P V P V * is a bijection with J s ̲ V the inverse map. ☐
Proposition 1.
Let s ̲ R > 0 r . For ξ P V * , one has:
Δ s ̲ ( J s ̲ ( ξ ) ) 1 = ( k = 1 r s k s k ) δ s ̲ ( ξ ) .
Proof. 
We prove the statement by induction on the rank. When r = 1 , the equality (45) is verified directly. Indeed, the left-hand side of (45) is computed as ( s 1 ξ 11 ) s 1 = s 1 s 1 ξ 11 s 1 .
When r > 1 , assume that (45) holds for a system of rank r 1 . We deduce from (31), (33), (43), (44) and the induction hypothesis that (45) holds for ξ P V * Z V 0 . Therefore, (45) holds for all ξ P V * by (30), (34) and (42). ☐
In general, for a non-zero function f, the function 1 f ( log f ) 1 is called the multiplicative Legendre transform of f. Thanks to Theorem 5 and Proposition 1, we see that the multiplicative Legendre transform of Δ s ̲ ( x ) is equal to δ s ̲ ( ξ ) on P V * up to constant multiple. As a corollary, we arrive at the following result.
Theorem 6.
The Fenchel–Legendre transform of the convex function log Δ s ̲ on P V is equal to the function log δ s ̲ ( ξ ) of ξ P * up to constant addition.

6. Application to Statistics and Optimization

Take s ̲ R r for which s k > q k / 2 ( k = 1 , , r ) . We define a measure ρ s ̲ V on P V * by:
ρ s ̲ V ( d ξ ) : = C V 1 Γ V ( s ̲ ) 1 δ s ̲ V ( ξ ) φ V ( ξ ) d ξ ( ξ P V * ) .
Theorem 4 states that:
P V * e ( x | ξ ) ρ s ̲ V ( d ξ ) = Δ s ̲ V ( x ) ( x P V ) .
Then, we obtain the natural exponential family generated by ρ s ̲ V , that is a family { μ s ̲ , x V } x P V of probability measures on P V * given by:
μ s ̲ , x V ( d ξ ) : = Δ s ̲ V ( x ) e ( x | ξ ) ρ s ̲ V ( d ξ ) .
In particular, when s ̲ = ( n 1 α , n 2 α , , n r α ) for sufficiently large α, we have μ s ̲ , x V ( d ξ ) = ( det x ) α e ( x | ξ ) ρ s ̲ V ( d ξ ) . We call μ s ̲ , x V the Wishart distributions on P V * in general.
From a sample ξ 0 P V * , let us estimate the parameter x P V in such a way that the likelihood function Δ s ̲ V ( x ) e ( x | ξ ) attains its maximum at the estimator x 0 . Then, we have the likelihood equation ξ 0 = I s ̲ V ( x 0 ) , whereas Theorem 5 gives a unique solution by x 0 = J s ̲ V ( ξ 0 ) .
The same argument leads us to the following result in semidefinite programming. For a fixed ξ 0 P V * and α > 0 , a unique solution x 0 of the minimization problem of ( x | ξ 0 ) α log det x subject to x P V = Z V P n is given by x 0 = J s ̲ V ( ξ 0 ) , where s ̲ = ( n 1 α , , n r α ) . Note that J s ̲ V is a rational map because δ s ̲ V is a product of powers of rational functions.

7. Special Cases

7.1. Matrix Realization of Homogeneous Cones

Let us assume that the system V = { V l k } 1 k < l r satisfies not only the conditions (V1) and (V2), but also the following:
(V3) A V l k , B V k j A B V l j ( 1 j < k < l r ) .
Then, the set H V of lower triangular matrices T of the form:
T = T 11 T 21 T 22 T r 1 T r 2 T r r
becomes a linear Lie group, and H V acts on the space Z V by ρ ( T ) x : = T x t T ( T H V , x Z V ) . The group H V acts on the cone P V simply transitively by this action ρ, so that P V is a homogeneous cone. Moreover, it is shown in [15] that every homogeneous cone is linearly isomorphic to such P V (see also [18]).
Let V 0 = { V l k 0 } 1 k < l 3 be the system given by V 21 0 = { 0 } and V l k 0 = R ( ( l , k ) ( 2 , 1 ) ) . Then:
Z V 0 = x 1 0 x 4 0 x 2 x 5 x 4 x 5 x 3 ; x 1 , , x 5 R ,
and P V 0 : = Z V 0 P 3 is homogeneous because (V1)–(V3) are satisfied in this case. On the other hand, let V 1 = { V l k 1 } 1 k < l 3 be the system given by V 31 1 = { 0 } and V l k 1 = R ( ( l , k ) ( 3 , 1 ) ) . Then:
Z V 1 = x 1 x 4 0 x 4 x 2 x 5 0 x 5 x 3 ; x 1 , , x 5 R .
Note that V 1 satisfies only (V1) and (V2), but P V 1 is homogeneous because P V 1 is isomorphic to the homogeneous cone P V 0 via the map:
P V 1 x 1 x 4 0 x 4 x 2 x 5 0 x 5 x 3 1 0 0 0 0 1 0 1 0 x 1 x 4 0 x 4 x 2 x 5 0 x 5 x 3 1 0 0 0 0 1 0 1 0 = x 1 0 x 4 0 x 3 x 5 x 4 x 5 x 2 P V 0 .
This example tells us that our matrix realization of a convex cone is not unique and that the condition (V3) is merely a sufficient condition for the homogeneity of the cone.
Many ideas in this work are inspired by the theory of homogeneous cones. The notion of generalized power functions, as well as the Γ-type integral formulas are due to Gindikin [8] (see also [23]). The Wishart distributions for homogeneous cones are studied in [17,21,24,25].

7.2. Cones Associated with Chordal Graphs

If n 1 = n 2 = = n r = 1 , then V l k equals either R or { 0 } . In this case, Z V is the space of symmetric matrices with prescribed zero components. Such a space is described by using an undirected graph in the graphical model theory.
Let us recall some notion in the graph theory. Let G be a graph and V G the set of vertices of G. We assume that G has no multiple edge, that is, for any two vertices i , j V G , either there is one edge connecting them or there is no edge between them. These relations of the vertices i and j are denoted by i j and i j , respectively. Assume further that G has no loop, which means that i i for i V G . We define the edge set E G V G × V G by:
E G : = ( i , j ) V G × V G ; i j .
Since V G and E G have all of the information of G, the graph G is often identified with the pair ( V G , E G ) . For a non-empty subset V of V G , put E : = E G ( V × V ) . The graph G : = ( V , E ) is called an induced subgraph of G. The graph G is said to be chordal or decomposable if G contains no cycle of length greater than three as an induced subgraph, and said to be A 4 -free if G contains no A 4 graph as an induced subgraph. Let ⪯ be a total order on the vertex set V G , and for i V G , put V G [ i ] : = j V G ; i j  and  i j V G . Then, ⪯ is said to be an eliminating order on the graph G if the induced subgraph with the vertex set V G [ i ] is complete for each i V G . It is known that there exists an eliminating order on G if and only if the graph G is chordal.
Let us identify the vertex set V G with { 1 , 2 , , r } . Let Z G be the space of symmetric matrices x = ( x i j ) Sym ( r , R ) , such that, if i j and i j , then x i j = 0 . Define P G : = Z G P r . We can show ([11] (Theorem 2.2), [26]) that the cone P G is homogeneous if and only if the graph G is chordal and A 4 -free. On the other hand, it is known in the graphical model theory as well as the sparse matrix linear algebra that even though P G is not homogeneous, various formulas still hold for P G if G is chordal.
The cone P G is expressed as P V with n 1 = n 2 = = n r = 1 and:
V j i = R ( j i ) , { 0 } ( j i ) .
Then, the condition (V2) means exactly that the order ≤ is an eliminating order on G. Therefore, any cone P G with chordal G can be treated as P V in our framework. Most of the integral formulas for P G in [11,27] can be deduced from Theorems 3 and 4, while the Wishart distribution is a central object in the theory of graphical model. In [28], the analysis for generalized power functions associated with all eliminating orders is discussed for a specific graph A n : by direct computations.

Acknowledgments

The author would like to express sincere gratitude to Piotr Graczyk and Yoshihiko Konno for stimulating discussions about this subject. He is also grateful to Frédéric Barbaresco for his interest in and encouragement of this work. He thanks to anonymous referees for valuable comments, which were helpful for the improvement of the present paper. This work was supported by JSPS KAKENHI Grant Number 16K05174.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koszul, J.L. Ouverts convexes homogènes des espaces affines. Math. Z. 1962, 79, 254–259. [Google Scholar] [CrossRef]
  2. Vinberg, E.B. The theory of convex homogeneous cones. Trans. Moscow Math. Soc. 1963, 12, 340–403. [Google Scholar]
  3. Vey, J. Sur les automorphismes affines des ouverts convexes saillants. Annali della Scuola Normale Superiore di Pisa 1970, 24, 641–665. [Google Scholar]
  4. Shima, H. The Geometry of Hessian Structures; World Scientific: Hackensack, NJ, USA, 2007. [Google Scholar]
  5. Nesterov, Y.; Nemirovskii, A. Interior-Point Polynomial Algorithms in Convex Programming; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1994. [Google Scholar]
  6. Barbaresco, F. Koszul information geometry and Souriau geometric temperature/capacity of Lie group thermodynamics. Entropy 2014, 16, 4521–4565. [Google Scholar] [CrossRef]
  7. Barbaresco, F. Symplectic structure of information geometry: Fisher etric and Euler-Poincaré equation of Souriau Lie group thermodynamics. In Geometric Science of Information; Nielsen, F., Barbaresco, F., Eds.; (Lecture Notes in Computer Science); Springer International Publishing: Basel, Switzerland, 2015; Volume 9389, pp. 529–540. [Google Scholar]
  8. Gindikin, S.G. Analysis in homogeneous domains. Russ. Math. Surv. 1964, 19, 1–89. [Google Scholar] [CrossRef]
  9. Truong, V.A.; Tunçel, L. Geometry of homogeneous convex cones, duality mapping, and optimal self-concordant barriers. Math. Program. 2004, 100, 295–316. [Google Scholar] [CrossRef]
  10. Tunçel, L.; Xu, S. On homogeneous convex cones, the Caratheodory number, and the duality mapping. Math. Oper. Res. 2001, 26, 234–247. [Google Scholar] [CrossRef]
  11. Letac, G.; Massam, H. Wishart distributions for decomposable graphs. Ann. Stat. 2007, 35, 1278–1323. [Google Scholar] [CrossRef]
  12. Rothaus, O.S. The construction of homogeneous convex cones. Ann. Math. 1966, 83, 358–376. [Google Scholar] [CrossRef]
  13. Xu, Y.-C. Theory of Complex Homogeneous Bounded Domains; Kluwer: Dordrecht, The Netherlands, 2005. [Google Scholar]
  14. Chua, C.B. Relating homogeneous cones and positive definite cones via T-algebras. SIAM J. Optim. 2003, 14, 500–506. [Google Scholar] [CrossRef]
  15. Ishi, H. On symplectic representations of normal j-algebras and their application to Xu’s realizations of Siegel domains. Differ. Geom. Appl. 2006, 24, 588–612. [Google Scholar] [CrossRef]
  16. Yamasaki, T.; Nomura, T. Realization of homogeneous cones through oriented graphs. Kyushu J. Math. 2015, 69, 11–48. [Google Scholar] [CrossRef]
  17. Ishi, H. Homogeneous cones and their applications to statistics. In Modern Methods of Multivariate Statistics; Graczyk, P., Hassairi, A., Eds.; Hermann: Paris, France, 2014; pp. 135–154. [Google Scholar]
  18. Ishi, H. Matrix realization of homogeneous cones. In Geometric Science of Information; Nielsen, F., Barbaresco, F., Eds.; (Lecture Notes in Computer Science); Springer International Publishing: Basel, Switzerland, 2015; Volume 9389, pp. 248–256. [Google Scholar]
  19. Lauritzen, S.L. Graphical Models; Clarendon Press: Oxford, UK, 1996. [Google Scholar]
  20. Faraut, J.; Korányi, A. Analysis on Symmetric Cones; Clarendon Press: Oxford, UK, 1994. [Google Scholar]
  21. Graczyk, P.; Ishi, H. Riesz measures and Wishart laws associated with quadratic maps. J. Math. Soc. Jpn. 2014, 66, 317–348. [Google Scholar] [CrossRef]
  22. Shima, H. Homogeneous Hessian manifolds. Ann. Inst. Fourier 1980, 30, 91–128. [Google Scholar] [CrossRef]
  23. Güler, O.; Tunçel, L. Characterization of the barrier parameter of homogeneous convex cones. Math. Program. A 1988, 81, 55–76. [Google Scholar] [CrossRef]
  24. Andersson, S.A.; Wojnar, G.G. Wishart distributions on homogeneous cones. J. Theor. Probab. 2004, 17, 781–818. [Google Scholar] [CrossRef]
  25. Graczyk, P.; Ishi, H.; Kołodziejek, B. Wishart exponential families and variance function on homogeneous cones. HAL, 2016; submitted for publication. [Google Scholar]
  26. Ishi, H. On a class of homogeneous cones consisting of real symmetric matrices. Josai Math. Monogr. 2013, 6, 71–80. [Google Scholar]
  27. Roverato, A. Cholesky decomposition of a hyper inverse Wishart matrix. Biometrika 2000, 87, 99–112. [Google Scholar] [CrossRef]
  28. Graczyk, P.; Ishi, H.; Mamane, S. Wishart exponential families on cones related to An graphs. HAL, 2016; submitted for publication. [Google Scholar]

Share and Cite

MDPI and ACS Style

Ishi, H. Explicit Formula of Koszul–Vinberg Characteristic Functions for a Wide Class of Regular Convex Cones. Entropy 2016, 18, 383. https://doi.org/10.3390/e18110383

AMA Style

Ishi H. Explicit Formula of Koszul–Vinberg Characteristic Functions for a Wide Class of Regular Convex Cones. Entropy. 2016; 18(11):383. https://doi.org/10.3390/e18110383

Chicago/Turabian Style

Ishi, Hideyuki. 2016. "Explicit Formula of Koszul–Vinberg Characteristic Functions for a Wide Class of Regular Convex Cones" Entropy 18, no. 11: 383. https://doi.org/10.3390/e18110383

APA Style

Ishi, H. (2016). Explicit Formula of Koszul–Vinberg Characteristic Functions for a Wide Class of Regular Convex Cones. Entropy, 18(11), 383. https://doi.org/10.3390/e18110383

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop