Asymptotic distributions for the intervals estimators of the extremal index and the cluster-size probabilities

Under appropriate long range dependence conditions, the point process of exceedances of a stationary sequence weakly converges to a homogeneous compound Poisson point process. This limiting point process can be characterized by the extremal index and the cluster-size probabilities. In this paper we address the problem of estimating these quantities and we consider the intervals estimators introduced in Ferro and Segers [2003. Inference for clusters of extreme values. J. Roy. Statist. Soc. Ser. B 545–556] and in Ferro [2004. Statistical methods for clusters of extreme values. Ph.D. Thesis, Lancaster University]. We establish asymptotic weak convergence to Gaussian random variables and we give their asymptotic variance.     

Asymptotic Location of Largest Values of a Stationary Random Field

The study of the relationship between extreme values of dependent random fields and their locations has important practical applications, for instance, when dealing with censored data.

In this article we study the asymptotic behavior of the joint locations of the largest order statistics generated by a stationary random field with extremal index as well as the joint limiting distribution of the location of a high level exceedance nearest of the origin and the location of the maximum.     

Regression Modelling of Disease Risk in Relation to Point Sources

We describe a class of models for the investigation of possible raised risk of disease around putative sources of environmental pollution. An adaptation of the point process method suggested by Diggle and Rowlingson is presented, allowing the use of routinely available aggregated data and incorporating the more general distance–risk model suggested by Elliott and co-workers. An application to data on cancers of the stomach around municipal solid waste incinerators is presented.     

Mixed-spectra analysis for stationary random fields

We consider a, discrete time, weakly stationary bidimensional process, for which the spectral measure is the sum of an absolutely continuous measure, a discrete measure of finite order and a finite number of absolutely continuous measures on several lines. In this paper we are interested in estimating the spectral density of the absolutely continuous measure and of the density on the lines. For this aim, by using the double kernel method, we construct consistent estimators of these densities and we study their asymptotic behaviors in term of the mean squared error with rate.     

Tests for Unit-Root versus Threshold Specification With an Application to the Purchasing Power Parity Relationship

We consider modeling the real exchange rate by a stationary three-regime self-exciting threshold autoregressive (SETAR) model with possibly a unit root in the middle regime. This representation is consistent with purchasing power parity in the presence of trading costs. Our main contribution is to provide statistical tools for testing unit root versus a SETAR. First, we show that a SETAR with a unit root in the middle regime is stationary and mixing under reasonable assumptions. Second, we derive analytically the asymptotic distribution of our unit-root test under the null. Using monthly real exchange rate data, our test rejects the null of unit-root against a threshold process for five European series.     

Re-colouring the Intensity-Based Bootstrap for Point Processes

Abstract

A method for obtaining bootstrapping replicates for one-dimensional point processes is presented. The method involves estimating the conditional intensity of the process and computing residuals. The residuals are bootstrapped using a block bootstrap and used, together with the conditional intensity, to define the bootstrap realizations. The method is applied to the estimation of the cross-intensity function for data arising from a reaction time experiment.     

A random cluster algorithm for some continuous Markov random fields

Based on a random cluster representation, the Swendsen–Wang algorithm for the Ising and Potts distributions is extended to a class of continuous Markov random fields. The algorithm can be described briefly as follows. A given configuration is decomposed into clusters. Probabilities for flipping the values of the random variables in each cluster are calculated. According to these probabilities, values of all the random variables in each cluster will be either updated or kept unchanged and this is done independently across the clusters. A new configuration is then obtained. We will show through a simulation study that, like the Swendsen–Wang algorithm in the case of Ising and Potts distributions, the cluster algorithm here also outperforms the Gibbs sampler in beating the critical slowing down for some strongly correlated Markov random fields.     

Extremal analysis of processes sampled at different frequencies

The observed extremes of a discrete time process depend on the process itself and the sampling frequency. We develop theoretical results which show how to account for the effect of sampling frequency on extreme values, thus enabling us to analyse systematically extremal data from series with different sampling rates. We present statistical methodology based on these results which we illustrate though simulations and by applications to sea-waves and rainfall data.     

Line transects of two-dimensional random fields: Estimation and design

In line-transect sampling data are collected continuously along the path travelled. These data are then used to estimate some characteristic of an entire region. We model the characteristic of the region as a linear functional of Z(), a realization of a random field, and our data as the values of Z(x) for x along the path. We study how to estimate the characteristic of interest and how to choose the path to accurately estimate it.     

Analyzing point processes subjected to random deletions

Certain nonstationary point process data are viewed as having arisen through time dependent random deletions of a stationary point process. Initially the probability Of deletion is assumed known and estimates of the rate and autointensity function of the inherent stationary process are constructed. Then an estimate of the deletion probability function is developed for the case of the function depending on a finite dimensional parameter. An estimate is provided for the variance of the autointensity estimate.     

Point and block prediction in log-Gaussian random fields: The non-constant mean case

This work considers the problems of point and block prediction in log-Gaussian random fields for the case when the mean of the log-process is not constant and depends linearly on unknown parameters. First, we propose a new point predictor that is optimal within a certain family of predictors, which extend a result in De Oliveira [2006. On optimal point and block prediction in log-Gaussian random fields. Scand. J. Statist. 33, 523–540.] that holds in the case when the mean of the log-process is constant. Second, we show that the results in De Oliveira [2006. On optimal point and block prediction in log-Gaussian random fields. Scand. J. Statist. 33, 523–540.] regarding optimal block prediction cannot be extended to the case when the mean of the log-process is not constant. Specifically, we show that the two families of block predictors considered by De Oliveira lack an optimal predictor. Finally, we numerically compare the predictive efficiency of the proposed point and block predictors.     

On the Multivariate Upcrossings Index

The multivariate extremal index function is a measure of the clustering among the extreme values of a multivariate stationary sequence. In this article, we introduce a measure of the degree of clustering of upcrossings in a multivariate stationary sequence, called multivariate upcrossings index, which is a multivariate generalization of the concept of upcrossings index. We derive the main properties of this function, namely the relations with the multivariate extremal index and the clustering of upcrossings.

Imposing general local and asymptotic dependence restrictions on the sequence or on its marginals we compute the multivariate upcrossings index from the marginal upcrossings indices and from the joint distribution of a finite number of variables. A couple of illustrative examples are exploited.     

Smoothed Extreme Value Estimators of Non-Uniform Point Processes Boundaries with Application to Star-Shaped Supports Estimation

We address the problem of estimating the edge of a bounded set in ? ^d given a random set of points drawn from the interior. Our method is based on a transformation of estimators dedicated to uniform point processes and obtained by smoothing some of its bias corrected extreme points. An application to the estimation of star-shaped supports is presented.     

Asymptotic results for random coefficient bifurcating autoregressive processes

The purpose of this paper is to study the asymptotic behaviour of the weighted least-squares estimators of the unknown parameters of random coefficient bifurcating autoregressive processes. Under suitable assumptions on the immigration and the inheritance, we establish the almost sure convergence of our estimators, as well as a quadratic strong law and central limit theorems. Our study mostly relies on limit theorems for vector-valued martingales.     

Estimating capability index cpk for processes with asymmetric tolerances

Pearn and Chen (1996) considered the process capability index C_pk, and investigated the statistical properties of its natural estimator under various process conditions. Their investigation, however, was restricted to processes with symmetric tolerances. Recently, Pearn and Chen (1998) considered a generalization of C_pk, referred to as C^? _pk, to cover processes with asymmetric tolerances. They investigated the statistical properties of the natural estimator of C^? _pk, and obtained the exact formulae for the expected value and variance. In this paper, we consider a new estimator of C^? _pk, assuming the knowledge on P(LI > T) = p is available, where 0 > p > 1, which can be obtained from historical information of a stable process. We obtain the exact distribution of the new estimator assuming the process characteristic follows the normal distribution. We show that the new estimator is consistent, asymptotically unbiased, which converges to a mixture of two normal distributions. We also show that by adding suitable correction factors to the new estimator, we may obtain the UMVUE and the MLE of the generalization C^? _pk.     

CARMA(p,q) generalized random processes

There is now a vast literature on the theory and applications of generalized random processes, pioneered by Itô (1953), Gel’fand (1955) and Yaglom (1957). In this note we make use of the theory of generalized random processes as defined in the book of Gel’fand and Vilenkin (1964) to extend the definition of continuous-time ARMA(p,q  ) processes to allow q≥p$q\ge p$, in which case the processes do not exist in the classical sense. The resulting CARMA generalized random processes provide a framework within which it is possible to study derivatives of CARMA processes of arbitrarily high order.     

Modeling material stress using integrated Gaussian Markov random fields

The equations of a physical constitutive model for material stress within tantalum grains were solved numerically using a tetrahedrally meshed volume. The resulting output included a scalar vonMises stress for each of the more than 94,000 tetrahedra within the finite element discretization. In this paper, we define an intricate statistical model for the spatial field of vonMises stress which uses the given grain geometry in a fundamental way. Our model relates the three-dimensional field to integrals of latent stochastic processes defined on the vertices of the one- and two-dimensional grain boundaries. An intuitive neighborhood structure of the said boundary nodes suggested the use of a latent Gaussian Markov random field (GMRF). However, despite the potential for computational gains afforded by GMRFs, the integral nature of our model and the sheer number of data points pose substantial challenges for a full Bayesian analysis. To overcome these problems and encourage efficient exploration of the posterior distribution, a number of techniques are now combined: parallel computing, sparse matrix methods, and a modification of a block update strategy within the sampling routine. In addition, we use an auxiliary variables approach to accommodate the presence of outliers in the data.     

Bayes estimation of the trendparameter in random fields

Based on continuous observations of a random field and on the first two prior moments of the trendparameter is derived. essentially by use of the solution of a WIENER-HOFF-type integral equation.     

Statistical inference for some Volterra population processes in a random environment

We study the problem of maximum-likelihood estimation in some random-environment population models, defined through nonlinear ltǒ stochastic differential equations. It is shown that a criticality parameter can be estimated consistently. The asymptotic behavior of the estimators is analyzed, and a goodness-of-fit test is proposed.