Non-stationarity in bivariate time series of counts may be induced by a number of time-varying covariates affecting the bivariate responses due to which the innovation terms of the individual series as well as the bivariate dependence structure becomes non-stationary. So far, in the existing models, the innovation terms of individual INAR(1) series and the dependence structure are assumed to be constant even though the individual time series are non-stationary. Under this assumption, the reliability of the regression and correlation estimates is questionable. Besides, the existing estimation methodologies such as the conditional maximum likelihood (CMLE) and the composite likelihood estimation are computationally intensive. To address these issues, this paper proposes a BINAR(1) model where the innovation series follow a bivariate Poisson distribution under some non-stationary distributional assumptions. The method of generalized quasi-likelihood (GQL) is used to estimate the regression effects while the serial and bivariate correlations are estimated using a robust moment estimation technique. The application of model and estimation method is made in the simulated data. The GQL method is also compared with the CMLE, generalized method of moments (GMM) and generalized estimating equation (GEE) approaches where through simulation studies, it is shown that GQL yields more efficient estimates than GMM and equally or slightly more efficient estimates than CMLE and GEE.

This paper introduces a first-order integer-valued autoregressive process with a new innovation distribution, shortly INARPQX(1) process. A new innovation distribution is obtained by mixing Poisson distribution with quasi-xgamma distribution. The statistical properties and estimation procedure of a new distribution are studied in detail. The parameter estimation of INARPQX(1) process is discussed with two estimation methods: conditional maximum likelihood and Yule-Walker. The proposed INARPQX(1) process is applied to time series of the monthly counts of earthquakes. The empirical results show that INARPQX(1) process is an important process to model over-dispersed time series of counts and can be used to predict the number of earthquakes with a magnitude greater than four.

This article proposes a novel non-stationary BINMA time series model by extending two INMA processes where their innovation series follow the bivariate Poisson under time-varying moment assumptions. This article also demonstrates, through simulation studies, the use and superiority of the generalized quasi-likelihood (GQL) approach to estimate the regression effects, which is computationally less complicated as compared to conditional maximum likelihood estimation (CMLE) and the feasible generalized least squares (FGLS). The serial and bivariate dependence correlations are estimated by a robust method of moments.

This article proposes a bivariate integer-valued autoregressive time-series model of order 1 (BINAR(1) with COM–Poisson marginals to analyze a pair of non stationary time series of counts. The interrelation between the series is induced by the correlated innovations, while the non stationarity is captured through a common set of time-dependent covariates that influence the count responses. The regression and dependence effects are estimated using generalized quasi-likelihood (GQL) approach. Simulation experiments are performed to assess the performance of the estimation algorithms. The proposed BINAR(1) process is applied to analyze a real-life series of day and night accidents in Mauritius.

Mauritius stands as one of the few countries in the world to have controlled the current pandemic, the novel coronavirus 2019 (COVID-19) to a significant extent in a relatively short lapse of time. Owing to uncertainties and crisis amid the pandemic, as an emergency announcement, the World Health Organization (WHO) solicits the help of health authorities, especially, researchers to conduct in-depth research on the evolution and treatment of COVID-19. This paper proposes an integer-valued time series model to analyze the series of COVID-19 cases in Mauritius wherein the corresponding innovation term accommodates for covariate specification. In this set-up, sanitary curfew followed by sanitization and sensitization campaigns, time factor and safe shopping guidelines have been tested as the most significant variables, unlike climatic conditions. The over-dispersion estimates and the serial auto-correlation parameter are also statistically significant. This study also confirms the presence of some unobservable effects like the pathological genesis of the novel coronavirus and environmental factors which contribute to rapid propagation of the zoonotic virus in the community. Based on the proposed COM-Poisson mixture models, we could predict the number of COVID-19 cases in Mauritius. The forecasting results provide satisfactory mean squared errors. Such findings will subsequently encourage the policymakers to implement strict precautionary measures in terms of constant upgrading of the current health care and wellness system and re-enforcement of sanitary obligations.