Fractional approaches for the distribution of innovation sequence of INAR(1) processes

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(1)] 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 two fractional approaches to obtain the distribution of innovation processes of the INAR(1) model and show that the distribution of the innovations sequence has geometric-type distribution. These approaches are discussed in detail and illustrated through a few examples.     

Process capability analysis for serially dependent processes of Poisson counts

Process capability indices evaluate the actual compliance of a process with given external specifications in a single number. For the case of a process of independent and identically distributed Poisson counts, two types of index have been proposed and investigated in the literature. The assumption of serial independence, however, is quite unrealistic for practice. We consider the case of an underlying Poisson INAR(1) process which has an AR(1)-like autocorrelation structure. We show that the performance of the estimated indices is degraded heavily if serial dependence is ignored. Therefore, we develop approaches for estimating the process capability (both for the observation and innovation process), which explicitly consider the observed degree of autocorrelation. For this purpose, we introduce a new unbiased estimator of the innovations’ mean of a Poisson INAR(1) process and derive its exact as well as asymptotic stochastic properties. In this context, we also present new explicit expressions for the third- and fourth-order moments of a Poisson INAR(1) process. Then the capability indices and the performance of their estimators are analysed and recommendations for practice are given.     

Fully observed INAR(1) processes

The innovations of an INAR(1) process (integer-valued autoregressive) are usually assumed to be unobservable. There are, however, situations in practice, where also the innovations can be uncovered, i.e. where we are concerned with a fully observed INAR(1) process. We analyze stochastic properties of such a fully observed INAR(1) process and explore the relation between the INAR(1) model and certain metapopulation models. We show how the additional knowledge about the innovations can be used for parameter estimation, for model diagnostics, and for forecasting. Our findings are illustrated with two real-data examples.     

Flexible Bivariate INAR(1) Processes Using Copulas

Multivariate count time series data occur in many different disciplines. The class of INteger-valued AutoRegressive (INAR) processes has the great advantage to consider explicitly both the discreteness and autocorrelation characterizing this type of data. Moreover, extensions of the simple INAR(1) model to the multi-dimensional space make it possible to model more than one series simultaneously. However, existing models do not offer great flexibility for dependence modelling, allowing only for positive correlation. In this work, we consider a bivariate INAR(1) (BINAR(1)) process where cross-correlation is introduced through the use of copulas for the specification of the joint distribution of the innovations. We mainly emphasize on the parametric case that arises under the assumption of Poisson marginals. Other marginal distributions are also considered. A short application on a bivariate financial count series illustrates the model.     

Efficient estimation of auto-regression parameters and innovation distributions for semiparametric integer-valued AR(p) models

Summary.  Integer-valued auto-regressive (INAR) processes have been introduced to model non-negative integer-valued phenomena that evolve over time. The distribution of an INAR( p ) process is essentially described by two parameters: a vector of auto-regression coefficients and a probability distribution on the non-negative integers, called an immigration or innovation distribution. Traditionally, parametric models are considered where the innovation distribution is assumed to belong to a parametric family. The paper instead considers a more realistic semiparametric INAR( p ) model where there are essentially no restrictions on the innovation distribution. We provide an (semiparametrically) efficient estimator of both the auto-regression parameters and the innovation distribution.     

Estimating monotonic change in the rate and dependence parameters of INAR(1) process (Case study: IP counts data)

In this article, we consider a first-order integer-valued autoregressive (INAR(1)) model. Then, we propose change point estimators for the rate and dependence parameters in INAR(1) model using maximum likelihood estimation method when the type of change belongs to a family of monotonic changes. To monitor the process, a combined EWMA and c control chart is considered. The results show that the proposed change point estimators provide efficient estimates of the change time. At the end, to illustrate the application of the proposed estimators, a real case related to IP counts data is investigated.     

Modelling of low count heavy tailed time series data consisting large number of zeros and ones

In this paper, we construct a new mixture of geometric INAR(1) process for modeling over-dispersed count time series data, in particular data consisting of large number of zeros and ones. For some real data sets, the existing INAR(1) processes do not fit well, e.g., the geometric INAR(1) process overestimates the number of zero observations and underestimates the one observations, whereas Poisson INAR(1) process underestimates the zero observations and overestimates the one observations. Furthermore, for heavy tails, the PINAR(1) process performs poorly in the tail part. The existing zero-inflated Poisson INAR(1) and compound Poisson INAR(1) processes have the same kind of limitations. In order to remove this problem of under-fitting at one point and over-fitting at others points, we add some extra probability at one in the geometric INAR(1) process and build a new mixture of geometric INAR(1) process. Surprisingly, for some real data sets, it removes the problem of under and over-fitting over all the observations up to a significant extent. We then study the stationarity and ergodicity of the proposed process. Different methods of parameter estimation, namely the Yule-Walker and the quasi-maximum likelihood estimation procedures are discussed and illustrated using some simulation experiments. Furthermore, we discuss the future prediction along with some different forecasting accuracy measures. Two real data sets are analyzed to illustrate the effective use of the proposed model.     

Risk aggregation with dependence and overdispersion based on the compound Poisson INAR(1) process

Abstract

This paper considers an extension of the classical discrete time risk model for which the claim numbers are assumed to be temporal dependence and overdispersion. The risk model proposed is based on the first-order integer-valued autoregressive (INAR(1)) process with discrete compound Poisson distributed innovations. The explicit expression for the moment generating function of the discounted aggregate claim amount is derived. Some numerical examples are provided to illustrate the impacts of dependence and overdispersion on related quantities such as the stop-loss premium, the value at risk and the tail value at risk.     

Monitoring correlated processes with binomial marginals

Few approaches for monitoring autocorrelated attribute data have been proposed in the literature. If the marginal process distribution is binomial, then the binomial AR(1) model as a realistic and well-interpretable process model may be adequate. Based on known and newly derived statistical properties of this model, we shall develop approaches to monitor a binomial AR(1) process, and investigate their performance in a simulation study. A case study demonstrates the applicability of the binomial AR(1) model and of the proposed control charts to problems from statistical process control.     

Self-Starting control charts for monitoring logistic regression profiles

In some applications, quality engineers cannot monitor the processes at the beginning of the production process. Because the process parameters are unknown and there are not enough initial samples to estimate the process parameters. Self-starting control charts are applied to monitor processes at the start-up stages with no enough initial samples. In this paper, we propose three self-starting control charts to monitor a logistic regression profile which models the relationship between a binomial response variable and explanatory variables. Also, we compare the proposed control charts with each other through simulation studies in terms of average run length (ARL) criterion.     

Control charts for fraction nonconforming in a bivariate binomial process

Many multivariate quality control techniques are used for multivariate variable processes, but few work for multivariate attribute processes. To monitor multivariate attributes, controlling the false alarms (type I errors) and considering the correlation between attributes are two important issues. By taking into account these two issues, a new control chart is presented to monitor a bivariate binomial process. An example is illustrated for the proposed method. To evaluate the performance of the proposed method, a simulation study is conducted to compare the results with those using both the multivariate np chart and skewness reduction approaches. The results show that the correlation is taken into account in the designed chart and the overall false alarm is controlled at the nominal value. Moreover, the process shift can be quickly detected and the variable that is responsible for a signal can be determined.     

The predictive distributions of thinning‐based count processes

This paper shows that the term structure of conditional (i.e. predictive) distributions allows for closed form expression in a large family of (possibly higher order or infinite order) thinning‐based count processes such as INAR(p), INARCH(p), NBAR(p), and INGARCH(1,1). Such predictive distributions are currently often deemed intractable by the literature and existing approximation methods are usually time consuming and induce approximation errors. In this paper, we propose a Taylor's expansion algorithm for these predictive distributions, which is both exact and fast. Through extensive simulation exercises, we demonstrate its advantages with respect to existing methods in terms of the computational gain and/or precision.     

Bootstrap-based bias corrections for INAR count time series

Integer-valued autoregressive (INAR) processes form a very useful class of processes suitable to model time series of counts. Several practically relevant estimators based on INAR data are known to be systematically biased away from their population values, e.g. sample autocovariances, sample autocorrelations, or the dispersion index. We propose to do bias correction for such estimators by using a recently proposed INAR-type bootstrap scheme that is tailor-made for INAR processes, and which has been proven to be asymptotically consistent under general conditions. This INAR bootstrap allows an implementation with and without parametrically specifying the innovations' distribution. To judge the potential of corresponding bias correction, we compare these bootstraps in simulations to several competitors that include the AR bootstrap and block bootstrap. Finally, we conclude with an illustrative data application.     

Modeling spatial overdispersion via the generalized Waring point process

Modeling spatial overdispersion requires point process models with finite‐dimensional distributions that are overdisperse relative to the Poisson distribution. Fitting such models usually heavily relies on the properties of stationarity, ergodicity, and orderliness. In addition, although processes based on negative binomial finite‐dimensional distributions have been widely considered, they typically fail to simultaneously satisfy the three required properties for fitting. Indeed, it has been conjectured by Diggle and Milne that no negative binomial model can satisfy all three properties. In light of this, we change perspective and construct a new process based on a different overdisperse count model, namely, the generalized Waring (GW) distribution. While comparably tractable and flexible to negative binomial processes, the GW process is shown to possess all required properties and additionally span the negative binomial and Poisson processes as limiting cases. In this sense, the GW process provides an approximate resolution to the conundrum highlighted by Diggle and Milne.     

A geometric time series model with a new dependent Bernoulli counting series

ABSTRACT

New generalized binomial thinning operator with dependent counting series is introduced. An integer valued time series model with geometric marginals based on this thinning operator is constructed. Main features of the process are analyzed and determined. Estimation of the parameters are presented and some asymptotic properties of the obtained estimators are discussed. Behavior of the estimators is described through the numerical results. Also, model is applied on the real data set and compared to some relevant INAR(1) models.     

Goodness-of-fit tests for binomial AR(1) processes

The binomial AR(1) model describes a nonlinear process with a first-order autoregressive (AR(1)) structure and a binomial marginal distribution. To develop goodness-of-fit tests for the binomial AR(1) model, we investigate the observed marginal distribution of the binomial AR(1) process, and we tackle its autocorrelation structure. Motivated by the family of power-divergence statistics for handling discrete multivariate data, we derive the asymptotic distribution of certain categorized power-divergence statistics for the case of a binomial AR(1) process. Then we consider Bartlett's formula, which is widely used in time series analysis to provide estimates of the asymptotic covariance between sample autocorrelations, but which is not applicable when the underlying process is nonlinear. Hence, we derive a novel Bartlett-type formula for the asymptotic distribution of the sample autocorrelations of a binomial AR(1) process, which is then applied to develop tests concerning the autocorrelation structure. Simulation studies are carried out to evaluate the size and power of the proposed tests under diverse alternative process models. Several real examples are used to illustrate our methods and findings.     

Time series with discrete semistable marginals

We propose several stationary integer-valued first-order autoregressive [INAR(1)] models with discrete semistable marginals and related distributions. The corresponding first-order moving average processes are also presented.     

Convolution of Binomial and Negative Binomial Variables

This paper considers a distribution formed by convolution of binomial and negative binomial variables. The distribution has the flexibility to adapt to the model under, equi, and over dispersion. Some properties of the proposed distribution are discussed, including characterization. Three stochastic processes leading to the distribution are also considered: (1) a three-dimensional random walk; (2) a birth, death, and immigration process; and (3) a thinned stochastic process.     

EWNA charts for monitoring the mean and the autocovariances of stationary processes

In this paper various types of EWMA control charts are introduced for the simultaneous monitoring of the mean and the autocovariances. The target process is assumed to be a stationary process up to fourth-order or an ARMA process with heavy tailed innovations. The case of a Gaussian process is included in our results as well. The charts are compared within a simulation study. As a measure of the performance the average run length is taken. The target process is an ARMA (1,1) process with Student-t distributed innovations. The behavior of the charts is analyzed with respect to several out-of-control models. The best design parameters are determined for each chart. Our comparisons show that the multivariate EWMA chart applied to the residuals has the best overall performance.     

Extended Poisson INAR(1) processes with equidispersion,underdispersion and overdispersion

Real count data time series often show the phenomenon of the underdispersion and overdispersion. In this paper, we develop two extensions of the first-order integer-valued autoregressive process with Poisson innovations, based on binomial thinning, for modeling integer-valued time series with equidispersion, underdispersion, and overdispersion. The main properties of the models are derived. The methods of conditional maximum likelihood, Yule–Walker, and conditional least squares are used for estimating the parameters, and their asymptotic properties are established. We also use a test based on our processes for checking if the count time series considered is overdispersed or underdispersed. The proposed models are fitted to time series of the weekly number of syphilis cases and monthly counts of family violence illustrating its capabilities in challenging the overdispersed and underdispersed count data.