1. Introduction
Since we will be dealing with scalar, vector and matrix variables in the real and complex domains; both mathematical and random variables; and real scalar, vector and matrix parameters in the real and complex domains, we will need a multiplicity of symbols. In order to simplify matters, we will be using the following simplified notations: Real scalar mathematical as well as random variables will be denoted by the same lower-case letters, such as . Real mathematical as well as random vector/matrix variables, whether square or rectangular matrices, will be denoted by capital letters such as . , etc., and , etc., will be used to denote scalar and vector/matrix constants, respectively. Variables, whether scalar vector or matrix in the complex domain will be denoted with a tilde placed over them. No tilde will be used on constants. When symbols or Greek letters are used, the notation will be explained then and there. Let be a matrix, where the elements, s, are functionally independent (distinct) real scalar variables. Then, the wedge product of differentials will be denoted by . A matrix in the complex domain can be written as , where real matrices. Then, . will denote the transpose of X and will denote the conjugate transpose of . When a matrix , then is Hermitian. When is Hermitian and if real, then (symmetric) and ( skew symmetric). The determinant of a matrix X will be denoted as or and the absolute value of the determinant of a square matrix in the complex domain will be denoted as and|.
The need for introducing an extended zeta function will be illustrated by considering a simple model building situation [
1]. Let
x be a real scalar variable. A popular model in the real scalar variable case in different fields, such as mathematical sciences, physical sciences, biological and social sciences and related fields, is the Gaussian model. But a real scalar variable logistic model behaves like a Gaussian model but with a thicker tail and, hence, a logistic model is preferred in many industrial applications. A generalized logistic model was introduced in [
2] as an exponentiation of a basic type-2 beta model. A basic type-2 beta model is the following:
and
elsewhere. Let
. Then,
where
and then the function corresponding to
y, denoted by
, is the following:
where
denotes the real part of
. It may be observed that
and when
, then (
4) reduces to
and
is the basic logistic model, which is the most important model in industrial applications as a replacement for the standard Gaussian model.
Let us look into some extensions of the logistic model in (
3) to vector/matrix cases in the real and complex domains. Consider a
real vector
X,
, where
are distinct real scalar variables. Then,
. If
is in the complex domain, then
real. Consider the evaluation of the following integral in the real domain:
for
and
, where, for example,
and
is the Pochhammer symbol. The binomial expansion is valid here since
. Note that the basic zeta function
and the generalized zeta function
are the following [
3]:
Note that (
5) is connected to a generalized zeta series but it is not a zeta series. Hence, we need to have an extension of the generalized zeta series. In the above integral, if
is replaced by
and if there is another multiplicative factor
, then the integral becomes complicated, but it can be evaluated by using some results of Jacobians of matrix transformation. For the discussions in
Section 2 and
Section 3, we will need some Jacobians of matrix transformations.
The organization of this paper is as follows. In
Section 2, an extended zeta function is introduced and then several extended logistic models are considered. All these models end up in extended zeta functions. In
Section 3, some Bayesian procedures are considered, some of which will also end up in extended zeta functions.
2. An Extension of the Logistic Model
We will introduce several models here; each one eventually will end up with the extended zeta function introduced here.
Definition 1 (Extended zeta function).
Let and . An extended zeta function will be defined by the following series, denoted bywhere is the Pochhammer symbol and it is assumed that or . Now, we examine the evaluations of a number of integrals. Let
, where
X is a
real vector;
;
;
, with the expected value or mean value
, where
denotes the expected value and
A is a constant
real positive definite matrix;
; and
is a
parameter vector. Consider the evaluation of the following integral:
where
is a constant. Note that
since
is positive definite. In order to bring this to a sum of squares, one can make a transformation
, then
, where
is the positive definite square root of the positive definite matrix
. Then, we can convert
into
. This will enable us to evaluate the integral in (
8). Two results of Jacobians of matrix transformations are needed here. These will be listed here as lemmas without proofs. For the proofs and for other Jacobians, see [
4].
Lemma 1. Let the matrix be in the real domain, where the elements, s, are functionally independent real scalar variables, and let A be a and B be a nonsingular constant matrices. Then, Let the matrix be in the complex domain and, where A and B are and nonsingular constant matrices, respectively, in the real or complex domain, thenwhere denotes the absolute value of the determinant of . Lemma 2. Let the matrix X of rank p be in the real domain with distinct elements s. Let the matrix , which is positive definite. Then, going through a transformation involving a lower triangular matrix with positive diagonal elements and a semi-orthonormal matrix and after integrating out the differential element corresponding to the semi-orthonormal matrix, we will have the following connection between and . See the details from [4]:where, for example, is the real matrix-variate gamma function given bywhere means the trace of the square matrix , and is the real part of . Since is associated with a real matrix-variate gamma integral of (12), we call a real matrix-variate gamma function. This is also known by different names in the literature. When the matrix of rank p, with distinct elements, is in the complex domain and letting , which is and Hermitian positive definite, then, going through a transformation involving a lower triangular matrix with real and positive diagonal elements and a semi-unitary matrix, we can establish the following connection between and :where, for example, is the complex matrix-variate gamma function given by We call the complex matrix-variate gamma because it is associated with a matrix-variate gamma integral in the complex domain.
Now, we can evaluate the integral in (
8). Make the transformation
from Lemma 1. Let
from Lemma 2, observing that
here is a
matrix. Note that since
, one can expand
by using a binomial expansion as
. Then,
where
is the case for
and
of (7) for
.
Hence, when
is a statistical density of
X, then
will be the normalizing constant there and
for
and
is defined in (7). We will be constructing our models as non-negative integrable functions with normalizing constants so that they can also be used as statistical densities. Then, we have the first extended logistic model as the following, which is also a statistical density.
2.1. Extended Logistic Model 1 in the Real Domain
where
X is a
vector in the real domain,
and
. The results in the complex domain, corresponding to the results in the real domain, will be denoted with a letter
a in the subscript as well as in the section number of the equation numbers. For example, the result in the complex domain corresponding to
will be denoted as
and the equation number will be denoted by (
18) so as to avoid the occurrence of too many equation numbers. Let
be a
vector in the complex domain with distinct elements. Then, the
j-th element in
is
, where
real and
. Let
denote the conjugate transpose of
. Then, the model in the complex domain, corresponding to
, is available by replacing
with
and
. Let
. Note that, here,
is real. Hence, we can apply Lemma 2 in the real case for
real variables or
p complex variables. Both will lead to the same results. For
real variables, we have
, which is the same result when applying the complex version of the lemma. Then, we have the equation for the density in the complex domain, shown in
Section 2.2.
2.2. Extended Logistic Model 1 in the Complex Domain
for
and
,
is the absolute value of the determinant of
A and
is the extended zeta function defined in (7).
Note 1. In logistic model 1, we have taken in the real case, where X is with distinct real scalar variables as elements and is a real positive definite constant matrix. Particular cases herein are (1): ; (2): , where I is the identity matrix; and (3): , both in the real and complex domains. In the complex domain, μ is denoted as . One can consider a generalization of model 1 by taking in the exponent in the real case and in the exponent in the complex case. Then, in the real and complex cases, the exponential parts will be the following, respectively:and the first factor will be the same, that is, and , respectively, in the real and complex domains. In this case, the normalizing constants will be multiplied by δ, both in the real and complex cases. In the real case, is to be replaced by . In the complex case, is to be replaced by . These will be the changes and, hence, this generalization of model 1 will not be listed here. Note 2. In the PhD thesis written by [5], details of the applications of various types of statistical models, especially in the complex domain, are given. Most of the models are Gaussian- and Wishart-based. Ref. [5] is mainly concerned with the analysis of PolSAR (Polarimetric Synthetic Aperture Radar) data on single-look and multi-look radar return signals. It is found that non-Gaussian models give better representations in certain regions such as forests, sea surface, urban areas, etc.; see, for example, Refs. [6,7,8,9]. In the real scalar case, it is found that in many of the industrial applications a logistic model is better than a Gaussian model. A logistic model behaves like a standard Gaussian model but with thicker tail. This is the main reason why a logistic model became very popular in industrial applications. Hence, in modeling multi-look return data from radar and sonar and in other communication problems, logistic-based multivariate and matrix-variate models in the complex domain are likely to give better representations where, currently, Gaussian- and Wishart-based models are used. Hence, it is hoped that the results obtained in this paper and the various models presented here will be highly useful in physics, engineering, communication theory, statistics and related areas. Note that the h-th moment of and , for an arbitrary h, are available from the normalizing constants in extended logistic model 1 by replacing with and then taking the ratios of the respective normalizing constants.
The next model will involve a positive definite matrix variable. Let
be a
symmetric positive definite matrix in the real domain. Due to symmetry, there are only
distinct real scalar variable
s in
. When symmetric matrices are involved, the transformations will not be based on Lemmas 1 and 2. We need one more result on Jacobians involving a symmetric matrix. This will be listed here as Lemma 3. For the proof and other details, see [
4].
Lemma 3. Let the matrix X be real symmetric, , with distinct real scalar variable s, and let A be a nonsingular matrix. Then, Let the matrix in the complex domain be Hermitian. Let A be a nonsingular constant matrix in the real or complex domain. Then, For
, consider the following integral:
for
and
are defined in (7). Hence, we have the model that begins
Section 2.3.
2.3. Extended Logistic Model 2 in the Real Domain
For
, a
real positive definite matrix,
,
If the
matrix
is Hermitian positive definite, then
is to be replaced with
, the rest of the procedure in the real case goes through and the model will be as follows in
Section 2.4.
2.4. Extended Logistic Model 2 in the Complex Domain
For the
Hermitian positive definite matrix
,
and
, as defined in Lemma 2. The model in the complex domain is the following:
Note that one can generalize model 2 by replacing in the exponent with or with in the numerator and with in the denominator. Such generalized integrals can be evaluated, but the process will be very lengthy and, hence, such generalizations involving a determinant and power of trace in the exponent will not be considered here. Now, we will look into a model involving a rectangular matrix in the real and complex domains.
Consider a real
matrix
X of rank
p. Then, we know that
(real symmetric positive definite). From Lemma 2,
Now, we examine the evaluation of the following integral:
for
nd
is defined in (7). The binomial expansion is valid since
. The integral is evaluated by making use of the real matrix-variate gamma integral of Lemma 2. Therefore, we have the model shown in
Section 2.5.
2.5. Extended Logistic Model 3 in the Real Domain
For
a
real matrix of rank
and then the
matrix
,
for
and
.
is defined in (7).
Now, we look at the complex domain. Let the
matrix
in the complex domain be of rank
. Then, the following are the replacement to obtain model 3 corresponding to model 3 in the real domain: Replace
with
. In the transformation
, we have
, where
is defined in Lemma 2. The rest of the procedure is parallel to that in the real case. Hence, model 3 in the complex domain is shown in
Section 2.6.
2.6. Extended Logistic Model 3 in the Complex Domain
for
and
.
Note 3. We can generalize model 3 in the real case by replacing with , where is a constant matrix; and are and constant positive definite matrices, respectively; and and are, respectively, the positive definite square roots of and B. Then, the normalizing constant is to be multiplied by in the real case and in the complex case the normalizing constant is multiplied by , that is, the absolute values of the determinants raised to q and p, respectively. in the real case is replaced by , and in the complex domain. Therefore, no separate listing will be there for this generalization. Another generalization of model 3 in the real domain is to replace with , that is, an arbitrary power is taken for the trace in the exponent, either the same for both the numerator and denominator exponents or different and and the corresponding replacements in the complex domain. Even with the same , the evaluation of the integral will become very lengthy and when the δ’s are different; it will be more difficult also. The same procedures will also work for . These types of generalizations will not be considered here.
Now, we consider a case where the multiplicative factor determinant is replaced by a trace and the exponential trace is given an arbitrary power. Let the
real matrix
X be of rank
. Consider
. Since
is the sum of squares of
real scalar variables, one can apply Lemma 2 by taking a row vector of
real scalar variables. Let
. Then, let us evaluate the following integral:
for and ; is defined in (7). From this integral, we can construct the model in Section 2.7. 2.7. Extended Logistic Model 4 in the Real Domain
for
and
. The following are the replacements to obtain model 4 in the complex domain. Replace
with
. Then,
is the sum of squares of the absolute values of
complex variables or the sum of squares of
real scalar variables. Therefore,
in the real case will go to
. Then, using steps parallel to those in the real case, one has model 4 in the complex domain, shown in
Section 2.8.
2.8. Extended Logistic Model 4 in the Complex Domain
For
and
.
is given in (7),
Note 4. One can generalize model 4 as follows: Take the δ in the numerator exponent as and that in the denominator exponent as . In this case, the evaluation of the integral becomes very complicated. One can introduce a factor containing a determinant such as as an additional factor into model 4 in the real case and in the complex case. Then, the evaluation of the integral will become very lengthy; see [10] and the references therein. Hence, these cases will not be discussed here. Another generalization possible is to replace with , where the constant matrix , the constant matrix , and is the positive definite square root of in the real case and the corresponding replacements in the complex domain with . This generalization brings only minor changes. In the real case, the normalizing constant is multiplied by , and in the complex case, the normalizing constant is multiplied by . The particular cases here will be (1): ; (2) ; (3): and . This generalized model will not be listed as a different model here. 3. Some Bayesian Models
Our aim here will be to look into a few Bayesian procedures which will end up in the extended zeta function. Bayesian models are very important in many applied areas. According to [
1] “Quantum Mechanics is just the Bayesian theory generalised to the complex Hilbert space”. The models that we construct in this section will hold in multivariate and matrix-variate complex domains, and will hold in complex Hilbert spaces as well; hence, the models are expected to be highly useful in many areas, especially in Quantum Physics.
Consider the standard Gaussian function
in the real case and
in the complex case, where
Z is
,
and
is
,
, and the functions are the following, denoted by
in the real case and
in the complex case, where
We will include the normalizing constants in all our models so that they will also be statistical densities. Consider a constant
matrix
H of rank
p. Then, in the real case,
and
X is called a scale mixture. Then,
(real positive definite). Let
. Then, the function or density of
, at a given
A, denoted by
, is the following:
for
and
In the complex case,
and the function or density of
in the conditional space, denoted by
, is the following:
for
and
. We will consider several cases of prior functions or prior distributions for the parameter matrix
A.
3.1. Bayesian Models
In the real case, let
be the prior function or prior density of
A, where
where
is defined in Lemma 2. The corresponding function in the complex domain is denoted again by
, where
A and
C in the complex case are Hermitian positive definite.
where
is defined in Lemma 2. Note that
C in (
30) may be complex or real. It is assumed that the
constant positive definite matrices
C in the real and complex cases are known. If
C has its own prior function or prior density, then we have one more step in the Bayesian procedure. The unconditional function or unconditional density of
X and
will be called the Bayesian models in this paper, and they will be denoted by
and
, respectively. The Bayesian procedure involves one more step of computing the posterior function or posterior density, denoted by
in the real case and
in the complex case, and then studying the properties of
A in the real case and
in the complex case in the conditional space of
and
, respectively. In the above real case,
is as follows in
Section 3.2.
3.2. Bayesian Model 1 in the Real Case
where
are the eigenvalues of
and
.
3.3. Bayesian Model 1 in the Complex Domain
Let the unconditional function or unconditional density of
in the complex domain be denoted by
. Then, going through parallel steps to those in the real case, we have
as the following:
for
and
, where
and
are the eigenvalues of
. We have used the same notations for the eigenvalues in the complex case as in the real case, but in the complex case the eigenvalues will be real because the matrix here is Hermitian positive definite, but the eigenvalues need not be the same as those in the real case.
The next step in the Bayesian procedure will be illustrated by computing the expected value of
in the conditional space of
A, given
X, in the real case here. Then, in the real case, the conditional density of
A, given
X, denoted by
and the Bayes’ estimate of
, given
X, is
in
and it is the following, where
denotes the expected value of
:
For the next model also, we will start with the conditional function as a Gaussian one as in (
27) and (
28). For the prior function or prior density, we will consider a matrix-variate logistic-type function. Let the prior function or prior density for
A in the real case be the following, again denoted by
:
for
and
, where the normalizing constant
c is already evaluated in
Section 2. The corresponding prior function or prior density in the complex domain, denoted by
, will be the following:
for
and
Then, the unconditional function or density of
X in the real case, again denoted by
, is as follows in
Section 3.4.
3.4. Bayesian Model 2 in the Real Domain
The corresponding
in the complex domain is the following:
By using the properties of the determinant of a partitioned matrix, we have the following equivalent relationships in the real case:
Then,
of (
33), that is, in the real case, is the following:
for
and
X is
real. By using steps parallel to those in the real case, we can obtain the corresponding unconditional function or density in the complex domain as follows in
Section 3.5.
3.5. Bayesian Model 2 in the Complex Domain
for
and
is
.
For the next model, we will also start with a real multivariate Gaussian function or density with the parameter matrix being the covariance matrix therein, denoted by a
real positive definite matrix
. The Gaussian function in the real domain is the following, again denoted by
:
for
and
. In the complex domain, the corresponding function or density is the following:
for
and
. The two models that we considered so far had
having a prior function or density. But, when
B has a prior function or density, then the problem becomes complicated. Then, the integral to be evaluated will belong to matrix-variate Bessel type or Krätzel type. Let the prior function or density for
B be an extended logistic function or density. We will denote the prior function in the real and complex domains by
and
, respectively, where in the real case
and in the complex case
. In the real case, let
for
and
. The corresponding prior function or density in the complex domain can be obtained by replacing
with
,
with
p,
by
and the condition
with
. Hence, it will not be listed here separately in order to save space. Then, the unconditional density of
X in the real case, again denoted by
, is as follows in
Section 3.6.
3.6. Bayesian Model 3 in the Real Domain
for
and
X is
real. Direct evaluation of the integral in (
38) is difficult because it is a real matrix-variate Bessel type or Krätzel type integral involving
B and
in the exponent. Ref. [
11] evaluated the general form of this integral for the real scalar variable case, in connection with the generalized reaction-rate probability integral in nuclear reaction rate theory. We can consider the Fourier transform of
in (
38) by taking an expected value of a conditional expectation. Let
T be a
parameter vector. The Fourier transform of
, denoted by
, is
. This we can evaluate by taking a conditional expectation at a given
B and then taking an expectation over
B. That is,
But from real Gaussian density, we have the conditional expectation of
as
, and in the complex case it is
, where
is the conjugate transpose of
T; see, for example, [
12]. Observe that in the real case
is now replaced by
. In the complex case,
is now replaced by
. Now, take the integral over
B to complete the integration in the real case:
for
. By using the properties of the determinant of a partitioned matrix, we have the following:
The corresponding quantity in the complex domain is
. The unconditional function or density of
X in the real case is now available from the inverse Fourier transform of
. In order to evaluate it, we will use the technique used in [
12]. We can replace
with an equivalent integral, namely,
where
y is a real positive scalar. The Fourier inverse of
is
. Hence, the Fourier inverse coming from (
39) is the following:
For evaluating the integral in (
41), we will use the following lemma:
Lemma 4. For where is the Bessel hypergeometric series. Proof. Let us take the Mellin transform of
with Mellin parameter
s. Then,
From the inverse Mellin transform,
For
, the poles of the integrand are simple and then the integral can be evaluated as the sum of the residues at the poles of
and
. Poles coming from
give the sum of residues, denoted by
, as the following:
where
is a Bessel hypergeometric series, the properties used are
and
whenever
and
is defined. The sum of the residues coming from
, denoted by
, is the following:
Note that proves the lemma.
Then, the integral part in (
41) is the following:
for
and
. In the complex case, the corresponding integral will be the following:
for
and
. □
Note 5. The integral in Lemma 4 can also be evaluated in terms of a modified Bessel function of the second kind:for and , where is the modified Bessel function of the second kind of order ν; see, for example, page 368 of [13]. Hence, from (41) and (42), we have the unconditional density of X in the real case as the following:for and . In the complex domain, the unconditional density is available from (44) by replacing with and with . That is,for and . 3.7. Bayesian Model 4 in the Real Case
Now, we consider a rectangular matrix. Let the
real matrix
X of rank
p have the following density:
for
and
, where
A is a
and
B is a
constant positive definite matrix. The normalizing constant
c can be evaluated by making the transformations
and
, and then integrate out by using a real matrix-variate gamma integral to obtain
The density in the complex case, denoted by
, is available from the real case structure by replacing
with
and
with
, and then the normalizing constant in the complex case, denoted by
, is the following:
If matrix
A has a prior function or density, then
in the real case is
. Let the prior density of
A in the real case be again denoted by
where,
where the constant positive definite matrix
is assumed to be known. Then, the unconditional density of
X in the real case, denoted by
, is the following:
One can have the following simplification by using the properties of the partitioned determinant:
where
are the eigenvalues of the positive definite matrix
From (
50) and (
51), the unconditional density can be written in terms of an extended zeta function. Observe that
is cancelled in the normalizing constant and the balance of the exponent is
. Then,
, the unconditional density in the real domain, is the following:
for
and
, where the
’s are defined in (
51).
In the complex case, the prior density of
, corresponding to the real case
, denoted by
, is the following:
The unconditional density of
, denoted by
, is the following:
for
and
, where the
s appearing here are the real positive eigenvalues of the Hermitian positive definite matrix
, and we have used the same notation,
, as in the real case, and the
s in the complex case need not be equal to those in the real case.
4. Concluding Remarks
Instead of
A in (
52) and (
54) having a prior distribution, if
B has a prior distribution or if both
A and
B have compatible prior distributions, then one can also evaluate the unconditional density, but the procedure will be a little more lengthy. Note that we have only mainly considered two cases of the unconditional function or density for
X, namely, one case where the
vector
X has a Gaussian density with a null mean value and
X has a
rectangular matrix-variate gamma distribution. Then, we considered different prior distributions for a parameter matrix, mainly in the class of logistic-type distributions. We did not consider the situations where a determinant appears as a multiplicative factor when the trace in the exponent has an arbitrary power, or a determinant and a trace appear as multiplicative factors, when the trace in the exponent has an arbitrary power, because these cases will involve lengthy steps and too much space would have been taken for the discussion. We could have also looked into Mathai’s pathway-type models for the conditional and prior functions or densities, which could be considered as extensions of multivariate-, matrix-variate Gaussian-, matrix-variate gamma- or Wishart-type models. Also, in order to save space, we did not study various properties of the models introduced here. In order to limit the size of the manuscript, we stop the discussion here.
For further studies, apart from the items mentioned above, one can explore the pathway idea by combining matrix-variate type-1 beta, type-2 beta and gamma densities, observing that particular cases of the matrix-variate gamma will produce matrix-variate Gaussian as well as Wishart densities. One can study various properties of the models introduced here. Such studies are likely to bring in more parameters from the set
s and
s. One can explore the applications of the models here, especially the ones in the complex domain, in quantum physics, communication theory and related areas in light of the evidence in [
1,
5].
This research received no external funding. The author declares no conflict of interest. The author acknowledges with thanks the comments and suggestions made by the five or six reviewers of the original version of this paper, which enabled him to improve the presentation of this paper.