Arxiv:2008.06200V4 [Math.PR] 4 Jun 2021 Characterizing the Zeta

arXiv:2008.06200v4 [math.PR] 4 Jun 2021 admvralssmtmsaeepoe omdlheavy-tai model to employed are arti sometimes present variables the t random for derived basis authors the the form which 2021), distribution, Xu, Zeta and Powers, Huang, Dai, (see spr fa netgto fhaytie iceedistri discrete heavy-tailed of investigation an of part As Introduction 1 n cwrmnClee snhaUiest,Biig Chi Beijing, University, Tsinghua College, Schwarzman and [email protected]. * † ‡ § 1 iadun netetMngmn;eal daiball@yeah. email: Management; Investment WizardQuant ii aia aaeet mi:[email protected] email: Management; Capital Sixie orsodn uhr 8CWiu ulig eatetof Department Building, Weilun 386C author; Corresponding eateto iac,Sho fEooisadManagement, and Economics of School Finance, of Department y“ev-ald,w enarno variable random a mean we “heavy-tailed”, By rcal otnosmxueo oso itiuin.I distributions. r Poisson of mixture continuous tractable a eofrtonvlcaatrztoso h eadistributi Zeta sha fixed the (with distributions of Binomial Negative of characterizations mixtures novel two offer We Keywords: values). negative some assumes function density itrs identifiability. mixtures; itr a eascae ihauiu iigdistribution mixing unique a with for associated be can mixture inhn Dai, Jiansheng ∈ r [1 ∈ , ∞ (0 hrceiigteZt itiuinvia Distribution Zeta the Characterizing ) , 1) n h oso ae h eutn eadsrbtosaei are distributions Zeta resulting the case, Poisson the and eadsrbto;Ngtv ioildsrbto;Poiss distribution; Binomial Negative distribution; Zeta h iigdsrbtosaeqaidsrbtos(o wh (for quasi-distributions are distributions mixing the , * iegHuang, Ziheng otnosMixtures Continuous ue4 2021 4, June X Abstract † ∼ ihe .Powers, R. Michael f 1 X a108;eal [email protected]. email: 100084; na ( x ) o which for net. m. butions iac,Sho fEooisadManagement, and Economics of School Finance, ohteNgtv ioilcs for case Binomial Negative the both n l.Wti h nuac otx,Zeta context, insurance the Within cle. eparameter, pe n rt statbecontinuous tractable as first, on: nteNgtv ioilcase Binomial Negative the In . E e osfeunis(.. onsof counts (i.e., frequencies loss led 1 X nisrneadatailscience actuarial and insurance in onwcaatrztoso the of characterizations new wo [ X ndsrbto;continuous distribution; on α c h quasi-probability the ich etfibebcueeach because dentifiable ] snhaUiest;email: University; Tsinghua ∞ → ‡ > r n ixnXu Jiaxin and 0 o some for ,adscn,as second, and ), α ∈ (0 , § ∞ ) . event occurrences, claim submissions, indemnity payments, etc.; see, e.g., Doray and Arsenault, 2002). However, their use is even more common in other fields of study, in which they serve as empirical models for a wide range of discrete processes. Examples include: the number of word occurrences in a text; various measures of communication and influence, such as numbers of telephone calls, emails, and website hits; and physical intensities, such as numbers of earthquakes and solar flares within specified discrete categories (see, e.g., Newman, 2005). The Zeta distribution also plays important roles in analytic number theory, especially with regard to the distribution of prime numbers and the Riemann Hypothesis (see Lin and Hu, 2001; and Aoyama and Nakamura, 2012).2 The novel formulations of the Zeta distribution presented in this article are likely to be of interest to researchers in various fields (as indicated above), especially those exploring social or physical mechanisms leading to heavy-tailed behavior. The first characterization, provided in Sec- tion 2, shows that Zeta random variables can be expressed as continuous mixtures of Negative Binomial counts with a fixed shape parameter, r > 0. This is accomplished via tractable mixing distributions that are well behaved for r ∈ [1, ∞), but consist of quasi-distributions (for which the quasi-probability density function assumes some negative values) for r ∈ (0, 1). The second characterization, given in Section 3, converts the mixtures of Negative Binomial counts from Section 2 into a tractable continuous mixture of Poisson counts by first expressing each Negative Binomial component as a continuous mixture of Poisson components. In both the Negative Binomial case for r ∈ [1, ∞) and the Poisson case, the resulting Zeta distributions are identifiable because the mixing distributions are unique. In Section 4, we conclude with some final observations.

2 Zeta as a Mixture of Negative Binomial Counts

In this section, we will show how Zeta random variables can be constructed as continuous mixtures of Negative Binomial counts with a ﬁxed shape parameter, r. To this end, let X|s ∼ (x+1)−s Zeta (s) have probability mass function (PMF) fX|s (x) = ζ(s) for x ∈ {0, 1, 2,...} and ∞ −s s ∈ (1, ∞), where ζ (s) = x=0(x + 1) denotes the Riemann zeta function, and let X | r, p ∼ 2It is important to note that thisP literature uses the term “Riemann Zeta” to describe random variables formed by taking the negative of the natural logarithm of conventional “Zeta” random variables (deﬁned on the sample space x ∈{1, 2, 3,...}).

2 Γ(r+x) x r Negative Binomial (r, p) have PMF fX|r,p (x | r, p) = Γ(r)Γ(x+1) p (1 − p) for x ∈ {0, 1, 2,...}, ﬁxed r ∈ (0, ∞), and p ∈ (0, 1).3 Our basic objective is to identify a mixing random variable, p|r, s ∼ fp|r,s (p) for p ∈ (0, 1), such that

fX|s (x)= fX|r,p (x) ∧ fp|r,s (p) . (1) p

Since the form of the mixing probability density function (PDF), fp|r,s (p), differs by the domain of the Negative Binomial shape parameter, r, we must consider the cases of r =1, r ∈ (1, ∞), and r ∈ (0, 1), respectively, in the following three subsections.

2.1 TheCaseof r = 1 (Geometric Distribution)

If the shape parameter equals 1, then the NegativeBinomial (r, p) distribution simpliﬁes to

x Geometric (p), with PMF fX|r=1,p (x)= p (1 − p). Rewriting (1) as

x 1 −1 E [px]=1 − (i + 1)−s. (2) p|r=1,s ζ (s) i=0 X 3The Zeta (s) and Negative Binomial (r, p) distributions often are deﬁned on the sample space x ∈ {1, 2, 3,...} rather than x ∈{0, 1, 2,...}. However, we have chosen the latter characterization both because it matches the sample space of the Poisson (λ) distribution and because it is the more commonly used formulation in insurance applications (where it is convenient for loss frequencies to admit the possibility of x =0).

3 The system (2) then can be used to express the moment-generating function of p|r =1,s as

tE [p] t2E [p2] t3E [p3] E etp =1+ p|r=1,s + p|r=1,s + p|r=1,s + ··· p|r=1,s 1! 2! 3! t 1 t2 1+2−s t3 1+2−s +3−s =1+ 1 − + 1 − + 1 − + ··· (3) 1! ζ (s) 2! ζ (s) 3! ζ (s) n 1 ∞ tn −1 = et − (i + 1)−s ζ (s) n! n=1 i=0 ! X X ∞ tn H = et − n,s , (4) n! ζ (s) n=1 X n−1 −s th where Hn,s = i=0 (i + 1) is the n generalized harmonic number. Alternatively, (3) may be written as P 1 ∞ tn ∞ (i + 1)−s ζ (s) n! n=0 i=n ! X X ∞ tn ζ (s, n + 1) = , (5) n! ζ (s) n=0 X ∞ −s where ζ (s, n +1)= i=n (i + 1) is the Hurwitz zeta function of order n +1.

Hn,s ζ(s,n+1) The series (4) andP (5) clearly converge for all t ∈ (0, ∞) because both ζ(s) and ζ(s) are bounded above as n → ∞ for all s ∈ (1, ∞). Although these series are not reducible to simpler expressions, it is possible to write the PDF fp|r=1,s (p) analytically, as shown in the following result.

Theorem 1: If X|s ∼ Zeta (s) and X|r =1,p ∼ Negative Binomial (r =1,p), then there exists a unique mixing random variable, p|r =1,s, with PDF

(− ln (p))s−1 f (p)= p|r=1,s ζ (s)Γ(s)(1 − p) for p ∈ (0, 1), such that

fX|s (x)= fX|r=1,p (x) ∧ fp|r=1,s (p) . p

Proof: First, consider

1 1 (− ln (p))s−1 f (x) f (p) dp = px (1 − p) dp X|r=1,p p|r=1,s ζ (s)Γ(s)(1 − p) Z0 Z0

4 1 1 = px (− ln (p))s−1 dp. (6) ζ (s)Γ(s) Z0 Then, using the substitution y = − ln (p) in the above integral, (6) can be rewritten as

0 1 x e−y ys−1 −e−y dy ζ (s)Γ(s) Z∞ ∞ 1 x = e−y +1 ys−1dy ζ (s)Γ(s) Z0 (x + 1)−s ∞ (x + 1)s ys−1e−(x+1)y = dy ζ (s) Γ(s) Z0 (x + 1)−s = . ζ (s)

The uniqueness of fp|r=1,s (p) – and equivalently, the identifiability of fX|s (x) – follows from the identifiability of Negative Binomial mixtures with fixed shape parameter r (see Theorem 2.1 of Sapatinas, 1995).

One fairly obvious, yet interesting, aspect of the above result is that it provides a natural connec- tion between two of the simplest convergent series in mathematical analysis: the inﬁnite geometric series, γ S =1+ γ−1 + γ−2 + γ−3 + ... = , G γ − 1 with γ ∈ (1, ∞); and the zeta function,

−s −s −s SZ =1+2 +3 +4 ... = ζ (s) ,

th with s ∈ (1, ∞). Letting τG (n) and τZ (n) denote the n terms of these two series, respectively 1 (for n ∈{1, 2, 3,...}), and treating γ as a random variable deﬁned by the transformation γ = p , it follows from Theorem 1 that

(ln (γ))s−1 f (γ)= γ|s ζ (s)Γ(s) γ (γ − 1) for γ ∈ (1, ∞) and τ (n) 1 Z = pn−1 (1 − p) f (p) dp S p|r=1,s Z Z0

5 1 (− ln (p))s−1 = pn−1 (1 − p) dp ζ (s)Γ(s)(1 − p) Z0 1 1 n−1 1 (− ln (1/γ))s−1 = − 1 − γ−2dγ γ γ 1 Z∞ ζ (s)Γ(s) 1 − γ ∞ γ−(n−1) (ln (γ))s−1 = γ dγ 1 ζ (s)Γ(s) γ (γ − 1) Z γ − 1 ∞ τ (n) = G f (γ) dγ S γ|s Z1 G τ (n) = E G . γ|s S G

2.2 The Case of Fixed r ∈ (1, ∞)

When the ﬁxed Negative Binomial shape parameter differs from 1, the mixing PDF, fp|r,s (p), gen- erally becomes more complex, but remains reasonably tractable. The following result addresses the case of r ∈ (1, ∞).

Theorem 2: If X|s ∼ Zeta (s) and X|r > 1,p ∼ Negative Binomial (r > 1,p), then there exists a unique mixing random variable, p|r> 1,s, with PDF

r − 1 1 (ω − p)r−2 (− ln (ω))s−1 f (p)= dω p|r>1,s ζ (s)Γ(s)(1 − p)r ωr−1 Zp for p ∈ (0, 1), such that

fX|s (x)= fX|r>1,p (x) ∧ fp|r>1,s (p) . p

Proof: Consider 1 fX|r>1,p (x) fp|r>1,s (p) dp Z0 1 Γ(r + x) r − 1 1 (ω − p)r−2 (− ln (ω))s−1 = px (1 − p)r dω dp Γ(r)Γ(x + 1) ζ (s)Γ(s)(1 − p)r ωr−1 Z0 " Zp # 1 Γ(r + x) 1 px 1 (ω − p)r−2 (− ln (ω))s−1 = dω dp. (7) ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) ωr−1 Z0 "Zp # Using the substitution t = − ln (ω) in the inside integral, followed by y = − ln (p) in the outside

6 integral, (7) can be rewritten as

1 Γ(r + x) 0 (e−y)x 0 (e−t − e−y)r−2 ts−1 −e−t dt −e−y dy −t r−1 ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) y (e ) Z∞ "Z # ∞ y 1 Γ(r + x) 1 r = e−(x+1)y 1 − et−y −2 ts−1dtdy. (8) ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) Z0 Z0 Now interchange the order of integration, and let ξ = y − t in the new inside integral, so (8) becomes

∞ ∞ 1 Γ(r + x) 1 r−2 e−(x+1)(ξ+t) 1 − e−ξ ts−1dξdt ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) Z0 Z0 ∞ ∞ 1 Γ(r + x) 1 r−2 = ts−1e−(x+1)tdt e−(x+1)ξ 1 − e−ξ dξ ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) Z0 Z0 −s ∞ s ∞ (x + 1) Γ(r + x) (x + 1) r−2 = ts−1e−(x+1)tdt e−(x+1)ξ 1 − e−ξ dξ ζ (s) Γ(r − 1)Γ(x + 1) Γ(s) Z0 Z0 −s ∞ (x + 1) Γ(r + x) r−2 = e−(x+1)ξ 1 − e−ξ dξ. ζ (s) Γ(r − 1)Γ(x + 1) Z0 Substituting q = e−ξ into the above integral then yields

(x + 1)−s Γ(r + x) 0 1 qx+1 (1 − q)r−2 − dq ζ (s) Γ(r − 1)Γ(x + 1) q Z1

(x + 1)−s 1 Γ(x + r) = qx (1 − q)r−2 dq ζ (s) Γ(x +1)Γ(r − 1) Z0 (x + 1)−s = . ζ (s)

The uniqueness of fp|r>1,s (p) follows in the same way as the uniqueness of fp|r=1,s (p) in the proof of Theorem 1.

It is worth noting that, for the special case of r =2, the PDF fp|r>1,s (p) simpliﬁes to

1 1 (− ln (ω))s−1 fp|r=2,s (p)= 2 dω ζ (s)Γ(s)(1 − p) p ω Z

7 1 (− ln (ω))s 1 = − 2 s ζ (s)Γ(s)(1 − p) " p#

s (− ln (p)) = , ζ (s)Γ(s + 1)(1 − p)2 an analytic form quite similar to fp|r=1,s (p).

2.3 The Case of Fixed r ∈ (0, 1)

For ﬁxed r ∈ (0, 1), the analysis is similar to that for r ∈ (1, ∞), with one major difference: the mixing distribution has a quasi-PDF, fp|r<1,s (p), that assumes some negative values. We provide the details in the following result.

Theorem 3: If X|s ∼ Zeta (s) and X|r < 1,p ∼ Negative Binomial (r < 1,p), then there exists a mixing quasi-random variable, p|r< 1,s, with quasi-PDF

1 1 (ω − p)r−1 (− ln (ω))s−2 f (p)= (s − 1) dω p|r<1,s ζ (s)Γ(s)(1 − p)r ωr " Zp

1 (ω − p)r−1 (− ln (ω))s−1 +(r − 1) dω (9) ωr Zp # for p ∈ (0, 1), such that

fX|s (x)= fX|r<1,p (x) ∧ fp|r<1,s (p) . p

Proof: After writing 1 fX|r<1,p (x) fp|r<1,s (p) dp Z0 1 Γ(r + x) 1 1 (ω − p)r−1 (− ln (ω))s−2 = px (1 − p)r (s − 1) dω Γ(r)Γ(x + 1) ζ (s)Γ(s)(1 − p)r ωr Z0 " Zp 1 (ω − p)r−1 (− ln (ω))s−1 +(r − 1) dω dp ωr Zp # 1 Γ(r + x) 1 1 1 (ω − p)r−1 (− ln (ω))s−2 = px dω ζ (s) Γ(r)Γ(x + 1) Γ(s − 1) ωr Z0 " Zp

8 (r − 1) 1 (ω − p)r−1 (− ln (ω))s−1 + dω dp, (10) Γ(s) ωr Zp # (x+1)−s (x+1)−s the proof that (10) equals ζ(s) is entirely analogous to the proof that (7) equals ζ(s) for Theorem 2.

To demonstrate fp|r<1,s (p) < 0 for some p ∈ (0, 1) for all s ∈ (1, ∞), it sufﬁces to insert p =0 into the right-hand side of (9), revealing lim fp|r<1,s (p)= −∞. This implies there exists an p→0+ interval (0,ε), for some ε> 0, such that p ∈ (0,ε)=⇒ fp|r<1,s (p) < 0.

The reason fp|r<1,s (p) < 0 in some neighborhood of zero is quite intuitive. Essentially, when r ∈ (0, 1), the Negative Binomial PMF becomes very steep at the lower end of its sample space

fX|r,p(0) 1 (in the sense that sup f (1) = lim rp = ∞), and this steepness is aggravated for values r X|r,p r→0+ of p close to zero. In this region of the sample space, it is impossible to construct the much ﬂatter

fX|s(0) s Zeta PMF (for which inf f = lim 2 = 2) as a convex combination of Negative Binomial s X|s(1) s→1+ PMFs. However, if one can assign negative weight to those Negative Binomial PMFs for which p is very small, then one can mitigate the impact of small r by offsetting it with negative contributions from small p.

2.4 The Case of Random r

If the Negative Binomial shape parameter, r, is not ﬁxed, but rather is a continuous random variable on (0, ∞), then it is possible to express X|s ∼ Zeta (s) as a mixture

fX|s (x)= fX|r,p (x) ∧ fr,p|s (r, p) , r,p

for some joint mixing PDF, fr,p|s (r, p). Unfortunately, this joint PDF is not unique, and the resulting Zeta distribution therefore is not identiﬁable. In fact, if fr,p|s (r, p) = fp|r,s (p) fr|s (r), where fp|r,s (p) denotes the mixing PDF given by Theorems 1 or 2 above, then

∞ 1 fX|s (x)= fX|r,p (x) fr,p|s (r, p) drdp Z0 Z0

∞ 1 = fX|r,p (x) fp|r,s (p) fr|s (r) drdp Z0 Z0

9 ∞ 1 = fX|r,p (x) fp|r,s (p) dp fr|s (r) dr Z0 Z0 ∞ = fX|s (x) fr|s (r) dr, Z0 so that fr|s (r) can be any well-deﬁned PDF on (0, ∞). Without the identiﬁability property, it is impossible to probe the statistical processes generating the mixed distribution of interest (in our case, fX|s (x)). In particular, one cannot estimate the parameters of the mixing distribution from observations of the mixed random variable (see, e.g., Karlis and Xekalaki, 2005).

3 Zeta as a Mixture of Poisson Counts

For any choice of ﬁxed r ∈ (0, ∞) and p ∈ (0, 1), the Negative Binomial (r, p) random variable can 1−p be expressed as a continuous mixture of Poisson (λ) counts, using a Gamma r, β = p mixing r 1−p λr−1 1−p distribution with PDF f 1−p (λ)= exp − λ for λ ∈ (0, ∞), r ∈ (1, ∞), λ r,β p Γ(r) p | = p 1−p and β = p ∈ (0, ∞). This mixture may be written as

fX|r,p (x)= fX|λ (x) ∧ f 1−p (λ) . λ r,β λ | = p

Substituting the right-hand side of the above equation into (1) yields

fX|s (x)= fX|λ (x) ∧ f 1−p (λ) ∧ fp|r,s (p) , λ λ|r,β= p p which, by the associative property of distribution mixing, is equivalent to

fX|s (x)= fX|λ (x) ∧ f 1−p (λ) ∧ fp|r,s (p) λ λ|r,β= p p

= fX|λ (x) ∧ fλ|r,s (λ) . λ

It thus follows that Zeta random variables can be expressed as continuous mixtures of Poisson counts with the mixing PDF,

fλ|r,s (λ)= f 1−p (λ) ∧ fp|r,s (p) λ r,β | = p p

10 1 − p r λr−1 1 p 1 − p = exp − λ f (p) dp. (11) Γ(r) p p|r,s Z0 Although the right-hand side of (11) appears to depend on the Negative Binomial shape parameter (r), that actually is not the case. As shown by Feller (1943), any mixed-Poisson random variable must possess the identiﬁability property, and thus a unique mixing distribution. Conse- quently, fλ|r,s (λ) must be invariant over r, and may be expressed using the simplest form of the above integral (i.e., by inserting r =1). Therefore,

fλ|r,s (λ)= fλ|s (λ)

1 1 − p 1 − p (− ln (p))s−1 = exp − λ dp p p ζ (s)Γ(s)(1 − p) Z0 1 1 1 1 − p = (− ln (p))s−1 exp − λ dp. (12) ζ (s)Γ(s) p p Z0 1−p Applying the substitution y = p to the above integral, (12) can be rewritten as

1 0 −1 (y + 1) (ln (y + 1))s−1 exp (−λy) dy ζ (s)Γ(s) (y + 1)2 Z∞ 1 ∞ 1 = (ln (y + 1))s−1 exp (−λy) dy, ζ (s)Γ(s) y +1 Z0 which is not further reducible.

4 Conclusions

The present article provided two novel characterizations of X|s ∼ Zeta (s). First, we showed that these random variables can be expressed as tractable continuous mixtures of X|r, p ∼ Negative Binomial (r, p) with ﬁxed shape parameter (r); that is,

fX|s (x)= fX|r,p (x) ∧ fp|r,s (p) , p

11 where

1 (ω − p)r−1 (− ln (ω))s−2 (s − 1) dω ωr  Zp  1 (ω − p)r−1 (− ln (ω))s−1  +(r − 1) dω for r ∈ (0, 1) 1  r  p ω fp|r,s (p)= r × Z , ζ (s)Γ(s)(1 − p)  s−1 (− ln (p)) for r =1 1 (ω − p)r−2 (− ln (ω))s−1 (r − 1) dω for r ∈ (1, ∞)  r−1  p ω  Z   fp|r,s (p)= fp|r=1,s (p) and fp|r,s (p)= fp|r>1,s (p) are unique PDFs, and fp|r,s (p)= fp|r<1,s (p) is a quasi-PDF (i.e., with some negative values). Next, based on the fact that Negative Binomial (r, p) 1−p random variables can be constructed as mixtures of X|λ ∼ Poisson (λ), witha Gamma r, β = p mixing distribution, we showed that Zeta random variables also can be expressed as a unique and tractable continuous mixture of Poisson counts; that is,

fX|s (x)= fX|λ (x) ∧ fλ|s (λ) , λ where 1 ∞ 1 f (λ)= (ln (y + 1))s−1 exp (−λy) dy. λ|s ζ (s)Γ(s) y +1 Z0 The appearance of quasi-PDFs in the Negative Binomial case for r ∈ (0, 1) was somewhat unexpected, but – as argued in Subsection 2.3 – has a fairly intuitive explanation. Therefore, it is natural to consider whether or not other heavy-tailed discrete random variables formed as mixtures of Negative Binomial counts also involve quasi-distributions. One obvious family to bΓ(b+1)Γ(x+1) consider is X|b ∼ Yule (b), with PMF fX|b (x) = Γ(x+b+2) for x ∈ {0, 1, 2,...} and b ∈ (0, ∞), which approximates X|s ∼ Zeta (s) for b = s − 1. In fact, Yule (b) random variables may be expressed as mixtures of Negative Binomial (r =1,p) (i.e., Geometric (p)) counts, using fp|a=1,b (p) ∼ Beta (a =1, b) as the mixing distribution; that is,

fX|b (x)= fX|r=1,p (x) ∧ fp|a=1,b (p) . p

fX|b(0) For the Yule distribution, inf f = lim (b +2) = 2, which is identical to the corre- b X|b(1) b→0+

12 sponding result for the Zeta distribution presented at the end of Subsection 2.3. Consequently, one might anticipate that constructing Yule random variables as mixtures of Negative Binomial counts with r ∈ (0, 1) would require quasi-distributions, as in the Zeta case. In Dai, Huang, Powers, and Xu (2021), we not only show that this is indeed true, but also provide complete results for the Yule distribution analogous to those of Theorems 1-3 and Section 3 of the present article.

References

[1] Aoyama, T. and Nakamura, T., 2012, “Zeros of Zeta Functions and Zeta Distributions on Rd”, RIMS Kokyˆ urokuˆ Bessatsu, 39-48.

[2] Dai, J., Huang, Z., Powers, M.R., and Xu, J., 2021, “Heavy-Tailed Loss Frequencies from Mix- tures of Negative Binomial and Poisson Counts”, working paper, Tsinghua University School of Economics and Management.

[3] Doray, L. G. and Arsenault, M., 2002, “Estimators of the Regression Parameters of the Zeta Distribution”, Insurance: Mathematics and Economics, 30, 3, 439-450.

[4] Feller, W., 1943, “On a Generalized Class of Contagious Distributions”, Annals of Mathemati- cal Statistics, 14, 389-400.

[5] Karlis, D. and Xekalaki, E., 2005, “Mixed Poisson Distributions”, International Statistical Review, 73, 1, 35-58.

[6] Lin, G. D. and Hu, C.-Y., 2001, “The Riemann Zeta Distribution”, Bernoulli, 7, 5, 817-828.

[7] Newman, M. E. J., 2005, “Power Laws, Pareto Distributions, and Zipf’s Law”, Contemporary Physics, 46, 5, 323-351.

[8] Sapatinas, T., 1995, Identiﬁability of Mixtures of Power-Series Distributions and Related Char- acterizations”, Annals of the Institute of Statistical Mathematics, 47, 447–459.