Next Article in Journal
A Tool for the Analysis and Characterization of School Mathematical Models
Next Article in Special Issue
Extreme Value Analysis for Mixture Models with Heavy-Tailed Impurity
Previous Article in Journal
Graph Theory for Primary School Students with High Skills in Mathematics
Previous Article in Special Issue
A Variant of the Necessary Condition for the Absolute Continuity of Symmetric Multivariate Mixture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Notion of Reproducibility and Its Full Implementation to Natural Exponential Families

Faculty of Industrial Engineering and Technology Management, Holon Institute of Technology, Holon 5810201, Israel
Mathematics 2021, 9(13), 1568; https://doi.org/10.3390/math9131568
Submission received: 26 May 2021 / Revised: 28 June 2021 / Accepted: 29 June 2021 / Published: 3 July 2021

Abstract

:
Let F = F θ : θ Θ R be a family of probability distributions indexed by a parameter θ and let X 1 , , X n be i.i.d. r.v.’s with L ( X 1 ) = F θ F . Then, F is said to be reproducible if for all θ Θ and n N , there exists a sequence ( α n ) n 1 and a mapping g n : Θ Θ , θ g n ( θ ) such that L ( α n i = 1 n X i ) = F g n ( θ ) F . In this paper, we prove that a natural exponential family F is reproducible iff it possesses a variance function which is a power function of its mean. Such a result generalizes that of Bar-Lev and Enis (1986, The Annals of Statistics) who proved a similar but partial statement under the assumption that F is steep as and under rather restricted constraints on the forms of α n and g n ( θ ) . We show that such restrictions are not required. In addition, we examine various aspects of reproducibility, both theoretically and practically, and discuss the relationship between reproducibility, convolution and infinite divisibility. We suggest new avenues for characterizing other classes of families of distributions with respect to their reproducibility and convolution properties.

1. Introduction

The notion of a distribution function, F, which is reproductive with respect to a parameter θ Θ , a nonempty set of R , was first introduced by Wilks [1] as follows: Let X 1 and X 2 be independent r.v.s with distributions F ( . ; θ 1 ) and F ( . ; θ 2 ) , respectively. If for any θ 1 , θ 2 Θ , the distribution of Y = X 1 + X 2 is F ( . ; θ 1 + θ 2 ) , then the distribution of F ( . ; θ ) is said to be reproductive with respect to θ . In this respect, Wilks mentioned the Poisson, chi-square and normal distributions.
After Wilks introduced this definition, there were early attempts to characterize such families, c.f.,Bolger and Harkness [2], Baringhaus, Davies, and Plachky [3] and Bar-Lev and Enis [4], but all ended in characterizing the Poisson distribution. A broader definition of Wilk’s concept was given by Bar-Lev and Enis [5], and then applied to one-parameter natural exponential families (NEFs) in-depth analysis. (Well known preliminaries on NEFs will be presented in Section 3). The definition of reproducibility given by Bar-Lev and Enis [5] is the following.
Definition 1.
Let F = F θ : θ Θ R be a family of distributions indexed by a parameter θ Θ , where Θ has a nonempty interior. Let X 1 , , X n be i.i.d. r.v.s with L ( X 1 ) = F θ F (where L ( X ) stands for the law of X 1 ). Then, F is said to be reproducible if for all θ Θ and n N , there exists a sequence ( α n ) n 1 (called the stabilizing constants), which is either identically equal to 1 or it is a one-to-one mapping α n : N R and there exists a mapping g n : Θ Θ , θ g n ( θ ) such that
L ( α n i = 1 n X i ) = F g n ( θ ) F .
Or equivalently, in terms of the characteristic functions, the reproducibility property is equivalent to requiring that
f ( t , g n ( θ ) ) = f n ( α n t , θ ) , n N , t R and θ Θ ,
where f is the characteristic function of X 1 .
Note that for α n 1 , the reproducibility property means that the convolution belongs to the same family, while for α n = 1 / n , it means that the distribution of the sample mean belongs to the same family as each of the X i ’s.
Bar-Lev and Enis [5] implemented their Definition 1 to the class of natural exponential families (NEFs) by imposing the following restrictions:
( i )   g n   is   an   onto   mapping , ( ii )   α n   has   the   form   α n = n β , β R , ( iii )   The   NEF   F   is   steep ,
where steepness of F means that the mean parameter space of F coincides with the interior of its common convex support. This definition of steepness will be further elaborated in the sequel. Their implementation of the notion of reproducibility of NEFs led to the characterization of the class of steep NEFs having power variance functions (VFs).
The analysis of the notion of reproducibility from a different angle than that formulated in Definition 1 and Equation (1) was given in Bar-Lev and Cassalis [6,7]. They defined an NEF F to be reproducible in the broad sense, as follows. Let F be an NEF which is associated with a Laplace transform L. If there exists a triple ( α , β , λ ) ( R , R , R + ) such that f α , β ( F ) = F λ , where f α , β ( F ) is the image of F under the affine transformation f α , β : x α x + β and F λ is the NEF associated with the Laplace transform L λ , then F is said to be reproducible in the broad sense. They obtained a complete classification of the class of reproducible NEF’s in the broad sense and showed that this class is composed of NEFs having power VFs, exponential VF, as well as NEF’s constituting discrete versions of the latter NEF’s. In a completely different aspect and not at all related to the notion of reproducibility, the class of NEFs having power VF was introduced by Tweedie [8] and won publicity only after being described in Jorgensen [9].
In this paper, we adhere to the original definition of reproducibility presented in Definition 1. In Section 2, we will present a number of issues related to reproducibility for general classes of parametric families of distributions and discuss some practical aspects. Section 3 is devoted to reproducibility in NEFs. Indeed, we will prove that the three requirements in (2) assumed by Bar-Lev and Enis [5] are not needed as one can obtain their characterization results and more without assuming that (2) holds. Our method of proof enable to delineate all NEFs having power VFs including those missed by Bar-Lev and Enis [5] due to their steepness restriction. Naturally, the latter assumption did not allow them to characterize the subclass of non-steep NEFs having VFs with negative power. Such a subclass of non-steep NEFs is huge as it contains all NEFs generated by stable distributions with a stable index in ( 1 , 2 ) . In the same context, one should note that Tweedie [8] was wrong when claimed that such a subclass does not exist (see Jorgensen [9] and Bar-Lev [10] for more details). In fact, we will prove in Section 3 that the three assumptions in (2) follow trivially from the proof presented in this paper.

2. Some General Noteworthy Comments on the Notion of Reproducibility

A number of noteworthy comments on the notion of reproducibility will now follow:
  • Restriction on the support of F
It would appear pretentious to think that one can characterize the property of reproducibility for all parametric families. It, therefore, seems necessary to limit ourselves in a number of aspects. One of them, for example, is to require that all members of F possess a common support. The second is inherently obvious as to require that the family F = F θ : θ Θ R is identifiable, i.e., if θ 1 , θ 2 Θ with θ 1 θ 2 then F θ 1 F θ 2 on at least one Borel set. Accordingly, in order to aim efficiently at characterizing reproducible parametric families, one should impose the following assumption.
Assumption A1.
(i)
For all θ Θ , F θ possesses a common support S . Let ( c , d ) denote the interior of the convex-hull of S, where c < d .
(ii)
F is identifiable.
Some immediate conclusions which follow from Assumption A1 are presented in the following lemma.
Lemma 1.
Suppose that F = F θ : θ Θ R is reproducible in the sense of Definition 1 and that it satisfies Assumption A1. Then,
(i)
c and d do not depend on θ ,
(ii)
for any fixed n N , the mapping g n is one-to-one from Θ into itself,
(iii)
if α n 1 / n , then either c = and d = 0 or ; or c = 0 and d = .
(iv)
α n > 0 , for all n N .
Proof .
(i)
Simple, as F depends on θ only and is assumed to have a common support S .
(ii)
If θ 1 θ 2 , then g n ( θ 1 ) g n ( θ 2 ) , otherwise a contradiction to the identifiably of F will follow.
(iii)
If α n 1 / n , then relation (1) implies that
n α n ( c , d ) = ( n α n c , n α n d ) = ( c , d )
and thus, the desired result is obtained.
(iv)
This follows from (3).
The situation where α n = 1 / n , which is excluded in Lemma 1, is satisfied by several reproducible families that are not NEFs. The most famous example is the Cauchy distribution (with an unknown scale parameter θ and a fixed location parameter δ ) having a density
f θ ( x ) = 1 π θ 1 + ( x δ θ ) 2 , < x < , θ R .
For this example, L ( X Θ n ) = L ( X 1 ) , c = , d = and g n ( θ ) θ for all n N .
2.
The restriction thatcandddo not depend on θ and that α n is a one-to-one mapping
We restricted the analysis to a situation where both c and d do not depend on θ , as otherwise, not placing such a restriction, would have made the task of the reproducibility characterization almost impossible to execute. Below are two examples (out of many) of reproducible families of distributions for which either c or d depend on θ .
Example 1.
Let F =   B i n ( θ , 1 / 2 ) ; θ N , α n 1 and g n ( θ ) = n θ , then α n Y n B i n ( g n ( θ ) , 1 / 2 ) F . Note that here d depends on θ and α n is not one-to-one.
Example 2.
Let U 1 , U n be n i.i.d. U ( 0 , θ ) r.v.s. For k 1 < θ k , denote by f θ the density of i = 1 k U k . Then, this density has the reproducibility property. Note that S depends on θ and α n is not one-to-one.
3.
Reproducibility and infinite divisibility
At first intuitive glance, it was conjectured that reproducibility of F will entail its infinite divisibility. Recall that a probability distribution is called infinitely divisible if it can be expressed as the probability distribution of the sum of an arbitrary number of independent and identically distributed random variables. Accordingly, an infinite divisibility of a probability distribution entails that the convex-hull of its support cannot be bounded (c.f., Chapter 17 of Feller [11]).
The two examples above show that such a conjecture is false as the corresponding families have bounded supports and thus are not infinitely divisible. However, then the question arose as to whether reproducible families with unbounded support are infinitely divisible. The following example (provided to the author via a personal communication by Gérard Letac - Institut de Mathématiques de Toulouse, Université Paul Sabatier, France) shows that the related conjecture is also false.
Example 3.
This is an example in which S is unbounded and does not depend on θ . Let a > 0 , n be a positive integer, Z N ( 0 , 1 ) , B n B ( n , 1 / 2 ) , where Z and B n are independent. Define
Y = n a Z + B n ,
then, the distribution P a , n of Y depends on two parameters n and a. However, if one defines the parameter
θ = 1 1 + a + n , θ > 1 ,
then, P θ = P a , n is parameterized by θ and has the reproducibility property, although it is not infinitely divisible.
4.
Reproducibility and convolutions
Let F be reproducible in the sense of (1), i.e.,
L ( α n i = 1 n X i ) = F g n ( θ ) F ,
from which it follows that the convolution of the X i s, i = 1 , , n , is
P i = 1 n X i x = F g n ( θ ) ( x α n ) ,
where F g n ( θ ) is the distribution of X 1 with parameter g n ( θ ) Θ , i.e., the distribution of a sum of i.i.d. r.v.s has the same distribution type as any of its components.
Convolutions of n i.i.d. random variables play an important role both in statistical inference, probability, stochastic processes, risk theory and insurance. Some immediate examples are (i) For NEFs in statistical modeling, Y n = i = 1 n X i is the minimal sufficient statistic for θ (see below), the distribution of which is needed for all statistical inference aspects on θ ; (ii) for an M / G / 1 queuing system with “first come, first served” queue discipline and service time distribution F, the n-fold convolution of F with itself is needed for computing the distribution of the length of the busy period; (iii) in insurance risk, the aggregated claim is a random sum of i.i.d. r.v.s with common distribution F and thus the n-fold convolution of F is needed for various computations (e.g., Shushi and Yaob [12], Bahnemann [13]). In general, however, the derivation of a convolution is usually rather complex and cumbersome. Thus, if F belongs to a reproducible family, then the convolution has the same distribution type as any of its components (as in (4)) and thus no further complex computations are needed at all. Indeed, the reproducibility property has been employed in Bar-Lev and Ridder [14] for computing the insurance risk of aggregated claim data and then implemented to some real car insurance claim data.
5.
An extension of the reproducibility notion to the multi-parameter case
Naturally, the reproducibility property can be easily extended in the sense of Definition 1 to the multi-parameter case. Three examples of reproducible two-parameter ( θ 1 , θ 2 ) NEFs appear in Bar-Lev and Reiser [15]. These are the two-parameter gamma, inverse Gaussian and normal NEFs. For these three NEFs, the distribution of the sample mean X ¯ n is of the same type as that of X 1 but with parameters ( n θ 1 , n θ 2 ) (see Bar-Lev and Reiser [15] for the definition of the parameters ( θ 1 , θ 2 ) for each of these NEFs).
Obviously, from both theoretical and practical (convolutions) aspects it would be interesting to characterize families of multi-parameter distributions by the reproducibility property.

3. NEFs—Preliminaries and Characterization by the Reproducibility Property

In this section, we first present some required preliminaries on NEFs and their associated VFs. These preliminaries are taken from the fundamental paper of Letac and Mora [16]. We will then introduce the main result and characterize reproducible NEFs in the sense of Definition 1. In particular, we will show that the sequence α n has the form n β for some β R , and that an NEF is reproducible iff it possesses a power VF. Afterwards, we present a table including the class reproducible NEFs with power VFs along with their associated α n , g n ( θ ) , S , and Θ .

3.1. Some Preliminaries on NEFs

Let ν be a positive Radon measure on R with convex support C ν . Consider the set
D ν θ R : L ν ( θ ) R exp ( θ x ) ν ( d x ) < ,
and assume that Θ ν int D ν is nonempty. Then, the NEF F ( ν ) generated by ν is defined by the set of probability distributions
F ( ν ) F θ , ν ( d x ) = exp θ x k ν ( θ ) ν ( d x ) : θ Θ ν ,
where k ν ( θ ) log L ν ( θ ) is the cumulant transform of ν and k ν is strictly convex and real analytic on Θ ν . Moreover, k ν ( θ ) and k ν ( θ ) , θ Θ ν , are the respective mean and variance corresponding to F ( θ , ν ) , and the open interval M ν k ν ( Θ ν ) is called the mean domain of F ( ν ) .
An important observation is that measure ν is not unique for F ( ν ) . Let M be the set of Radon measures ν on R for which L ν ( θ ) < on domain Θ ν . Consider two measures ν , ν M , and suppose that ν is an exponential shift of ν ; i.e., ν ( d x ) = e a + b x ν ( d x ) for some real a , b . Then, a simple calculation shows that F ( ν ) = F ( ν ) . This also holds in reverse, if F ( ν ) = F ( ν ) for two measures ν , ν M , then one is an exponential shift of the other and, obviously, they are equivalent in the sense that each is absolutely continuous with respect to the other. A VF ( V , M ) of NEF F determines the NEF uniquely within the class of NEFs in the following sense: If F 1 and F 2 are NEFs with respective VFs ( V 1 , M 1 ) and ( V 2 , M 2 ) such that V 1 = V 2 on J = V 1 V 2 ϕ , then F 1 = F 2 , c.f., Mora [17,18]. This would imply the following interesting observation: if ( V , M ) is a VF of an NEF F , then
The   mean   parameter   space   M   is   the   largest   open   interval      on   which   V   is   positive   real   analytic .
Consequently, we may denote the NEF by F = F ( ν ) and the mean domain M = M ν to stress the fact that these do not depend on ν .
Since the function k ν : Θ ν M is one-to-one, its inverse function k ν 1 : M Θ ν is well defined. Since k ν is continuous and increasing, then M, the image of k ν ( Θ ) is an interval, say ( a , b ) ,on the real line. Now, when we compute the variance V ν ( θ ) k ν ( θ ) of the distribution F ( θ , ν ) as a function of the mean m M , i.e.,
V ν ( m ) = k ν ( k ν ) 1 ( m ) ,
then, it also does not depend on ν , and we can suppress the dependence on ν and write V ( m ) instead of V ν ( m ) . Furthermore, when no confusion is made, we shall suppress the dependence on ν of all other notations introduced above and write C , D , Θ ,..., for C ν , D ν , Θ ν
Finally, F is called steep iff M = int C (or equivalently, if k ν is essentially a smooth convex function on Θ —for a definition, see Barndorff-Nielsen [19]).
Now, F is said to have a power, VF, iff it has the form
V ( m ) = δ m γ , δ > 0 , γ R .
We are now ready to present the reproducibility characterization of NEFs in the following proposition. This proposition states that F is reproducible (in the sense of Definition 1) iff it possesses a power VF. Accordingly, for brevity, the details of a classification of all NEFs having a power VFs is given in Section 3.2 and not in the proof of the proposition.
In addition, in a manner similar to Part (iii) of Lemma 1, it follows that
n α n M = ( n α n a , n α n b ) = ( a , b )
and thus M is either ( 0 , ) or ( , ) (for more details, see Bar-Lev and Enis [5]).
Proposition 1.
Let F be an NEF generated by a basis ν and associated by a VF ( V , M ) . Then, F is reproducible in the sense of Definition 1 iff V is a power VF of the form (9). The values of γ for which (9) is a VF of an NEF along with the corresponding α n , g n ( θ ) , k p ( θ ) , Θ and M are presented in Table 1 below.
Proof. 
If F is an NEF with power VF of the form (9), then trivially F is reproducible (see also Table 1 below). The reverse implication is made as follows. Assume now that F is a reproducible NEF. We shall henceforth exclude the case where α n 1 (i.e., convolution) as this case was treated by Bar-Lev and Enis [4]—without assuming (2)—and led to the characterization of the Poisson NEF. We shall also exclude the case α n = 1 / n as by Lemma 2.1 in Bar-Lev and Enis [5], no one-parameter reproducible NEF exists for this case. Accordingly, we shall assume throughout the sequel that
α n 1 and α n 1 / n .
Let ν 0 M be a specific chosen basis of F . Then, by (1), the probabilities of α n i = 1 n X i are
P ( ν 0 ( d x ) ; θ ) = exp ( g n ( θ ) x k g n ( θ ) ν 0 ( d x ) , θ Θ .
Let ν 0 n be the n-th fold convolution of the chosen basis ν 0 of F . Then, the probabilities of i = 1 n X i are given by
exp ( θ x n k θ ) ν 0 n ( d x ) , θ Θ ,
and thus from (12), it follows that the probabilities of α n i = 1 n X i are given by
exp ( θ x α n n k θ ) 1 α n ν 0 n ( d x ) , θ Θ .
Consequently, by comparing (11) and (13) we obtain that
l ν 0 = 1 α n ν 0 n ,   almost   everywhere   ν 0 , g n ( θ ) = θ α n , θ Θ , n N ,   and k g n ( θ ) = n k θ ) , θ Θ , n N .
Hence, by comparing (11) and (13), we obtain
g n ( θ ) = θ α n
and
n k ( θ ) = k g n ( θ ) .
Substitute (15) into (16), we have
n k ( θ ) = k θ α n .
Denote α n α ( n ) , n N , and let y = y ( n ) 1 / α ( n ) . Since by assumption α ( n ) is a one-to-one mapping from N into R (see part (ii) of Lemma 1), its inverse
n = α 1 ( 1 y ) h ( y )
is well defined and thus by using (18) in (17), we obtain
k ( θ y ) = k ( θ ) h ( y ) , θ Θ ,
By Aczél ([20], Theorem 4, pp. 144–145), the general solution of the functional equation
f ( x y ) = g ( x ) h ( y )
with positive x and y and f continuous at a point is
f ( t ) = a b t c , g ( t ) = a t c , h ( t ) = b t c
(supplemented with a trivial solution, which is irrelevant in our situation). In order to apply the solution (21) of (20) to (19), we need to show that the premises of Theorem 4 in Aczél [20] are met. Indeed, k is real analytic on Θ and thus is continuous there. In addition, by Lemma 1, y = 1 / α ( n ) is positive. Now, if Θ contains an open interval ( a 1 , b 1 ) R + , then we will confine any further analysis of the problem to this interval. However, if otherwise Θ R , we define x = θ and r ( x ) = k ( x ) and apply the solution (21) of (20), to r. Therefore, for simplicity and without any loss of generality, we may assume that Θ contains an open interval ( a 1 , b 1 ) R + and thus we will continue with k instead of r. Hence, Theorem 4 of Aczél [20] can be applied for k = f = g . In which case, we obtain
b = 1 ,
k ( t ) = a t c , t ( a 1 , b 1 ) R +
and
n = h ( y ) = y c = α ( n ) c , n N ,
which implies that
α ( n ) = n 1 c and g n ( θ ) = θ n 1 c , where by ( 10 ) c 1 .
We now exclude the case c = 2 as it will be treated separately in the sequel and thus we assume that c 2 . Let m = k ( t ) = a c t c 1 and V ( t ) = k ( t ) = a c ( c 1 ) t c 2 be the mean and variance of F defined on ( a 1 , b 1 ) and let M 1 = k ( ( a 1 , b 1 ) ) be the image of ( a 1 , b 1 ) under k . Since V ( t ) > 0 on t ( a 1 , b 1 ) , k is strictly increasing on ( a 1 , b 1 ) and thus its inverse k ν 1 : M 1 ( a 1 , b 1 ) is well defined and is given by k ν ( m ) 1 = m a c 1 / ( c 1 ) , a c > 0 , m M 1 . Hence, V expressed in terms of m M 1 has the form
V ( m ) = a c ( c 1 ) a c c 1 c 2 m c 1 c 2 , a c > 0 , c 1 > 0 , m M 1 .
Let
δ = a c ( c 1 ) a c c 1 c 2 and γ = c 1 c 2 , c 1 , 2
then, V has the form
V ( m ) = δ m γ , δ > 0 , γ < 0 or γ > 0 , γ 1 , m M 1 = ( a 1 , b 1 ) R + .
As by (7), the mean parameter space M of V is the largest open interval on which V is positive real analytic; it follows that the mean parameter space corresponding to V in (24) is M = R + . Accordingly, if F is reproducible in the sense of Definition 1, then its VF is necessarily of the form
( V , M ) = ( δ m γ , R + ) , δ > 0 , γ < 0 or γ > 0 , γ 0 , 1 .
In order to conclude the proof we need to consider the two remaining cases: γ = 0 ( c = 2 ) and γ = 1 ( c = 1 ) . Note that for γ = 0 , the VF is, identically, a positive constant δ > 0 and thus the corresponding mean parameter space is M = R . Indeed, this case with γ = 0 was treated in Theorem 2.1 of Bar-Lev and Enis [5]. It was shown there that F is reproducible with stabilizing constants α n = n 1 / 2 iff F is the family of normal distributions with constant variance δ (i.e., F is an NEF having a power VF with power parameter γ = 0 ) . The case γ = 1 was analyzed in Bar-Lev and Enis [4]. It was shown that F is reproducible with stabilizing constants α n 1 iff F is the family of Poisson distributions (i.e., F is an NEF with power VF and power parameter γ = 1 ) .
Thus, we have proven that an NEF F is reproducible in the sense of Definition 1 iff it possesses a power VF. The appropriate forms of the stabilizing constants α n and of g n ( θ ) are presented in Table 1 below and this concludes the proof. Furthermore, for less common used distributions, the corresponding probability densities will be specified. □

3.2. A Classification of NEFs with Power VFs and Their Associated α n , g n ( θ ) , k γ ( θ ) , Θ and M

The permissible values of γ for which (9) is a VF of an NEF along with the corresponding α n , g n ( θ ) , k p ( θ ) , Θ and M are taken from Jorgensen [9], Bar-Lev and Enis [5] and Bar-Lev and Cassalis [7]. Only if γ ( 0 , 1 ) , then no NEF exists. For all other values of γ , we have the following:
  • For γ ( , 0 ) , the NEF F is generated by an extreme stable distribution with stable index 1 < τ < 2 , where τ = ( γ 2 ) / ( γ 1 ) , in which case C = R and M = R + , i.e., F is not steep. Here, α n = n 1 γ γ 2 and g n ( θ ) = θ n γ 1 γ 2 . (As already mentioned above, Bar-Lev and Enis [5] showed that non-steep NEFs exist if γ < 0 , whereas Tweedie [8] claimed that such NEFs do not exist by utilizing an incorrect claim). Here, the associated absolutely continuous probability density is quite cumbersome as it depends on several parameters. Consequently, we do not present it here and the interested reader is referred to Chapters 6 and 7 of Lukacs [21].
  • For γ = 0 , the corresponding NEF F is the normal one with variance equaling the constant δ . F is steep with C = M = R + , α n = n 1 / 2 and g n ( θ ) = n 1 / 2 θ .
  • For γ ( 0 , 1 ) , no NEF exists with VF in the form (9).
  • For γ = 1 , the corresponding NEF F is Poisson. F is steep with C = M = R + , α n 1 and g n ( θ ) = n 1 / 2 θ .
  • For γ ( 1 , 2 ) , the corresponding NEF F is a compound Poisson NEF generated by gamma distributions. F is steep with C = M = R + , α n = n 1 γ γ 2 and g n ( θ ) = θ n γ 1 γ 2 . Here, the corresponding cumulative distribution function is given by
    P ( X x ) = e p E ( x ) + e p 0 x e θ t n = 1 b n t n ρ 1 n ! Γ ( n ρ ) d t ,
    where E ( x ) is a cumulative distribution function degenerated at 0 and
    ρ = 2 γ 1 γ , p = ρ 1 δ ρ δ θ ρ 1 ρ , b = ρ 1 1 ρ δ 1 ρ , ρ < 0 , θ < 0 and δ > 0 .
  • For γ = 2 , the corresponding NEF is gamma one with shape parameter δ 1 . However, as was shown in Bar-Lev and Enis (1986) it is not reproducible when considered as a one-parameter NEF. It is reproducible when considered as a two-parameter NEF (see part 5 of the previous section).
  • For γ ( 2 , ) , the corresponding NEF is generated by a positive stable distribution with stable index 0 < α < 1 , where α = ( γ 2 ) / ( γ 1 ) . F is steep with C = M = R + , α n = n 1 γ γ 2 and g n ( θ ) = θ n γ 1 γ 2 . Here, the associated absolutely continuous stable probability density is ν ρ ( d x ) = h ρ ( x ) d x (c.f., Bar-Lev and Enis [5]) where
    h ρ ( x ) = 1 π k = 0 ( 1 ) k k ! sin ( π ρ x ) 1 ρ k ( 1 ρ ) Γ ( ρ k + 1 ) ρ k δ k ( 1 ρ ) x ρ k + 1 , x > 0 , θ < 0 , α > 0 , 0 < ρ < 1 .
The above NEFs are displayed in the following table along with their corresponding α n , g n ( θ ) , k γ ( θ ) , Θ and M .

4. Conclusions and Topics for Further Research

Convolutions of n i.i.d. random variables play an important role in statistical inference, probability and stochastic processes. However, they are typically very complex, intricate and cumbersome to calculate. Reproducible families of distributions in the sense of Definition 1 are, therefore, very useful in allowing a simple computation of such convolutions by employing (4). The convolution of i.i.d. r.v.s has the same distribution type as any of its components up to a dilation. Various applications of reproducibility and convolution are presented in part 5 of Section 2.
In this study, we have classified all one-parameter reproducible NEFs in the sense of Definition 1 and showed that their corresponding VFs are a power of their mean, as in (9). However, there are still various research avenues concerning the notion of reproducibility that need to be explored. Below are some examples:
1.
Reproducibility for the multi-parameter case
When F is a multi-parameter family of distributions depending, say, on k parameters ( θ 1 , , θ k ) , then under which conditions is it reproducible in the sense of Definition 1? As we have already seen from the one-parameter case, such conditions cannot be found for general families. We will, therefore, need to limit the analysis to families of distributions with certain characteristics such as multi-parameter NEFs. Indeed, in Part 5 of Section 2, we already indicated that for a sub-class of two-parameter NEFs, which include the gamma, normal and inverse Gaussian NEFs with parameters ( θ 1 , θ 2 ) ; the distribution of X ¯ n is also gamma, normal and inverse Gaussian but with parameters ( n θ 1 , n θ 2 ) . For these three cases, the sequence of stabilizing constants is α n = n 1 with g n 1 ( θ 1 ) , g n 2 ( θ 1 ) = ( n θ 1 , n θ 2 ) .
2.
Reproducibility for non-NEFs families
We have seen in Examples 1 and 2 that one of the terminal points of the support of the reproducible families of distributions presented depends on the parameter θ . Accordingly, one might consider the reproducibility property for huge classes of distributions for which at least one of the terminal points of their respective support depends on a parameter (as the uniform distribution on the interval ( θ 1 , θ 2 ) or the shifted exponential distribution with scale and location parameters θ 1 and θ 2 , respectively). One class of this type was introduced by Hogg and Craig [22]. Another class of truncated exponential families was introduced by Bar-Lev [23] and further elaborated by numerous authors in various fields of statistical inference, c.f., Vancak, Goldberg, Bar-Lev, and Boukai, B. [24] and the referees cited therein.
3.
Reproducibility and infinite divisibility
We have seen that all one-parameter reproducible NEFs are also infinitely divisible. On the other hand, Examples 1–3 demonstrate that there are reproducible families that are not infinitely divisible. Accordingly, one can pose the following challenging problem: Under which necessary and/or sufficient conditions are reproducible families also infinitely divisible and vice versa?

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

I am thankful to two reviewers for helpful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Wilks, S. Mathematical Statistics; Wiley: New York, NY, USA, 1963. [Google Scholar]
  2. Bolger, E.M.; Harkness, W.L. Some characterizations of exponential-type distributions. Pacific J. Math. 1966, 16, 5–11. [Google Scholar] [CrossRef] [Green Version]
  3. Baringhaus, L.; Davies, P.L.; Plachky, D. A characterization of the Poisson distribution by convolution properties of one-parameter exponential families. Z. Angew. Math. Mech. 1976, 56, T333–T334. [Google Scholar] [CrossRef]
  4. Bar-Lev, S.K.; Enis, P. Reproducibility in the one-parameter exponential family. Metika 1985, 32, 391–394. [Google Scholar] [CrossRef]
  5. Bar-Lev, S.K.; Enis, P. Reproducibility and natural exponential families with power variance functions. Ann. Stat. 1986, 14, 1507–1522. [Google Scholar] [CrossRef]
  6. Bar-Lev, S.K.; Casalis, M. Les familles exponentielles naturelles reproduisantes. C. R. Acad. Sci. Paris Ser. I 1994, 319, 1323–1326. [Google Scholar]
  7. Bar-Lev, S.K.; Casalis, M.A. Classification of reducible natural exponential families in the broad sense. J. Theor. Probab. 2003, 16, 175–196. [Google Scholar] [CrossRef]
  8. Tweedie, M.C.K. An index which distinguishes between some important exponential families. In Statistics: Applications and New Directions, Proceedings of the Indian Institute Golden Jubilee International Conference, Calcutta, India, 27 September–1 October 1984; Ghosh, J.K., Roy, J., Eds.; Indian Statistical Institute: Calcutta, India, 1984; pp. 579–604. [Google Scholar]
  9. Jørgensen, B. Exponential dispersion models (with discussion). J. R. Stat. Soc. Ser. B 1987, 49, 127–162. [Google Scholar]
  10. Bar-Lev, S.K. Independent, tough Identical results: The class of Tweedie on power variance functions and the Class of Bar-Lev and Enis on reproducible natural exponential families. Int. J. Stat. Probab. 2020, 9, 30–35. [Google Scholar] [CrossRef]
  11. Feller, W. An Introduction to Probability Theory and Its Applications 2; Wiley: New York, NY, USA, 1966. [Google Scholar]
  12. Shushi, T.; Yao, J. Multivariate risk measures based on conditional expectation and systemic risk for Exponential Dispersion Models. Insur. Math. Econ. 2020, 93, 178–186. [Google Scholar] [CrossRef]
  13. Bahnemann, D. Distributions for Actuaries; Casualty Actuarial Society: Arlington, VA, USA, 2015. [Google Scholar]
  14. Bar-Lev, S.K.; Carlo, A.R.M. Methods for insurance risk computations. Int. J. Stat. Probab. 2019, 8, 54–74. [Google Scholar] [CrossRef]
  15. Bar-Lev, S.K.; Reiser, B. An exponential subfamily which admits UMPU tests based on a single test statistic. Ann. Stat. 1982, 10, 979–989. [Google Scholar] [CrossRef]
  16. Letac, G.; Mora, M. Natural real exponential families with cubic variance functions. Ann. Stat. 1990, 18, 1–37. [Google Scholar] [CrossRef]
  17. Mora, M. Classification des fontions-variance cubiques des families exponentielles sur ℝ. C. R. Acad. Sci. Paris Ser. 1986, 116, 582–591. [Google Scholar]
  18. Mora, M. La convergence des fonctions variance des familles exponentielles naturelles. Ann. Fac. Sci. Toulouse 1990, 11, 105–120. [Google Scholar] [CrossRef]
  19. Barndorff-Nielsen, O. Information and Exponential Families in Statistical Theory; Wiley: New York, NY, USA, 1978. [Google Scholar]
  20. Aczél, J. Functional Equations and Their Applications; Academic Press: New York, NY, USA, 1966. [Google Scholar]
  21. Lukacs, E. Characteristic Functions, 2nd ed.; Hafner: New York, NY, USA, 1970. [Google Scholar]
  22. Hogg, R.V.; Craig, A.T. Sufficient statistics in elementary distribution theory. Sankhya 1956, 17, 209–216. [Google Scholar]
  23. Bar Lev, S.K. Large sample properties of the MLE and MCLE for the natural parameter of a truncated exponential family. Ann. Inst. Stat. Math. Part A 1984, 36, 217–222. [Google Scholar] [CrossRef]
  24. Vancak, V.; Goldberg, Y.; Bar-Lev, S.K.; Boukai, B. Continuous statistical models: With or without truncation parameter? Math. Methods Stat. 2015, 24, 55–73. [Google Scholar] [CrossRef] [Green Version]
Table 1. NEFs having Power VFs.
Table 1. NEFs having Power VFs.
γ NEF Type Θ MC k γ ( θ ) α ( n ) g n ( θ )
( , 0 ) Extreme stable R + R + R 1 a ( 2 γ ) a ( 1 γ ) θ γ 2 γ 1 n 1 γ γ 2 θ n γ 1 γ 2
0Normal R R R 1 2 a θ 2 ± n 1 / 2 ± n 1 / 2 θ
1Poisson R R + R 0 + 1 a e a θ 1 θ + 1 a ln n
1 , 2 Compound Poisson R R + R 0 + 1 a ( 2 γ ) a ( 1 γ ) θ γ 2 γ 1 n 1 γ γ 2 θ n γ 1 γ 2
2gamma R R + R + 1 a ln ( 1 a θ )
( 2 , ) positive stable R R + R + 1 a ( 2 γ ) a ( 1 γ ) θ γ 2 γ 1 n 1 γ γ 2 θ n γ 1 γ 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bar-Lev, S.K. On the Notion of Reproducibility and Its Full Implementation to Natural Exponential Families. Mathematics 2021, 9, 1568. https://doi.org/10.3390/math9131568

AMA Style

Bar-Lev SK. On the Notion of Reproducibility and Its Full Implementation to Natural Exponential Families. Mathematics. 2021; 9(13):1568. https://doi.org/10.3390/math9131568

Chicago/Turabian Style

Bar-Lev, Shaul K. 2021. "On the Notion of Reproducibility and Its Full Implementation to Natural Exponential Families" Mathematics 9, no. 13: 1568. https://doi.org/10.3390/math9131568

APA Style

Bar-Lev, S. K. (2021). On the Notion of Reproducibility and Its Full Implementation to Natural Exponential Families. Mathematics, 9(13), 1568. https://doi.org/10.3390/math9131568

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