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.     

Local likelihood estimation for nonstationary random fields

We develop a weighted local likelihood estimate for the parameters that govern the local spatial dependency of a locally stationary random field. The advantage of this local likelihood estimate is that it smoothly downweights the influence of faraway observations, works for irregular sampling locations, and when designed appropriately, can trade bias and variance for reducing estimation error. This paper starts with an exposition of our technique on the problem of estimating an unknown positive function when multiplied by a stationary random field. This example gives concrete evidence of the benefits of our local likelihood as compared to unweighted local likelihoods. We then discuss the difficult problem of estimating a bandwidth parameter that controls the amount of influence from distant observations. Finally we present a simulation experiment for estimating the local smoothness of a local Matérn random field when observing the field at random sampling locations in [0,1]². The local Matérn is a fully nonstationary random field, has a closed form covariance, can attain any degree of differentiability or Hölder smoothness and behaves locally like a stationary Matérn. We include an appendix that proves the positive definiteness of this covariance function.     

Hazard rate estimation on random fields

Consider observations (representing lifelengths) taken on a random field indexed by lattice points. Our purpose is to estimate the hazard rate r(x), which is the rate of failure at time x for the survivors up to time x. We estimate r(x) by the nonparametric estimator constructed in terms of a kernel-type estimator for f(x) and the natural estimator for . Under some general mixing assumptions, the limiting distribution of the estimator at multiple points is shown to be multivariate normal. The result is useful in establishing confidence bands for r(x) with x in an interval.     

Some inequalities for strong mixing random variables with applications to density estimation

In this paper, we establish an inequality of the characteristic functions for strongly mixing random vectors, by which, an upper bound is provided for the supremum of the absolute value of the difference of two multivariate probability density functions based on strongly mixing random vectors. As its application, we consider the consistency and asymptotic normality of a kernel estimate of a density function under strong mixing. Our results generalize some known results in the literature.     

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.     

A central limit theorem for stationary random fields

This paper establishes a central limit theorem and an invariance principle for a wide class of stationary random fields under natural and easily verifiable conditions. More precisely, we deal with random fields of the form _Xk=g(_εk−s,s∈Z^d)$X_k=g({\epsilon }_{k-s},s\in {\mathbb{Z}}^d)$, k∈Z^d$k\in {\mathbb{Z}}^d$, where _(εi)i∈Zd${\left({\epsilon }_i\right)}_{i\in {\mathbb{Z}}^d}$ are iid random variables and g$g$ is a measurable function. Such kind of spatial processes provides a general framework for stationary ergodic random fields. Under a short-range dependence condition, we show that the central limit theorem holds without any assumption on the underlying domain on which the process is observed. A limit theorem for the sample auto-covariance function is also established.     

On the mode-change problem for random measures

The classical change-point problem in modern terms, i.e., the mode-change problem, is stated for sufficiently general set-indexed random processes, namely for random measures. A method is shown for solving this problem both in the general form and for the intensity of compound Poisson random measures. The results obtained are novel for the change-point problem, too.     

Drift parameter estimation in fractional diffusions driven by perturbed random walks

We estimate the drift parameter in a simple linear model driven by fractional Brownian motion. We propose maximum likelihood estimators (MLE) for the drift parameter construct by using a random walk approximation of the fractional Brownian motion.     

Properties of nonparametric estimators of autocovariance for stationary random fields

Summary We introduce nonparametric estimators of the autocovariance of a stationary random field. One of our estimators has the property that it is itself an autocovatiance. This feature enables the estimator to be used as the basis of simulation studies such as those which are necessary when constructing bootstrap confidence intervals for unknown parameters. Unlike estimators proposed recently by other authors, our own do not require assumptions such as isotropy or monotonicity. Indeed, like nonparametric function estimators considered more widely in the context of curve estimation, our approach demands only smoothness and tail conditions on the underlying curve or surface (here, the autocovariance), and moment and mixing conditions on the random field. We show that by imposing the condition that the estimator be a covariance function we actually reduce the numerical value of integrated squared error.     

Nonparametric estimation for stochastic differential equations with random effects

We consider N$N$ independent stochastic processes (_Xj(t),t∈[0,T])$(X_j\left(t\right),t\in [0,T])$, j=1,…,N$j=1,\dots ,N$, defined by a one-dimensional stochastic differential equation with coefficients depending on a random variable _?j${?}_j$ and study the nonparametric estimation of the density of the random effect _?j${?}_j$ in two kinds of mixed models. A multiplicative random effect and an additive random effect are successively considered. In each case, we build kernel and deconvolution estimators and study their L²$L^2$-risk. Asymptotic properties are evaluated as N$N$ tends to infinity for fixed T$T$ or for T=T(N)$T=T\left(N\right)$ tending to infinity with N$N$. For T(N)=N²$T\left(N\right)=N^2$, adaptive estimators are built. Estimators are implemented on simulated data for several examples.     

A normal inverse Gaussian model for a risky asset with dependence

We present a new construction of the normal inverse Gaussian (NIG) fractal activity time model for a risky asset. The construction uses superpositions of diffusion processes and allows for specified exact NIG marginal distributions of the returns and flexible and tractable dependence structure including short or long range dependence. In the case of finite superposition, the fractal activity time is asymptotically self-similar, which is a desired feature seen in practice. The support for the distributional and dependence features of the risky asset model is provided by the data of currency exchange rates.     

Martingale transforms with non-atomic limits and stochastic approximation

Summary A sufficient condition for the limit of a martingale transform to posses a continuous distribution is given. The result is used to show that for a stochastic approximation procedure, if the adjustment rate is too small then it would not converge to the target value a.s. Furthermore, if the adjustment rate is taken to be 1/n as usual but the derivative of the regression function at the target value is 0, then the convergence rate is shown to be logn instead of , the rate obtained when the derivative is non-zero.Research supported by the National Science Council, R.O.C.     

Local Whittle estimator for anisotropic random fields

A local Whittle estimator is developed to simultaneously estimate the long memory parameters for stationary anisotropic scalar random fields. It is shown that these estimators are consistent and asymptotically normal, under some weak technical conditions. A brief simulation study illustrates a practical application of the estimator.     

Consistent estimation of the intensity function of a cyclic Poisson process

We construct and investigate a consistent kernel-type nonparametric estimator of the intensity function of a cyclic Poisson process when the period is unknown. We do not assume any particular parametric form for the intensity function, nor do we even assume that it is continuous. Moreover, we consider the situation when only a single realization of the Poisson process is available, and only in a bounded window. We prove, in particular, that the proposed estimator is consistent when the size of the window indefinitely expands. We also obtain complete convergence of the estimator.     

Density estimation for compound Poisson processes from discrete data

In this article we investigate the nonparametric estimation of the jump density of a compound Poisson process from the discrete observation of one trajectory over [0,T]$[0,T]$. We consider the case where the sampling rate Δ=_ΔT→0$\Delta ={\Delta }_T\to 0$ as T→∞$T\to \infty$. We propose an adaptive wavelet threshold density estimator and study its performance for _Lp$L_p$ losses, p≥1$p\ge 1$, over Besov spaces. The main novelty is that we achieve minimax rates of convergence for sampling rates _ΔT${\Delta }_T$ that vanish slowly. The estimation procedure is based on the explicit inversion of the operator giving the law of the increments as a nonlinear transformation of the jump density.     

Wiener processes with random effects for degradation data

This article studies the maximum likelihood inference on a class of Wiener processes with random effects for degradation data. Degradation data are special case of functional data with monotone trend. The setting for degradation data is one on which n independent subjects, each with a Wiener process with random drift and diffusion parameters, are observed at possible different times. Unit-to-unit variability is incorporated into the model by these random effects. EM algorithm is used to obtain the maximum likelihood estimators of the unknown parameters. Asymptotic properties such as consistency and convergence rate are established. Bootstrap method is used for assessing the uncertainties of the estimators. Simulations are used to validate the method. The model is fitted to bridge beam data and corresponding goodness-of-fit tests are carried out. Failure time distributions in terms of degradation level passages are calculated and illustrated.     

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.     

A modified functional delta method and its application to the estimation of risk functionals

The classical functional delta method (FDM) provides a convenient tool for deriving the asymptotic distribution of statistical functionals from the weak convergence of the respective empirical processes. However, for many interesting functionals depending on the tails of the underlying distribution this FDM cannot be applied since the method typically relies on Hadamard differentiability w.r.t. the uniform sup-norm. In this article, we present a version of the FDM which is suitable also for nonuniform sup-norms, with the outcome that the range of application of the FDM enlarges essentially. On one hand, our FDM, which we shall call the modified FDM, works for functionals that are “differentiable” in a weaker sense than Hadamard differentiability. On the other hand, it requires weak convergence of the empirical process w.r.t. a nonuniform sup-norm. The latter is not problematic since there exist strong respective results on weighted empirical processes obtained by Shorack and Wellner (1986) [25], Shao and Yu (1996) [23], Wu (2008) [32], and others. We illustrate the gain of the modified FDM by deriving the asymptotic distribution of plug-in estimates of popular risk measures that cannot be treated with the classical FDM.     

The point process approach for fractionally differentiated random walks under heavy traffic

We prove some heavy-traffic limit theorems for some nonstationary linear processes which encompass the fractionally differentiated random walk as well as some FARIMA processes, when the innovations are in the domain of attraction of a non-Gaussian stable distribution. The results are based on an extension of the point process methodology to linear processes with nonsummable coefficients and make use of a new maximal type inequality.     

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.