arxiv: v1 [stat.me] 29 Nov 2017

Transcription

1 Fractional approaches for the distribution of innovation sequence of INAR() processes arxiv:70645v [statme] 29 Nov 207 Josemar Rodrigues, Marcelo Bourguignon 2, Manoel Santos-Neto 3,4 N Balakrishnan 5 Departamento de Estatística, Universidade de São Paulo, São Carlos, Brazil 2 Departamento de Estatística, Universidade Federal do Rio Gre do Norte, Natal, Brazil 3 Departamento de Estatística, Universidade Federal de Campina Gre, Campina Gre, Brasil 4 Departamento de Estatística, Universidade Federal de São Carlos, São Carlos, Brasil 5 Department of Mathematics Statistics, McMaster University, Ontario, Canada L8S4KI Abstract In this paper, we present a fractional decomposition of the probability generating function of the innovation process of the first-order non-negative integer-valued autoregressive [INAR()] process to obtain the corresponding probability mass function We also provide a comprehensive review of integer-valued time series models, based on the concept of thinning operators, with geometric-type marginals In particular, we develop four fractional approaches to obtain the distribution of innovation processes of the INAR() model show that the distribution of the innovations sequence has geometric-type distribution These approaches are discussed in detail illustrated through a few examples Finally, using the methods presented here, we develop four new first-order non-negative integer-valued autoregressive process for autocorrelated counts with overdispersion with known marginals, derive some properties of these models Keywords: Time series Introduction Count data, Fractional decomposition, Geometric-type distribution, Innovation process, Models for count time series with overdispersion, based on thinning operators, have been discussed by many authors The vast majority of the proposed integer-valued autoregressive processes are based on geometrictype distribution There is an exhaustive list of INAR() processes with geometric-type marginal distributions, geometric-type count variables or geometric-type innovations McKenzie (986) Al-Osh & Aly (992) proposed INAR() processes with geometric negative binomial distributions as marginals Similarly, Alzaid & Al-Osh (988) discussed the geometric INAR() process [GINAR()] Ristić et al (2009) introduced the geometric first-order integer-valued autoregressive [NG- INAR()] model with geometric marginal distribution Integer-valued time series, with geometric marginal distribution, generated by mixtures of binomial negative binomial thinning operators have been considered by Nastić & Ristić (202), Nastić et al (202) Ristić & Nastić (202) Recently, Nastić et al (206a) constructed a new stationary time series model with geometric marginals, based on thinning operator, which is a mixture of Bernoulli geometric distributed rom variables Of course, a simple approach is to change the distribution of innovations In this context, Jazi et al (202) proposed the geometric INAR() model with geometric innovations, while Bourguignon (208) introduced a first-order non-negative integer-valued autoregressive process with zero-modified geometric innovations based on binomial thinning Corresponding author: Josemar Rodrigues josemar@icmcuspbr

2 In context of dependent Bernoulli counting variables, Ristić et al (203), Miletić Ilić (206) Nastić et al (207) introduced an integer-valued time series model with geometric marginals based on dependent Bernoulli count variables Recently, Miletić Ilić et al (207) introduced an INAR() model based on a mixed dependent independent count series with geometric marginals Nastić et al (206b) introduced anr-states rom environment non-stationary INAR() which, by its different values, represents the marginals selection mechanism from a family of different geometric distributions Borges et al (206) Borges et al (207) introduced geometric first-order integer-valued autoregressive processes with geometric marginal distribution based on ρ-binomial thinning operator ρ-geometric thinning operator, respectively Other works that have recently appeared in the literature dealing with geometric-type INAR() processes are those of Nastić (202) (shifted geometric INAR() process), Barreto-Souza (205) (INAR() process with zeromodified geometric marginals) Yang et al (206) (threshold INAR() process) Final mention should be made to the work of Ristić et al (202), Popović et al (206) Popović (206) proposed bivariate INAR() time series models with geometric marginals In many recent works, it is not clear how the distribution of the innovation processes was derived In this context, motivated by Feller (see 2008, p 276), we formulate here some procedures to obtain the probability mass function of the innovation process of an INAR() process by rewriting the probability generating function (pgf) of the innovation process as a quotient of two polynomial functions of enequal degrees In particular, we present four fractional approaches to obtain the probability mass function of innovation processes of the INAR() model when the distribution of the innovations sequence has geometric-type distribution Furthermore, we put forward four new first-order non-negative integer-valued autoregressive processes (Examples 6, 7, 8 9) with inflated-parameter Bernoulli inflated-parameter geometric marginals (Kolev et al, 2000) The rest of this paper is organized as follows In Section 2, we develop four fractional approaches to obtain the distribution of innovation processes of the INAR() model show that the distribution of the innovations sequence has geometric-type distribution In addition, some illustrative examples new models are presented Finally, some concluding remarks are made in Section 3 2 Main results Let Z, N R denote the set of integers, positive integers real numbers, respectively All rom variables will be defined on a common probability space (Ω, A, P) A discrete-time stationary non-negative integervalued stochastic process {X t } t Z is said to be a first-order integer-valued autoregressive [INAR()] process if it satisfies the equation X t = α X t +ǫ t = X t j= N j +ǫ t, t Z, where α is a thinning operator, α (0,), {ǫ t } t Z is an innovation sequence of independent identically distributed non-negative integer-valued rom variables, not depending on past values of{x t } t Z, mean µ ǫ (< ) variance σ 2 ǫ(< ) It is also assumed that the N j variables that define α X t are independent of the variables from which other values of the series are calculated, are such that E(N j ) = α Moreover, we assume that all N j variables defining the thinning operations are independent of the innovation sequence {ǫ t } t Z The autocorrelation function of {X t } t Z is of the same form as in the case of the usual AR() processes The stationary marginal distribution of{x t } t Z can be determined from the equation ϕ X (s) = ϕ X (ϕ N (s)) ϕ ǫ (s), () 2

3 where ϕ X (s) := E(s X ), ϕ N (s) := E(s N ) ϕ ǫ (s) := E(s ǫ ) denote the pgf s of {X t } t Z, {N j } X t j= {ǫ t } t Z, respectively Equation () can be used to obtain the distribution of the innovations sequence if the marginal distribution of the observable INAR() process is known By deriving the pgf s in (), using the stationarity of {X t } t Z, it is easily shown that the pgf of {ǫ t } t Z is given by ϕ ǫ (s) = ϕ X(s) ϕ X (ϕ N (s)) Next, we present some methods to obtain the probability mass function of the innovation sequence of an INAR() process when we can rewrite the pgf of the innovation process as a quotient of two polynomial functions of enequal degrees 2 Method : Fractional decomposition of the innovation distribution As in Feller (2008), let us assume that the pgf of the innovations process is a rational function given by ϕ ǫ (s) = U p(s) V q (s), where U p (s) V q (s) are two polynomials of degrees p q (p < q), respectively, that the equation V q (s) = 0 has q distinct roots s,s 2,,s q with s < min{s,s 2,,s q } Then, we readily have that V q (s) = q i= (s s i) The probability mass function (pmf) of the innovation process can then be expressed (see Feller, 2008, p 276) as q where ρ i = Up(s i) dvq(s) Pr[ǫ t = m] = i= ρ i s m+ i, (2) If s is smaller in absolute value than all other roots, the above pmf can be approxi- ds s=si as m Assuming that p = q + r, r 0, using some mated (see Feller, 2008, p 277) by ρ s m+ well-known algebraic manipulations, it is possible to rewrite the pgf of the innovation process (see Feller, 2008, p 277) as ϕ ǫ (s) = U r (s)+ U k(s), k < q (3) V q (s) Then, the results above can be applied to the rational function U k (s)/v q (s) Remark 2 Note that taking s i = ( + µ i )/u i in (2), the innovation process is a mixture of geometric distributions with parameters µ i /(+µ i ) expressed as Pr[ǫ t = m] = q i= ( ) m ( ) µi c i, (4) +µ i +µ i where c i = ρ i µ i, i =,,q We can find many different INAR() models proposed in the literature whose innovation processes are given in (4) These models will be called fractional integer-valued autoregressive process denoted here by FINAR() processes Remark 22 Feller (2008) has mentioned in p 276 that it is a hard work to find the exact mixture distribution in (4) a simple approximation could be of practical interest in order to get satisfactory solutions for inferential problems In fact, if we suppose that µ is smaller than all the others, then as m increases, we see that Pr[ǫ t = m] can be approximated by the geometric distribution with parameter µ /(+µ ) 3

4 Given the marginal pmf of INAR(0,) process, the next remark presents in detail the proof of a recursive formula to obtain the innovation pmf This alternative recursive procedure could be of interest if the pmf of the thinning operator is available Remark 23 (Alternative approach) The innovation pmf can be written by the recursive formula Proof Suppose Pr[ǫ t = l] = l b 0 i=l q c i b l i, l = q +,q +2, ϕ X (s) = U p (s) = p q a i s i ϕ α X (s) = V p (s) = b j s j with p < q, i=0 j=0 where a i s b i s are some coefficients such that a 0 +a + +a p = b 0 +b + +b q = Then, we readily have q p ϕ ǫ (s) b j s j = a i s i, s (5) j=0 i=0 Now, let ϕ ǫ (s) = c j s l, where c l s are some coefficients such that c l = Then, upon substituting this in (5), we obtain l=0 l=0 q c j s l b j s j = l=0 j=0 p a i s i, s (6) i=0 Upon comparing coefficients ofs 0,s,,s p on both sides of (6), we obtain the following equations: c 0 b 0 = a 0, c 0 b +c b 0 = a, c 0 b 2 +c b +c 2 b 0 = a 2, c 0 b p +c b p +c 2 b p 2 + +c p b 0 = a p, which can be expressed equivalently as b c 0 a 0 b b c a b 2 b b 0 0 c 2 = a 2 b p b p b p 2 b 0 c p a p 4

5 This gives c 0 c c 2 c p = B p a 0 a a 2 a p b b b 0 0 0, where B p = b 2 b b 0 0, b p b p b p 2 b 0 which is clearly invertible since b 0 0(b 0 = V q (0) > 0) Next, upon comparing coefficients of s p+,s p+2,,s q on both sides of (6), we get the following equations: c 0 b p+ +c b p +c 2 b p + +c p+ b 0 = 0, c 0 b p+2 +c b p+ +c 2 b p + +c p+2 b 0 = 0, c 0 b q +c b q +c 2 b q 2 + +c q b 0 = 0 This system of equations can be readily rewritten as follows c p+ b 0 = (c 0 b p+ +c b p +c 2 b p + +c p b ), c p+ b +c p+2 b 0 = (c 0 b p+2 +c b p+ +c 2 b p + +c p b 2 ), c p+ b q p + +c q b 0 = (c 0 b q +c b q +c 2 b q 2 + +c p b q p ) which can be expressed equivalently as b b b b 2 b b 0 0 b q p b q p 2 b q p 3 b 0 This readily yields c p+ c p+2 c p+3 c q = B q p B Bp a 0 a a 2 a p c p+ c p+2 c p+3 c q b p+ b p b p b c 0 b p+2 b p+ b p b 2 c = b p+3 b p+2 b p+ b 3 c 2 b q b q b q 2 b q p c p b p+ b p b p b b p+2 b p+ b p b 2, where B = b p+3 b p+2 b p+ b 3 b q b q b q 2 b q p 5

6 Next, upon comparing coefficients of s q+,s q+2,, on both sides of (6), we get the following equations c b q +c 2 b q +c 3 b q 2 + +c q+ b 0 = 0, (s q+ ) c 2 b q +c 3 b q +c 4 b q 2 + +c q+2 b 0 = 0, (s q+2 ) c 3 b q +c 4 b q +c 5 b q 2 + +c q+3 b 0 = 0, (s q+3 ) which readily yields the following solutions: c q+ = q c i b q i+, b 0 i= c q+2 = q+ c i b q i+2, b 0 i=2 c q+3 = q+2 c i b q i+3, b 0 i=3 in general, we have c l = l b 0 i=l q c i b l i, l = q +,q +2, Example (Geometric INAR() process) As in McKenzie (985), let us consider the marginal distribution of {X t } t Z to be a geometric distribution with parameter θ (0,) with pmf Pr[X t = m] = ( θ) m θ, m = 0,,2, The GINAR() process is based on the binomial thinning operator (Steutel Van Harn 979), α X t := Xt j= Y j, where the so-called counting series{y j } j is a sequence of independent identically distributed Bernoulli rom variables with Pr(Y j = ) = Pr(Y j = 0) = α [0,), satisfies the equation X t = α X t +ǫ t, t Z Thus, the pgf of the innovation processe can be written as ϕ ǫ (s) = ϕ X t (s) ϕ α Xt (s) = U(s), s < /( θ), V(s) where U(s) = θ +( θ)α ( θ)αs V(s) = ( θ)s Since we have here two linear functions with p = q =, upon using a simple algebraic manipulation 6

7 suggested in (3), the pgf of the innovation process can be rewritten as ϕ ǫ (s) = U(s) V(s) = α+( α) θ ( θ)s In this example, we have shown that the GINAR() process is a FINAR() process with the polynomial functions U(s) V(s) having the same degrees p = q = By using the rational decomposition in (2) for the fractional function θ( α) ( θ)s, it can be seen that the innovation process follows a zero-inflated geometric distribution given by Pr[ǫ t = m] = α {0} (m)+( α)θ( θ) m, m = 0,,2, Example 2 (New geometric INAR() process) In this example, we consider the stationary NGINAR() process with negative binomial thinning operator, introduced by Ristić et al (2009), given by X t = α X t +ǫ t, t Z, where α X t = X t i= W i, α [0,), {W i } t Z ( is a ) sequence of iid rom variables with geometric m distribution, ie, with pmf given by Pr[W i = m] = α +α +α, {X t} t Z is a stationary process with geometric(µ/( +µ)) marginals with pmf given by Pr[X t = m] = ( µ +µ )m +µ,µ > 0 Using these assumptions of the NGINAR() process, the pgf of the innovation sequence is given by ϕ ǫ (s) = U(s), s <, V(s) ( where U(s) = +α(+µ) α(+µ)s V(s) = αµ s +u µ ) (s ) +α α Since, NGINAR() process is a fractional INAR() process with p = q = 2, the goal is to apply the fraction decomposition described before to obtain the pmf of the innovation processe {ǫ t } t Z The decomposition proceeds as follows: i) The roots of the quadratic function V(s) are s = +µ µ s 2 = +α α ; ii) The coefficients ρ ρ 2 are ρ = U(s ) V (s ) = µ ( αµ ) µ α ρ 2 = U(s 2) V (s 2 ) = µ µ α Thus, the pmf of the innovation sequence is obtained from (2) is given by ( Pr[ǫ t = m] = αµ )( ) µ m µ α +µ +µ + αµ ( ) α m µ α +α +α, m = 0,,2, The pmf of the innovation sequence is well defined with the condition α [0,µ/( + µ)] for µ > 0 being necessary to guarantee that all probabilities are non-negative, is undefined for values outside this range Note that this result was obtained earlier in Ristić et al (2009) 7

8 Example 3 (Dependent counting INAR() process) Ristić et al (203) proposed a geometrically distributed time series generated by dependent Bernoulli count series, called DCINAR() process The DCINAR() process is based on new generalized binomial thinning operator α θ satisfies the following equation: X t = α θ X t +ǫ t, t Z, where the operator α θ is defined as α θ X = X t i= U i, α,θ [0,], {U i } i N is a sequence of Bernoulli(α) variables with E[U i ] = α Var[U i ] = α( α) Moreover, these rom variables are dependent, since Corr[U i,u j ] = θ 2, θ 0 i j For more details, see Ristić et al (203) The pgf of the innovation sequence can be expressed into partial fractions as where ρ = A µ, ρ 2 = ϕ ǫ (s) = α( θ)(α+θ αθ) α+θ 2αθ + ρ s s + ρ 2, s <, (7) s 2 s A 2 (α+θ 2αθ)µ s = +µ µ, s 2 = +(α+θ 2αθ)µ α+θ 2αθ)µ The expressions for A A 2 are given in Ristić et al (203) Taking µ = µ µ 2 = (α +θ 2αθ)µ, the mixture innovation distribution in Ristić et al (203) follows from (4) (7) as ( ) m ( ) µ m µ2 Pr[ǫ t = m] = A 0 {0} (m)+a +A 2, (8) +µ +µ +µ 2 +µ 2 where A 0 = (A +A 2 ) = α( θ)(α+θ αθ),a 0 +A +A 2 = A i 0, i = 0,,2 α+θ 2αθ We note that the probability distribution in (8) is a zero-inflated model involving two geometric distributions This result is not cleary stated in Theorem 2 of Ristić et al (203) 22 Method 2: Linear rational probability generation function In this section, we extend the innovation sequence of the geometric INAR() process discussed in Example We consider here the a linear fractional probability generating function with four-parameters given by ϕ ǫ (s) = a+bs c+ds, s [0,] ϕ ǫ(0) <, (9) with the following parametric restrictions: a < c a + b = c + d After a simple algebraic manipulation, the pgf in (9) can be decomposed into partial fractions as ϕ ǫ (s) = b d + ρ s s, where s = c/d ρ = bc/d 2 a/d Following the same procedure as in Feller (2008, p 276), we obtain an exact expression for the pmf of the innovation sequence as Pr[ǫ t = m] = b d {0}(m)+ ρ 8 s m+, m = 0,,2,

9 Taking s = /( θ), 0 θ <, we have the pmf of the three-parameter innovation sequence Pr[ǫ t = m] = b ( d {0}(m)+ b ) ( θ) m θ, m = 0,,2, (0) d Note that the two-parameter innovation distribution in Example is deduced by taking b = ( θ)α d = ( θ) The mean the variance of the innovation process in (0) are given by ( )( θ E[ǫ t ] = b ), θ d ( )( θ Var[ǫ t ] = b )[ ( + b ) b ], θ d θ d d respectively The Fisher index of the dispersion of innovation process is given by I ǫ = Var[ǫ t] = ( + b ) b E[ǫ t ] θ d d This Fisher index indicates equidispersion of the innovation process if b = d, underdispersion if b > d overdispersion if b < d Note that for b = d, we have the mean to equal the variance, but the pmf in (0) is not the classical Poisson distribution Example 4 (Two-parameter innovation sequence) Let us now consider the following two-parameter pgf discussed by Aly & Bouzar (207) ϕ ǫ (s) = m s +r( s), 0 s, where r 0 m r + It can be easily seen that the above pgf is a fraction linear function with a = r m, b = m r, c = +r d = r For this particular case we have s = (r +)/r b/d = m/r The corresponding two-parameter innovation process is given by ( Pr[ǫ t = m] = m ) ( m ) ( )( ) r m r {0} (m)+, m = 0,,2, r +r +r Example 5 (Zero-modified geometric innovation process) Johnson et al (2005) have formulated in Section 823 the following pgf ϕ ǫ (s) = k+k, s <, () +µ( s) with the parametric restriction /µ k It is straightforward to see that the pgf in () has a linear representation given by a = +kµ, b = kµ, c = +µ d = µ Note thats = c/d = (+µ)/µ the parametric restriction 0 a c implies that /µ k Also, we have b/d = k in view of (0), withθ = /(+µ), that the pmf of {ǫ t } t Z is given by ( )( ) µ m Pr[ǫ t = m] = k {0} (m)+( k) m = 0,, +µ +µ 9

10 This pmf is the so-called zero-modified geometric distribution which presents equidispersion when k = ; underdispersion when µ (0, ) k [ /µ, ); overdispersion when k (, ) For 0 < k <, it is known as zero-inflated distribution zero-deflated distribution for /µ < k < 0 The zero-modified geometric distribution was considered in a recent paper by Bourguignon (208) to formulate a new INAR() process with zero-modified geometric innovations Also, if we take k = π/µ, 0 π, the zero-modified geometric distribution is the Bernoulli-Geometric distribution with parameters µ π introduced by Bourguignon & Weiß (207) This can be easily seen by taking the linear representation (see Bourguignon & Weiß, 207) a = π, b = π, c = +µ d = µ 23 Method 3: Quadratic rational probability generation function In this section, we consider the pgf of the innovation process as a quadratic rational function expressed as follows: ϕ ǫ (s) = as2 +bs+c as 2 + bs+ c, s [0,] ϕ ǫ(0) <, (2) where a 0, b+c = b+ c c c Initially, to facilitate the partial fraction decomposition, we have to reduce the degree of the polynomial function in the numerator of (2) by using a simple algebraic manipulation Thus, we obtain a new expression for the pgf as ϕ ǫ (s) = + (b b)s+c c as 2 + bs+ c = + U (s) (s s )(s s 2 ), (3) with U (s) = (b b)(s ), a where s s 2 are two different roots of the polynomial function of the denominator in (2) such that as 2 + bs+ c = a(s s )(s s 2 ) In order to obtain the innovation distribution, from here on, we assume that s 2 s > The rational U function (s) (s s )(s s 2 ) in (3) can now be decomposed into partial fractions (see Feller, 2008) the pgf of the innovation variable can be rewritten as where ϕ ǫ (s) = + ρ s s + ρ 2 s 2 s, ρ = U (s ) s s 2 ρ 2 = U (s 2 ) s 2 s Using the same geometric expansion for s s sequence is obtained as s 2 s as in Feller (2008), the pmf of the innovation { π, if m = 0, Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m 2 ], if m, (4) where p = /s, p 2 = /s 2 (p 2 p ) π = c/ c Note that w ( p )p m + w 2 ( p 2 )p m 2 is a 0

11 mixture of two geometric distributions, wherew = p /(p p 2 ) w 2 = p 2 /(p 2 p ) are the weights with w +w 2 = If π p 2 0, then ǫ t Geo( p ) if p 2 0, then ǫ t has a hurdle geometric distribution with parameter p Also, if π p 0, then ǫ t Geo( p 2 ) if p 0, ǫ t has a hurdle geometric distribution with parameter p 2 Since in (4) the zeros come from the Bernoulli process with parameter π the nonzeros come from a different process characterized by a mixture of two geometric distributions, it may be referred to as hurdlefractional geometric innovation distribution denoted by HFG The mean the variance of the innovation process in (4) are given by respectively E[ǫ t ] = ( π) ( p )( p 2 ), Var[ǫ t ] = ( π)[ p 2 p 2 3p p 2 +p +p 2 +π ] ( p ) 2 ( p 2 ) 2, Example 6 (Inflated-parameter geometric INAR() process based on binomial thinning operator) In this example, we present a new stationary first-order non-negative integer valued autoregressive process, {X t } t Z, with inflated-parameter geometric (Kolev et al, 2000) marginals The proposed process is based on the binomial thinning operator satisfies the equation: X t = α X t +ǫ t, t Z, withα [0,) {X t } t Z being a stationary process with inflated-parameter geometric marginals, ie, with pmf given by Pr[X t = m] = ρ ( µ+ρ {0}(m)+ ρ )( ) µ+ρ m ( ) ρ, m = 0,,2,, µ+ρ +µ +µ where ρ [0, ) µ > 0 The inflated-parameter geometric distribution is well-known as ρ-geometric distribution (see Kolev et al, 2000) considered in a recent paper by Borges et al (207) to formulate a new thinning operator The corresponding pgf is given by ϕ X (s) = ρs, s < sρ+µ( s) The mean variance of the process {X t } t Z are given, respectively, by E[X t ] = µ ρ Var[X t ] = µ(+µ+ρ) ( ρ) 2, so the dispersion index, denoted by I X, becomes I X = +α+ρ ρ = + 2ρ+µ ρ +µ >, that is, the proposed inflated-parameter geometric INAR() process accommodates overdispersion If ρ = 0, then the process{x t } t Z reduces to the NGINAR() process introduced by Ristić et al (2009); see Example 2 The additional parameter ρ has a natural interpretation in terms of both the zero-inflated proportion the correlation coefficient For more details about the inflated-parameter geometric distribution, see Kolev et al

12 (2000) The pgf of the innovation sequence is given by ϕ ǫ (s) = ϕ X(s) ϕ α X (s) = ( ρs){ ρ[ α( s)]+αµ( s)}, s < (5) [ sρ+µ( s)]{ ρ[ α( s)]} The main goal here is to obtain the innovation sequence by applying the fractional decomposition method for the quadratic rational function in (5) After a simple algebraic manipulation, the pgf in (5) can be rewritten as a quotient of two quadratic polynomial functions, where a = αρ(µ+ρ),b = [ρ( ρ)+α(µ+ρ)(+ρ)],c = ρ+α(µ+ρ), b = {(µ+ρ)[ ρ( α)]+ρα(+µ)} with ρ+α(µ+ρ) = c c = (+µ)[ ρ( α)] π = ρ+α(µ+ρ) function in the denominator of the pgf in (5) are given by c = (+µ)[ ρ( α)], (+µ)[ ρ( α)] The roots of the quadratic s = +µ ρ+µ s 2 = ρ( α), ρα withs 2 s >, consequently p = ρ+µ +µ, p 2 = αρ ρ( α), w = (ρ+µ)[ ρ( α)] (ρ+µ)[ ρ( α)] αρ(+µ) αρ(+µ) w 2 = (ρ+µ)[ ρ( α)] αρ(+µ) The pmf of the innovation sequence is then obtained from (4) is given by { π, if m = 0; Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m 2 ], if m The mean variance of{ǫ t } t Z are given by E[ǫ t ] = µ( α) ρ Var[ǫ t ] = µ( α)[ µ 2 (α ρ+)+µ ( αρ+α ρ 2 ρ+2 ) +αρ ( ρ 2 2ρ+2 ) ρ 2 + ] (+µ)( ρ) 3 Example 7 (Hurdle geometric INAR() process based on binomial thinning operator) In this example, we propose a new stationary first-order non-negative integer valued autoregressive process, {X t } t Z, with hurdle geometric marginals The proposed process is based on the binomial thinning operator satisfies the equation X t = α X t +ǫ t, t Z, (6) with α [0,) {X t } t Z being a stationary process with hurdle geometric marginals, ie, with pmf given 2

13 by { µ, m = 0, Pr(X t = m) = ( ) m ( ) µ ρ +ρ +ρ, m =,2,, where µ, ρ (0, ) its corresponding pgf is given by (7) Consequently, the pgf of α X t is given by ϕ Xt (s) = ( s)(µ+µρ ρ) +ρ ρs ϕ α Xt (s) = ϕ X (ϕ Y (s)) = α( s)(µ+µρ ρ), +αρ αρs where Pr(Y = ) = Pr(Y = 0) = α [0,) The hurdle geometric distribution is well-known as inflated-parameter Bernoulli distribution ρ- Bernoulli distribution (see Kolev et al, 2000) has also been discussed in a recent paper by Borges et al (206) to formulate a new thinning operator If ρ = 0 in (7), the hurdle geometric distribution coincides with the usual Bernoulli distribution with parameter µ, taking values 0 l, while if we put µ = ρ/( + ρ), the pmf of the geometric distribution with mean ρ is obtained The mean variance of{x t } t Z are given, respectively, by E[X t ] = α(+ρ) Var[X t ] = µ(+ρ)[ρ+(+ρ)( µ)], so the dispersion index is given by I X = ρ+(+ρ)( µ), that is, the hurdle geometric distribution has equidispersion/overdispersion/underdispersion according to the following conditions (for fixed 0 < µ < ): ρ = µ 2 µ equidispersion, µ 2 µ < ρ < overdispersion, 0 ρ < µ 2 µ underdispersion Now, let us derive the distribution of the rom variable {ǫ t } t Z Let ϕ X (s), ϕ Y (s) ϕ ǫ (s) be the pgf s of the rom variables X t,y j ǫ t, respectively Since the process {X t } t Z is a stationary process, it follows from (6) that ϕ ǫ (s) = ϕ X t (s) ϕ α Xt (s) = ( s)(µ+µρ ρ) +αρ αρs, s < (8) +ρ ρs α( s)(µ+µρ ρ) 3

14 After a simple algebra, we can identify the polynomial functions with corresponding coefficients a = αρ[ρ µ(+ρ)], with a > 0 if µ < ρ/(+ρ), b = {2αρ[ρ µ(+ρ)]+[ρ µ(+ρ)]+αρ}, c = +αρ[ρ µ(+ρ)]+[ρ µ(+ρ)]+αρ = (+ρ)(+α[ρ µ(+ρ)] ( α)µ), b = {2αρ[ρ µ(+ρ)]+α[ρ µ(+ρ)]+ρ}, c = +αρ[ρ µ(+ρ)]+α[ρ µ(+ρ)]+ρ = (+ρ){+α[ρ µ(+ρ)]} Note thatb+c = b+ c = αρ[ρ µ(+ρ)], c c, the roots of polynomial function in the denominator of (8) are s = +ρ ρ s 2 = +α[ρ µ(+ρ)], α[ρ µ(+ρ)] withs 2 s > From (4), the pmf of the innovation sequence is then obtained as { π, if m = 0; Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m 2 ], if m, µ with π = ( α) +α[ρ µ(+ρ)], p = ρ +ρ, p 2 = α[ρ µ(+ρ)] +α[ρ µ(+ρ)], w = ρ(+α[ρ µ(+ρ)]) ρ α[ρ µ(+ρ)] w 2 = α(+ρ)[ρ µ(+ρ)] ρ α[ρ µ(+ρ)] The mean variance of{ǫ t } t Z are given by E[ǫ t ] = µ( α)(+ρ), Var[ǫ t ] = µ( α) [ αρ 3 µ(+ρ) [ α ( ρ 2 +ρ+ ) +ρ+ ] +2(+α)ρ 2 +(2α+3)ρ+ ] 24 Method 4: A flexible quadratic rational probability generation function In this subsection, we extend the Method 3 to include a new family of innovation distributions by allowing the coefficients of the quadratic functions to be different In this general case, the pgf will be of the form ϕ ǫ (s) = as2 +bs+c ās 2 + bs+ c, s [0,] ϕ ǫ t (0) <, with the restrictions c c a + b + c = ā + b + c Since the rational pgf above involves two quadratic polynomial functions we need to reduce the degree of the quadratic function of the numerator which, after a simple algebra,can be rewritten as ϕ ǫ (s) = ā a + U (s) V 2 (s), where U (s) = [( b a b ) ( s+ c a c ) ] V 2 (s) = ās 2 + bs+ c = ā(s s )(s s 2 ) ā ā ā 4

15 Note that we are assuming that the quadratic function V 2 (s) has two distinct roots s s 2 Then, the rational function U (s) V 2 (s) can be decomposed into partial fractions a U (s) V 2 (s) = ρ (s s) + ρ 2 (s 2 s), where ρ = U (s ) s s 2 ρ 2 = U (s 2 ) s 2 s Using the same fractional expansion as the one the one in Method 3, the innovation distribution referred to as Modified hurdle-fractional geometric distribution (MHFG) is given by { π, if m = 0, Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m (9) 2 ], if m, where p = s, p 2 = s 2, w = ρ p 2 ( p )( π) w 2 = ρ 2 p 2 2 ( p 2 )( π) The MHFG innovation distribution is a two-stage distribution where the first stage is characterized by a Bernoulli variable with parameter π = c/ c the second stage is a mixture of two geometric distributions It is not difficult to see that if a = ā, the MHFG becomes the HFG obtained by Method 3 The mean the variance of the innovation process in (4) are given by Var[ǫ t ] = ( π) respectively E[ǫ t ] = ( π) { w { w ( p ) +w 2 [ ( ) ] 2 p µ Z + p ( p ) 2 } = ( π)µ Z ( p 2 ) +w 2 [ ( ) ]} 2 p 2 µ Z + p 2 ( p 2 ) 2 +π( π)µ 2 Z {, w ( w )(p p 2 ) 2 +w (p p 2 )+p 2 ( p ) 2 } = ( π) ( p ) 2 ( p 2 ) 2 +πµ 2 Z, Example 8 (Inflated-parameter geometric INAR() process based on NB thinning operator) A discretetime non-negative integer-valued stochastic process {X t } t Z is said to be an INAR() process with inflatedparameter geometric based on NB thinning operator if it satisfies the equation (0 α < ) X t = α X t +ǫ t, t Z, where * is the negative binomial thinning operator defined in Example 2 We assume that the rom variables of the process are identically distributed If ρ = 0, then the process {X t } t Z is reduced to the NGINAR() process introduced by Ristić et al (2009); see Example 2 Thus, the pgf of the innovation sequence is given by ϕ ǫ (s) = ϕ X t (s) ϕ α X (s) = sρ sρ+µ( s) +α( s) ρ+αµ( s), s < (20) +α( s) ρ After a simple algebraic manipulation, the pgf in (20) can be rewritten as a quotient of two quadratic polynomials wherea = αρ(+µ), b = [α(+µ)ρ+α(+µ)+( ρ)ρ], c = (+µ)[ ρ+α] 5

16 µ( ρ), ā = α(ρ+µ), b = [ρ+µ+α+αρ+2αµ µρ ρ 2 ], c = (+µ)[ ρ+α] Also, note thata+b+c = ā+ b+ c = ( ρ) 2,(+µ)[ ρ+α] µ( ρ) = c c = (+µ)[ ρ+α] the roots of the polynomial function in the denominator of (20) are s = +µ ρ+µ s 2 = ρ+α α withs 2 s > From (9), the pmf of the innovation sequence is then given by { π, if m = 0, Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m 2 ], if m, where π = µ( ρ) (+µ)[ ρ+α], p = ρ+µ +µ, p 2 = The mean variance of {ǫ t } t Z are given by α ρ+α, w = (α ρ+)(αµ+α µ ρ) (ρ )( α+µ+ρ) w 2 = α(+µ)(α ρ) (ρ )(α µ ρ) E[ǫ t ] = µ( α) ρ, Var[ǫ t ] = µ[ α 2 (+µ) 2 (ρ 2)+α ( µ 2 µ ( ρ 2 4ρ+ ) +2ρ ) +(+µ)( ρ)(+µ+ρ) ] (+µ)( ρ) 3 Example 9 (Hurdle geometric INAR() process based on NB thinning operator) Finally, in this example, we consider a new stationary INAR() process with hurdle geometric marginals based on NB thinning operator The proposed process satisfies the equation X t = α X t +ǫ t, t Z, where 0 α < The pgf of the innovation sequence is given by ϕ ǫ (s) = ( s)[µ(+ρ) ρ] +ρ( s) +α( s)+ρα( s) +α( s) α( s)[µ(+ρ) ρ] (2) Then, after a simple algebraic manipulation, the pgf in (2) can be rewritten as a quotient of two quadratic polynomials where: a = α(+ρ)[(+ρ)( µ) ], b = [2 α(+ρ)[ρ µ(+ρ)]+α(+ρ) µ(+ρ) +ρ], c = ( + ρ)[α( + ρ)( µ) + µ], ā = αρ( + ρ)( µ), b = [α( + ρ)( µ)(2ρ + ) +ρ] c = (+ρ)[α(+ρ)( µ)+] Note that a+b+c = ā+ b+ c = c c Moreover, the roots of the polynomial function in the denominator of (2) are s = +ρ ρ s 2 = +α(+ρ)( µ) α(+ρ)( µ) 6

17 withs 2 s > From (9), the pmf of the innovation sequence is then given by { π, if m = 0, Pr[ǫ t = m] = ( π)[w ( p )p m +w 2 ( p 2 )p m 2 ], if m, where π = α( µ)(ρ+) µ+, p = ρ α( µ)(ρ+)+ +ρ, p 2 = α(+ρ)( µ) +α(+ρ)( µ), w = (αρ+α ρ)[α(µ )(ρ+) ] α(µ )(ρ+)+ρ The mean variance of {ǫ t } t Z are given by w 2 = α(ρ+)[α(µ )(ρ+) µ(ρ+)+ρ] α(µ )(ρ+)+ρ E[ǫ t ] = µ( α)(+ρ) Var[ǫ t ] = µ(+ρ) [ α 2 (+ρ)(µρ+µ ρ 2)+α [ ( µ)ρ 2 ] µ(+ρ)+2ρ+ ] 3 Concluding remarks In the recent literature concerning INAR() models, many papers have assumed a known marginal distribution but usually it is difficult to underst how the innovation distribution was obtained In this paper, a new technique based on the fractional approach developed in Feller (see 2008, p 276), is formulated to find the innovation process by using a simple algebraic manipulations This fractional procedure has been discussed in detail for the linear quadratic probability generating functions illustrated with many recent INAR() models We also see that some of these examples reproduce some new innovation processes which could be of interest in this area For example, a new innovation process can be obtained if we assume an inflated-parameter geometric marginal distribution for the INAR() model as described in Method 4 As part of our future research study, we plan to study the four newly proposed processes in detail their properties, the recursive method associated inferential issues References Al-Osh, MA, & Aly, EEAA 992 First order autoregressive time series with negative binomial geometric marginals Communications in Statistics - Theory Methods, 2(9), Aly, EEAA, & Bouzar, N 207 Expectation thinning operators based on linear fractional probability generating functions arxiv Alzaid, A, & Al-Osh, M 988 First-order integer-valued autoregressive [INAR ()] process: distributional regression properties Statistica Neerlica, 42(), 53 6 Barreto-Souza, W 205 Zero-modified geometric INAR() process for modelling count time series with deflation or inflation of zeros Journal of Time Series Analysis, 36(6), Borges, P, Molinares, FF, & Bourguignon, M 206 A geometric time series model with inflated-parameter Bernoulli counting series Statistics & Probability Letters, 9(Supplement C),

18 Borges, P, Bourguignon, M, & Molinares, FF 207 A generalised NGINAR() process with inflatedparameter geometric counting series Australian & New Zeal Journal of Statistics, 59(), Bourguignon, M 208 Modelling time series of counts with deflation or inflation of zeros Manuscript submitted for publication Bourguignon, M, & Weiß, CH 207 An INAR() process for modeling count time series with equidispersion, underdispersion overdispersion TEST, 26(4), Feller, W 2008 An introduction to probability theory Its applications 3 edn Vol Hoboker, New Jersey: Jonh Wiley & Sons Jazi, MA, Jones, G, & Lai, C-D 202 Integer valued AR() with geometric innovations Journal of the Iranian Statistical Society, (2), Johnson, NL, Kemp, AW, & Kotz, S 2005 Univariate discrete distribution 3rd edition edn Hoboker, New Jersey: Jonh Wiley & Sons Kolev, N, Minkova, L, & Neytchev, P 2000 Inflated-parameter family of generalized power series distributions their application in analysis of overdispersed insurance data ARCH research clearing ARCH Research Clearing, 2, McKenzie, E 985 Some simple models for discrete variate time series Journal of the American Water Resources Association, 2(4), McKenzie, E 986 Autoregressive moving-average processes with negative-binomial geometric marginal distributions Advances in Applied Probability, 8(3), Miletić Ilić, AV 206 A geometric time series model with a new dependent Bernoulli counting series Communications in Statistics - Theory Methods, 45(2), Miletić Ilić, AV, Ristić, MM, Nastić, AS, & Bakouch, HS 207 An INAR() model based on a mixed dependent independent counting series Journal of Statistical Computation Simulation, 0, 5 Nastić, AS 202 On shifted geometric INAR() models based on geometric counting series Communications in Statistics - Theory Methods, 4(23), Nastić, AS, & Ristić, MM 202 Some geometric mixed integer-valued autoregressive (INAR) models Statistics & Probability Letters, 82(4), Nastić, AS, Ristić, MM, & Bakouch, HS 202 A combined geometric INAR(p) model based on negative binomial thinning Mathematical Computer Modelling, 55(5), Nastić, AS, Ristić, MM, & Janjić, AD 206a A mixed thinning based geometric INAR() model Filomat, to appear Nastić, AS, Laketa, PN, & Ristić, MM 206b Rom environment integer-valued autoregressive process Journal of Time Series Analysis, 37(2), Nastić, AS, Ristić, MM, & Miletić Ilić, AV 207 A geometric time-series model with an alternative dependent Bernoulli counting series Communications in Statistics - Theory Methods, 46(2),

19 Popović, PM 206 A bivariate INAR() model with different thinning parameters Statistical Papers, 57(2), Popović, PM, Ristić, MM, & Nastić, AS 206 A geometric bivariate time series with different marginal parameters Statistical Papers, 57(3), Ristić, MM, & Nastić, AS 202 A mixed INAR(p) model Journal of Time Series Analysis, 33(6), Ristić, MM, Bakouch, HS, & Nastić, AS 2009 A new geometric first-order integer-valued autoregressive [NGINAR()] process Journal of Statistical Planning Inference, 39(7), Ristić, MM, Nastić, AS, Jayakumar, K, & Bakouch, HS 202 A bivariate INAR() time series model with geometric marginals Applied Mathematics Letters, 25(3), Ristić, MM, Nastić, AS, & Miletić Ilić, AV 203 A geometric time series model with dependent Bernoulli counting series Journal of Time Series Analysis, 34(4), Yang, K, Wang, D, Jia, B, & Li, H 206 An integer-valued threshold autoregressive process based on negative binomial thinning Statistical Papers, to appear, 30 9