Next Article in Journal
Study on Queue Length in the Whole Process of a Traffic Accident in an Extra-Long Tunnel
Previous Article in Journal
Formal Group Fairness and Accuracy in Automated Decision Making
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Zero-and-One Integer-Valued AR(1) Time Series with Power Series Innovations and Probability Generating Function Estimation Approach

by
Vladica S. Stojanović
1,
Hassan S. Bakouch
2,3,
Eugen Ljajko
4 and
Najla Qarmalah
5,*
1
Department of Informatics & Computer Sciences, University of Criminal Investigation and Police Studies, 11060 Belgrade, Serbia
2
Department of Mathematics, College of Science, Qassim University, Buraydah 51452, Saudi Arabia
3
Department of Mathematics, Faculty of Science, Tanta University, Tanta 31111, Egypt
4
Department of Mathematics, Faculty of Sciences & Mathematics, University of Kosovska Mitrovica, 38220 Kosovska Mitrovica, Serbia
5
Department of Mathematical Sciences, Princess Nourah bint Abdulrahman University, Riyadh 11671, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(8), 1772; https://doi.org/10.3390/math11081772
Submission received: 20 March 2023 / Revised: 1 April 2023 / Accepted: 5 April 2023 / Published: 7 April 2023

Abstract

:
Zero-and-one inflated count time series have only recently become the subject of more extensive interest and research. One of the possible approaches is represented by first-order, non-negative, integer-valued autoregressive processes with zero-and-one inflated innovations, abbr. ZOINAR(1) processes, introduced recently, around the year 2020 to the present. This manuscript presents a generalization of ZOINAR processes, given by introducing the zero-and-one inflated power series (ZOIPS) distributions. Thus, the obtained process, named the ZOIPS-INAR(1) process, has been investigated in terms of its basic stochastic properties (e.g., moments, correlation structure and distributional properties). To estimate the parameters of the ZOIPS-INAR(1) model, in addition to the conditional least-squares (CLS) method, a recent estimation technique based on probability-generating functions (PGFs) is discussed. The asymptotic properties of the obtained estimators are also examined, as well as their Monte Carlo simulation study. Finally, as an application of the ZOIPS-INAR(1) model, a dynamic analysis of the number of deaths from the disease COVID-19 in Serbia is considered.

1. Introduction

Integer-valued time series modeling is a very popular research topic, with various applications (cf. [1,2,3,4,5]). One of the most popular approaches in modeling the dynamics of count data is provided by non-negative integer-valued autoregressive (INAR) time series models. This approach started from the famous work by Al-Osh and Alzaid [6], which first introduced the so-called INAR(1) process, and since then many results related to these models have been obtained (cf. [7,8,9,10,11,12,13,14,15,16,17]). One of the recently frequent problems in count data modeling is the presence of inflated zero-and-one values in the data, which can appear in various areas of human activity (e.g., the number of requests for issuing policies, breakdowns in the production process, injury in traffic accidents, etc.). To investigate this and similar problems, Saito et al. [18] and Zhang et al. [19] considered a modification (and generalization) of the traditional Poisson distribution, i.e., the so-called zero-and-one inflated Poisson (ZOIP) distribution. As an example of the application of the ZOIP distribution, the frequency of visits to the dentist in Swedish cities was considered in both mentioned manuscripts. Subsequently, Zhang et al. [20] introduced the multivariate ZOIP distribution, with applications in respect of healthcare demand data in Australia and car portfolio data in France.
Using the ZOIP distributed innovations, Qi et al. [21] introduced the first zero-and-one inflated INAR-based model, named the first-order zero-and-one INAR (ZOINAR(1)) process. Another class of ZOINAR time series, named the ZOIPLINAR process, has recently been introduced by Mohammadi et al. [22], where a Poisson–Lindley distribution of innovations inflated by zero and one is considered. Our main motivation is to introduce a more general form of the ZOINAR process, where the power series (PS) distribution with zero-and-one inflation is observed as its innovations. It should be noted that PS distributions represent a wide family of stochastic distributions, based on which many known integer-valued distributions can be obtained. In this way, the first-order zero-and-one inflated power series INAR process (abbr. ZOIPS-INAR(1) process) is proposed here, and it can be seen as a generalization of the previous ZOINAR models. The definition and the basic stochastic characteristics of this process are described in Section 2 and Section 3.
In order to estimate the parameters of INAR-based processes, various techniques have been developed. Conditional least-squares (CLS) estimation is a commonly utilized method, proposed in [23], as well as the method of conditional maximum likelihood (CML) [24]. However, in [25,26], among other methods, moment-based estimation procedures, i.e., Yule–Walker (YW) equations, are discussed. Nevertheless, to apply any of the aforementioned methods, there is usually an assumption that the estimation functions are given in a closed form and also bounded on some parameter space. As will be seen below, such conditions are not fully met in the case of our ZOIPS-INAR(1) process. Therefore, an alternative and a more contemporary approach, named the probability generating function (PGF) method, is proposed here. PGF estimation was theoretically described by Esquivel [27], first practically applied in Stojanović et al. [28], and recently examined in its general form by Stojanović et al. [29]. In order to apply the PGF method in the parameter estimation of the ZOIPS-INAR(1) process, some basic facts about this estimation method are given in Section 4. In addition, the asymptotic properties and efficiency of the PGF estimators of the ZOIPS-INAR(1) process, under some regulatory conditions, are considered here.
In Section 5, the PGF estimators for some specific ZOIPS innovations are analyzed. As typical members of the PS distributions family, but also for some practical reasons, the Poisson and geometric zero-and-one inflated distributions are considered here. For both of them, Monte Carlo simulations of PGF estimates were calculated and compared with the corresponding CLS estimates, which were taken as initial values for the PGF procedure. The asymptotic properties of both types of estimators are also examined here. The application of the ZOIPS-INAR(1) process in modeling the distribution of the number of deaths from the disease COVID-19 in the Republic of Serbia is presented in Section 6. In addition, by comparing the ZOIPS-INAR(1) model with the standard INAR(1) model, it is shown that the observed actual series has pronounced zero-and-one inflation, and that the proposed model has better efficiency and predictive accuracy. Finally, some concluding remarks are given in Section 7.

2. Structure of the ZOIPS-INAR(1) Process

In this section, similarly as in [28,29,30], we firstly introduce the independent identically distributed (IID) time series with the so-called power series (PS) distribution.
Definition 1.
The IID time series ( ε t ) , t Z is PS-distributed if its probability mass distribution (PMF) is as follows:
p ε ( x ; θ ) : = P { ε t = x } = a ( x ) θ x f ( θ ) , x S .
Here, S Z + = 0 , 1 , 2 , is the discrete set of values that the series ( ε t ) can take, and:
(i) 
a ( x ) 0 is a function defined on the set S ;
(ii) 
θ > 0 is the (unknown) one-dimensional parameter;
(iii) 
f ( θ ) : = x S a ( x ) θ x is the function that depends (only) on θ, and such that it is 0 < f ( θ ) < + , when 0 < θ < R .
Equation (1) can, for particular choices of a ( x ) , f ( θ ) and θ , give some of the most well-known types of discrete distributions (see Table 1 below). Nevertheless, we assume that the condition 0 , 1 S is fulfilled, as is usual in zero-and-one inflated distributions. Additionally, note that according to ( i i i ) , the power series f ( θ ) converges in fact on ( R , R ) . Nevertheless, the assumption θ > 0 is common for the PMF of PS-distributed series ( ε t ) , and we observe the convergence of f ( θ ) only on the positive interval ( 0 , R ) . Moreover, in this interval, the function f ( θ ) has positive, increasing values, as well as positive derivatives, as follows:
f ( n ) ( θ ) : = d n d θ n x = 0 a ( x ) θ x = x = n k = 0 n 1 ( x k ) a ( x ) θ x n , n N .
Equality (2) can be useful for determining the moments μ n ( ε ) : = E ( ε t n ) of series ( ε t ) . For this purpose, we have applied a similar procedure as in Stojanović et al. [28], based on the calculation of the moment-generating function (MGF):
M ε ( u ) = E e u ε t = 1 f ( θ ) x S a ( x ) ( θ e u ) x = f ( θ e u ) f ( θ ) .
Using (2) and the properties of the MGFs, after some calculations, one obtains:
μ n ( ε ) = E ε t n = d n M ( u ) d u n | u = 0 = 1 f ( θ ) k = 1 n b n ; k θ e u k f ( k ) θ e u | u = 0 = 1 f ( θ ) k = 1 n b n ; k θ k f ( k ) ( θ ) .
In doing so, the coefficients b n ; k , for each n N , are calculated recursively:
b n ; 1 = 1 , b n ; k + 1 = ( k + 1 ) n 1 k ! j = 1 k b n ; j ( k j + 1 ) ! , k = 1 , 2 , , n 1 .
Using the first two moments, obtained by Equation (3), the mathematical expectation and the variance of the random variables (RVs) ( ε t ) are obtained as follows:
μ ε : = E [ ε t ] = μ 1 ( ε ) = θ g ( θ ) σ ε 2 : = V a r [ ε t ] = μ 2 ( ε ) μ 1 ( ε ) 2 = μ ε + θ 2 g ( θ ) ,
where g ( θ ) = log f ( θ ) . If D ε ( θ ) : = σ ε 2 μ ε = θ 2 g ( θ ) is the so-called over-dispersion index, then the series ( ε t ) is over-dispersed; that is, D ε ( θ ) > 0 , if and only if the inequality g ( θ ) > 0 holds, for any θ ( 0 , R ) . Hence, the convexity of g ( θ ) indicates an overdispersion of the series ( ε t ) , as can be seen in Table 1. Moreover, for an arbitrary u 1 , 1 , we can introduce the first-order PGF of RVs ( ε t ) in the following way:
Ψ ε ( u ; θ ) : = E u ε t = 1 f ( θ ) x S a ( x ) ( θ u ) x = f ( θ u ) f ( θ ) .
The sum obtained above obviously converges on the interval u ( 0 , R / θ ) . In addition, the expression in (5) gives the possibility of the simple calculation of first-order PGFs for PS distributions, which are also given in Table 1.
In the following, we define a zero-and-one inflated distribution for an arbitrary PS-distributed time series ( ε t ) .
Definition 2.
Let ( ε t ) , t Z be the IID time series with the PS distribution, given by Equation (1). The series ( η t ) , t Z has a zero-and-one inflated power series (ZOIPS) distribution if for some ϕ 0 , ϕ 1 > 0 , such that 0 < ϕ 0 + ϕ 1 1 , its PMF is given as follows:
p η ( x ; λ η ) : = P { η t = x } = ϕ 0 + ( 1 ϕ 0 ϕ 1 ) p ε ( x ; θ ) , x = 0 ϕ 1 + ( 1 ϕ 0 ϕ 1 ) p ε ( x ; θ ) , x = 1 ( 1 ϕ 0 ϕ 1 ) p ε ( x ; θ ) , x 2 ,
where λ η = ( θ , ϕ 0 , ϕ 1 ) is the vector of (unknown) parameters.
Note that the ZOIPS distribution is a mixture of three distributions: I 0 ( x ) concentrated in zero, I 1 ( x ) concentrated in one, and the PS distribution of the series ( ε t ) . Thus, the PMF of RVs ( η t ) can be written as
p η ( x ; λ η ) = ϕ 0 I 0 ( x ) + ϕ 1 I 1 ( x ) + ϕ 2 p ε ( x ; θ ) ,
where we set ϕ 2 = 1 ϕ 0 ϕ 1 [ 0 , 1 ) , or equivalently, ϕ 0 + ϕ 1 + ϕ 2 = 1 . It is obvious that when ϕ 0 = ϕ 1 = 0 , the ZOIPS distribution is reduced to the previous PS distribution. For these reasons, we assume that ϕ 0 > 0 and ϕ 1 > 0 , so that these coefficients represent, respectively, the additional proportions of zeros and ones compared to those allowed by the PS distribution of the series ( ε t ) . Using the previous facts, and similar to Qi et al. [21], for the n-th moments μ n ( η ) : = E [ η t n ] of series ( η t ) one obtains:
μ n ( η ) = E η t n = ϕ 1 + ϕ 2 μ n ( ε ) .
According to (4) and (8), the mean and variance of the series ( η t ) are, respectively,
μ η : = E [ η t ] = ϕ 1 + ϕ 2 μ ε = ϕ 1 + ϕ 2 θ g ( θ ) σ η 2 : = V a r [ η t ] = μ 2 ( η ) μ η 2 = ϕ 1 + ϕ 2 μ 2 ( ε ) μ η 2 = ϕ 1 + ϕ 2 μ ε 2 + σ ε 2 μ η 2 = μ η μ η 2 + ϕ 2 μ ε 2 + θ 2 g ( θ ) = μ η μ η 2 + ϕ 2 θ 2 g ( θ ) 2 + g ( θ ) .
Using Equations (4) and (9), similarly as with the PS series ( ε t ) , one can obtain the necessary and sufficient conditions for the over-dispersion of the ZOIPS series ( η t ) . According to
D η ( λ η ) : = σ η 2 μ η = ϕ 2 μ ε 2 + σ ε 2 μ ε ϕ 1 + ϕ 2 μ ε 2 ,
the series ( η t ) will be over-dispersed if and only if ϕ 2 ( μ ε 2 + D ε ( θ ) ) > ( ϕ 1 + ϕ 2 μ ε ) 2 . This condition is more flexible than the “ordinary“ over-dispersion of the series ( ε t ) . For instance, if ( ε t ) is equally-dispersed, i.e., D ε ( θ ) = 0 holds, the inequality D η ( λ η ) > 0 is fulfilled when μ ε > ϕ 1 / ϕ 2 ϕ 2 . Then, the ZOIPS series ( η t ) will be overdispersed, which is the same result as for the Poisson ZOINAR process introduced in Qi et al. [21].
Finally, using Equations (1), (5) and (6), the first-order PGF of the RVs ( η t ) can be easily obtained. After some simple computations, for an arbitrary u 1 , 1 , it follows that:
Ψ η ( u ; λ η ) : = E u η t = k = 0 u k p η ( k ; λ η ) = ϕ 0 + ϕ 1 u + ϕ 2 k = 0 u k p ε ( k ; θ ) = ϕ 0 + ϕ 1 u + ϕ 2 f ( θ u ) f ( θ ) .
Assuming that the aforementioned notations are valid, we now introduce the INAR-based time series with the ZOIPS-distributed innovations ( η t ) .
Definition 3.
The time series ( X t ) , t Z , represents an INAR(1) process with ZOIPS innovations or, simply, a ZOIPS-INAR(1) process, if it fulfills the recurrence relation:
X t = α X t 1 + η t .
Here, ( η t ) is the ZOIPS-distributed time series with the PMF given by (6), α ( 0 , 1 ) is an unknown parameter, and
α X : = j = 1 X B j ( α )
is the binomial thinning operator. More precisely, for an arbitrary non-negative integer-valued RV X, the RVs B j = B j ( α ) have Bernoulli’s distribution P { B j = 1 } = 1 P { B j = 0 } = α . In addition, RVs B j are mutually independent (and also independent of X).
As an illustration, Figure 1 shows the realizations of the ZOIPS series and ZOIPS-INAR(1) process, where a Poisson distribution with the parameter θ = 5 is taken as the PS distribution. As can be easily seen, although the value of the parameter θ is significantly greater than zero (and one), both time series have emphasized zero-and-one inflation.

3. Stochastic Properties of the ZOIPS-INAR Process

Based on the mentioned properties of the ZOIPS series ( η t ) , some special properties of the ZOIPS-INAR(1) process can be shown. First, we compute the k-step conditional measures of X t + k on X t . Using some well-known properties of binomial thinning (cf. [31,32]) and Equation (11), for the first-step conditional mean one obtains:
E X t + 1 | X t = α X t + μ η = α X t + ϕ 1 + ϕ 2 θ g ( θ ) ,
and for the conditional variance:
V a r X t + 1 | X t = α ( 1 α ) X t + σ η 2 = α ( 1 α ) X t + μ η μ η 2 + ϕ 2 θ 2 g ( θ ) 2 + g ( θ ) .
In the general case, by the method of induction and after some computation, conditional measures of the k-degree can be computed for each k N . Hence, it follows that:
Theorem 1.
Let ( X t ) be a ZOIPS-INAR(1) process defined by Equation (11). Then, for each k = 1 , 2 , the k-step conditional mean and variance for the series ( X t ) are, respectively,
E X t + k | X t = α k X t + 1 α k 1 α μ η , Var X t + k | X t = α k ( 1 α k ) X t + ( 1 α k ) ( α α k ) 1 α 2 μ η + 1 α 2 k 1 α 2 σ η 2 ,
and the autocorrelation function (ACF) at lag k is ρ ( k ) = α k .
According to Equalities (13), when k , the unconditional mean and the variance of RVs ( X t ) can be obtained as follows:
μ X : = E [ X t ] = lim k E X t + k | X t = μ η 1 α = ϕ 1 + ϕ 2 θ g ( θ ) 1 α , σ X 2 : = V a r [ X t ] = lim k V a r X t + k | X t = α μ η + σ η 2 1 α 2 = μ X + ϕ 2 μ ε 2 + σ ε 2 μ ε μ η 2 1 α 2 .
Remark 1.
It can easily be seen that differences D X ( θ ) : = σ X 2 μ X and D η ( θ ) = σ η 2 μ η justify equality D X ( θ ) = D η ( θ ) / ( 1 α 2 ) . Thus, similarly to other INAR processes, the series ( X t ) and ( η t ) have the equivalent over-dispersed properties, i.e., they are both at the same time over-, equal- or under-dispersed.
In the following, we examine some characteristics of the distribution of the ZOIPS-INAR(1) process. First of all, we conduct an integer-valued moving average (INMA) representation of infinite order for the series ( X t ) .
Theorem 2.
Let us assume that PS series ( ε t ) , defined by Equation (1), has finite moments up to the order two, uniformly bounded for any θ ( 0 , R ) . Then, for any α ( 0 , 1 ) , ZOIPS-INAR(1) series ( X t ) , defined by Equation (11), has an INMA ( ) representation:
X t = d k = 0 α k η t k .
In addition, the sum in (14) converges in mean-square sense and almost surely.
Proof. 
Using the assumptions of the theorem, one can find a constant C > 0 such that
0 < μ ε ( θ ) C < + , θ ( 0 , R ) .
According to this, it follows that:
k = 1 P { ε t k } = k = 1 k P ε t = k = μ ε ( θ ) C < + ,
and using the definition of the ZOIPS series ( η t ) , given by Equation (6), one obtains:
k = 1 P { η t k } = ϕ 1 + ϕ 2 k = 1 P { ε t k } = ϕ 1 + ϕ 2 μ ε ( θ ) ϕ 1 + ϕ 2 C < + .
Hence, the above sum converges uniformly on θ ( 0 , R ) . Further, the sequence 1 / k , k = 1 , 2 , is monotone and bounded, so Abel’s convergence criterion for infinite sums gives
k = 1 1 k P { η t k } = ϕ 1 + ϕ 2 k = 1 1 k P { ε t k } < + ,
where the convergence above is uniformly on θ ( 0 , R ) . According to Theorem 2.1 in Alzaid and Al-Osh [33], condition (15) is sufficient for the equality
Ψ X ( u ; λ ) = k = 0 Ψ η 1 + α k ( u 1 ) ; λ η .
Here, Ψ η ( u ; λ η ) : = E [ u η t ] and Ψ X ( u ; λ ) : = E [ u X t ] are, respectively, the PGFs of ZOIPS series ( η t ) and ZOIPS-INAR(1) process ( X t ) , and λ : = ( θ , α , ϕ 0 , ϕ 1 ) is the vector of (unknown) parameters. Moreover, at least on u [ 1 , 1 ] , the product above converges absolutely. Hence, according to the one-to-one correspondence between the PMFs and PGFs of discrete RVs, it follows that the INMA ( ) representation in (14) is equivalent to (16).
To prove the second part of the theorem, note that, using Equation (11), for any n N it holds that:
X t = d α n X t n + j = 0 n 1 α j η t j .
According to this, as well as the well-known properties of the binomial thinning operator [31,32], it follows that:
E X t j = 0 n 1 α j η t j 2 = E α n X t n 2 = α 2 n E [ X t n 2 ] + α n ( 1 α n ) E [ X t n ] 0 , n .
Therefore, the mean-square convergence of the sum in (14) is valid.
To prove the almost certain convergence in (14), we define the event
A : = lim n j = 0 n 1 α j η t j = X t .
According to (17) and the definition of the limit value of real functions, we can write the event A as
A : = lim n α n X t n = 0 = δ > 0 k = 1 n = k 0 α n X t n < δ = k = 1 A k ,
where A k : = n = k α n X t n = 0 . Again using (17), for each (fixed) k = 1 , 2 , , and m = 1 , 2 , , one obtains:
α k X t k = d α m + k X t m k + n = k m + k 1 α n η t n ,
where the expression on the right is the sum of the uncorrelated RVs. Thus, applying the continuity of probability and the definition of the thinning operator (12), for events A k we have:
P ( A k ) = lim m P α m + k X t m k = 0 × n = k m + k 1 P α n η t n = 0 = lim m j = 0 1 α m + k j P X t m k = j × lim m n = k m + k 1 j = 0 1 α n j P η t n = j = lim m Ψ X ( 1 α m + k ; λ ) × lim m n = k m + k 1 Ψ η ( 1 α n ; λ η ) = Ψ X ( 1 ; λ ) × n = k Ψ η ( 1 α n ; λ η ) = n = k Ψ η ( 1 α n ; λ η ) .
By re-applying the continuity property of the probability, as well as the convergence of the product in (16), it follows that P ( A ) = lim k P ( A k ) = 1 , i.e., the sum in (14) converges almost surely.    □
Remark 2.
Using a similar procedure as in the previous theorem and some general results about the PGFs of non-negative discrete-valued stationary time series (cf. Stojanović et al. [29]), an explicit expression for PGFs of the ZOIPS-INAR(1) process can be obtained. According to Equations (10) and (16), for arbitrary α ( 0 , 1 ) , the series ( X t ) has a PGF of the first order:
Ψ X ( u ; λ ) = k = 0 ϕ 0 + ϕ 1 1 + α k ( u 1 ) + ϕ 2 f 1 + α k ( u 1 ) θ f θ .
Furthermore, suppose that u = ( u 1 , , u r ) R r , r 2 and X t ( r ) : = ( X t , , X t + r 1 ) , t Z are the so-called overlapping blocks of series ( X t ) . Putting k = 1 , , r 1 into Equation (17), after some calculations, the explicit expression of the r-dimensional PGF random vector X t ( r ) can be obtained:
Ψ X ( r ) ( u ; λ ) : = E u 1 X t u r X t + r 1 = Ψ X k = 0 r 1 1 + α k ( u k + 1 1 ) ; λ = 2 r Ψ η k = 0 r 1 + α k ( u k + 1 ) ; λ η .
It can be noted that the PGF Ψ X ( r ) ( u ; λ ) of the order r = 2 will be used in the estimation of parameters of the ZOIPS-INAR(1) process (see Section 4 below).
Let us now consider the Markov properties and marginal distribution of our model.
Theorem 3.
Let ( η t ) be the ZOIPS series, with the PMF given by (6). Then, the process ( X t ) , defined by (11), is a homogeneous Markovian process with the first-step transition probabilities:
p i j : = P X t = j | X t 1 = i = ϕ 0 i j α j ( 1 α ) i j I j i + ϕ 1 i j 1 α j 1 ( 1 α ) i j + 1 I 1 j i + 1 + ϕ 2 x = 0 m i x α x ( 1 α ) i x p ε ( j x ; θ ) ,
where i , j S and m = min { i , j } .
Proof. 
According to Equation (11) and the definition of binomial thinning, the conditional distribution of X t at a given X t 1 can be expressed as follows:
p i j : = P α X t 1 + η t = j | X t 1 = i = x = 0 m P α X t 1 = x | X t 1 = i P { η t = j x } = x = 0 m i x α x ( 1 α ) i x p η ( j x ; λ η ) = x = 0 m p i ( x ; α ) p η ( j x ; λ η ) ,
where
p i ( x ; α ) : = i x α x ( 1 α ) i x , x = 0 , 1 , , i
is the PMF of the binomial B ( i ; α ) distribution. Using the definition of the ZOIPS distribution, that is, the PMF of RVs ( η t ) given by Equation (7), one obtains:
p i j = x = 0 m p i ( x ; α ) ϕ 0 p 0 ( j x ) + ϕ 1 p 1 ( j x ) + ϕ 2 p ε ( j x ; θ ) ,
where p 0 ( k ) and p 1 ( k ) are, respectively, the PMFs of I 0 = a s 0 and I 1 = a s 1 . Based on these, we obtain:
p i j = ϕ 0 x = 0 m p i ( x ; α ) I ( j = x ) + ϕ 1 x = 0 m p i ( x ; α ) I ( j = x + 1 ) + ϕ 2 x = 0 m p i ( x ; α ) p ε ( j x ; θ ) = ϕ 0 p i ( k ; α ) I ( j m ) + ϕ 1 p i ( j 1 ; α ) I ( 1 j m + 1 ) + ϕ 2 x = 0 m p i ( x ; α ) p ε ( j x ; θ ) .
It is easy to see that the obtained equality is equivalent to (20), which proves the theorem.    □
Remark 3.
Let us recall once again that it is usual to assume that for the zero-and-one distribution, the condition 0 , 1 S is usually taken. It follows that the first singular part in Equation (20) exists if and only if ( X t ) passes from the state X t 1 = i to the non-increasing state X t = j i . Similarly, the second singular part exists if and only if 1 X t = j i + 1 . Finally, the transition probabilities (20), as well as the use of the conditional probability, give the marginal PMF of the ZOIPS-INAR(1) process:
p X ( x ; λ ) : = P X t = x = k S P X t = x | X t 1 = k P X t 1 = k = ϕ 0 α 1 α x k = 0 k x ( 1 α ) k I ( x k ) p X ( k ; λ ) + ϕ 1 α 1 α x 1 k = 0 k x 1 ( 1 α ) k I ( 1 x k + 1 ) p X ( k ; λ ) + ϕ 2 k = 0 p X ( k ; λ ) j = 0 m p k ( j ; α ) p ε ( x j ; θ ) .
Thus, the ZOIPS-INAR(1) process ( X t ) is a strictly stationary and ergodic time series.
At the end of this section, similar to Qi et al. [21] and Mohammadi et al. [22], we examine the distribution of the zero-and-one lengths of the ZOIPS-INAR(1) process. Starting with the basic results of the zero-inflated INAR processes (cf. Jazi et al. [34], Wang et al. [35]), we observe the distribution of the “runs”of zeros (resp. ones) in the ZOIPS-INAR(1) process. Thereby, the “runs“ are defined as the number of zeros (resp. ones) between two different non-zeros, i.e., non-one values, respectively. In the following statements, we give the expected lengths and proportions of zeros and ones in our model.
Theorem 4.
The expected lengths of the runs of zeros (resp. ones) for a ZOIPS-INAR(1) process are, respectively, given by:
L 0 = ϕ 0 f ( θ ) + ϕ 2 a ( 0 ) ( 1 ϕ 0 ) f ( θ ) ϕ 2 a ( 0 ) , L 1 = α ϕ 0 f ( θ ) + ϕ 2 a ( 0 ) + ( 1 α ) ϕ 1 f ( θ ) + ϕ 2 a ( 1 ) θ α ( 1 ϕ 0 ) f ( θ ) ϕ 2 a ( 0 ) + ( 1 α ) ( 1 ϕ 1 ) f ( θ ) ϕ 2 a ( 1 ) θ ,
where a ( x ) and f ( θ ) are the functions introduced in Definition 1 for the PS distributions.
Proof. 
According to Equation (20), the transition probabilities from zero to zero, and from zero to non-zero values are, respectively,
p 0 * : = P X t = 0 | X t 1 = 0 = ϕ 0 + ϕ 2 p ε ( 0 ; θ ) , 1 p 0 * = 1 ϕ 0 ϕ 2 p ε ( 0 ; θ ) = ϕ 1 + ϕ 2 1 p ε ( 0 ; θ ) .
Since the zero run length is defined as the number of zeros between two non-zero values, it can easily be seen to follow a geometric distribution with parameter 1 p 0 * . Therefore, the expected length of zero is
L 0 = p 0 * 1 p 0 * = ϕ 0 + ϕ 2 p ε ( 0 ; θ ) ϕ 1 + ϕ 2 1 p ε ( 0 ; θ ) ,
where p ε ( 0 ; θ ) = a ( 0 ) / f ( θ ) , and the first equality in (21) immediately follows. Similarly, the transition probabilities from one to one and from one to non-one values are, respectively,
p 1 * : = P X t = 1 | X t 1 = 1 = α ϕ 0 + ϕ 2 p ε ( 0 ; θ ) + ( 1 α ) ϕ 1 + ϕ 2 p ε ( 1 ; θ ) , 1 p 1 * = α 1 ϕ 0 ϕ 2 p ε ( 0 ; θ ) + ( 1 α ) 1 ϕ 1 ϕ 2 p ε ( 1 ; θ ) .
Applying the same procedure as before, the second equality in (21) is easily obtained.    □
Let us point out that, in the same way as in similar INAR-based models, the expected length of the zero runs in the ZOIPS-INAR(1) process is independent of the parameter α .
Theorem 5.
The proportions of zeros and ones in the ZOIPS-INAR(1) process are, respectively,
P X t = 0 = k = 0 ϕ 0 + ϕ 1 1 α k + ϕ 2 f ( 1 α k ) θ f θ , P X t = 1 = j = 0 α j ϕ 1 + ϕ 2 θ f ( 1 α j ) θ f θ k j ϕ 0 + ϕ 1 1 α k + ϕ 2 f ( 1 α k ) θ f θ .
Proof. 
Using the well-known properties of the PGFs of discrete-valued RVs ( X t ) , for an arbitrary k = 0 , 1 , the PMF of ( X t ) can be expressed as follows:
P X t = k = 1 k ! k Ψ X ( u ; λ ) u k | u = 0 .
From here, using expression (18) for the first-order PGF of the ZOIPS-INAR process, and replacing k = 0 , 1 in Equation (22), the statement of the theorem immediately follows.    □

4. Parameter Estimation Procedure

Due to the specific structure, the parameter estimation procedure of the ZOIPS-INAR(1) model is more complex than for the ordinary INAR model. The main reason is that the basic INAR models with ‘ordinary’ PS innovations ( ε t ) have (only) two unknown parameters θ , λ . However, in the non-trivial case, the ZOIPS-INAR(1) process has two additional parameters ϕ 0 , ϕ 1 0 , 1 . Therefore, the structure of this process affects the fact that even some simpler types of estimators, such as YW estimators, cannot be obtained by simple calculation. In previous works on ZOINAR processes, Qi et al. [21] proposed the CML estimation method, while Mohammadi et al. [22] additionally described the CLS estimates. In the following, CLS estimators will also be used as initial values for the PGF estimation procedure, which will now be given more attention.
We emphasize once again that the general aspect of the PGF method was recently described by Stojanović et al. [29]. Accordingly, a specific PGF estimation procedure is examined here, in the case of the ZOIPS-INAR(1) process. The basic idea of the PGF method is close to the empirical characteristic function (ECF) estimation method introduced in [36,37]. It is based on minimizing the ‘distance’ between the theoretical PGF of order r N , defined by Equations (18) and (19), as well as the corresponding empirical PGF:
Ψ ˜ T ( r ) ( u ) : = 1 T r + 1 t = 1 T r + 1 u 1 X t u r X t + r 1 ,
where { X 1 , X 2 , , X T } is some finite realization of the ZOIPS-INAR series ( X t ) . Since the ZOIPS-INAR(1) series ( X t ) is stationary and ergodic, hence it follows that:
E Ψ ˜ T ( r ) ( u ) = Ψ X ( r ) ( u ; λ 0 ) ,
where λ 0 R 4 is the (unknown) true value of parameter λ = ( θ , α , ϕ 0 , ϕ 1 ) . Thus, Ψ ˜ T ( r ) ( u ) is an unbiased estimator of Ψ X ( r ) ( u ; λ 0 ) . Further, as the PGF Ψ X ( r ) ( u ; λ ) is well-defined at least for all u [ 1 , 1 ] r , the objective (minimization) function can be defined as follows:
S T ( r ) ( λ ) : = 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ ) Ψ ˜ T ( r ) ( u ) 2 d u ,
where d u : = d u 1 d u r and ω : R r R + is a weight function, integrable on [ 1 , 1 ] r .
The PGF estimators are then obtained by minimizing the objective function (24) with respect to λ . In other words, they represent the solutions of the following minimization equation:
λ ^ T = a r g min λ Λ S T ( r ) ( λ ) ,
where Λ = ( 0 , R ) × 0 , 1 3 R 4 is a regular parameter space of the ZOIPS-INAR(1) process. To solve Equation (25), numerical integration procedures have been used, which are discussed in Section 5. In the following, under some regularity conditions, we examine the consistency and asymptotic normality (AN) of the PGF estimators.
Theorem 6.
Let λ 0 be the exact value of the parameter λ, and λ ^ T , for arbitrary T = 1 , 2 , , the solutions of Equation (25). In addition, assume that the following regularity conditions are fulfilled:
( A 1 )
λ 0 Λ and λ ^ T Λ for large enough T.
( A 2 )
At the point λ = λ 0 function
S 0 ( r ) ( λ ) : = 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ ) Ψ X ( r ) ( u ; λ 0 ) 2 d u
has a unique minimum S 0 ( r ) ( λ 0 ) = 0 .
( A 3 )
2 S T ( r ) ( λ 0 ) λ λ is a regular matrix.
( A 4 )
Ψ X ( r ) ( u ; λ 0 ) λ Ψ X ( r ) ( u ; λ 0 ) λ is a non-zero matrix uniformly bounded by a positive ω-integrable function W : R r R + .
Then, the estimator λ ^ T ( r ) is strictly consistent and AN for the parameter λ.
Proof. 
To prove the consistency of the estimator λ ^ T , we firstly check the sufficient conditions for the consistency of the extremum estimators (cf. Newey and McFadden [38]). As it was previously shown that the ZOIPS-INAR(1) series ( X t ) is ergodic, applying Equation (23) and the strong law of large numbers (SLLN) follows:
Ψ ˜ T ( r ) ( u ) as Ψ X ( r ) ( u ; λ 0 ) , T ,
where “as” denotes the almost sure convergence. Further, under assumption ( A 1 ) , the closed set Λ ¯ = [ 0 , R ] × [ 0 , 1 ] 3 is compact, and λ 0 belongs to its interior. Hence, the functions Ψ X ( r ) ( u ; λ ) and Ψ ˜ T ( r ) ( u ) are continuous on the compacts [ 1 , 1 ] r × Λ ¯ and [ 1 , 1 ] r , respectively, and therefore, there exist constants C 1 , C 2 > 0 such that
max ( u ; λ ) [ 1 , 1 ] r × Λ ¯ | Ψ X ( r ) ( u ; λ ) | C 1 < + , max u [ 1 , 1 ] r | Ψ ˜ T ( r ) ( u ) | C 2 < + .
According to these, similarly as in Stojanović et al. [29], one obtains
| S T ( r ) ( λ ) S 0 ( r ) ( λ ) | ( 3 C 1 + C 2 ) 1 1 1 1 ω ( u ) | Ψ ˜ T ( r ) ( u ) Ψ X ( r ) ( u ; λ 0 ) | d u ,
so the last inequality and Equation (26) imply
sup λ Λ ¯ | S T ( r ) ( λ ) S 0 ( r ) ( λ 0 ) | as 0 , T + .
Thus, the function S T ( r ) ( λ ) converges uniformly and almost surely to S 0 ( r ) ( λ ) . According to these, as well as assumption ( A 2 ) and Theorem 2.1 in Newey and McFadden [38], it follows that:
λ ^ T ( r ) λ 0 as 0 , T + ,
that is, the estimator λ ^ T is strictly consistent for λ .
To prove the AN property of the estimator λ ^ T , note that the first two orders of partial derivatives of the function S T ( r ) ( λ ) are continuous functions. Therefore, using the Taylor expansion for S T ( r ) ( λ ) / λ at the point λ = λ 0 , one obtains:
S T ( r ) ( λ ) λ = S T ( r ) ( λ 0 ) λ + 2 S T ( r ) ( λ 0 ) λ λ ( λ λ 0 ) + o ( λ λ 0 ) .
By replacing λ with λ ^ T , for sufficiently large T, under assumption ( A 3 ) and the fact that S T ( r ) ( λ ^ T ) / λ = 0 , we have:
λ ^ T ( r ) λ 0 = 2 S T ( r ) ( λ 0 ) λ λ 1 S T ( r ) ( λ 0 ) λ + o λ ^ T λ 0 .
Furthermore, according to the properties mentioned above, the function S T ( r ) ( λ ) can be differentiated under the sign of the integral, such as:
S T ( r ) ( λ ) λ = 2 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ ) Ψ ˜ T ( u ) Ψ X ( r ) ( u ; λ ) λ d u , 2 S T ( r ) ( λ ) λ λ = 2 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ ) λ Ψ X ( r ) ( u ; λ ) λ
+ Ψ X ( r ) ( u ; λ ) Ψ ˜ T ( u ) 2 Ψ X ( r ) ( u ; λ ) λ λ d u .
By taking the mathematical expectations in Equations (28) and (29), one obtains:
E S T ( r ) ( λ 0 ) λ = 0 r × 1 , E 2 S T ( r ) ( λ 0 ) λ λ = 2 V ,
where, according to Equation (23), it follows that:
V = 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ 0 ) λ Ψ X ( r ) ( u ; λ 0 ) λ d u .
According to assumption ( A 4 ) , there exists a ω -integrable function V : R r R + such that
0 < Ψ X ( r ) ( u ; λ 0 ) λ Ψ X ( r ) ( u ; λ 0 ) λ W ( u ) , u [ 1 , 1 ] r ,
where · is the matrix norm on R r × R r consistent with the Euclidean vector norm. Hence, inequalities
0 < V 1 1 1 1 ω ( u ) W ( u ) d u < +
hold, so Equation (30) and SLLN give
S T ( r ) ( λ 0 ) λ , 2 S T ( r ) ( λ 0 ) λ λ a s ( 0 , 2 V ) , T + .
Further, note that Equation (28) for the gradient of function S T ( r ) ( λ ) can be written as
S T ( r ) ( λ ) λ = 2 T r + 1 t = 1 T r + 1 K t ( r ) ( λ ) ,
where
K t ( λ ) : = 1 1 1 1 ω ( u ) Ψ X ( r ) ( u ; λ ) u 1 X t u r X t + r 1 Ψ X ( r ) ( u ; λ ) λ d u .
Thereby, according to Equation (23), the equality E [ K t ( λ 0 ) ] = 0 holds. It can also be shown (cf. Stojanović et al. [39,40]) that
W 2 : = lim T Var 1 T r + 1 t = 1 T r + 1 K t ( λ 0 ) = lim T 1 T r + 1 E t = 1 T r + 1 K t ( λ 0 ) 2 = lim T 1 ( T r + 1 ) 2 t = 1 T r + 1 s = 1 T r + 1 Cov K t ( λ 0 ) , K s ( λ 0 )
is the finite non-zero limit value if the covariance function γ ( k ) : = Cov ( X t , X t + k ) , when k = 0 , ± 1 , ± 2 , , is absolutely summable. In the case of the ZOIPS-INAR ( 1 ) process, using Theorem 1 we obtain:
C : = k = + γ ( k ) = γ ( 0 ) k = + ρ ( k ) = σ X 2 2 k = 0 + α k 1 = σ X 2 1 + α 1 α .
Thus, 0 < C < + hold for any λ Λ . By applying the central limit theorem for stationary processes (cf. Billingsley [41]), the convergence that is proved and Equations (31)–(33) give:
T r + 1 S T ( r ) ( λ 0 ) λ d N ( 0 r × 1 , 4 W 2 ) , T + ,
where “d” denotes the convergence in the distribution. Finally, according to Equations (27), (31) and (34) it follows that:
T r + 1 λ ^ T ( r ) λ 0 d N ( 0 r × 1 , V 1 W 2 V 1 ) , T + ,
which completes the proof of the theorem.    □
Remark 4.
Using similar considerations as in ECF estimates (cf. Knight & Yu [36]), the PGF procedure for estimating the true values of the parameter λ = λ 0 is based on the realization of the two-dimensional random vector X t ( 2 ) : = ( X t , X t + 1 ) . Then, the objective function S T ( 2 ) represents a double integral with respect to the weight function ω : R 2 R + , and can be numerically approximated by some cubature formulas. For that purpose, it is necessary to determine the two-dimensional PGF (as well as EPGF) of the ZOIPS-INAR(1) series ( X t ) . By replacing r = 2 in Equation (19), and using Equation (18), the two-dimensional theoretical PGF can be obtained as follows:
Ψ X ( 2 ) ( u 1 , u 2 ; λ ) = Ψ X u ; λ Ψ ε u 2 ; λ η = k = 0 ϕ 0 + ϕ 1 1 + α k ( u 1 ) + ϕ 2 f 1 + α k ( u 1 ) θ f θ × ϕ 0 + ϕ 1 u 2 + ϕ 2 f u 2 θ f θ ,
where u = u 1 1 + α ( u 2 1 ) . As an illustration, Figure 2 presents the theoretical and empirical PGF of the ZOIPS-INAR(1) process with geometric PS innovations, which were obtained using the ‘R’ function “persp()”.

5. Numerical Simulations

In this section, numerical simulations of the proposed PGF procedure for estimating the unknown parameters λ = ( θ , α , ϕ 0 , ϕ 1 ) of the ZOIPS-INAR(1) process are performed. For this purpose, as previously noted, different PS-distributed series ( ε t ) can be observed. As an illustration, but also for practical reasons that are stated in the next section, two different distributions of the series ( ε t ) are observed. First it is assumed that RVs ( ε t ) have a Poisson distribution, and then they are assumed to have a geometric distribution. For both of these distributions, we generated samples of length T = 1000 , whose size is close to the length of the COVID-19 count series data that will be analyzed in Section 6. These samples were generated through 500 independent Monte Carlo simulations of the PS series ( ε t ) as well as the ZOIPS series ( η t ) . After that, according to Equations (11) and (12), the corresponding realizations { X 1 , , X T } of the ZOIPS-INAR ( X t ) series were obtained.
Using a similar procedure as in Mohammadi et al. [22], we firstly computed the CLS estimates by minimizing the objective function:
Q T ( λ ) : = t = 2 T X t E X t | X t 1 2 = t = 2 T X t α X t 1 μ η 2 = t = 2 T X t α X t 1 ϕ 1 ϕ 2 θ g ( θ ) 2 .
By applying the usual procedure, that is, by solving coupled equations
Q T ( λ ) α = Q T ( λ ) μ η = 0 ,
parameter estimators α and μ η can be easily obtained as follows:
α ˜ : = ρ ^ X ( 1 ) = ( T 1 ) t = 2 T X t X t 1 t = 2 T X t X t 1 t = 2 T ( T 1 ) t = 2 T X t 2 t = 2 T X t 2 , μ ˜ η = t = 2 T X t α ˜ t = 2 T X t 1 T .
In the next step, estimates of the parameters θ , ϕ 0 , ϕ 1 can be obtained by minimizing the objective function:
Q T * ( λ ) : = t = 2 T X t E X t | X t 1 2 V a r X t | X t 1 2 = t = 2 T X t α X t 1 + μ η 2 α ( 1 α ) μ η ( 1 μ η ) ϕ 2 θ 2 g ( θ ) 2 + g ( θ ) 2 ,
where α and μ η were replaced by their CLS estimators α ˜ and μ ˜ η , respectively. Minimization of the function Q T * ( λ ) was conducted using a numerical procedure based on the R-function "nlminb", where the initial values of the parameters were taken randomly from the uniform distribution U ( 0 , 1 ) . The asymptotic properties of the obtained CLS estimates can be proven by applying some basic results of the CLS theory [42,43], in the same way as this was achieved in Mohammadi et al. [22]. The results of these simulations are given in the left part of Table 2, where the minimums (Min.), mean values (Mean), maximums (Max.) and mean squared estimated errors (MSEE ) of the CLS estimates are shown.
Next, the PGF method was applied, with initial values obtained from the previous CLS procedure. The parameter estimates λ = ( θ , α , ϕ 0 , ϕ 1 ) were calculated based on the minimization of the double integral
S T ( 2 ) ( λ ) = [ 1 , 1 ] 2 ω ( u 1 , u 2 ) Ψ X ( 2 ) ( u 1 , u 2 ; λ ) Ψ ˜ T ( u 1 , u 2 ) 2 d u 1 d u 2 ,
where ω : [ 1 , 1 ] 2 R + is the weight function, and Ψ X ( 2 ) ( u 1 , u 2 ; λ ) is the two-dimensional PGF of the ZOIPS-INAR(1) process ( X t ) . Using Equation (35), in the case of the series ( X t ) with ZOIPS for the Poisson distribution case, this PGF can be obtained as follows:
Ψ X ( 2 ) ( u 1 , u 2 ; λ ) = k = 0 ϕ 0 + ϕ 1 1 + α k ( u 1 ) + ϕ 2 exp θ α k ( u 1 ) × ϕ 0 + ϕ 1 u 2 + ϕ 2 exp θ ( u 2 1 ) .
Similarly, for the appropriate PGF of the series ( X t ) with ZOIPS for the geometric innovations, one obtains:
Ψ X ( 2 ) ( u 1 , u 2 ; λ ) = k = 0 ϕ 0 + ϕ 1 1 + α k ( u 1 ) + ϕ 2 1 θ 1 1 + α k ( u 1 ) θ × ϕ 0 + ϕ 1 u 2 + ϕ 2 1 θ 1 θ u 2 .
It can be observed that these PGFs are not in closed form, but they can be approximated by finite k-term products with an arbitrary precession.
Thereafter, the integral in (37) can be numerically approximated using some of the N-point cubature formulas of the form:
I ( φ ; ω ) : = [ 1 , 1 ] 2 ω ( u 1 , u 2 ) φ ( u 1 , u 2 ) d u 1 d u 2 j = 1 N ω j φ ( u 1 j , u 2 j ) .
Here, ( u 1 j , u 2 j ) are the cubature nodes, and ω j denotes the appropriate weight coefficients. In this simulation study, we used cubature formulas with N = 36 nodes, based on Gauss–Legendre orthogonal polynomials and weight function ω ( u 1 , u 2 ) 1 . The numerical construction of these formulas was carried out using the package “Orthogonal polynomials” within the software Mathematica, authored by Cvetković and Milovanović [44]. Next, the objective function (24) was minimized using the “R” procedure for linearly constrained minimization, based on the Nelder–Mead optimization method [45]. Summary statistics of the PGF estimates, which were obtained using the aforementioned estimation procedure and with additional values of the objective function S T ( 2 ) , are shown in the right part of Table 2.
Comparison of the CLS and PGF estimated values indicates that the mean values of CLS estimates are somewhat closer to the true parameter values (only) for the parameter α . This is expected, because α represents the first correlation of the ZOIPS-INAR(1) series ( X t ) . Since the CLS-estimate α ˜ , given by the first of Equation (36), is a sampled correlation, it is the most efficient estimate for α . However, it is obvious that the other parameter estimates are more efficient in the case when the PGF estimation procedure is applied, as well as that they have smaller mean squared estimation errors (MSEEs).
In addition, the AN test results are shown in Table 2, where the Anderson–Darling normality test was conducted. The test statistic, denoted AD, along with the corresponding p-values, were calculated using the procedure from the R-package “nortest”, authored by Gross [46]. According to the obtained values, it can be seen that the AN property is confirmed for most of the PGF estimates of parameters λ = ( θ , α , ϕ 0 , ϕ 1 ) . However, the CLS estimates of the parameter θ do not have the AN property, at the significance level of p = 0.05 . Therefore, this would be another advantage of the PGF estimates. Certain confirmation of these facts can be observed visually in Figure 3 and Figure 4.

6. Application of the Model

Here we point out some possibilities of practical application of the ZOIPS-INAR process in real-world data modeling. In this regard, it is worth noting that the COVID-19 pandemic has received a great deal of attention since its appearance. From a mathematical point of view, various theoretical models have been proposed to investigate this still current problem (for more recent examples, see [47,48,49,50,51,52,53]). To that end, here we explore some additional possibilities in modeling the dynamics of COVID data.
More precisely, we observed an actual count data set, which represents the dynamics of the number of deaths due to COVID-19 in the Republic of Serbia, based on the data of the World Health Organization (WHO) [54], over a period from 1 January 2020 to 31 December 2022. In this way, time series of counting data of length T = 1094 were observed, the dynamics part of which is shown in the diagram in Figure 5. In addition, the autocorrelation function (ACF) and partial autocorrelation function (PACF) of this time series can be seen. Depending on the lag k = 0 , 1 , 2 , , it is clear that there is an exponentially decreasing autocorrelation, as well as that partial autocorrelation indicates suitable modeling with the INAR process of the order one (or two). Therefore, it could be expected that INAR-based processes can be adequate stochastic models to describe these dynamics.
Furthermore, based on the summary statistics of the data of this time series, shown in Table 3, it can be observed that COVID-19 data have a significant number of zero and one values. Therefore, it could be assumed that the ZOIPS-INAR(1) process can be taken as a suitable stochastic model. In doing so, as members of the PS family, we consider innovations with the Poisson distribution (and the ‘small’ parameter θ ), and with the geometric distribution. In the following, we first assume Poisson PS-innovations, since a zero-one inflation testing procedure can be applied. In this regard, we point out that the testing procedure here is different than one used in Qi et al. [21]. We form a null hypothesis that a certain INAR time series is not of the ZOI-type; that is, H 0 : ϕ 0 = ϕ 1 = 0 . In that case, according to Theorem 5, that is, Equation (22), the proportions of the occurrences of the zero and one values can be, respectively, expressed as follows:
p 0 : = P X t = 0 = Ψ X ( 0 ; λ ) = k = 0 exp θ 1 + α k ( u 1 ) exp ( θ ) | u = 0 = exp θ ( u 1 ) 1 α | u = 0 = exp θ 1 α ,
p 1 : = P X t = 1 = Ψ X ( 0 ; λ ) u = θ 1 α exp θ ( u 1 ) 1 α | u = 0 = θ 1 α exp θ 1 α .
Therefore, in the INAR(1) process with ‘ordinary’ Poisson innovations, the zero-and-one proportions are exponentially related to the mean μ X = θ / ( 1 α ) . If we define the so-called sample zero-and-one proportions:
p ^ 0 = 1 T t = 1 T I { X t = 0 } , p ^ 1 = 1 T t = 1 T I { X t = 1 } ,
then, according Equations (38) and (39), the so-called zero-and-one test statistics can be taken as follows:
I ^ 0 = p ^ 0 e X ¯ T σ ^ 0 , I ^ 1 = p ^ 1 X ¯ T e X ¯ T σ ^ 1 .
Here, X ¯ T = T 1 t = 1 T X t is the sample mean, and σ ^ 0 , σ ^ 1 are sample deviations of the statistics p ^ 0 e X ¯ T , p ^ 1 X ¯ T e X ¯ T , respectively. By applying some general asymptotic results related to Poisson INAR processes (cf. Weiß et al. [55]), it is shown that the central limit theorem (CLS) holds for these statistics; that is,
I ^ 0 , I ^ 1 d N ( 0 , 1 ) , T + .
Thus, by applying a simple testing procedure based on the standard Gaussian distribution, hypothesis H 0 can be verified. The test results are shown in the lower part of Table 3, from which it is evident that in both cases hypothesis H 0 was rejected. Thus, the observed count data series can be modeled using the ZOIPS-INAR(1) process.
In a similar way, for the appropriate INAR(1) process with geometric PS-innovations, one obtains:
p 0 : = Ψ X ( 0 ; λ ) = k = 0 1 θ 1 θ 1 + α k ( u 1 ) | u = 0 = exp k = 0 ln 1 + θ 1 + θ α k = exp j = 1 1 j ( 1 α j ) θ 1 θ j , p 1 : = P X t = 1 = Ψ X ( 0 ; λ ) u = j = 1 1 1 α j θ 1 θ j exp j = 1 1 j ( 1 α j ) θ 1 θ j .
It can also be seen that the earlier expressions are not in closed form. However, they can be computed approximately with an arbitrary precision.
The estimated values of parameters for both of those PS innovations (Poisson and geometric) are shown in the upper parts of Table 4 and Table 5. In order to compare the ZOIPS-INAR(1) model with the ordinary INAR(1) model, the same estimation procedures were used to fit the observed data, assuming ϕ 0 = ϕ 1 = 0 . Note that the INAR(1) model has only two parameters ( θ , α ) whose estimated values are also shown in Table 4 and Table 5. Thereby, the estimated values of the parameter α obtained using the CLS method, as the first sample correlation, are the same for both models. As in the previous simulation study, the PGF estimates were calculated using the previously described two-step estimation procedure. It is worth noting that the estimated parameter values obtained by applying both estimation methods are quite close. At the same time, the relatively small estimated values of the parameter θ are a consequence of the other small emphasized frequencies of the observed time series, which are not zero and one (and can also be seen in the following Figure 6). However, the estimated values of the parameter α are close to 1. This can also be explained by the previously described properties of the observed series, i.e., its significantly high first correlation. Finally, note that by using Equation (7), and similarly as in Mohammadi et al. [22], innovations of the ZOIPS-INAR(1) process can be represented as follows:
η t = d Y t Z t + ( 1 Y t ) ε t .
Here, Y t : B ( α ) , Z t : B ( β ) are the IID Bernoulli time series, also mutually independent of the RVs ( ε t ) , and ϕ 0 = ( 1 α ) ( 1 β ) , ϕ 1 = ( 1 α ) β . Thus, in the observed series, the estimates of ϕ 0 are related to the zero-proportion p 0 , and the estimates of ϕ 1 to the one-proportion p 1 .
In addition, we analyzed the efficiency of the fit for both INAR-based models and both estimation procedures. To this end, using the estimated parameters, we generated 500 independent simulations of the INAR(1) and ZOIPS-INAR(1) time series. To check the effectiveness of the fit to real-life data, two typical goodness-of-fit statistics were calculated: the mean squared error of estimation (MSEE) and the Akaike information criterion (AIC). Average values of these statistics, obtained on the basis of the previously described simulations, are shown in the middle parts of Table 4 and Table 5. In both cases, it is noticeable that the MSEE and AIC statistics have relatively close and small estimated values, i.e., the CLS and PGF estimates apparently have similar efficiency. Nevertheless, it is certainly noticeable that the PGF estimates have slightly smaller fitting errors. The goodness-of-fit statistics are also significantly lower in the case when the ZOIPS-INAR(1) model is applied. This means that it is a more suitable stochastic model for fitting the observed real-life time series. Finally, in the case of Poisson PS-innovations, the error statistics have slightly smaller values, so they seem more suitable for fitting with the ZOPIS-INAR(1) model. Some of the mentioned facts can also be seen in Figure 6, where the empirical and fitted frequencies of both INAR-based models are shown, whereby parameter estimates obtained by the PGF method were used.
Finally, for both of INAR-based models, as well as the both estimation procedures, the forecast accuracy analysis of the obtained models was checked. To that end, the time interval from 1 January 2023 to 28 February 2023 was taken as the forecast horizon of length h = 59 . The testing procedure was conducted using the one-sided Diebold–Mariano prediction accuracy test [56]. The null hypothesis was that the time series fitted with the INAR and ZOIPS-INAR models had the same forecast accuracy, while the alternative was that the ZOIPS-INAR model had better accuracy. The test statistic, denoted as DM, along with the appropriate p-values, was computed using the R-package “forecast”, authored by Hyndman [57]. They both are shown in the lower parts of Table 4 and Table 5, and indicate that ZOIPS-INAR models have better forecast accuracy; that is, the alternative hypothesis is valid. This is particularly evident in the case of the Poisson PS-innovations, where for both CLS and PGF estimation procedures the obtained p-values are less than the significance level of p = 0.01 .

7. Conclusions

In this paper, a generalized zero-and-one inflated ZOIPS-INAR(1) process, based on power series (PS) distributed innovations, was presented. As already noted, zero-and-one inflationary data are of particular interest in contemporary real-world data research. Let us emphasize once again that there are only two contributions in this direction: the ZOINAR model with Poisson [21] and another model with Poisson–Lindley innovations [22]. The ZOIPS-INAR(1) process proposed here, according to the general form of its PS distributions, can be viewed as a generalization of the previous models. Thus, for instance, the Poisson innovations represent only a special case of the PS-distributed innovation series. In a similar way, the ZOIPS-INAR(1) process with any other PS innovations (such as, for example, geometrically distributed innovations) may be suitable for estimating and fitting different kinds of time series. However, it should be noted that there are distributions of discrete types that do not belong to the PS family, which can be a certain limitation of this model.
The stochastic properties of the ZOIPS-INAR(1) process were considered in detail, and two parameter estimation methods were proposed. As a more contemporary method, we once again emphasize the estimation procedure based on probability generating functions (PGF method). It was shown here that the PGF estimators have slightly better asymptotic properties as well as efficiency compared to the widely used conditional least-squares (CLS) estimators. Additionally, the application of the ZOIPS-INAR(1) model to fitting real-life data on the number of deaths from the COVID-19 pandemic in the Republic of Serbia was presented. In order to determine the effectiveness of the proposed ZOIPS-INAR(1) model, it was compared with the standard INAR(1) model. In doing so, it was shown that the ZOIPS-INAR(1) model provides a better fit with the considered data; that is, it has greater efficiency and predictive accuracy compared to the standard INAR(1) model. Finally, we emphasize that two particular cases of the ZOIPS-INAR(1) model were obtained, and for both of them the fitting procedures were checked with different aspects, along with the predicting accuracy. The results obtained in this way indicate the appropriateness of the proposed model in fitting zero-and-one processes, both from a theoretical and a practical point of view. This can also be a motivation for some further research.

Author Contributions

Conceptualization, V.S.S. and H.S.B.; methodology, V.S.S. and H.S.B.; software, V.S.S. and H.S.B.; validation, E.L. and N.Q.; formal analysis, V.S.S. and H.S.B.; data curation, V.S.S. and H.S.B.; writing—original draft preparation, V.S.S., H.S.B. and E.L.; writing—review and editing, E.L. and N.Q.; visualization, V.S.S. and H.S.B.; supervision, V.S.S. and E.L.; project administration, E.L. and N.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R376), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Data Availability Statement

Official data on deaths from the disease COVID-19 can be found on the website of the World Health Organization (https://covid19.who.int/data, accessed on 19 March 2023).

Acknowledgments

The authors gratefully acknowledge Princess Nourah bint Abdulrahman University Researchers’ Supporting Project number (PNURSP2023R376), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia, for the financial support for this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Johnson, N.L.; Kemp, A.W.; Kotz, S. Univariate Discrete Distributions, 3rd ed.; Wiley Series in Probability and Statistics; John Wiley & Sons: Washington, DC, USA, 2005. [Google Scholar]
  2. Xu, H.-Y.; Xie, M.; Goh, T.N.; Fu, X. A Model for Integer–Valued Time Series with Conditional Overdispersion. Comput. Stat. Data Anal. 2012, 56, 4229–4242. [Google Scholar] [CrossRef]
  3. Weiß, C.H.; Pollet, P.K. Binomial Autoregressive Processes with Density-Dependent Thinning. J. Time Ser. Anal. 2014, 35, 115–132. [Google Scholar] [CrossRef]
  4. Martin, V.L.; Tremayne, A.R.; Jung, R.C. Efficient Method of Moments Estimators for Integer Time Series Models. J. Time Ser. Anal. 2014, 35, 491–516. [Google Scholar] [CrossRef]
  5. Graziadei, H.; Lijoi, A.; Lopes, H.F.; Marques, F.P.C.; Prünster, I. Prior Sensitivity Analysis in a Semi-Parametric Integer-Valued Time Series Model. Entropy 2020, 22, 69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Al-Osh, M.A.; Alzaid, A.A. First-Order Integer-Valued Autoregressive (INAR(1)) Process. J. Time Ser. Anal. 1987, 8, 261–275. [Google Scholar] [CrossRef]
  7. Silva, M.E.; Oliveira, V.L. Difference Equations for the Higher-Order Moments and Cumulants of the INAR(1) Model. J. Time Ser. Anal. 2004, 25, 317–333. [Google Scholar] [CrossRef]
  8. Jung, R.C.; Ronning, G.; Tremayne, A.R. Estimation in Conditional First Order Autoregression with Discrete Support. Stat. Pap. 2005, 46, 195–224. [Google Scholar] [CrossRef]
  9. Ristić, M.; Bakouch, H.S.; Nastić, A. A New Geometric First-Order Integer-Valued Autoregressive (NGINAR(1)) Process. J. Stat. Plan. Infer. 2009, 139, 2218–2226. [Google Scholar] [CrossRef]
  10. Schweer, S.; Weiß, C.H. Compound Poisson INAR(1) Processes: Stochastic Properties and Testing for Over-Dispersion. Comput. Stat. Data Anal. 2014, 77, 267–284. [Google Scholar] [CrossRef]
  11. Asgharzadeh, A.; Bakouch, H.S.; Nadarajah, S.; Sharafi, F. A New Weighted Lindley Distribution with Application. Braz. J. Probab. Stat. 2016, 30, 1–27. [Google Scholar]
  12. Mohammadpour, M.; Bakouch, H.S.; Shirozhan, M. Poisson–Lindley INAR(1) Model with Applications. Braz. J. Probab. Stat. 2018, 32, 262–280. [Google Scholar] [CrossRef]
  13. Bermúdez, L.; Karlis, D. Multivariate INAR(1) Regression Models Based on the Sarmanov Distribution. Mathematics 2021, 9, 505. [Google Scholar] [CrossRef]
  14. Li, Q.; Chen, H.; Liu, X. A New Bivariate Random Coefficient INAR(1) Model with Applications. Symmetry 2022, 14, 39. [Google Scholar] [CrossRef]
  15. Maya, R.; Chesneau, C.; Krishna, A.; Irshad, M.R. Poisson Extended Exponential Distribution with Associated INAR(1) Process and Applications. Stats 2022, 5, 755–772. [Google Scholar] [CrossRef]
  16. Maya, R.; Irshad, M.R.; Chesneau, C.; Nitin, S.L.; Shibu, D.S. On Discrete Poisson–Mirra Distribution: Regression, INAR(1) Process and Applications. Axioms 2022, 11, 193. [Google Scholar] [CrossRef]
  17. Khoo, W.C.; Ong, S.H.; Atanu, B. Coherent Forecasting for a Mixed Integer-Valued Time Series Model. Mathematics 2022, 10, 2961. [Google Scholar] [CrossRef]
  18. Saito, M.Y.; Rodrigues, J. A Bayesian Analysis of Zero and One Inflated Distributions. Rev. Mat. Estatíst. 2005, 23, 47–57. [Google Scholar]
  19. Zhang, C.; Tian, G.; Ng, K. Properties of the Zero-and-One Inflated Poisson Distribution and Likelihood-Based Inference Methods. Stat. Interface 2016, 9, 11–32. [Google Scholar] [CrossRef]
  20. Zhang, C.; Tian, G.-L.; Yuen, K.C.; Wu, Q.; Li, T. Multivariate Zero-and-One Inflated Poisson Model with Applications. J. Comput. Appl. Math. 2020, 365, 112356. [Google Scholar] [CrossRef]
  21. Qi, X.; Li, Q.; Zhu, F. Modeling Time Series of Count with Excess Zeros and Ones Based on INAR(1) Model with Zero-and-One Inflated Poisson Innovations. J. Comput. Appl. Math. 2019, 346, 572–590. [Google Scholar] [CrossRef]
  22. Mohammadi, Z.; Sajjadnia, Z.; Bakouch, H.S.; Sharafi, M. Zero-and-One Inflated Poisson–Lindley INAR(1) Process for Modelling Count Time Series with Extra Zeros and Ones. J. Stat. Comput. Simulat. 2022, 92, 2018–2040. [Google Scholar] [CrossRef]
  23. Franke, J.; Seligmann, T.H. Conditional Maximum Likelihood Estimates for INAR(1) Processes and Their Application to Modelling Epileptic Seizure Counts; Developments in Time Series; Chapman & Hall: London, UK, 1993; pp. 310–330. [Google Scholar]
  24. Du, J.-G.; Li, Y. The Integer-Valued Autoregressive (INAR(p)) Model. J. Time Ser. Anal. 1991, 12, 129–142. [Google Scholar]
  25. Latour, A. Existence and Stochastic Structure of a Non-negative Integer-valued Autoregressive Process. J. Time Ser. Anal. 1998, 19, 439–455. [Google Scholar] [CrossRef]
  26. Silva, I.; Silva, M.E. Asymptotic Distribution of the Yule-Walker Estimator for INAR(p) Processes. Stat. Prob. Lett. 2006, 76, 1655–1663. [Google Scholar] [CrossRef]
  27. Esquìvel, M.L. Some Applications of Probability Generating Function Based Methods to Statistical Estimation. Discuss. Math. 2009, 29, 131–153. [Google Scholar] [CrossRef]
  28. Stojanović, V.; Randjelović, D.; Kuk, K. Noise-Indicator Non-negative Integer-Valued Autoregressive Time Series of the First Order. Braz. J. Probab. Stat. 2018, 32, 147–171. [Google Scholar] [CrossRef]
  29. Stojanović, V.; Ljajko, E.; Tošić, M. Parameters Estimation in Non-Negative Integer-Valued Time Series: Approach Based on Probability Generating Functions. Axioms 2023, 12, 112. [Google Scholar] [CrossRef]
  30. Bourguignon, M.; Vasconcellos, K.L.P. First Order Non-Negative Integer Valued Autoregressive Processes with Power Series Innovations. Braz. J. Probab. Stat. 2015, 29, 71–93. [Google Scholar] [CrossRef]
  31. Weiß, C.H. Thinning Operations for Modelling Time Series of Counts—A Survey. Adv. Statist. Anal. 2008, 92, 319–341. [Google Scholar] [CrossRef]
  32. Kella, O.; Löpker, A. On Binomial Thinning and Mixing. Indag. Math. 2022, in press. [Google Scholar] [CrossRef]
  33. Alzaid, A.A.; Al-Osh, M.A. An Integer-Valued pth-order Autoregressive Structure (INAR(p)) Process. J. App. Prob. 1990, 27, 314–324. [Google Scholar] [CrossRef]
  34. Jazi, M.A.; Jones, G.; Lai, C.D. First-Order Integer-Valued Ar Processes with Zero-Inflated Poisson Innovations. J. Time Ser. Anal. 2012, 33, 954–963. [Google Scholar] [CrossRef]
  35. Li, C.; Wang, D.; Zhang, H. First-Order Mixed Integer-Valued Autoregressive Processes with Zero-Inflated Generalized Power Series Innovations. J. Korean Stat. Soc. 2015, 44, 232–246. [Google Scholar] [CrossRef]
  36. Knight, J.L.; Yu, J. Empirical Characteristic Function in Time Series Estimation. Econ. Theory 2002, 18, 691–721. [Google Scholar] [CrossRef] [Green Version]
  37. Yu, J. Empirical Characteristic Function Estimation and Its Applications. Econ. Rev. 2004, 23, 93–123. [Google Scholar] [CrossRef]
  38. Newey, W.K.; McFadden, D. Large Sample Estimation and Hypothesis Testing. In Handbook of Econometrics; Elsevier: Amsterdam, The Netherlands, 1994; Volume 4. [Google Scholar]
  39. Stojanović, V.; Milovanović, G.V.; Jelić, G. Distributional Properties and Parameters Estimation of GSB Process: An Approach Based on Characteristic Functions. ALEA Lat. Am. J. Probab. Math. Stat. 2016, 13, 835–861. [Google Scholar] [CrossRef]
  40. Stojanović, V.; Popović, B.Č.; Milovanović, G.V. The Split-SV Model. Comput. Statist. Data Anal. 2016, 100, 560–568. [Google Scholar] [CrossRef]
  41. Billingsley, P. Probability and Measure; John Wiley & Sons: New York, NY, USA, 1995. [Google Scholar]
  42. Klimko, L.A.; Nelson, P.I. On Conditional Least Squares Estimation for Stochastic Processes. Ann. Stat. 1978, 6, 629–642. [Google Scholar] [CrossRef]
  43. Tjøstheim, D. Estimation in Non-linear Time Series Models. Stochastic Process. Appl. 1986, 21, 251–273. [Google Scholar] [CrossRef] [Green Version]
  44. Cvetković, A.S.; Milovanović, G.V. The Mathematica Package “Orthogonal Polynomials”. Facta Univ. Ser. Math. Inform. 2004, 19, 17–36. [Google Scholar]
  45. Lange, K. Numerical Analysis for Statisticians; Statistics and Computing; Springer: New York, NY, USA, 2001. [Google Scholar]
  46. Gross, L. Tests for Normality. R Package Version 1.0-2. 2013. Available online: http://CRAN.R-project.org/package=nortest (accessed on 28 February 2023).
  47. Vaz, S.; Torres, D.F.M. A Discrete-Time Compartmental Epidemiological Model for COVID-19 with a Case Study for Portugal. Axioms 2021, 10, 314. [Google Scholar] [CrossRef]
  48. Ghosh, S.; Volpert, V.; Banerjee, M. An Epidemic Model with Time Delay Determined by the Disease Duration. Mathematics 2022, 10, 2561. [Google Scholar] [CrossRef]
  49. Sivakumar, B.; Deepthi, B. Complexity of COVID-19 Dynamics. Entropy 2022, 24, 50. [Google Scholar] [CrossRef] [PubMed]
  50. Hassan, S.M.; Riveros Gavilanes, J.M. First to React Is the Last to Forgive: Evidence from the Stock Market Impact of COVID 19. J. Risk Financ. Manag. 2021, 14, 26. [Google Scholar] [CrossRef]
  51. Zakharov, V.; Balykina, Y.; Ilin, I.; Tick, A. Forecasting a New Type of Virus Spread: A Case Study of COVID-19 with Stochastic Parameters. Mathematics 2022, 10, 3725. [Google Scholar] [CrossRef]
  52. Jovanović, M.; Stojanović, V.; Kuk, K.; Popović, B.; Čisar, P. Asymptotic Properties and Application of GSB Process: A Case Study of the COVID-19 Dynamics in Serbia. Mathematics 2022, 10, 3849. [Google Scholar] [CrossRef]
  53. Stapper, M. Count Data Time Series Modelling in Julia—The CountTimeSeries.jl Package and Applications. Entropy 2021, 23, 666. [Google Scholar] [CrossRef]
  54. World Health Organization. Available online: https://covid19.who.int/data (accessed on 3 March 2023).
  55. Weiß, C.H.; Homburg, A.; Puig, P. Testing for Zero Inflation and Overdispersion in INAR(1) Models. Stat. Pap. 2019, 60, 823–848. [Google Scholar] [CrossRef]
  56. Diebold, F.X.; Mariano, R.S. Comparing Predictive Accuracy. J. Bus. Econ. Stat. 1995, 13, 253–263. [Google Scholar]
  57. Hyndman, R. Forecasting Functions for Time Series and Linear Models. R Package Version 7.1. 2016. Available online: http://CRAN.R-project.org/package=forecast (accessed on 3 March 2023).
Figure 1. (a) Dynamics of the ZOIPS series and ZOIPS-INAR(1) process. (b) Empirical frequency distributions of both time series. (Parameters values are: α = 0.5 , θ = 5 , ϕ 0 = ϕ 1 = 0.45 ).
Figure 1. (a) Dynamics of the ZOIPS series and ZOIPS-INAR(1) process. (b) Empirical frequency distributions of both time series. (Parameters values are: α = 0.5 , θ = 5 , ϕ 0 = ϕ 1 = 0.45 ).
Mathematics 11 01772 g001
Figure 2. 3D plots of the two-dimensional PGF (a) and the corresponding EPGF (b) innovations are the ZOIPS series with a geometric PS distribution. (Parameters values are θ = α = 0.5 , ϕ 0 = ϕ 1 = 0.4 ).
Figure 2. 3D plots of the two-dimensional PGF (a) and the corresponding EPGF (b) innovations are the ZOIPS series with a geometric PS distribution. (Parameters values are θ = α = 0.5 , ϕ 0 = ϕ 1 = 0.4 ).
Mathematics 11 01772 g002
Figure 3. Histograms of empirical distributions of the PGF estimated parameters of the ZOIPS-INAR(1) process with zero-and-one inflated Poisson PS innovations. (Parameters are the same as in Table 2).
Figure 3. Histograms of empirical distributions of the PGF estimated parameters of the ZOIPS-INAR(1) process with zero-and-one inflated Poisson PS innovations. (Parameters are the same as in Table 2).
Mathematics 11 01772 g003
Figure 4. Histograms of empirical distributions of the estimated parameters of the ZOIPS-INAR(1) process with zero-and-one inflated geometric PS innovations. (Parameters are the same as in Table 2).
Figure 4. Histograms of empirical distributions of the estimated parameters of the ZOIPS-INAR(1) process with zero-and-one inflated geometric PS innovations. (Parameters are the same as in Table 2).
Mathematics 11 01772 g004
Figure 5. Graph above: dynamics of the number of deaths caused by the disease COVID-19. Graphs below: autocorrelation (a) and partial autocorrelation (b) of the observed time series.
Figure 5. Graph above: dynamics of the number of deaths caused by the disease COVID-19. Graphs below: autocorrelation (a) and partial autocorrelation (b) of the observed time series.
Mathematics 11 01772 g005
Figure 6. Frequency distributions of the observed time series and data fitted by INAR(1) and ZOIPS-INAR(1) models: (a) Poisson-distributed innovations; (b) geometric-distributed innovations.
Figure 6. Frequency distributions of the observed time series and data fitted by INAR(1) and ZOIPS-INAR(1) models: (a) Poisson-distributed innovations; (b) geometric-distributed innovations.
Mathematics 11 01772 g006
Table 1. Some specific PS distributions, along with their over-dispersion indices and PGFs.
Table 1. Some specific PS distributions, along with their over-dispersion indices and PGFs.
Distributions S a ( x ) ( 0 , R ) f ( θ ) μ ε D ε ( θ ) Ψ ε ( u ; θ ) R / θ
1. Bernoulli 0 , 1 1 ( 0 , ) 1 + θ θ 1 + θ θ 2 ( 1 + θ ) 2 1 + θ u 1 + θ
2. Binomial 0 , , n n x ( 0 , ) ( 1 + θ ) n n θ 1 + θ n θ 2 ( 1 + θ ) 2 1 + θ u 1 + θ n
3. Poisson 0 , , 1 x ! ( 0 , ) exp ( θ ) θ 0 exp θ ( u 1 )
4. Geometric 0 , , 1 ( 0 , 1 ) 1 1 θ θ 1 θ θ 2 ( 1 θ ) 2 1 θ 1 θ u 1 / θ
5. Negative binomial 0 , , Γ ( x + n ) x ! Γ ( n ) ( 0 , 1 ) 1 ( 1 θ ) n n θ ( 1 θ ) n n θ 2 ( 1 θ ) 2 1 θ 1 θ u n 1 / θ
6. Pascal n , , x 1 n 1 ( 0 , 1 ) θ 1 θ n n 1 θ n ( 2 θ 1 ) ( 1 θ ) 2 ( 1 θ ) u 1 θ u n 1 / θ
Table 2. Estimated values of parameters of the ZOIPS-INAR ( 1 ) process. (True parameters are θ = α = 0.5 , ϕ 0 = ϕ 1 = 0.35 ).
Table 2. Estimated values of parameters of the ZOIPS-INAR ( 1 ) process. (True parameters are θ = α = 0.5 , ϕ 0 = ϕ 1 = 0.35 ).
SampleCLS EstimatesPGF Estimates S T ( 2 )
θ ˜ α ˜ ϕ ˜ 0 ϕ ˜ 1 θ ^ α ^ ϕ ^ 0 ϕ ^ 1
Min.0.21370.41160.12090.15150.29260.33030.09750.13943.45 × 10 6
Mean0.48460.49800.33800.34120.50540.50400.35910.35741.18 × 10 4
Max.0.76240.58410.72860.71710.71550.63520.75150.60841.84 × 10 3
MSEE0.00910.00120.00820.00800.00570.00490.00600.0022
Poisson A D 2.0393 **0.34350.51780.9376 *0.57250.35940.7822 *0.4727
(p-value)(3.34 × 10 5 )(0.2670)(0.1870)(0.0173)(0.1361)(0.4474)(0.0418)(0.2414) –
Min.0.15220.39550.10280.12710.30430.29730.15150.17972.52 × 10 7
Mean0.51310.49660.35810.33870.49520.49630.35460.34834.04 × 10 6
Max.0.64870.57640.69380.67020.69730.67680.59580.52941.28 × 10 5
MSEE0.01370.00110.01750.01740.00740.00530.00860.0064
Geometric A D 1.0235 *0.28590.70600.51080.66280.18670.70690.4753
(p-value)(0.0106)(0.6227)(0.0646)(0.1945)(0.0826)(0.9038)(0.0643)(0.2379) –
* p < 0.05; ** p < 0.01.
Table 3. Summary statistics of the number of deaths per day in the Republic of Serbia, with zero-and-one inflation testing.
Table 3. Summary statistics of the number of deaths per day in the Republic of Serbia, with zero-and-one inflation testing.
StatisticsEstimated Values
Sample size1094
Minimum0
Mode0
Median9
Mean16.01
Maximum79
St. deviation18.21
Variance331.57
Skewness1.269
Kurtosis3.4367
Zero-count104
Zero-proportion0.0951
Zero-statistics ( I ^ 0 )5.601 **
(p-value)(2.10 × 10 6 )
Ones-count102
Ones-proportion0.0932
Ones-statistics ( I ^ 1 )5.492 **
(p-value)(3.54 × 10 6 )
** p < 0.01.
Table 4. Estimated fitting parameters of the considered data with the INAR(1) and ZOIPS-INAR(1) models, along with their goodness-of-fit and predictive statistics. (Poisson distribution is assumed for PS-distributed innovations).
Table 4. Estimated fitting parameters of the considered data with the INAR(1) and ZOIPS-INAR(1) models, along with their goodness-of-fit and predictive statistics. (Poisson distribution is assumed for PS-distributed innovations).
Parameters/StatisticsCLS EstimatesPGF Estimates
INARZOIPS-INARINARZOIPS-INAR
θ 0.17310.14200.17680.1592
α 0.98900.98900.94310.9902
ϕ 0 0.07090.0775
ϕ 1 0.14920.1416
S T 2 3.92 × 10 2 3.09 × 10 2 1.17 × 10 2 6.96 × 10 4
MSEE8.61 × 10 3 5.98 × 10 3 6.79 × 10 3 5.96 × 10 3
AIC4.83934.24954.82654.2489
D M 2.4269 * * 2.3312 * *
(p-value)(7.69 × 10 3 )(9.96 × 10 3 )
** p < 0.01.
Table 5. Estimated fitting parameters of the considered data with the INAR(1) and ZOIPS-INAR(1) models, along with their goodness-of-fit and predictive statistics. (Geometric distribution is assumed for PS-distributed innovations).
Table 5. Estimated fitting parameters of the considered data with the INAR(1) and ZOIPS-INAR(1) models, along with their goodness-of-fit and predictive statistics. (Geometric distribution is assumed for PS-distributed innovations).
Parameters/StatisticsCLS EstimatesPGF Estimates
INARZOIPS-INARINARZOIPS-INAR
θ 0.14610.09190.22360.1075
α 0.98900.98900.90990.9908
ϕ 0 0.10680.1169
ϕ 1 0.12310.0859
S T 2 2.90 × 10 2 2.78 × 10 2 8.99 × 10 3 9.91 × 10 5
MSEE1.11 × 10 2 6.81 × 10 3 6.84 × 10 3 6.78 × 10 3
AIC4.86864.83014.83554.8291
D M 1.7961 *2.3446 **
(p-value)(3.64 × 10 2 )(9.61 × 10 3 )
* p < 0.05; ** p < 0.01.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stojanović, V.S.; Bakouch, H.S.; Ljajko, E.; Qarmalah, N. Zero-and-One Integer-Valued AR(1) Time Series with Power Series Innovations and Probability Generating Function Estimation Approach. Mathematics 2023, 11, 1772. https://doi.org/10.3390/math11081772

AMA Style

Stojanović VS, Bakouch HS, Ljajko E, Qarmalah N. Zero-and-One Integer-Valued AR(1) Time Series with Power Series Innovations and Probability Generating Function Estimation Approach. Mathematics. 2023; 11(8):1772. https://doi.org/10.3390/math11081772

Chicago/Turabian Style

Stojanović, Vladica S., Hassan S. Bakouch, Eugen Ljajko, and Najla Qarmalah. 2023. "Zero-and-One Integer-Valued AR(1) Time Series with Power Series Innovations and Probability Generating Function Estimation Approach" Mathematics 11, no. 8: 1772. https://doi.org/10.3390/math11081772

APA Style

Stojanović, V. S., Bakouch, H. S., Ljajko, E., & Qarmalah, N. (2023). Zero-and-One Integer-Valued AR(1) Time Series with Power Series Innovations and Probability Generating Function Estimation Approach. Mathematics, 11(8), 1772. https://doi.org/10.3390/math11081772

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop