Semi-aligned rank tests

Summary The distribution-free test based on semi-aligned rankings for no treatment effects in a two-way layout, with unequal number of replications in each cell is considered. The asymptotic χ-square distribution of the test statistic under the null hypothesis is derived. The Pitman asymptotic relative efficiency of the test (i) based on semi-aligned rankings with respect to the test (ii) based on within-block rankings, is shown to be larger than one as the number of blocks tends to infinity. Also the asymptotic properties of linear rank statistics (i) and (ii) are investigated and the asymptotic relative efficiency of the test (i) with respect to the test (ii) is again shown to be larger than one.     

Exact asymptotic minimax constants for the estimation of analytical functions in L p

The L _p minimax risks (1≤p<∞) are studied for statistical estimation in the Gaussian white noise model. The asymptotic rate and constants are given, and the optimal estimator is proposed. This, together with the work of Golubev, Levit and Tsybakov (1996) establishes the classification of the L _p minimax constants on the classes of analytical functions. Received: 10 December 1996 / Revised version: 14 December 1997     

Nonparametric estimation of an extreme-value copula in arbitrary dimensions

Inference on an extreme-value copula usually proceeds via its Pickands dependence function, which is a convex function on the unit simplex satisfying certain inequality constraints. In the setting of an i.i.d. random sample from a multivariate distribution with known margins and an unknown extreme-value copula, an extension of the Capéraà-Fougères-Genest estimator was introduced by D. Zhang, M. T. Wells and L. Peng [Nonparametric estimation of the dependence function for a multivariate extreme-value distribution, Journal of Multivariate Analysis 99 (4) (2008) 577-588]. The joint asymptotic distribution of the estimator as a random function on the simplex was not provided. Moreover, implementation of the estimator requires the choice of a number of weight functions on the simplex, the issue of their optimal selection being left unresolved.A new, simplified representation of the CFG-estimator combined with standard empirical process theory provides the means to uncover its asymptotic distribution in the space of continuous, real-valued functions on the simplex. Moreover, the ordinary least-squares estimator of the intercept in a certain linear regression model provides an adaptive version of the CFG-estimator whose asymptotic behavior is the same as if the variance-minimizing weight functions were used. As illustrated in a simulation study, the gain in efficiency can be quite sizable.     

Some central limit theorems for ℓ∞-valued semimartingales and their applications

Summary. This paper is devoted to the generalization of central limit theorems for empirical processes to several types of ℓ^∞(Ψ)-valued continuous-time stochastic processes t⇝X _t ⁿ =(X _t ^{n ,ψ}|ψ∈Ψ), where Ψ is a non-empty set. We deal with three kinds of situations as follows. Each coordinate process t⇝X _t ^{n ,ψ} is: (i) a general semimartingale; (ii) a stochastic integral of a predictable function with respect to an integer-valued random measure; (iii) a continuous local martingale. Some applications to statistical inference problems are also presented. We prove the functional asymptotic normality of generalized Nelson-Aalen's estimator in the multiplicative intensity model for marked point processes. Its asymptotic efficiency in the sense of convolution theorem is also shown. The asymptotic behavior of log-likelihood ratio random fields of certain continuous semimartingales is derived. Received: 6 May 1996 / In revised form: 4 February 1997     

Godambe estimating functions and asymptotic optimal inference

Godambe (1985) introduced a class of optimum estimating functions which can be regarded as a generalization of quasilikelihood score functions. The “optimality” established by Godambe (1985) within a certain class is for estimating functions and it is based on finite samples. The question that arises naturally is what (if any) asymptotic optimality properties do the estimators and tests based on optimum estimating functions possess. In this paper, we establish, via presenting a convolution theorem, asymptotic optimality of estimators and tests obtained from Godambe optimum estimating functions. It is noted that we do not require the knowledge of the likelihood function.     

Asymptotic equivalence for nonparametric generalized linear models

Summary. We establish that a non-Gaussian nonparametric regression model is asymptotically equivalent to a regression model with Gaussian noise. The approximation is in the sense of Le Cam's deficiency distance Δ; the models are then asymptotically equivalent for all purposes of statistical decision with bounded loss. Our result concerns a sequence of independent but not identically distributed observations with each distribution in the same real-indexed exponential family. The canonical parameter is a value f(t _i ) of a regression function f at a grid point t _i (nonparametric GLM). When f is in a H?lder ball with exponent we establish global asymptotic equivalence to observations of a signal Γ(f(t)) in Gaussian white noise, where Γ is related to a variance stabilizing transformation in the exponential family. The result is a regression analog of the recently established Gaussian approximation for the i.i.d. model. The proof is based on a functional version of the Hungarian construction for the partial sum process. Received: 4 February 1997     

Asymptotics in Bayesian decision theory with applications to global robustness

We provide the rate of convergence of the Bayes action derived from non smooth loss functions involved in Bayesian robustness. Such loss functions are typically not twice differentiable but admit right and left second derivatives. The asymptotic limit of three measures of global robustness is given. These measures are the range of the Bayes actions set associated with a class of loss functions, the maximum regret of using a particular loss when the subjective loss belongs to a given class and the range of the posterior expected loss when the loss ranges over a given class. An application to prior robustness with density ratio classes is provided.     

Convergence analysis of a global optimization algorithm using stochastic differential equations

We establish the convergence of a stochastic global optimization algorithm for general non-convex, smooth functions. The algorithm follows the trajectory of an appropriately defined stochastic differential equation (SDE). In order to achieve feasibility of the trajectory we introduce information from the Lagrange multipliers into the SDE. The analysis is performed in two steps. We first give a characterization of a probability measure (Π) that is defined on the set of global minima of the problem. We then study the transition density associated with the augmented diffusion process and show that its weak limit is given by Π.     

Bias-corrected empirical likelihood in a multi-link semiparametric model

In this paper, we investigate the empirical likelihood for constructing a confidence region of the parameter of interest in a multi-link semiparametric model when an infinite-dimensional nuisance parameter exists. The new model covers the commonly used varying coefficient, generalized linear, single-index, multi-index, hazard regression models and their generalizations, as its special cases. Because of the existence of the infinite-dimensional nuisance parameter, the classical empirical likelihood with plug-in estimation cannot be asymptotically distribution-free, and the existing bias correction is not extendable to handle such a general model. We then propose a link-based correction approach to solve this problem. This approach gives a general rule of bias correction via an inner link, and consists of two parts. For the model whose estimating equation contains the score functions that are easy to estimate, we use a centering for the scores to correct the bias; for the model of which the score functions are of complex structure, a bias-correction procedure using simpler functions instead of the scores is given without loss of asymptotic efficiency. The resulting empirical likelihood shares the desired features: it has a chi-square limit and, under-smoothing technique, high order kernel and parameter estimation are not needed. Simulation studies are carried out to examine the performance of the new method.     

Parametric estimation of a bivariate stable Lévy process

We propose a parametric model for a bivariate stable Lévy process based on a Lévy copula as a dependence model. We estimate the parameters of the full bivariate model by maximum likelihood estimation. As an observation scheme we assume that we observe all jumps larger than some ε>0 and base our statistical analysis on the resulting compound Poisson process. We derive the Fisher information matrix and prove asymptotic normality of all estimates when the truncation point ε→0. A simulation study investigates the loss of efficiency because of the truncation.     

Statistical approach to some ill-posed problems for linear partial differential equations

For linear partial differential equations, some inverse source problems are treated statistically based on nonparametric estimation ideas. By observing the solution in a small Gaussian white noise, the kernel type of estimators is used to estimate the unknown source function and its partial derivatives.. It is proved that such estimators are consistent as the noise intensity tends to zero. Depending on the principal part of the differential operator, the optimal asymptotic rate of convergence is ascertained within a wide class of risk functions in a minimax sense. Received: 5 May 1997 / Revised version: 18 June 1998     

Dual divergence estimators and tests: Robustness results

The class of dual ?-divergence estimators (introduced in Broniatowski and Keziou (2009) [5]) is explored with respect to robustness through the influence function approach. For scale and location models, this class is investigated in terms of robustness and asymptotic relative efficiency. Some hypothesis tests based on dual divergence criteria are proposed and their robustness properties are studied. The empirical performances of these estimators and tests are illustrated by Monte Carlo simulation for both non-contaminated and contaminated data.     

Robust discrimination under a hierarchy on the scatter matrices

Under normality, Flury and Schmid [Quadratic discriminant functions with constraints on the covariances matrices: some asymptotic results, J. Multivariate Anal. 40 (1992) 244-261] investigated the asymptotic properties of the quadratic discrimination procedure under hierarchical models for the scatter matrices, that is: (i) arbitrary scatter matrices, (ii) common principal components, (iii) proportional scatter matrices and (iv) identical matrices. In this paper, we study the properties of robust quadratic discrimination rules based on robust estimates of the involved parameters. Our analysis is based on the partial influence functions of the functionals related to these parameters and allows to derive the asymptotic variances of the estimated coefficients under models (i)-(iv). From them, we conclude that the asymptotic variances verify the same order relations as those obtained by Flury and Schmid [Quadratic discriminant functions with constraints on the covariances matrices: some asymptotic results, J. Multivariate Anal. 40 (1992) 244-261] for the classical estimators. We also perform a Monte Carlo study for different sample sizes and different hierarchies which shows the advantage of using robust procedures over classical ones, when anomalous data are present. It also confirms that better rates of misclassification can be achieved if a more parsimonious model among all the correct ones is used instead of the standard quadratic discrimination.     

On Runge-Kutta-type methods for two-dimensional stochastic differential equations

In the multidimensional case, second-order weak Runge-Kutta methods for stochastic differential equation (SDE) need simulation of correlated random variables, unless the diffusion matrix of SDE satisfies the commutativity condition. In this paper, we show that this can be avoided for some types of diffusion matrices and test functions important for applications. Published in Lietuvos Matematikos Rinkinys, Vol. 46, No. 3, pp. 403–412, July–September, 2006.     

Asymptotics of Bayesian median loss estimation

We establish the consistency, asymptotic normality, and efficiency for estimators derived by minimizing the median of a loss function in a Bayesian context. We contrast this procedure with the behavior of two Frequentist procedures, the least median of squares (LMS) and the least trimmed squares (LTS) estimators, in regression problems. The LMS estimator is the Frequentist version of our estimator, and the LTS estimator approaches a median-based estimator as the trimming approaches 50% on each side. We argue that the Bayesian median-based method is a good tradeoff between the two Frequentist estimators.     

Modified Whittle estimation of multilateral models on a lattice

In the estimation of parametric models for stationary spatial or spatio-temporal data on a d-dimensional lattice, for d?2, the achievement of asymptotic efficiency under Gaussianity, and asymptotic normality more generally, with standard convergence rate, faces two obstacles. One is the “edge effect”, which worsens with increasing d. The other is the possible difficulty of computing a continuous-frequency form of Whittle estimate or a time domain Gaussian maximum likelihood estimate, due mainly to the Jacobian term. This is especially a problem in “multilateral” models, which are naturally expressed in terms of lagged values in both directions for one or more of the d dimensions. An extension of the discrete-frequency Whittle estimate from the time series literature deals conveniently with the computational problem, but when subjected to a standard device for avoiding the edge effect has disastrous asymptotic performance, along with finite sample numerical drawbacks, the objective function lacking a minimum-distance interpretation and losing any global convexity properties. We overcome these problems by first optimizing a standard, guaranteed non-negative, discrete-frequency, Whittle function, without edge-effect correction, providing an estimate with a slow convergence rate, then improving this by a sequence of computationally convenient approximate Newton iterations using a modified, almost-unbiased periodogram, the desired asymptotic properties being achieved after finitely many steps. The asymptotic regime allows increase in both directions of all d dimensions, with the central limit theorem established after re-ordering as a triangular array. However our work offers something new for “unilateral” models also. When the data are non-Gaussian, asymptotic variances of all parameter estimates may be affected, and we propose consistent, non-negative definite estimates of the asymptotic variance matrix.     

Superefficient estimation of the marginals by exploiting knowledge on the copula

We consider the problem of estimating the marginals in the case where there is knowledge on the copula. If the copula is smooth, it is known that it is possible to improve on the empirical distribution functions: optimal estimators still have a rate of convergence n^−1/2, but a smaller asymptotic variance. In this paper we show that for non-smooth copulas it is sometimes possible to construct superefficient estimators of the marginals: we construct both a copula and, exploiting the information our copula provides, estimators of the marginals with the rate of convergence logn/n.     

Purely sequential bounded-risk point estimation of the negative binomial mean under various loss functions: one-sample problem

A negative binomial (NB) distribution is useful to model over-dispersed count data arising from agriculture, health, and pest control. We design purely sequential bounded-risk methodologies to estimate an unknown NB mean \(\mu (>0)\) under different forms of loss functions including customary and modified Linex loss as well as squared error loss. We handle situations when the thatch parameter \(\tau (>0)\) may be assumed known or unknown. Our proposed methodologies are shown to satisfy properties including first-order asymptotic efficiency and first-order asymptotic risk efficiency. Summaries are provided from extensive sets of simulations showing encouraging performances of the proposed methodologies for small and moderate sample sizes. We follow with illustrations obtained by implementing estimation strategies using real data from statistical ecology: (1) weed count data of different species from a field in Netherlands and (2) count data of migrating woodlarks at the Hanko bird sanctuary in Finland.     

Asymptotically efficient two-sample rank tests for modal directions on spheres

A general class of optimal and distribution-free rank tests for the two-sample modal directions problem on (hyper-) spheres is proposed, along with an asymptotic distribution theory for such spherical rank tests. The asymptotic optimality of the spherical rank tests in terms of power-equivalence to the spherical likelihood ratio tests is studied, while the spherical Wilcoxon rank test, an important case for the class of spherical rank tests, is further investigated. A data set is reanalyzed and some errors made in previous studies are corrected. On the usual sphere, a lower bound on the asymptotic Pitman relative efficiency relative to Hotelling’s T²-type test is established, and a new distribution for which the spherical Wilcoxon rank test is optimal is also introduced.     

Determining <Emphasis Type="Italic">p</Emphasis>-values for systems cointegration tests with a prior adjustment for deterministic terms

Following Doornik (J Econ Surv 12:573–593, 1998) I present a procedure to approximate the asymptotic distributions of systems cointegration tests with a prior adjustment for deterministic terms suggested by Lütkepohl (Econometrica 72:647–662, 2004), Saikkonen and Lütkepohl (Econometric Theory 16:373–406, 2000a, J Business Econ Stat 18:451–464, 2000b, Time Series Anal 21:435–456, 2000c) and Saikkonen and Luukkonen (J Econ 81:93–126, 1997). These tests rely upon different assumptions as to the inclusion of deterministic components such as a constant, a linear trend or a level shift. The asymptotic distributions, which are functions of Brownian motions, are approximated by Gamma distributions. Only estimates of the mean and variance of the asymptotic test distributions are needed to fit the Gamma distributions. Such estimates are obtained from response surfaces. The required coefficients to compute the asymptotic moments are presented in this paper. Via the fitted Gamma distributions one can, then, easily derive p-values or arbitrary percentiles.