A spatiotemporal auto-regressive moving average model for solar radiation

Summary.  To investigate the variability in energy output from a network of photovoltaic cells, solar radiation was recorded at 10 sites every 10 min in the Pentland Hills to the south of Edinburgh. We identify spatiotemporal auto-regressive moving average models as the most appropriate to address this problem. Although previously considered computationally prohibitive to work with, we show that by approximating using toroidal space and fitting by matching auto-correlations, calculations can be substantially reduced. We find that a first-order spatiotemporal auto-regressive (STAR(1)) process with a first-order neighbourhood structure and a Matern noise process provide an adequate fit to the data, and we demonstrate its use in simulating realizations of energy output.     

A Note on the Innovation Distribution of a Gamma Distributed Autoregressive Process

A representation of the innovation random variable for a gamma distributed first-order autoregressive process was found by Lawrance (1982) in the form of a compound Poisson distribution, connected with a shot-noise process. In this note we simplify the representation of Lawrance by providing a direct representation in terms of density functions.     

The effect of autocorrelation on the retrospective X-chart

Quality control chart interpretation is usually based on the assumption that successive observations are independent over time. In this article we show the effect of autocorrelation on the retrospective Shewhart chart for individuals, often referred to as the X-chart, with the control limits based on moving ranges. It is shown that the presence of positive first lag autocorrelation results in an increased number of false alarms from the control chart. Negative first lag autocorrelation can result in unnecessarily wide control limits such that significant shifts in the process mean may go undetected. We use first-order autoregressive and first-order moving average models in our simulation of small samples of autocorrelated data.     

Information loss due to incomplete data from a spatial gaussian one-parameter first-order conditional process

Formulae are given for the Fisher information loss on parameters for the mean and the variance when some values of a Gaussian process are not observed. The special case of a one-parameter first-order conditional process on a rectangular lattice is considered in detail, and formulae are compared with numerical results.     

ARL-unbiased control charts with alarm and warning lines for monitoring Weibull percentiles using the first-order statistic

This article discusses methodology for constructing control charts to monitor the percentiles of a Weibull process with known shape parameter. Periodic samples are censored at the smallest observed value. Charts with alarm and warning limits are studied, and these limits are derived using theoretical results based on the first-order statistic. The performance of the proposed charts is evaluated and compared using average run lengths. A numerical application concerning life tests of an electronic product is presented to illustrate the methods.     

A Note on Whittle's Likelihood

The approximate likelihood function introduced by Whittle has been used to estimate the spectral density and certain parameters of a variety of time series models. In this note we attempt to empirically quantify the loss of efficiency of Whittle's method in nonstandard settings. A recently developed representation of some first-order non-Gaussian stationary autoregressive process allows a direct comparison of the true likelihood function with that of Whittle. The conclusion is that Whittle's likelihood can produce unreliable estimates in the non-Gaussian case, even for moderate sample sizes. Moreover, for small samples, and if the autocorrelation of the process is high, Whittle's approximation is not efficient even in the Gaussian case. While these facts are known to some extent, the present study sheds more light on the degree of efficiency loss incurred by using Whittle's likelihood, in both Gaussian and non-Gaussian cases.     

Inference on C pk for autocorrelated data in the presence of random measurement errors

The present paper examines the properties of the C _pk estimator when observations are autocorrelated and affected by measurement errors. The underlying reason for this choice of subject matter is that in industrial applications, process data are often autocorrelated, especially when sampling frequency is not particularly low, and even with the most advanced measuring instruments, gauge imprecision needs to be taken into consideration. In the case of a first-order stationary autoregressive process, we compare the statistical properties of the estimator in the error case with those of the estimator in the error-free case. Results indicate that the presence of gauge measurement errors leads the estimator to behave differently depending on the entity of error variability.     

Control charts based on dependent count data with deflation or inflation of zeros

The process of serially dependent counts with deflation or inflation of zeros is commonly observed in many applications. This paper investigates the monitoring of such a process, the first-order zero-modified geometric integer-valued autoregressive process (ZMGINAR(1)). In particular, two control charts, the upper-sided and lower-sided CUSUM charts, are developed to detect the shifts in the mean process of the ZMGINAR(1). Both the average run length performance and the standard deviation of the run length performance of these two charts are investigated by using Markov chain approaches. Also, an extensive simulation is conducted to assess the effectiveness or performance of the charts, and the presented methods are applied to two sets of real data arising from a study on the drug use.     

On first-order integer-valued autoregressive process with Katz family innovations

This paper considers the first-order integer-valued autoregressive (INAR) process with Katz family innovations. This family of INAR processes includes a broad class of INAR(1) processes with Poisson, negative binomial, and binomial innovations, respectively, featuring equi-, over-, and under-dispersion. Its probabilistic properties such as ergodicity and stationarity are investigated and the formula of the marginal mean and variance is provided. Further, a statistical process control procedure based on the cumulative sum control chart is considered to monitor autocorrelated count processes. A simulation and real data analysis are conducted for illustration.     

Multivariate Quality Control Chart for Autocorrelated Processes

Traditional multivariate statistical process control (SPC) techniques are based on the assumption that the successive observation vectors are independent. In recent years, due to automation of measurement and data collection systems, a process can be sampled at higher rates, which ultimately leads to autocorrelation. Consequently, when the autocorrelation is present in the data, it can have a serious impact on the performance of classical control charts. This paper considers the problem of monitoring the mean vector of a process in which observations can be modelled as a first-order vector autoregressive VAR (1) process. We propose a control chart called Z-chart which is based on the single step finite intersection test (Timm, 1996). An important feature of the proposed method is that it not only detects an out of control status but also helps in identifying variable(s) responsible for the out of control situation. The proposed method is illustrated with the help of suitable illustrations.     

A Nonparametric Test for the Parallelism of Two First-Order Autoregressive Processes

A nonparametric testing procedure for the parallelism of two first-order autoregressive processes is presented. This paper discuss the Mann–Whitney statistic, its natural competitor two-sample t -test, and the bootstrap method. It studies the asymptotic efficacies of the studentized Mann–Whitney statistic and the t -test statistic with their relative efficiency. Simulation results for comparing the powers of these test statistics are also presented.     

Generalized Poisson integer-valued autoregressive processes with structural changes

In this paper, we introduce a new first-order generalized Poisson integer-valued autoregressive process, for modeling integer-valued time series exhibiting a piecewise structure and overdispersion. Basic probabilistic and statistical properties of this model are discussed. Conditional least squares and conditional maximum likelihood estimators are derived. The asymptotic properties of the estimators are established. Moreover, two special cases of the process are discussed. Finally, some numerical results of the estimates and a real data example are presented.     

A new geometric first-order integer-valued autoregressive (NGINAR(1)) process

A new stationary first-order integer-valued autoregressive process with geometric marginal distributions is introduced based on negative binomial thinning. Some properties of the process are established. Estimators of the parameters of the process are obtained using the methods of conditional least squares, Yule–Walker and maximum likelihood. Also, the asymptotic properties of the estimators are derived involving their distributions. Some numerical results of the estimators are presented with a discussion to the obtained results. Real data are used and a possible application is discussed.     

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.     

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.     

Lindley first-order autoregressive model with applications

ABSTRACT

A new stationary first-order autoregressive process with Lindley marginal distribution, denoted as LAR(1) is introduced. We derive the probability function for the innovation process. We consider many properties of this process, involving spectral density, some multi-step ahead conditional measures, run probabilities, stationary solution, uniqueness and ergodicity. We estimate the unknown parameters of the process using three methods of estimation and investigate properties of the estimators with some numerical results to illustrate them. Some applications of the process are discussed to two real data sets and it is shown that the LAR(1) model fits better than other known non Gaussian AR(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.     

Robust first-order rotatable lifetime improvement experimental designs

Experimental designs are widely used in predicting the optimal operating conditions of the process parameters in lifetime improvement experiments. The most commonly observed lifetime distributions are log-normal, exponential, gamma and Weibull. In the present article, invariant robust first-order rotatable designs are derived for autocorrelated lifetime responses having log-normal, exponential, gamma and Weibull distributions. In the process, robust first-order D-optimal and rotatable conditions have been derived under these situations. For these lifetime distributions with correlated errors, it is shown that robust first-order D-optimal designs are always robust rotatable but the converse is not true. Moreover, it is observed that robust first-order D-optimal and rotatable designs depend on the respective error variance–covariance structure but are independent from these considered lifetime response distributions.     

Reliable error estimation for Sobol’ indices

In the field of sensitivity analysis, Sobol’ indices are sensitivity measures widely used to assess the importance of inputs of a model to its output. The estimation of these indices is often performed through Monte Carlo or quasi-Monte Carlo methods. A notable method is the replication procedure that estimates first-order indices at a reduced cost in terms of number of model evaluations. An inherent practical problem of this estimation is how to quantify the number of model evaluations needed to ensure that estimates satisfy a desired error tolerance. This article addresses this challenge by proposing a reliable error bound for first-order and total effect Sobol’ indices. Starting from the integral formula of the indices, the error bound is defined in terms of the discrete Walsh coefficients of the different integrands. We propose a sequential estimation procedure of Sobol’ indices using the error bound as a stopping criterion. The sequential procedure combines Sobol’ sequences with either Saltelli’s strategy to estimate both first-order and total effect indices, or the replication procedure to estimate only first-order indices.     

Inference for Random Coefficient INAR(1) Process Based on Frequency Domain Analysis

In this article, we present the explicit expressions for the higher-order moments and cumulants of the first-order random coefficient integer-valued autoregressive (RCINAR(1)) process. The spectral and bispectral density functions are also obtained, which can characterize the RCINAR(1) process in the frequency domain. We use a frequency domain approach which is named Whittle criterion to estimate the parameters of the process. We propose a test statistic which is based on the frequency domain approach for the hypothesis test, H₀: α = 0?H₁: 0 < α < 1, where α is the mean of the random coefficient in the process. The asymptotic distribution of the test statistic is obtained. We compare the proposed test statistic with other statistics that can test serial dependence in time series of count via a typically numerical simulation, which indicates that our proposed test statistic has a good power.