Integer-valued Autoregressive Model Research Articles

Bivariate Random Coefficient Integer-Valued Autoregressive Model Based on a ρ-Thinning Operator

While overdispersion is a common phenomenon in univariate count time series data, its exploration within bivariate contexts remains limited. To fill this gap, we propose a bivariate integer-valued autoregressive model. The model leverages a modified binomial thinning operator with a dispersion parameter ρ and integrates random coefficients. This approach combines characteristics from both binomial and negative binomial thinning operators, thereby offering a flexible framework capable of generating counting series exhibiting equidispersion, overdispersion, or underdispersion. Notably, our model includes two distinct classes of first-order bivariate geometric integer-valued autoregressive models: one class employs binomial thinning (BVGINAR(1)), and the other adopts negative binomial thinning (BVNGINAR(1)). We establish the stationarity and ergodicity of the model and estimate its parameters using a combination of the Yule–Walker (YW) and conditional maximum likelihood (CML) methods. Furthermore, Monte Carlo simulation experiments are conducted to evaluate the finite sample performances of the proposed estimators across various parameter configurations, and the Anderson-Darling (AD) test is employed to assess the asymptotic normality of the estimators under large sample sizes. Ultimately, we highlight the practical applicability of the examined model by analyzing two real-world datasets on crime counts in New South Wales (NSW) and comparing its performance with other popular overdispersed BINAR(1) models.

On Comparing and Assessing Robustness of Some Popular Non-Stationary BINAR(1) Models

Intra-day transactions of stocks from competing firms in the financial markets are known to exhibit significant volatility and over-dispersion. This paper proposes some bivariate integer-valued auto-regressive models of order 1 (BINAR(1)) that are useful to analyze such financial series. These models were constructed under both time-variant and time-invariant conditions to capture features such as over-dispersion and non-stationarity in time series of counts. However, the quest for the most robust BINAR(1) models is still on. This paper considers specifically the family of BINAR(1)s with a non-diagonal cross-correlation structure and with unpaired innovation series. These assumptions relax the number of parameters to be estimated. Simulation experiments are performed to assess both the consistency of the estimators and the robust behavior of the BINAR(1)s under mis-specified innovation distribution specifications. The proposed BINAR(1)s are applied to analyze the intra-day transaction series of AstraZeneca and Ericsson. Diagnostic measures such as the root mean square errors (RMSEs) and Akaike information criteria (AICs) are also considered. The paper concludes that the BINAR(1)s with negative binomial and COM–Poisson innovations are among the most suitable models to analyze over-dispersed intra-day transaction series of stocks.

The Circumstance-Driven Bivariate Integer-Valued Autoregressive Model.

The novel circumstance-driven bivariate integer-valued autoregressive (CuBINAR) model for non-stationary count time series is proposed. The non-stationarity of the bivariate count process is defined by a joint categorical sequence, which expresses the current state of the process. Additional cross-dependence can be generated via cross-dependent innovations. The model can also be equipped with a marginal bivariate Poisson distribution to make it suitable for low-count time series. Important stochastic properties of the new model are derived. The Yule-Walker and conditional maximum likelihood method are adopted to estimate the unknown parameters. The consistency of these estimators is established, and their finite-sample performance is investigated by a simulation study. The scope and application of the model are illustrated by a real-world data example on sales counts, where a soap product in different stores with a common circumstance factor is investigated.

Comparison of estimation and prediction methods for a zero-inflated geometric INAR(1) process with random coefficients

This study explores zero-inflated count time series models used to analyze data sets with characteristics such as overdispersion, excess zeros, and autocorrelation. Specifically, we investigate the ZIGINAR RC ( 1 ) process, a first-order stationary integer-valued autoregressive model with random coefficients and a zero-inflated geometric marginal distribution. Our focus is on examining various estimation and prediction techniques for this model. We employ estimation methods, including Whittle, Taper Spectral Whittle, Maximum Empirical Likelihood, and Sieve Bootstrap estimators for parameter estimation. Additionally, we propose forecasting approaches, such as median, Bayesian, and Sieve Bootstrap methods, to predict future values of the series. We assess the performance of these methods through simulation studies and real-world data analysis, finding that all methods perform well, providing 95% highest predicted probability intervals that encompass the observed data. While Bayesian and Bootstrap methods require more time for execution, their superior predictive accuracy justifies their use in forecasting.

A Two-Step Estimation Method for a Time-Varying INAR Model

This paper proposes a new time-varying integer-valued autoregressive (TV-INAR) model with a state vector following a logistic regression structure. Since the autoregressive coefficient in the model is time-dependent, the Kalman-smoothed method is applicable. Some statistical properties of the model are established. To estimate the parameters of the model, a two-step estimation method is proposed. In the first step, the Kalman-smoothed estimation method, which is suitable for handling time-dependent systems and nonstationary stochastic processes, is utilized to estimate the time-varying parameters. In the second step, conditional least squares is used to estimate the parameter in the error term. This proposed method allows estimating the parameters in the nonlinear model and deriving the analytical solutions. The performance of the estimation method is evaluated through simulation studies. The model is then validated using actual time series data.

New One-Parameter Over-Dispersed Discrete Distribution and Its Application to the Nonnegative Integer-Valued Autoregressive Model of Order One

Count data arise in inference, modeling, prediction, anomaly detection, monitoring, resource allocation, evaluation, and performance measurement. This paper focuses on a one-parameter discrete distribution obtained by compounding the Poisson and new X-Lindley distributions. The probability-generating function, moments, skewness, kurtosis, and other properties are derived in the closed form. The maximum likelihood method, method of moments, least squares method, and weighted least squares method are used for parameter estimation. A simulation study is carried out. The proposed distribution is applied as the innovation in an INAR(1) process. The importance of the proposed model is confirmed through the analysis of two real datasets.

Change-points analysis for generalized integer-valued autoregressive model via minimum description length principle

This article considers the problem of modeling a class of count time series with multiple change-points using segmented generalized integer-valued autoregressive (S-GINAR) processes. The minimum description length principle (MDL) is applied to study the statistical inference for the S-GINAR model, and the consistency results of the MDL model selection procedure are established respectively under the condition of known and unknown number of change-points. To find the “best” combination of the number of change-points, the locations of change-points, the order of each segment and its parameters, a genetic algorithm with simulated annealing is implemented to solve this difficult optimization problem. In particular, the simulated annealing process makes up for the precocious problem of the traditional genetic algorithm. Numerical results from simulation experiments and three examples of real data analysis show that the procedure has excellent empirical properties.

Artificial Neural Networks and Time Series of Counts: A Class of Nonlinear INGARCH Models

Abstract Time series of counts are frequently analyzed using generalized integer-valued autoregressive models with conditional heteroskedasticity (INGARCH). These models employ response functions to map a vector of past observations and past conditional expectations to the conditional expectation of the present observation. In this paper, it is shown how INGARCH models can be combined with artificial neural network (ANN) response functions to obtain a class of nonlinear INGARCH models. The ANN framework allows for the interpretation of many existing INGARCH models as a degenerate version of a corresponding neural model. Details on maximum likelihood estimation, marginal effects and confidence intervals are given. The empirical analysis of time series of bounded and unbounded counts reveals that the neural INGARCH models are able to outperform reasonable degenerate competitor models in terms of the information loss.

Robust estimation for general integer-valued autoregressive models based on the exponential-polynomial divergence

In this study, we develop a robust estimator for integer-valued one-parameter exponential family autoregressive models, named general integer-valued autoregressive models. This model accommodates a broad class of integer-valued time series models. In particular, we propose a robust estimation method that minimizes the exponential-polynomial divergence (EPD) belonging to the Brègman divergence family. EPD subsumes the density power divergence (DPD), which has been extensively studied by many authors for the past decades. Under regularity conditions, the minimum EPD estimator (MEPDE) is shown to be consistent and asymptotically normal. Comparing the performance of MEPDE with the minimum DPD estimator, we substantiate the validity of MEPDE through a simulation study and real data analysis.

Two-Threshold-Variable Integer-Valued Autoregressive Model

In the past, most threshold models considered a single threshold variable. However, for some practical applications, models with two threshold variables may be needed. In this paper, we propose a two-threshold-variable integer-valued autoregressive model based on the binomial thinning operator and discuss some of its basic properties, including the mean, variance, strict stationarity, and ergodicity. We consider the conditional least squares (CLS) estimation and discuss the asymptotic normality of the CLS estimator under the known and unknown threshold values. The performances of the CLS estimator are compared via simulation studies. In addition, two real data sets are considered to underline the superior performance of the proposed model.

Robust estimation for bivariate integer-valued autoregressive models based on minimum density power divergence

This study considers a robust estimation method for bivariate integer-valued autoregressive (BINAR) models using a minimum density power divergence (MDPD) scheme. We first address this problem within a general framework of time series models to verify the strong consistency and asymptotic normality of the MDPD estimator (MDPDE) under regularity conditions. Based on these results, we assert the same asymptotic properties of MDPDE for BINAR models and propose an outlier detection rule based on MDPDE. To assess the performance of the proposed method, we conduct a simulation study. Real data analysis is also carried out for illustration, using the number of monthly earthquakes in the United States.

New bivariate Poisson extended exponential distributions and associated BINAR(1) processes with applications

This study proposes two types of bivariate Poisson extended exponential distributions: the basic bivariate Poisson extended exponential distribution and the Sarmanov-based bivariate Poisson extended exponential distribution. The two bivariate Poisson extended exponential distributions are then introduced as joint innovation distributions in a bivariate first-order integer-valued autoregressive process based on binomial thinning. The model parameters are estimated using maximum likelihood in basic bivariate Poisson extended exponential and Sarmanov-based bivariate Poisson extended exponential distributions. Conditional maximum likelihood is applied to the bivariate first-order integer-valued autoregressive process. Simulation experiments are used to evaluate the performance of small and large samples. In addition, the newly developed bivariate first-order integer-valued autoregressive models are then applied to the Pittsburgh crime series data and candies data. We show that they fit better than existing bivariate first-order integer-valued autoregressive models described in the literature.

An Observation-Driven Random Parameter INAR(1) Model Based on the Poisson Thinning Operator.

This paper presents a first-order integer-valued autoregressive time series model featuring observation-driven parameters that may adhere to a particular random distribution. We derive the ergodicity of the model as well as the theoretical properties of point estimation, interval estimation, and parameter testing. The properties are verified through numerical simulations. Lastly, we demonstrate the application of this model using real-world datasets.

INAR approximation of bivariate linear birth and death process

In this paper, we propose a new type of univariate and bivariate Integer-valued autoregressive model of order one (INAR(1)) to approximate univariate and bivariate linear birth and death process with constant rates. Under a specific parametric setting, the dynamic of transition probabilities and probability generating function of INAR(1) will converge to that of birth and death process as the length of subintervals goes to 0. Due to the simplicity of Markov structure, maximum likelihood estimation is feasible for INAR(1) model, which is not the case for bivariate and multivariate birth and death process. This means that the statistical inference of bivariate birth and death process can be achieved via the maximum likelihood estimation of a bivariate INAR(1) model.

High-Order Self-excited Threshold Integer-Valued Autoregressive Model: Estimation and Testing

High-Order Self-excited Threshold Integer-Valued Autoregressive Model: Estimation and Testing

On Periodic Generalized Poisson INAR(1) Model

In this paper, we introduce a first-order Periodic Generalized Poisson Integer-Valued Autoregressive model which has been shown to be useful to describe overdispersion, equidispersion and underdispersion feature encountered in periodically correlated Integer-Valued time series. Some probabilistic and statistical properties are established, such as the periodically correlated stationarity conditions, in the first and the second moments are provided and the closed-forms of these moments are, under these conditions, derived. Moreover, the structure of the periodic autocovariance is obtained. The estimation problem is addressed through the Yule-Walker the Two-Stage Conditional Least Squares and the Conditional Maximum Likelihood methods. The performance of these methods is done through an intensive simulation study and an application on real data set is accomplished. Keywords and phrases: Periodic Generalized Poisson, Integer-Valued Autoregressive, Periodically correlated process, periodically stationary condition.

A new bivariate INAR(1) model with paired Poisson-weighted exponential distributed innovations

This paper proposes a novel Bivariate integer-valued auto-regressive model of order 1 with paired Poisson Weighted Exponential (PWE) distributed innovations which is denoted by INAR(1)-PWE with two Sarmanov and classical versions. The CML and CLS estimators of the parameters are obtained and the performance of the proposed models are assessed through some Monte Carlo simulation experiments. Also, the BINAR(1)-PWE is applied to the two real data sets and is compared with some bivariate INAR processes. The research findings commend the BINAR(1)-PWE as another suitable alternative to analyze bivariate series of counts and open the avenues to explore the Sarmanov-based bivariate models.

The balanced discrete Burr–Hatke model and mixing INAR(1) process: properties, estimation, forecasting and COVID-19 applications

The main concern of this paper is providing a flexible discrete model that captures every kind of dispersion (equi-, over- and under-dispersion). Based on the balanced discretization method, a new discrete version of Burr–Hatke distribution is introduced with the partial moment-preserving property. Some statistical properties of the new distribution are introduced, and the applicability of proposed model is evaluated by considering counting series. A new integer-valued autoregressive (INAR) process based on the mixing Pegram and binomial thinning operators with discrete Burr–Hatke innovations is introduced, which can model contagious data properly. The different estimation approaches of parameters of the new process are provided and compared through the Monte Carlo simulation scheme. The performance of the proposed process is evaluated by four data sets of the daily death counts of the COVID-19 in Austria, Switzerland, Nigeria and Slovenia in comparison with some competitor INAR(1) models, along with the Pearson residual analysis of the assessing model. The goodness of fit measures affirm the adequacy of the proposed process in modeling all COVID-19 data sets. The fundamental prediction procedures are considered for new process by classic, modified Sieve bootstrap and Bayesian forecasting methods for all COVID-19 data sets, which is concluded that the Bayesian forecasting approach provides more reliable results.

Integer-valued Autoregressive Model Based on Innovations with Discrete Exponential-Weibull Distribution

Integer-valued Autoregressive Model Based on Innovations with Discrete Exponential-Weibull Distribution

Bayesian Analysis of Spatial Model for Frequency of Tornadoes

Frequency analysis of tornadoes is very important for risk analysis and disaster control. In this paper, the annual frequency of tornadoes that occurred in the United States from 1967 to 2016 is analyzed. The simple analysis shows that frequencies of tornadoes of different sites are spatially correlated and over-dispersed. To explain the two characteristics of the data, the Bayesian hierarchical model is proposed. For comparison purposes, the Bayesian model with negative binomial distribution, Poisson distribution, Polya distribution, and first-order, non-negative, integer-valued autoregressive model with Bell innovations(BL-INAR(1)) are considered to fit the frequency of tornado. The distribution parameters of all sites are assumed to be spatially correlated, and the corresponding Bayesian hierarchical models were established. MCMC (Markov Chain Monte Carlo) method is applied to parameter estimations and relative statistical inference. By comparison of the analysis results, the negative binomial distribution is recommended to analyze the overdispersion and spatial correlation among the sites of the data. The comparison between the simulated frequencies based on the proposed model and the actual frequencies also verifies that the proposed method is a better model for the data.