Regular variation of order 1 nonlinear AR-ARCH models

We prove both geometric ergodicity and regular variation of the stationary distribution for a class of nonlinear stochastic recursions that includes nonlinear AR-ARCH models of order 1. The Lyapounov exponent for the model, the index of regular variation and the spectral measure for the regular variation all are characterized by a simple two-state Markov chain.     

A new class of interacting Markov chain Monte Carlo methods

We present a new class of interacting Markov chain Monte Carlo methods to approximate numerically discrete-time nonlinear measure-valued equations. These stochastic processes belong to the class of self-interacting Markov chains with respect to their occupation measures. We provide several convergence results for these new methods including exponential estimates and a uniform convergence theorem with respect to the time parameter, yielding what seems to be the first results of this kind for this type of self-interacting models. We illustrate these models in the context of Feynman–Kac distribution semigroups arising in physics, biology and in statistics.     

Estimating the inter-arrival time density of Markov renewal processes under structural assumptions on the transition distribution

We consider a stationary Markov renewal process whose inter-arrival time density depends multiplicatively on the distance between the past and present state of the embedded chain. This is appropriate when the jump size is governed by influences that accumulate over time. Then we can construct an estimator for the inter-arrival time density that has the parametric rate of convergence. The estimator is a local von Mises statistic. The result carries over to the corresponding semi-Markov process.     

Representations of the optimal filter in the context of nonlinear filtering of random fields with fractional noise

The problem of nonlinear filtering of multiparameter random fields, observed in the presence of a long-range dependent spatial noise, is considered. When the observation noise is modelled by a persistent fractional Wiener sheet, several pathwise representations of the optimal filter are derived. The representations involve series of multiple stochastic integrals of different types and are particularly important since the evolution equations, satisfied by the best mean-square estimate of the signal random field, have a complicated analytical structure and fail to be proper (measure-valued) stochastic partial differential equations. Several of the above optimal filter representations involve a new family of strong martingale transforms associated to the multiparameter fractional Brownian sheet; the latter martingale family is of independent interest in fractional stochastic calculus of multiparameter random fields.     

The rank of the covariance matrix of an evanescent field

Evanescent random fields arise as a component of the 2D Wold decomposition of homogeneous random fields. Besides their theoretical importance, evanescent random fields have a number of practical applications, such as in modeling the observed signal in the space-time adaptive processing (STAP) of airborne radar data. In this paper we derive an expression for the rank of the low-rank covariance matrix of a finite dimension sample from an evanescent random field. It is shown that the rank of this covariance matrix is completely determined by the evanescent field spectral support parameters, alone. Thus, the problem of estimating the rank lends itself to a solution that avoids the need to estimate the rank from the sample covariance matrix. We show that this result can be immediately applied to considerably simplify the estimation of the rank of the interference covariance matrix in the STAP problem.     

Neighborhood radius estimation for variable-neighborhood random fields

We consider random fields defined by finite-region conditional probabilities depending on a neighborhood of the region which changes with the boundary conditions. To predict the symbols within any finite region, it is necessary to inspect a random number of neighborhood symbols which might change according to the value of them. In analogy with the one-dimensional setting we call these neighborhood symbols the context associated to the region at hand. This framework is a natural extension, to d-dimensional fields, of the notion of variable length Markov chains introduced by Rissanen [24] in his classical paper. We define an algorithm to estimate the radius of the smallest ball containing the context based on a realization of the field. We prove the consistency of this estimator. Our proofs are constructive and yield explicit upper bounds for the probability of wrong estimation of the radius of the context.     

Unbiased estimates for gradients of stochastic network performance measures

Three classes of stochastic networks and their performance measures are considered. These performance measures are defined as the expected value of some random variables and cannot normally be obtained analytically as functions of network parameters in a closed form. We give similar representations for the random variables to provide a useful way of analytical study of these functions and their gradients. The representations are used to obtain sufficient conditions for the gradient estimates to be unbiased. The conditions are rather general and usually met in simulation study of the stochastic networks. Applications of the results are discussed and some practical algorithms of calculating unbiased estimates of the gradients are also presented.     

A class of nonlinear evolution equations subjected to nonlocal initial conditions

We prove the existence of C⁰-solutions for a class of nonlinear evolution equations subjected to nonlocal initial conditions, of the form:

     

Fractional-order regularization and wavelet approximation to the inverse estimation problem for random fields

The least-squares linear inverse estimation problem for random fields is studied in a fractional generalized framework. First, the second-order regularity properties of the random fields involved in this problem are analysed in terms of the fractional Sobolev norms. Second, the incorporation of prior information in the form of a fractional stochastic model, with covariance operator bicontinuous with respect to a certain fractional Sobolev norm, leads to a regularization of this problem. Third, a multiresolution approximation to the class of linear inverse problems considered is obtained from a wavelet-based orthogonal expansion of the input and output random models. The least-squares linear estimate of the input random field is then computed using these orthogonal wavelet decompositions. The results are applied to solving two important cases of linear inverse problems defined in terms of fractional integral operators.     

On testing equality of pairwise rank correlations in a multivariate random vector

Spearman’s rank-correlation coefficient (also called Spearman’s rho) represents one of the best-known measures to quantify the degree of dependence between two random variables. As a copula-based dependence measure, it is invariant with respect to the distribution’s univariate marginal distribution functions. In this paper, we consider statistical tests for the hypothesis that all pairwise Spearman’s rank correlation coefficients in a multivariate random vector are equal. The tests are nonparametric and their asymptotic distributions are derived based on the asymptotic behavior of the empirical copula process. Only weak assumptions on the distribution function, such as continuity of the marginal distributions and continuous partial differentiability of the copula, are required for obtaining the results. A nonparametric bootstrap method is suggested for either estimating unknown parameters of the test statistics or for determining the associated critical values. We present a simulation study in order to investigate the power of the proposed tests. The results are compared to a classical parametric test for equal pairwise Pearson’s correlation coefficients in a multivariate random vector. The general setting also allows the derivation of a test for stochastic independence based on Spearman’s rho.     

Local asymptotic normality in a stationary model for spatial extremes

De Haan and Pereira (2006) [6] provided models for spatial extremes in the case of stationarity, which depend on just one parameter β>0 measuring tail dependence, and they proposed different estimators for this parameter. We supplement this framework by establishing local asymptotic normality (LAN) of a corresponding point process of exceedances above a high multivariate threshold. Standard arguments from LAN theory then provide the asymptotic minimum variance within the class of regular estimators of β. It turns out that the relative frequency of exceedances is a regular estimator sequence with asymptotic minimum variance, if the underlying observations follow a multivariate extreme value distribution or a multivariate generalized Pareto distribution.     

Weak convergence in the functional autoregressive model

The functional autoregressive model is a Markov model taylored for data of functional nature. It revealed fruitful when attempting to model samples of dependent random curves and has been widely studied along the past few years. This article aims at completing the theoretical study of the model by addressing the issue of weak convergence for estimates from the model. The main difficulties stem from an underlying inverse problem as well as from dependence between the data. Traditional facts about weak convergence in non-parametric models appear: the normalizing sequence is not an , a bias term appears. Several original features of the functional framework are pointed out.     

Discrete-time stochastic processes on Riesz spaces

It has been recognised that order is closely linked with probability theory, with lattice theoretic approaches being used to study Markov processes but, to our knowledge, the complete theory of (sub, super) martingales and their stopping times has not been formulated on Riesz spaces. We generalize the concepts of stochastic processes, (sub, super) martingales and stopping times to Riesz spaces. In this paper we consider discrete time processes with bounded stopping times.     

Multivariate conditional versions of Spearman's rho and related measures of tail dependence

A new family of conditional-dependence measures based on Spearman's rho is introduced. The corresponding multidimensional versions are established. Asymptotic distributional results are derived for related estimators which are based on the empirical copula. Particular emphasis is placed on a new type of multidimensional tail-dependence measure and its relationship to other measures of tail dependence is shown. Multivariate tail dependence describes the limiting amount of dependence in the vertices of the copula's domain.     

Estimation of two ordered bivariate mean residual life functions

Situations occur frequently in which the mean residual life (mrl) functions of two populations must be ordered. For example, if a mechanical device is improved, the mrl function for the improved device should not be less than that of the original device. Also, mrl functions for medical patients should often be ordered depending on the status of concomitant variables. This paper proposes nonparametric estimators of the bivariate mrl function under a mrl ordering. The estimators are shown to be asymptotically unbiased, strongly uniformly consistent and weakly convergent to a bivariate Gaussian process. The estimators are shown to be the projections, in a sense to be made precise, of the empirical mrl function onto an appropriate convex set of mrl functions. In the one-sample problem, the new estimators dominate the empirical mrl function in terms of risk with respect to a wide class of loss functions.     

Nonlinear least-squares estimation

The paper uses empirical process techniques to study the asymptotics of the least-squares estimator (LSE) for the fitting of a nonlinear regression function. By combining and extending ideas of Wu and Van de Geer, it establishes new consistency and central limit theorems that hold under only second moment assumptions on the errors. An application to a delicate example of Wu's illustrates the use of the new theorems, leading to a normal approximation to the LSE with unusual logarithmic rescalings.     

Sequential tracking of a hidden Markov chain using point process observations

We study finite horizon optimal switching problems for hidden Markov chain models with point process observations. The controller possesses a finite range of strategies and attempts to track the state of the unobserved state variable using Bayesian updates over the discrete observations. Such a model has applications in economic policy making, staffing under variable demand levels and generalized Poisson disorder problems. We show regularity of the value function and explicitly characterize an optimal strategy. We also provide an efficient numerical scheme and illustrate our results with several computational examples.     

An analytic approach to purely nonlocal Bellman equations arising in models of stochastic control

Given a bounded domain Ω⊂Rd and two integro-differential operators L¹, L² of the form we study the fully nonlinear Bellman equation

(0.1)