The tempered discrete Linnik distribution
LUCIO BARABESI
Department of Economics & Statistics, University of Siena, Italy

arXiv:1605.02326v1 [math.ST] 8 May 2016

CAROLINA BECATTI
Department of Economics & Statistics, University of Siena, Italy
MARZIA MARCHESELLI
Department of Economics & Statistics, University of Siena, Italy
ABSTRACT. A tempered version of the discrete Linnik distribution is introduced in order to obtain integer-valued distribution
families connected to stable laws. The proposal constitutes a generalization of the well-known Poisson-Tweedie law, which is actually a tempered discrete stable law. The features of the new
tempered discrete Linnik distribution are explored by providing a
series of identities in law - which describe its genesis in terms of
mixture and compound Poisson law, as well as in terms of mixture
discrete stable law. A manageable expression of the corresponding
probability function is also provided and several special cases are
analysed.
Key words: Lévy-Khintchine representation; Positive Stable distribution;
Discrete Stable distribution; Mixture Poisson distribution; Compound Poisson distribution.
MSC 2010 Primary 62E10, secondary 62E15.

1. Introduction
In recent years, heavy-tailed models - in primis, stable distributions - have
been used in a variety of fields, such as statistical physics, mathematical finance and financial econometrics (see e.g. Rachev et al., 2011, and references
therein). However, these models may be partially appropriate to provide a
good fit to data, since their tails are too “fat” to describe empirical distributions, as remarked by Klebanov and& Slámová (2015). In order to overcome
this drawback, the so-called “tempered” versions of heavy-tailed distributions have been successfully introduced (see e.g. Rosı́nski, 2007). Indeed,
tempering allows for models that are similar to original distributions in some
central region, even if possess lighter - i.e. tempered - tails. Klebanov &

1

Slámová (2015) have suitably discussed these issues and have suggested various tempering techniques - emphasizing that tempering is not necessarily
unambiguous.
In the framework of integer-valued laws, the Poisson-Tweedie is a wellknown tempered distribution introduced independently by Gerber (1991) and
Hougaard et al. (1997). This law is de facto the tempered counterpart of
the Discrete Stable law originally suggested by Steutel & van Harn (1979).
The Poisson-Tweedie distribution encompasses classical families (such as the
Poisson), as well as large families (such as the Generalized Poisson Inverse
Gaussian and the Poisson-Pascal). Hence, this law may be very useful for
modelling data arising in a plethora of frameworks - for example, clinical
experiments (Hougaard, Lee & Whitmore, 1997), environmental studies (ElShaarawi, Zhu & Joe, 2011) and scientometric analysis (Baccini, Barabesi &
Stracqualursi, 2016).
Christoph & Schreiber (1998) emphasize that the Discrete Stable distribution may be seen as the special case - for the limiting value of a parameter
- of the so-called Discrete Linnik distribution introduced by Devroye (1993)
and Pakes (1995). Therefore, the Discrete Linnik distribution is more flexible
than the Discrete Stable distribution and it could be very useful to achieve
its tempered version. In the present paper, we preliminarily explore in detail
the method - roughly outlined by Barabesi & Pratelli (2014a) - for obtaining
integer-valued families of distributions linked to stable and tempered stable
laws. On the basis of these findings, after revising some properties of the
Discrete Linnik distribution, we introduce its tempered counterpart. The
new Tempered Discrete Linnik is analysed thoroughly and its properties are
given. In particular, some stochastic representations for this law are obtained
and the corresponding probability function is achieved as a manageable finite
sum.
The paper is organized as follows. In Sections 2 and 3, we revise and
expand the issues suggested by Barabesi & Pratelli (2014a) for introducing
families of distributions connected to stable and tempered stable laws, respectively. In Section 4, we survey the main features of the Discrete Linnik
distribution. Section 5 contains our proposal for the new Tempered Discrete Linnik distribution, while in Section 6 we consider its main properties.
Finally, some conclusions are given in Section 7.
2. Integer-valued distribution families linked to stable laws
Barabesi & Pratelli (2014a) have suggested an approach - based on the definition of subordinator - for devising integer-valued distribution families as
mixtures of stable laws, as well as tempered stable laws. In order to describe
2

and to develop at length their proposal, we first consider the absolutelycontinuous Positive Stable random variable (r.v.) - say XP S - with Laplace
transform given by
LXP S (t) = exp(−λtγ ) , Re(t) > 0,

(1)

where (γ, λ) ∈]0, 1] × R+ and R+ =]0, ∞[ (see e.g. Zolotarev, 1986, p.114).
As is well known, γ is the so-called characteristic exponent - i.e. a “tail”
index - while λ is actually a scale parameter. In order to emphasize the
dependence on γ and λ, we eventually adopt the notation XP S := XP S (γ, λ).
We have also to introduce the integer-valued counterpart of the Positive
Stable r.v., i.e. the Discrete Stable r.v. XDS proposed by Steutel & van
Harn (1979) with probability generating function (p.g.f.) given by
gXDS (s) = exp(−λ(1 − s)γ ) , s ∈ [0, 1],

(2)

where in turn (γ, λ) ∈]0, 1] × R+ . For a survey of the properties of this
distribution, see e.g. Marcheselli, Barabesi & Baccini (2008). Similarly to
the Positive Stable r.v., we also write XDS := XDS (γ, λ). The parameter γ
is again a “tail” index, while λ is a “scale” parameter.
In order to clarify the meaning of the scale parameter for the Discrete
Stable law, we remind the “thinning” operator
introduced by Steutel &
van Harn (1979). If X is an integer-valued r.v., the dot product α X is
defined as
X
X
α X=
Zi ,
i=1

where the Zi ’s are copies of Bernoulli r.v’s of parameter α ∈ [0, 1] independent
of X. Obviously, the p.g.f. of α X is given by gα X (s) = gX (1 − α + αs),
where gX is the p.g.f. of X. The dot product is also defined for α > 1,
whenever gα X is a proper p.g.f. (see Christoph & Schreiber, 2001). Hence,
λ is a “scale” parameter for the Discrete Stable law in the sense that
L

XDS (γ, λ) = λ1/γ

XDS (γ, 1).
R
Let ν be a measure on R+ in such a way that min(1, x)ν(dx) < ∞. From
the Lévy-Khintchine representation (see e.g. Sato, 1999, p.197) there exists
a positive r.v. Y with Laplace transform given by
LY (t) = exp(−ηψ(t)) , Re(t) > 0 ,
R
where η ∈ R+ and ψ(t) = (1−exp(−tx))ν(dx). Moreover, let XP := XP (λ)
represent a Poisson r.v. with parameter λ, i.e. the p.g.f. of XP (λ) is given
3

by gXP (s) = exp(−λ(1 − s)) with s ∈ [0, 1]. Hence, if the r.v.’s XP and Y
are independent, the Mixture Poisson r.v.
L

XM P = XP (Y )
displays the p.g.f. given by
gXM P (s) = LY (1 − s) = exp(−ηψ(1 − s)) , s ∈ [0, 1].

(3)

The Discrete Stable r.v. XDS (γ, λ) is obtained from expression (3) when a
stable subordinator is actually considered, i.e. ν(dx)/dx ∝ x−γ−1 IR+ (x) with
γ ∈]0, 1[ and where IB represents the indicator function of the set B. Indeed,
since the expression of ν gives rise to ψ(t) ∝ tγ , the p.g.f. (2) promptly
follows from expression (3) with η = λ and ψ(t) = tγ . Moreover, since
L
gXDS (s) = LXP S (1 − s), it follows that Y = XP S (γ, λ) and
L

XDS (γ, λ) = XP (XP S (γ, λ)),

(4)

which is actually equivalent to the identity in distribution emphasized by
Devroye (1993, Theorem in Section 1). On the basis of expression (4), a
general “scale” mixture of Discrete Stable r.v.’s, say XM DS , with a mixturing
absolutely-continuous positive r.v. V having Laplace transform LV , may be
achieved by considering the identity in distribution
L

L

XM DS = XDS (γ, V ) = XP (XP S (γ, V )),

(5)

where the r.v.’s involved in the right-hand side are independent. Obviously,
(4) is achieved from (5) by assuming a degenerate distribution for V , i.e.
P (V = λ) = 1. In addition, expression (5) may be also seen as the stochastic
“scaling” of the Discrete stable law in terms of the operator , i.e.
L

XM DS = V 1/γ

L

XDS (γ, 1) = V 1/γ

XP (XP S (γ, 1)) .

Moreover, from (5), it is apparent that the p.g.f. of the r.v. XM DS turns out
to be
gXM DS (s) = LV ((1 − s)γ ) , s ∈ [0, 1].
(6)
Hence, families of mixture of Discrete Stable r.v.’s can be generated on the
basis of (5) and (6) by suitably selecting the r.v. V .
We conclude with a final remark on the p.g.f. (6). Let XS be a Sibuya
r.v. (as named by Devroye, 1993) with p.g.f.
gXS (s) = 1 − (1 − s)γ , s ∈ [0, 1],
4

where γ ∈]0, 1] (for a recent survey of this law, see Huillet, 2016). Indeed, the
Sibuya distribution is a special case of the (shifted) Negative Binomial Beta
distribution introduced by Sibuya (1979) with parameters given by 1, γ and
(1 − γ). In the following, the Sibuya r.v. is also denoted by XS := XS (γ).
Therefore, expression (6) may be also interestingly reformulated as
gXM DS (s) = LV (1 − gXS (s)) , s ∈ [0, 1].

(7)

Thus, if LV displays a suitable structure, expression (7) eventually gives rise
to a representation of the r.v. XM DS in terms of a compound r.v. with a
compounding Sibuya r.v. As a quite easy example, the r.v. XDS may be also
expressed as a compound Poisson r.v. as
gXDS (s) = exp(−λ(1 − gXS (s))) , s ∈ [0, 1],
and hence
L

XDS (γ, λ) =

Z
X

Wi ,

i=1
L

L

where Z = XP (λ) and the Wi ’s are i.i.d. r.v.’s such that Wi = XS (γ) - which
are in turn independent of Z.
3. Integer-valued distribution families linked to tempered stable
laws
First, for subsequent use, we provide some issues on the so-called Tweedie
distribution (for a recent survey of this law, see Barabesi et al., 2016). The
Tweedie distribution is actually a Tempered Positive Stable distribution introduced by Hougaard (1986). Hence, for this reason, in the following we
denote the Tweedie r.v. as XT P S . With a slight change in the parameterization proposed by Hougaard (1986), the Laplace transform of the r.v. XT P S
is given by
LXT P S (t) = exp(sgn(γ)λ(θγ − (θ + t)γ )) , Re(t) > 0,

(8)

where (γ, λ, θ) ∈ {] − ∞, 1] × R+ × R+ } ∪ {]0, 1] × R+ × {0}}. The formulation
proposed in expression (8) is convenient, since it avoids to define the Laplace
transform for analytical continuity for γ = 0 as in the case of the parameterization considered by Hougaard (1986). Moreover, it is worth noting that θ
is actually the “tempering” parameter. This is at once apparent for γ ∈]0, 1]
from the following identity
LXT P S (t) =

LXP S (θ + t)
,
LXP S (θ)
5

which reveals the exponential “nature” of the tempering. By following the
usual route, we also write XT P S := XT P S (γ, λ, θ). Obviously, for γ ∈]0, 1] it
holds
L
XT P S (γ, λ, 0) = XP S (γ, λ).
It should be strongly remarked that the tempering extends the range of
parameter values (with respect to the Positive Stable distribution) for the
parameter γ - which may assume negative values, even if θ must be strictly
positive in such a case. This is an interesting feature, since for γ ∈ R−
where R− =] − ∞, 0[, it is immediate to reformulate the r.v. XT P S as a
compound Poisson of Gamma r.v.’s (see e.g. Barabesi et al., 2016). More
precisely, let the r.v. XG := XG (λ, δ) be distributed according to a Gamma
law with corresponding Laplace transform given by LXG (t) = (1 + λt)−δ
with Re(t) > 0 and (λ, δ) ∈ R+ × R+ (obviously, λ is the scale parameter
and δ is the shape parameter). Hence, on the basis of (8) and owing to the
reproductive property of the Gamma distribution with respect to the shape
parameter, the following identity in distribution holds for γ ∈ R−
L

XT P S (γ, λ, θ) = XG (1/θ, −γXP (λθγ )).

(9)

Hence, in this case the r.v. XT P S displays a mixed distribution, given by
a convex combination of a Dirac distribution (with mass at zero) and an
absolutely-continuous distribution (a very useful property for modelling data
with an excess of zeroes, see Barabesi et al., 2016).
In order to extend the issues of Section 2, we introduce a tempered stable
subordinator, i.e.
ν(dx)/dx ∝ exp(−θx)x−γ−1 IR+ (x),
where (γ, θ) ∈ {] − ∞, 1[×R+ } ∪ {]0, 1[×{0}}. In such a case, on the basis of
(3) the so-called Poisson-Tweedie distribution - which is actually a Tempered
Discrete Stable distribution - is achieved (for more about the Poisson-Tweedie
law, see Baccini, Barabesi & Stracqualursi, 2016, and El-Shaarawi, Zhu &
Joe, 2011). Indeed, if XT DS denotes a Tempered Discrete Stable r.v., since
the expression of ν provides ψ(t) ∝ (θ + t)γ − θγ , from (3) with η = λ and
ψ(t) = (θ + t)γ − θγ the following p.g.f. is obtained
gXT DS (s) = exp(sgn(γ)λ(θγ − (θ + 1 − s)γ )) , s ∈ [0, 1].
Hence - as expected - it holds gXT DS (s) = LXT P S (1 − s). However, in order to
be consistent with the existing literature and for practical convenience, we
prefer to reparameterize the previous p.g.f. by assuming that γ = a, λ = bca
and θ = 1/c − 1, in such a way that
gXT DS (s) = exp(sgn(a)b((1 − c)a − (1 − cs)a )), s ∈ [0, 1],
6

(10)

where (a, b, c) ∈ {] − ∞, 0] × R+ × [0, 1[} ∪ {]0, 1] × R+ × [0, 1]}. The p.g.f.
(10) is provided as a slight modification of the formulation suggested by
El-Shaarawi et al. (2011). Indeed, the parameterization considered in (10)
avoids to define the p.g.f. for analytical continuity for a = 0. In this special
case, XT DS degenerates at 0, i.e. P (XT DS = 0) = 1. As usual, we write
XT DS := XT DS (a, b, c).
It is worth noting that c actually represents the “tempering” parameter.
Indeed, for a ∈]0, 1] and by considering the r.v. XDS (a, b), the following
identity
gX (cs)
LXT DS (t) = DS
gXDS (c)
emphasizes the geometric “nature” of the tempering. In addition, it straightL
forwardly holds XDS (a, b) = XT DS (a, b, 1) for a ∈]0, 1]. In turn, similarly to
the Tweedie distribution, tempering extends the range of parameter values
(with respect to the Discrete Stable distribution) for the parameter a - which
may assume negative values, even if c must be strictly less than unity in such
L
a case. Finally, from (3) with Y = XT P S (a, bca , 1/c − 1), it holds
L

XT DS (a, b, c) = XP (XT P S (a, bca , 1/c − 1)),

(11)

which actually constitutes the identity in distribution remarked by Hougaard
et al. (1997) and which generalizes expression (4).
On the basis of expression (11), it is at once apparent that a general
“scale” mixture of Tempered Discrete Stable r.v.’s, say XM T DS , with a mixturing absolutely-continuous positive r.v. V having Laplace transform LV ,
may be achieved by considering the identity in distribution - which generalizes (5) - given by
L

L

XM T DS = XT DS (a, V, c) = XP (XT P S (a, ca V, 1/c − 1)),

(12)

where the r.v.’s involved in the right-hand side are independent. Obviously,
(11) is achieved from (12) by assuming a degenerate distribution for V , i.e.
P (V = b) = 1. Moreover, the corresponding p.g.f. turns out to be
gXM T DS (s) = LV (sgn(a)((1 − cs)a − (1 − c)a )) , s ∈ [0, 1].

(13)

Hence, families of mixture of Tempered Discrete Stable r.v.’s can be generated by means of (12) - and accordingly (13) - by suitably selecting the
r.v. V . In order to reformulate (13) similarly to (7) when a ∈]0, 1], let
XGDS := XGDS (γ, τ ) be the Geometric Down-weighting Sibuya r.v. (introduced by Zhu & Joe, 2009) with the p.g.f.
gXGDS (s) = 1 − gXS (τ ) + gXS (τ s) = 1 + (1 − τ )γ − (1 − τ s)γ ,
7

where (γ, τ ) ∈]0, 1]×]0, 1]. Hence, by considering the Geometric Downweighting Sibuya r.v. XGDS (a, c), for a ∈]0, 1] expression (13) may be also
interestingly rewritten as
gXM T DS (s) = LV (1 − gXGDS (s)) , s ∈ [0, 1].

(14)

Moreover, let XN B := XN B (π, δ) be a Negative Binomial r.v. with p.g.f.
given by

gXN B (s) =

1−π
1 − πs

δ

= (1 + λ(1 − s))−δ , s ∈ [0, 1],

where (π, δ) ∈]0, 1[×R+ , while λ = π/(1 − π) ∈ R+ . Thus, when a ∈ R− and
by considering the Negative Binomial r.v. XN B (c, −a), expression (13) may
be also rewritten as
gXM T DS (s) = LV ((1 − c)a (1 − gXN B (s))) , s ∈ [0, 1].

(15)

Therefore, if LV displays a suitable structure, expression (14) and (15) eventually gives rise to a representation of the r.v. XM DS in terms of a compound
r.v. with compounding Sibuya and Negative Binomial r.v.’s, respectively.
4. Some remarks on the genesis and properties of the discrete
Linnik distribution
By following the formulation adopted by Christoph & Schreiber (1998), the
p.g.f. of the r.v. XDL distributed according to the Discrete Linnik law is
given by
gXDL (s) = (1 + λ(1 − s)γ /δ)−δ , s ∈ [0, 1],
(16)
where (γ, λ, δ) ∈]0, 1] × R+ × R+ . For a detailed description of the main
features of the law, see Christoph & Schreiber (1998, 2001). It is apparent
that the p.g.f. of the Discrete Stable r.v. XDS is achieved as δ → ∞. In
addition, it should be remarked that the Discrete Linnik law is defined for
some negative δ in such a way that λ ≤ |δ|(1 − γ) and, in this case, the
distribution is also named Generalized Sibuya (Huillet, 2016). However, the
results of this Section are mainly given for positive δ. Thus, it is apparent that
the Discrete Linnik law is very flexible and may encompass a large variety
of distribution families ranging from light-tailed laws (e.g. the Binomial law
for suitable λ and negative integer δ when γ = 1) to heavy-tailed laws (e.g.
the Discrete Mittag-Leffler law for δ = 1).
In order to emphasize the dependence on the parameters, we also adopt
the notation XDL := XDL (γ, λ, δ). It should be remarked that the parameter
8

γ is in turn a “tail” index and λ is a “scale” parameter, while δ is a shape
parameter. Indeed, λ is a “scale” parameter for the Discrete Linnik law in
the sense that
L
XDL (γ, λ, δ) = λ1/γ XDL (γ, 1, δ).
For δ ∈ R+ and on the basis of the remarks provided in Section 2, it is
promptly proven that the Discrete Linnik r.v. is a “scale” mixture of Discrete
L
Stable r.v.’s obtained by selecting V = XG (λ/δ, δ) in (5). In this case, the
p.g.f. (16) is obtained by means of (6), while from (5) it also holds
L

L

XDL (γ, λ, δ) = XDS (γ, XG (λ/δ, δ)) = XP (XP S (γ, XG (λ/δ, δ))),

(17)

which is actually similar to the identity in distribution obtained by Devroye
(1993, Section 2). A further remark is in order, since expression (17) may
be also rephrased in terms of the absolutely-continuous Positive Linnik r.v.
XP L := XP L (γ, λ, δ) with Laplace transform given by
LXP L (t) = (1 + λtγ /δ)−δ , Re(t) > 0,
where (γ, λ, δ) ∈]0, 1] × R+ × R+ - according to the formulation provided by
Christoph & Schreiber (2001). The absolutely-continuous Positive Linnik has
been introduced by Pakes (1995) after that Linnik (1962, p.67) had identified
its symmetric version in a simple case. Obviously, the Laplace transform (1)
of the Positive Stable r.v. XP S (γ, λ) is achieved as δ → ∞. For more details
on the Positive Linnik distribution, see e.g. Jose et al. (2010). Since it easily
proven that
L
XP L (γ, λ, δ) = XP S (γ, XG (λ/δ, δ)),
where the r.v.’s involved in the right-hand side are independent, by means
of expression (17) it promptly follows that
L

XDL (γ, λ, δ) = XP (XP L (γ, λ, δ)),

(18)

which actually generalizes (4). Indeed, expression (4) is recovered from expression (18) as δ → ∞. Finally, we remark that the identity in distribution
(17) is very suitable for random variate generation. Indeed, many generators
for Poisson and Gamma variates are available in statistical literature, while
Positive Stable variates are readily obtained by means of the well-known
Kanter’s representation (Kanter, 1975). A similar remark applies to identity
(4).
For δ ∈ R+ , the r.v. XDL may be also expressed as a compound Negative
Binomial r.v. with a compounding Sibuya r.v. Indeed, from expression (7)
L
with V = XG (λ/δ, δ), it follows that
gXDL (s) = (1 + λ(1 − gXS (s))/δ)−δ , s ∈ [0, 1] ,
9

and hence
L

XDL (γ, λ, δ) =

Z
X

Wi ,

i=1
L

L

where Z = XN B (λ/(δ + λ), δ) and the Wi ’s are i.i.d. r.v.’s such that Wi =
XS (γ) - which are in turn independent of Z.
5. The tempered discrete Linnik distribution
On the basis of the issues discussed in Sections 2, 3 and 4, we are ready to
introduce a tempered version of the Discrete Linnik distribution. Indeed,
in a complete parallelism with Section 4, the Tempered Discrete Linnik r.v.
may be given as a “scale” mixture of the Tempered Discrete Stable r.v.’s
L
by selecting V = XG (bd, 1/d) in (12). In this case, if the Tempered Discrete
Linnik r.v. is denoted by XT DL , from expression (13) the corresponding p.g.f.
turns out to be
gXT DL (s) = (1 + sgn(a)bd((1 − cs)a − (1 − c)a ))−1/d , s ∈ [0, 1],

(19)

with (a, b, c, d)∈{] − ∞, 0] × R+ × [0, 1[×R+ } ∪ {]0, 1] × R+ × [0, 1] × R+ }.
Several comments on the parameterization are in order. First, similarly to the
Discrete Linnik law, the Tempered Discrete Linnik law is also defined for some
negative δ - a distribution that could be named Tempered Generalized Sibuya
by extending the definition by Huillet (2016) - even if we mainly deal with
positive d. Secondly, the adopted formulation allows to achieve the PoissonTweedie distribution as d approaches zero - and this could be preferable with
respect to a limit at infinity. Moreover, for a = 0 the r.v. XT DL degenerates
at 0. Thus, the p.g.f. (19) seems a natural generalization of the p.d.f. (10).
As usual, we also adopt the notation XT DL := XT DL (a, b, c, d). Obviously, it
L
promptly holds XDL (a, b, d) = XT DL (a, b, 1, d) for a ∈]0, 1]. Analogously to
the Poisson-Tweedie distribution, it is worth noting that tempering extends
the range of parameter values for the parameter a with respect to the Discrete
Linnik distribution.
For d ∈ R+ , similarly to (11) and by means of (12), the following identity
in distribution holds
L

XT DL (a, b, c, d) = XT DS (a, XG (bd, 1/d), c))

(20)

L

= XP (XT P S (a, XG (bdca , 1/d), 1/c − 1)
which actually generalizes (17). We remark that the identity (20) is suitable for random variate generation, if Tempered Positive Stable variates are
10

available. Such a generation is straightforward on the basis of (9) if a ∈ R− ,
since in this case
L

XT DL (a, b, c, d) = XP (XG (c/(1 − c), −aXP (XG (bd(1 − c)a , 1/d)))).

(21)

In contrast, for a ∈]0, 1] the algorithms proposed by Barabesi et al. (2016)
and Devroye (2009) for the generation of Tempered Positive Stable variates
may be considered. In any case, Barabesi & Pratelli (2014b, 2015) have suggested algorithms based on a rejection method linked to the p.g.f., which may
be suitable for the direct generation of Tempered Discrete Linnik variates.
Actually, Baccini, Barabesi & Stracqualursi (2016) have introduced a similar
proposal for the Tempered Discrete Stable law.
It should be remarked that - on the basis of (3) - an absolutely-continuous
Tempered Positive Linnik r.v. XT P L := XT P L (γ, λ, θ, δ) may be given with
Laplace transform given by
LXT P L (t) = (1 + sgn(γ)λ((θ + t)γ − θγ )/δ)−δ , Re(t) > 0,
where (γ, λ, θ, δ)∈{]−∞, 1]×R+ ×R+ ×R+ }∪{]0, 1]×R+ ×{0}×R+ }. In turn,
the Laplace transform (8) of the Tempered Positive Stable r.v. XT P S (γ, λ)
is achieved as δ → ∞. Since it holds that
L

XT P L (γ, λ, θ, δ) = XT P S (γ, XG (λ/δ, δ), θ),
where the r.v.’s involved in the right-hand side are independent, from expression (20) it promptly follows that
L

XT DL (a, b, c, d) = XP (XT P L (a, b, c, 1/d)),
which actually generalizes (18).
For d ∈ R+ , on the basis of (14) and (15), the r.v. XT DL may be also
expressed as a compound distribution. Indeed, when a ∈]0, 1] from (14) it
holds
gXT DL (s) = (1 + bd(1 − gXGDS (s)))−1/d , s ∈ [0, 1],
and hence
L

XT DL (a, b, c, d) =

Z
X

Wi ,

(22)

i=1
L

where Z = XN B (bd/(1 + bd), 1/d) and the Wi ’s are i.i.d. r.v.’s such that
L
Wi = XGDS (a, c) - which are in turn independent of Z. Moreover, when
a ∈ R− expression (15) may be rewritten as
gXT DL (s) = (1 + bd(1 − c)a (1 − gXN B (s)))−1/d , s ∈ [0, 1],
11

L

and hence (22) holds with Z = XN B (bd(1 − c)a /(1 + bd(1 − c)a ), 1/d) and
L
the Wi ’s are i.i.d. r.v.’s such that Wi = XN B (c, −a) - which are in turn independent of Z. Owing to the reproductive property of the Negative Binomial
law, expression (22) also provides
L

XT DL (a, b, c, d) = XN B (c, −aXN B (bd(1 − c)a /(1 + bd(1 − c)a ), 1/d)). (23)
It is worth remarking that identities (21) and (23) are straightforwardly
L
equivalent since XN B (π, δ) = XP (XG (π/(1 − π), δ)).
6.Further features of the tempered discrete Linnik distribution
First, in order to obtain a closed form for the probability function (p.f.) of
the Tempered Discrete Linnik r.v., we provide a general result on the p.f. of
r.v.’s displaying a general type of p.g.f. - which encompasses (13) as special
case. Indeed, let us suppose that the r.v. X has p.g.f. given by
gX (s) = ϕ(α + β(1 − φs)γ ) , s ∈ [0, 1],
where γ, α, β and φ are real constants, in such a way that φ ∈ [0, 1]. If ϕ is
analytic in a neighbourhood of (α + β), it follows that
1 dk
ϕ(α + β + β((1 − φs)γ − 1)|s=0
k! dsk
k
1 X β m dm
dk
=
ϕ(s)|
((1 − φs)γ − 1)m |s=0
s=α+β
k! m=0 m! dsm
dsk
  k
k
m
X
d
1 X β m dm
m−j m
=
ϕ(s)|s=α+β
(1 − φs)jγ
(−1)
m
k
k! m=0 m! ds
j
ds
j=0

P (X = k) =

s=0

k
X
(−β)m dm
ϕ(s)|s=α+β Cγ,m (k),
= (−φ)
m! dsm
m=0
k

where

m
X

  
m γj
Cγ,m (k) =
(−1)
.
j
k
j=0
j

Thus, the p.d.f. of X may be achieved as a finite sum. This result is very
useful for obtaining the p.f. of the r.v. XT DL by setting ϕ(x) = x−1/d and
the values γ = a, α = 1−sgn(a)bd(1 − c)a , β =sgn(a)bd and φ = c, i.e.

k 
X
−1/d
(−sgn(a)bd)m
k
Ca,m (k)
P (XT DL = k) = (−c)
(1 + sgn(a)bd(1 − (1 − c)a ))m+1/d
m
m=0
(24)
12

for k ∈ N. Similarly, the p.f. of Tempered Discrete Stable r.v. is achieved by
setting ϕ(x) = exp(x) and the values γ = a, α =sgn(a)b(1−c)a , β = −sgn(a)b
and φ = c, i.e.
k
X
(sgn(a)b)m
Ca,m (k)
P (XT DS = k) = exp(−sgn(a)b(1 − (1 − c) ))(−c)
m!
m=0
(25)
for k ∈ N. It is worth remarking that (25) is a ctually equivalent to the
expression provided by Baccini, Barabesi & Stracqualursi (2016).
On the basis of these findings, a series of comments on some particular values of the parameters of the Tempered Discrete Linnik distribution
for d ∈ R+ are worthwhile. Indeed, when a = 1, expression (19) provides
L
XT DL (1, b, c, d) = XN B (bcd/(1 + bcd), 1/d) and hence the r.v. XP (bc) is obtained as d → 0. It should be noticed that the Negative Binomial is also
achieved when a ∈]0, 1] and bd(1 − c)a = 1. Moreover, for a = 0, the r.v.
XT DL (0, b, c, d) has a degenerate distribution at zero, as previously remarked.
In addition, for d = 1 a tempered version of the Discrete Mittag-Leffler is
achieved.
For a = 1/2, the Tempered Discrete Linnik law is actually a generalization
of the Poisson Inverse Gaussian law (for a discussion of this law, see e.g.
Johnson et al., 2005, p.484), which is obtained as d → 0. Thus, this new
distribution could be suitably named Negative Binomial Inverse Gaussian.
Moreover, in this special case expression (24) remarkably reduce to a single
sum, since


−2k+m
m 2k − m
k 2
C1/2,m (k) = (−1)
2k − m
k
a

k

for k ∈ N+ , while C1/2,m (0) = 1. In addition, since it holds
C1/2,m (k + 1) = −

2k − m
C1/2,m (k),
4(2k − m + 2)

the probabilities in (24) may be sequentially evaluated - so that the actual
computational burden is small.
When a = −1, the Tempered Discrete Linnik is in turn a generalization
of the Polyá-Aeppli law (for more about this distribution, see e.g. Johnson
et al., 2005, p.410), which is obtained as d → 0. In this case also, expression
(24) and (25) remarkably reduce to a single sum, since


k+m k − 1
C−1,m (k) = (−1)
m−1
13

for k ∈ N+ , while C−1,m (0) = 1. Similarly to the case a = 1/2, the probabilities in (24) may be sequentially evaluated, since it holds
C−1,m (k + 1) = −

m
C−1,m (k).
k

As to the main descriptive indexes, on the basis of (19) and after tedious
algebra, it follows that the expectation and the variance of the r.v. XT DL
are respectively given by
µ = E[XT DL ] = sgn(a)abc(1 − c)a−1
and

(1 − ac)µ
.
1−c
It is worth noting that µ does not depend on the parameter d - i.e., the
Tempered Discrete Linnik r.v. and the Poisson-Tweedie r.v. actually display
the same expectation. However, since the dispersion index is given by
σ 2 = Var[XT DL ] = dµ2 +

D=

1 − ac
σ2
= dµ +
,
µ
1−c

the Poisson-Tweedie distribution may solely display over-dispersion (since
D ≥ 1 when d = 0), while the Tempered Discrete Linnik distribution may
accommodate for under-dispersion (for some admissible values d < 0), as
well as for over-dispersion (when d > 0). Hence, the Tempered Discrete
Linnik law substantially extend the range of the dispersion index with respect
to the Poisson-Tweedie law.
After further tedious algebra, it also follows
c(1 − a)µ
σ4
+ dµσ 2 +
m3 = E[(XT DL − µ) ] =
µ
(1 − c)2
3

and
m4 = E[(XT DL −µ)4 ] = 3(2d+1)σ 4 +

(4c(1 − a) + (1 − ac)2 )σ 2 c2 (1 − a2 )µ
+
,
(1 − c)2
(1 − c)3

from which the skewness and kurtosis indexes may be promptly expressed as
α3 =

m3
D dσ
c(1 − a)
=
+
+
3
σ
σ
D
(1 − c)2 σD

and
α4 =

m4
4c(1 − a) + (1 − ac)2
c2 (1 − a2 )
=
3(2d
+
1)
+
+
.
σ4
(1 − c)2 σ 2
(1 − c)3 σ 2 D
14

On the basis of these expressions, we just provide a sketch of the flexibility
of the Tempered Discrete Linnik distribution in Figures 1-3. Indeed, these
Figures display the parametric plot of (α3 , α4 ) as c and d vary (for fixed a
and b) and it is apparent that the Tempered Discrete Linnik distribution
may cover an extended region in the (α3 , α4 )-plane.
Figure 1 about here
Figure 2 about here
Figure 3 about here
Figure 4 about here
7. Conclusions
We have proposed a tempered version of the Discrete Linnik law. The new
distribution is introduced on the basis of several stochastic mixture and compound representations, which strongly justify its genesis. The Tempered Discrete Linnik law encompasses many well-known distributions as special cases
and displays appealing probabilistic features. In addition, we provide a manageable form for the corresponding probability function, which may enhance
the practical adoption of the Tempered Discrete Linnik law as a statistical
model for many different types of dataset.
Acknowledgements The authors would like to thank Prof. Luca Pratelli
for many useful advices.
References
Baccini, A., Barabesi, L., & Stracqualursi, L.(2016). Random variate generation and connected computational issues for the
Poisson-Tweedie distribution. Computational Statistics, in press, DOI:
10.1007/s00180-015-0623-5.
Barabesi, L., Cerasa, A., Perrotta, D., & Cerioli, A. (2016). Modelling international trade data with the Tweedie distribution for antifraud purposes. European Journal of Operational Research, in press.
Barabesi, L., & Pratelli, L. (2014a). Discussion of “On simulation and
properties of the Stable law” by L. Devroye and L. James, Statistical
Methods and Applications 23, 345–351.
Barabesi, L., & Pratelli, L. (2014b). A note on a universal random variate generator for integer-valued random variables. Statistics and Computing 24, 589–596.

15

Barabesi, L., & Pratelli, L. (2015). Universal methods for generating
random variables with a given characteristic function. Journal of Statistical Computation and Simulation 85, 1679–1691.
Christoph, G., & Schreiber, K. (1998). The generalized Linnik distributions, in Advances in Stochastic Models for Reliability, Quality and Safety
(W. Kahle, E. von Collani, J. Franz and U. Jensen, eds.), Birkhäuser,
Boston, 3-18.
Christoph, G., & Schreiber, K. (2001). Positive Linnik and discrete
Linnik distributions, in Asymptotic Methods in Probability and Statistics
with Applications (N. Balakrishnan, I.A. Ibragimov and V.B. Nevzorov,
eds.), Birkhäuser, Boston, 3–17.
Devroye, L. (1993). A triptych of discrete distribution related to the stable
law. Statistics and Probability Letters 18, 349–351.
Devroye, L. (2009). Random variate generation for exponentially and polynomially tilted stable distributions. ACM Transactions on Modeling and
Computer Simulation 19, Article 18.
El-Shaarawi, A.H., Zhu, R., & Joe, H. (2011). Modelling species abundance using the Poisson-Tweedie family. Environmetrics 22, 152–164.
Gerber, H.U. (1991). From the generalized gamma to the generalized negative binomial distribution. Insurance: Mathematics and Economics 10,
303–309.
Hougaard, P. (1986). Survival models for heterogeneous populations derived from stable distributions. Biometrika 73, 387–396.
Hougaard, P., Lee, M.T., & Whitmore, G.A. (1997). Analysis of
overdispersed count data by mixtures of Poisson variables and Poisson
processes. Biometrics 53, 1225–1238.
Huillet, T.E. (2016). On Mittag-Leffler distributions and related stochastic processes. Journal of Computational and Applied Mathematics 296,
181–211.
Johnson, N.L., Kemp, A.W., & Kotz, S. (2005). Univariate Discrete
Distributions, 3rd edn. Wiley, New York.
Jose, K.K., Uma, P., Lekshmi, V.S., & Haubold, H.J. (2010).
Generalized Mittag-Leffler distributions and processes for applications
16

in astrophysics and time series modeling, in Proceedings of the Third
UN/ESA/NASA Workshop on the International Heliophysical Year 2007
and Basic Space Science, H.J. Haubold and A.M. Mathai (eds.), Springer,
New York.
Kanter, M. (1975). Stable densities under change of scale and total variation inequalities, Annals of Probability 3, 697–707.
Klebanov, L.B. & Slámová, L. (2015). Tempered distributions: does
universal tempering procedure exist?, arXiv:1505.02068v1[math.PR].
Linnik, Y.V. (1962). Linear forms and statistical criteria II. Translations in
Mathematical Statistics and Probability 3, American Mathematical Society, Providence, 41–90.
Marcheselli, M., Baccini, A., & Barabesi, L. (2008). Parameter estimation for the discrete stable family. Communications in Statistics Theory and Methods 37, 815–830.
Pakes, A.G. (1995). Characterization of discrete laws via mixed sums and
Markov branching processes. Stochastic Processes and their Applications
55, 285–300.
Rachev, S.T., Kim, Y.S., Bianchi, M.L., & Fabozzi, F.J. (2011).
Financial Models with Lévy Processes and Volatility Clustering, Wiley,
New York.
Rosı́nski, J. (2007). Tempering stable processes. Stochastic Processes and
Their Applications 117, 677–707.
Sato, K. (1999). Lévy Processes and Infinitely Divisible Distributions, Cambridge University Press, Cambridge.
xSibuya, M. (1979). Generalized hypergeometric, digamma and trigamma
distributions. Annals of the Institute of Statistical Mathematics 31, 373–
390.
Steutel, F.W., & van Harn, K. (1979). Discrete analogues of selfdecomposability and stability, Annals of Probability 7, 893–899.
Zhu, R., & Joe, H. (2009). Modelling heavy-tailed count data using a generalized Poisson-inverse Gaussian family. Statistics & Probability Letters
79, 1695–1703.

17

Zolotarev, V.M. (1986). One-dimensional Stable Distributions. Translations of Mathematical Monographs 65, American Mathematical Society,
Providence.

Lucio Barabesi, Dipartimento di Economia Politica e Statistica, P.zza S.Francesco
7, 53100 Siena, Italy.
E-mail: lucio.barabesi@unisi.it
Carolina Becatti, IMT School for advanced studies, Piazza S.Francesco, 19,
55100 Lucca.
E-mail: carolina.becatti@imtlucca.it
Marzia Marcheselli, Dipartimento di Economia Politica e Statistica, P.zza
S.Francesco 7, 53100 Siena, Italy.
E-mail: marzia.marcheselli@unisi.it

18

40

35

30

4.2

4.4

4.6

4.8

5.0

5.2

5.4

Figure 1: Parametric plot of (α3 , α4 ) for a = 1/4, b = 1 and c ∈ [3/10, 7/10],
d ∈ [−1, 3].
30

25

20

15

2.5

3.0

3.5

4.0

Figure 2: Parametric plot of (α3 , α4 ) for a = 1/2, b = 1 and c ∈ [3/10, 7/10],
d ∈ [−1, 3].
19

25

20

15

10

1.5

2.0

2.5

3.0

3.5

4.0

Figure 3: Parametric plot of (α3 , α4 ) for a = 3/4, b = 1 and c ∈ [3/10, 7/10],
d ∈ [−1, 3].

30

25

20

15

10

5
1

2

3

4

Figure 4: Parametric plot of (α3 , α4 ) for a = −1, b = 1 and c ∈ [1/10, 9/10],
d ∈ [0, 3].
20

