The Wiener disorder problem with finite horizon

The Wiener disorder problem seeks to determine a stopping time which is as close as possible to the (unknown) time of ‘disorder’ when the drift of an observed Wiener process changes from one value to another. In this paper we present a solution of the Wiener disorder problem when the horizon is finite. The method of proof is based on reducing the initial problem to a parabolic free-boundary problem where the continuation region is determined by a continuous curved boundary. By means of the change-of-variable formula containing the local time of a diffusion process on curves we show that the optimal boundary can be characterized as a unique solution of the nonlinear integral equation.     

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.     

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.     

Cramér–Karhunen–Loève representation and harmonic principal component analysis of functional time series

We develop a doubly spectral representation of a stationary functional time series, and study the properties of its empirical version. The representation decomposes the time series into an integral of uncorrelated frequency components (Cramér representation), each of which is in turn expanded in a Karhunen–Loève series. The construction is based on the spectral density operator, the functional analogue of the spectral density matrix, whose eigenvalues and eigenfunctions at different frequencies provide the building blocks of the representation. By truncating the representation at a finite level, we obtain a harmonic principal component analysis of the time series, an optimal finite dimensional reduction of the time series that captures both the temporal dynamics of the process, as well as the within-curve dynamics. Empirical versions of the decompositions are introduced, and a rigorous analysis of their large-sample behaviour is provided, that does not require any prior structural assumptions such as linearity or Gaussianity of the functional time series, but rather hinges on Brillinger-type mixing conditions involving cumulants.     

Symmetry analysis of a model of stochastic volatility with time-dependent parameters

We provide the solutions for the Heston model of stochastic volatility when the parameters of the model are constant and when they are functions of time. In the former case, the solution follows immediately from the determination of the Lie point symmetries of the governing 1+1 evolution partial differential equation. This is not the situation in the latter case, but we are able to infer the essential structure of the required nonlocal symmetry from that of the autonomous problem and hence can present the solution to the nonautonomous problem. As in the case of the standard Black-Scholes problem the presence of time-dependent parameters is not a hindrance to the demonstration of a solution.     

A stochastic partially reversible investment problem on a finite time-horizon: Free-boundary analysis

We study a continuous-time, finite horizon, stochastic partially reversible investment problem for a firm producing a single good in a market with frictions. The production capacity is modeled as a one-dimensional, time-homogeneous, linear diffusion controlled by a bounded variation process which represents the cumulative investment–disinvestment strategy. We associate to the investment–disinvestment problem a zero-sum optimal stopping game and characterize its value function through a free-boundary problem with two moving boundaries. These are continuous, bounded and monotone curves that solve a system of non-linear integral equations of Volterra type. The optimal investment–disinvestment strategy is then shown to be a diffusion reflected at the two boundaries.     

Weak Dirichlet processes with a stochastic control perspective

The motivation of this paper is to prove verification theorems for stochastic optimal control of finite dimensional diffusion processes without control in the diffusion term, in the case where the value function is assumed to be continuous in time and once differentiable in the space variable (C^0,1$C^{0,1}$) instead of once differentiable in time and twice in space (C^1,2$C^{1,2}$), like in the classical results. For this purpose, the replacement tool of the Itô formula will be the Fukushima–Dirichlet decomposition for weak Dirichlet processes. Given a fixed filtration, a weak Dirichlet process is the sum of a local martingale M$M$ plus an adapted process A$A$ which is orthogonal, in the sense of covariation, to any continuous local martingale. The decomposition mentioned states that a C^0,1$C^{0,1}$ function of a weak Dirichlet process with finite quadratic variation is again a weak Dirichlet process. That result is established in this paper and it is applied to the strong solution of a Cauchy problem with final condition.     

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.     

Locally stationary long memory estimation

There exists a wide literature on parametrically or semi-parametrically modelling strongly dependent time series using a long-memory parameter d, including more recent work on wavelet estimation. As a generalization of these latter approaches, in this work we allow the long-memory parameter d to be varying over time. We adopt a semi-parametric approach in order to avoid fitting a time-varying parametric model, such as tvARFIMA, to the observed data. We study the asymptotic behavior of a local log-regression wavelet estimator of the time-dependent d. Both simulations and a real data example complete our work on providing a fairly general approach.     

On Nummelin splitting for continuous time Harris recurrent Markov processes and application to kernel estimation for multi-dimensional diffusions

We introduce a sequence of stopping times that allow us to study an analogue of a life-cycle decomposition for a continuous time Markov process, which is an extension of the well-known splitting technique of Nummelin to the continuous time case. As a consequence, we are able to give deterministic equivalents of additive functionals of the process and to state a generalisation of Chen’s inequality. We apply our results to the problem of non-parametric kernel estimation of the drift of multi-dimensional recurrent, but not necessarily ergodic, diffusion processes.     

Concentration of measure and isoperimetric inequalities in product spaces

The concentration of measure phenomenon in product spaces roughly states that, if a set A in a product Ω^N of probability spaces has measure at least one half, “most” of the points of Ωn are “close” to A. We proceed to a systematic exploration of this phenomenon. The meaning of the word “most” is made rigorous by isoperimetrictype inequalities that bound the measure of the exceptional sets. The meaning of the work “close” is defined in three main ways, each of them giving rise to related, but different inequalities. The inequalities are all proved through a common scheme of proof. Remarkably, this simple approach not only yields qualitatively optimal results, but, in many cases, captures near optimal numerical constants. A large number of applications are given, in particular to Percolation, Geometric Probability, Probability in Banach Spaces, to demonstrate in concrete situations the extremely wide range of application of the abstract tools. Dedicated to Vitali Milman     

Power variation for Gaussian processes with stationary increments

We develop the asymptotic theory for the realised power variation of the processes X=?•G$X=?•G$, where G$G$ is a Gaussian process with stationary increments. More specifically, under some mild assumptions on the variance function of the increments of G$G$ and certain regularity conditions on the path of the process ?$?$ we prove the convergence in probability for the properly normalised realised power variation. Moreover, under a further assumption on the Hölder index of the path of ?$?$, we show an associated stable central limit theorem. The main tool is a general central limit theorem, due essentially to Hu and Nualart [Y. Hu, D. Nualart, Renormalized self-intersection local time for fractional Brownian motion, Ann. Probab. (33) (2005) 948–983], Nualart and Peccati [D. Nualart, G. Peccati, Central limit theorems for sequences of multiple stochastic integrals, Ann. Probab. (33) (2005) 177–193] and Peccati and Tudor [G. Peccati, C.A. Tudor, Gaussian limits for vector-valued multiple stochastic integrals, in: M. Emery, M. Ledoux, M. Yor (Eds.), Seminaire de Probabilites XXXVIII, in: Lecture Notes in Math, vol. 1857, Springer-Verlag, Berlin, 2005, pp. 247–262], for sequences of random variables which admit a chaos representation.     

Vectors of two-parameter Poisson-Dirichlet processes

The definition of vectors of dependent random probability measures is a topic of interest in applications to Bayesian statistics. They represent dependent nonparametric prior distributions that are useful for modelling observables for which specific covariate values are known. In this paper we propose a vector of two-parameter Poisson-Dirichlet processes. It is well-known that each component can be obtained by resorting to a change of measure of a σ-stable process. Thus dependence is achieved by applying a Lévy copula to the marginal intensities. In a two-sample problem, we determine the corresponding partition probability function which turns out to be partially exchangeable. Moreover, we evaluate predictive and posterior distributions.     

A simulation approach to optimal stopping under partial information

We study the numerical solution of nonlinear partially observed optimal stopping problems. The system state is taken to be a multi-dimensional diffusion and drives the drift of the observation process, which is another multi-dimensional diffusion with correlated noise. Such models where the controller is not fully aware of her environment are of interest in applied probability and financial mathematics. We propose a new approximate numerical algorithm based on the particle filtering and regression Monte Carlo methods. The algorithm maintains a continuous state space and yields an integrated approach to the filtering and control sub-problems. Our approach is entirely simulation-based and therefore allows for a robust implementation with respect to model specification. We carry out the error analysis of our scheme and illustrate with several computational examples. An extension to discretely observed stochastic volatility models is also considered.     

On Toeplitz type quadratic functionals of stationary Gaussian processes

Summary A central limit theorem for Toeplitz type quadratic functionals of a stationary Gaussian processX(t),t, is proved, generalizing the result of Avram [1] for discrete time processes. The result is applied to the problem of nonparametric estimation of linear functionals of an unknown spectral density function. We give some upper bounds for the minimax mean square risk of the nonparametric estimators, similar to those by Ibragimov and Has'minskii [12] for a probability density function.     

The hyper-Dirichlet process and its discrete approximations: The butterfly model

The aim of this paper is the study of some random probability distributions, called hyper-Dirichlet processes. In the simplest situation considered in the paper these distributions charge the product of three sample spaces, with the property that the first and the last component are independent conditional to the middle one. The law of the marginals on the first two and on the last two components are specified to be Dirichlet processes with the same marginal parameter measure on the common second component. The joint law is then obtained as the hyper-Markov combination, introduced in [A.P. Dawid, S.L. Lauritzen, Hyper-Markov laws in the statistical analysis of decomposable graphical models, Ann. Statist. 21 (3) (1993) 1272-1317], of these two Dirichlet processes. The processes constructed in this way in fact are in fact generalizations of the hyper-Dirichlet laws on contingency tables considered in the above paper. Our main result is the convergence to the hyper-Dirichlet process of the sequence of hyper-Dirichlet laws associated to finer and finer “discretizations” of the two parameter measures, which is proved by means of a suitable coupling construction.     

Sequential testing of simple hypotheses about compound Poisson processes

One of two simple hypotheses for the unknown arrival rate and jump distribution of a compound Poisson process is correct. We start observing the process, and the problem is to decide on the correct hypothesis as soon as possible and with the smallest probability of wrong decision. We find a Bayes-optimal sequential decision rule and describe completely how to calculate its parameters without any restrictions on the arrival rate and the jump distribution.     

Stochastic wave equation of pure jumps: Existence,uniqueness and invariant measures

In this paper, we investigate a class of nonlinear damped stochastic hyperbolic equations with jumps. The jump component considered here is described as a Poisson point process. This paper is divided into two parts. The first part deals with existence and uniqueness of global weak and strong solutions to this type of equations, based on the energy approach. The second part devotes to the existence and support of invariant measures corresponding to the weak solution semi-group, based on Markov property of the solution.