Four general multivariate stationary extremal Markovian processes

Two general multivariate stationary Markovian process with maximization structure (denoted by Max-AR(1) and MaxI-AR(1)) are developed respectively. Max-AR(1) is a subclass of MaxI-AR(1). The characterization of the Max-AR(1) and MaxI-AR(1) to be stationary is studied. Some properties of the two maximization processes are derived. Two more related general multivariate stochastic Markovian process with minification structure are analogously constructed (denoted by Min-AR(1) and MinI-AR(1)). Some well known maximization and minification processes are special cases of these four extermal Markovian processes. Two of them are simulated and some point estimations are provided as an illustration of the wide application of these four processes.     

Some Related Minima Stability and Minima Infinite Divisibility of the General Multivariate Pareto Distributions

This article studies the minima stable property of the general multivariate Pareto distributions MP^(k)(I), MP^(k)(II), MP^(k)(III), MP^(k)(IV) which can be applied to characterize the MP^(k) distribution via its weighted ordered coordinates minima and marginal distribution. Also, the multivariate semi-Pareto distribution (denoted by MSP) is discerned in the class of geometric minima infinite divisible and geometric minima stable distributions. If the exponent measure is satisfied by some functional equation, then the geometric minima stable property can be used to characterize the MSP distribution. Finally, the finite sample minima infinite divisible property of the MP^(k)(I), (II), and (IV) distributions is also discussed.     

Bridging the Gap: A Generalized Stochastic Process for Count Data

The Bernoulli and Poisson processes are two popular discrete count processes; however, both rely on strict assumptions. We instead propose a generalized homogenous count process (which we name the Conway–Maxwell–Poisson or COM-Poisson process) that not only includes the Bernoulli and Poisson processes as special cases, but also serves as a flexible mechanism to describe count processes that approximate data with over- or under-dispersion. We introduce the process and an associated generalized waiting time distribution with several real-data applications to illustrate its flexibility for a variety of data structures. We consider model estimation under different scenarios of data availability, and assess performance through simulated and real datasets. This new generalized process will enable analysts to better model count processes where data dispersion exists in a more accommodating and flexible manner.     

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.     

Non‐stationary Cross‐Covariance Models for Multivariate Processes on a Globe

Abstract. In geophysical and environmental problems, it is common to have multiple variables of interest measured at the same location and time. These multiple variables typically have dependence over space (and/or time). As a consequence, there is a growing interest in developing models for multivariate spatial processes, in particular, the cross‐covariance models. On the other hand, many data sets these days cover a large portion of the Earth such as satellite data, which require valid covariance models on a globe. We present a class of parametric covariance models for multivariate processes on a globe. The covariance models are flexible in capturing non‐stationarity in the data yet computationally feasible and require moderate numbers of parameters. We apply our covariance model to surface temperature and precipitation data from an NCAR climate model output. We compare our model to the multivariate version of the Matérn cross‐covariance function and models based on coregionalization and demonstrate the superior performance of our model in terms of AIC (and/or maximum loglikelihood values) and predictive skill. We also present some challenges in modelling the cross‐covariance structure of the temperature and precipitation data. Based on the fitted results using full data, we give the estimated cross‐correlation structure between the two variables.     

New control charts for simultaneous monitoring of the mean vector and covariance matrix of multivariate multiple linear profiles

Simultaneous monitoring of the mean vector and covariance matrix in multivariate processes allows practitioners to avoid the inflated false alarm rate that results from using two independent control charts. In this paper, we extend exponentially weighted moving average semicircle and generally weighted moving average semicircle control charts to monitor the mean vector and covariance matrix of multivariate multiple linear regression profiles in Phase II simultaneously. These new control charts are compared with the existing control charts in the literature in terms of the average run length criterion. Finally, a case is considered to show the application of the proposed charts.     

One-Class Classification-Based Control Charts for Monitoring Autocorrelated Multivariate Processes

In recent years, statistical process control (SPC) of multivariate and autocorrelated processes has received a great deal of attention. Modern manufacturing/service systems with more advanced technology and higher production rates can generate complex processes in which consecutive observations are dependent and each variable is correlated. These processes obviously violate the assumption of the independence of each observation that underlies traditional SPC and thus deteriorate the performance of its traditional tools. The popular way to address this issue is to monitor the residuals—the difference between the actual value and the fitted value—with the traditional SPC approach. However, this residuals-based approach requires two steps: (1) finding the residuals; and (2) monitoring the process. Also, an accurate prediction model is necessary to obtain the uncorrelated residuals. Furthermore, these residuals are not the original values of the observations and consequently may have lost some useful information about the targeted process. The main purpose of this article is to examine the feasibility of using one-class classification-based control charts to handle multivariate and autocorrelated processes. The article uses simulated data to present an analysis and comparison of one-class classification-based control charts and the traditional Hotelling's T ² chart.     

AN ADAPTIVE TRIMMED LIKELIHOOD ALGORITHM FOR IDENTIFICATION OF MULTIVARIATE OUTLIERS

This article describes an algorithm for the identification of outliers in multivariate data based on the asymptotic theory for location estimation as described typically for the trimmed likelihood estimator and in particular for the minimum covariance determinant estimator. The strategy is to choose a subset of the data which minimizes an appropriate measure of the asymptotic variance of the multivariate location estimator. Observations not belonging to this subset are considered potential outliers which should be trimmed. For α less than about 0.5, the correct trimming proportion is taken to be that α > 0 for which the minimum of any minima of this measure of the asymptotic variance occurs. If no minima occur for an α > 0 then the data set will be considered outlier free.     

On improvement of prediction in arma processes

Necessary and sufficient conditions are derived in the paper that enable to decide whether an additional multivariate process will improve the prediction in a given multivariate discrete stationary process. The both processes are assumed to form together a process ARMAm n Further it was investigated wnen one can asser t that the both processes are uncorrelated provided the additional process did not improve the prediction in the original process, Some hints for the actual construction of predictors in a multivariate ARMA. (m n) process can be found in the paper.     

Continuous Time Wishart Process for Stochastic Risk

Risks are usually represented and measured by volatility–covolatility matrices. Wishart processes are models for a dynamic analysis of multivariate risk and describe the evolution of stochastic volatility–covolatility matrices, constrained to be symmetric positive definite. The autoregressive Wishart process (WAR) is the multivariate extension of the Cox, Ingersoll, Ross (CIR) process introduced for scalar stochastic volatility. As a CIR process it allows for closed-form solutions for a number of financial problems, such as term structure of T-bonds and corporate bonds, derivative pricing in a multivariate stochastic volatility model, and the structural model for credit risk. Moreover, the Wishart dynamics are very flexible and are serious competitors for less structural multivariate ARCH models.     

Phase-II monitoring and diagnosing of multivariate categorical processes using generalized linear test-based control charts

In this paper, two control charts based on the generalized linear test (GLT) and contingency table are proposed for Phase-II monitoring of multivariate categorical processes. The performances of the proposed methods are compared with the exponentially weighted moving average-generalized likelihood ratio test (EWMA-GLRT) control chart proposed in the literature. The results show the better performance of the proposed control charts under moderate and large shifts. Moreover, a new scheme is proposed to identify the parameter responsible for an out-of-control signal. The performance of the proposed diagnosing procedure is evaluated through some simulation experiments.     

EVVMA and cusum control charts in the presence of correlation

The statistical properties of control charts are usually evaluated under the assumption that the observations from the process are independent. For many processes however, observations which are closely spaced in time will be correlated. This paper considers EWMA and CUSUM control charts for the process mean when the observations are from an AR(1) process with additional random error. This simple model may be a reasonable model for many processes encountered in practice. The ARL and steady state ARL of the EWMA and CUSUM charts are evaluated numerically using an integral equation approach and a Markov chain approach. The numerical results show that correlation can have a significant effect on the properties of these charts. Tables are given to aid in the design of these charts when the observations follow the assumed model.     

Dependence Estimation for High‐frequency Sampled Multivariate CARMA Models

The paper considers high‐frequency sampled multivariate continuous‐time autoregressive moving average (MCARMA) models and derives the asymptotic behaviour of the sample autocovariance function to a normal random matrix. Moreover, we obtain the asymptotic behaviour of the cross‐covariances between different components of the model. We will see that the limit distribution of the sample autocovariance function has a similar structure in the continuous‐time and in the discrete‐time model. As a special case, we consider a CARMA (one‐dimensional MCARMA) process. For a CARMA process, we prove Bartlett's formula for the sample autocorrelation function. Bartlett's formula has the same form in both models; only the sums in the discrete‐time model are exchanged by integrals in the continuous‐time model. Finally, we present limit results for multivariate MA processes as well, which are not known in this generality in the multivariate setting yet.     

SPC for short-run multivariate autocorrelated processes

This paper discusses the development of a multivariate control charting technique for short-run autocorrelated data manufacturing environment. The proposed approach is a combination of the multivariate residual charts for autocorrelated data and the multivariate transformation technique for i.i.d. process observations of short lengths. The proposed approach consists in fitting adequate multivariate time-series model of various process outputs and computes the residuals, transforming them into standard normal N(0, 1) data and then using standardized data as inputs to plot conventional univariate i.i.d. control charts. The objective for applying multivariate finite horizon techniques for autocorrelated processes is to allow continuous process monitoring, since all process outputs are controlled trough the use of a single control chart with constant control limits. Throughout simulated examples, it is shown that the proposed short-run process monitoring technique provides approximately similar shifts detection properties as VAR residual charts.     

Testing the Equality of Multivariate Distributions Using the Bootstrap and Integrated Empirical Processes

We introduce some projected integrated empirical processes for testing the equality of two multivariate distributions. The bootstrap is used for determining the approximate critical values. We show that the bootstrap test is consistent. A number-theoretic method is used for efficient computation of the bootstrap critical values. Some simulation results are also given.     

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.     

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.     

Some principles for surveillance adopted for multivariate processes with a common change point

The surveillance of multivariate processes has received growing attention during the last decade. Several generalizations of well-known methods such as Shewhart, CUSUM and EWMA charts have been proposed. Many of these multivariate procedures are based on a univariate summarized statistic of the multivariate observations, usually the likelihood ratio statistic. In this paper we consider the surveillance of multivariate observation processes for a shift between two fully specified alternatives. The effect of the dimension reduction using likelihood ratio statistics are discussed in the context of sufficiency properties. Also, an example of the loss of efficiency when not using the univariate sufficient statistic is given. Furthermore, a likelihood ratio method, the LR method, for constructing surveillance procedures is suggested for multivariate surveillance situations. It is shown to produce univariate surveillance procedures based on the sufficient likelihood ratios. As the LR procedure has several optimality properties in the univariate, it is also used here as a benchmark for comparisons between multivariate surveillance procedures     

Parameter estimation for multivariate diffusion processes with the time inhomogeneously positive semidefinite diffusion matrix

Statistical inference for the diffusion coefficients of multivariate diffusion processes has been well established in recent years; however, it is not the case for the drift coefficients. Furthermore, most existing estimation methods for the drift coefficients are proposed under the assumption that the diffusion matrix is positive definite and time homogeneous. In this article, we put forward two estimation approaches for estimating the drift coefficients of the multivariate diffusion models with the time inhomogeneously positive semidefinite diffusion matrix. They are maximum likelihood estimation methods based on both the martingale representation theorem and conditional characteristic functions and the generalized method of moments based on conditional characteristic functions, respectively. Consistency and asymptotic normality of the generalized method of moments estimation are also proved in this article. Simulation results demonstrate that these methods work well.     

Multivariate convex risk statistics with scenario analysis

Abstract

In this article, we introduce three new classes of multivariate risk statistics, which can be considered as data-based versions of multivariate risk measures. These new classes are multivariate convex risk statistics, multivariate comonotonic convex risk statistics and multivariate empirical-law-invariant convex risk statistics, respectively. Representation results are provided. The arguments of proofs are mainly developed by ourselves. It turns out that all the relevant existing results in the literature are special cases of those obtained in this article.