L p –Estimators as Estimates of a Parameter of Location for a Sharp–pointed Symmetric Density

We study the asymptotics of L _p estimators, p > 0, over a sample having a symmetric density with a sharp–point at the centre of symmetry of the distribution. The rates of convergence of the L _p estimators in this situation depend on p and on the shape of the density. To obtain some of the limit distributions, we present new results in the asymptotics of M–estimators. We extend the delta method to the case when the Euclidean norm of the conveniently normalized M–estimators converge to a power of the Euclidean norm of a (possibly Gaussian) stable distribution.     

Sequential monitoring of high-dimensional time series

In the paper we derive new types of multivariate exponentially weighted moving average (EWMA) control charts which are based on the Euclidean distance and on the distance defined by using the inverse of the diagonal matrix consisting of the variances. The design of the proposed control schemes does not involve the computation of the inverse covariance matrix and, thus, it can be used in the high-dimensional setting. The distributional properties of the control statistics are obtained and are used in the determination of the new control procedures. Within an extensive simulation study, the new approaches are compared with the multivariate EWMA control charts which are based on the Mahalanobis distance.     

Riesz inverse Gaussian distributions

The notion of generalized power of a positive definite symmetric matrix and a related notion of generalized Bessel function are used to introduce an extension of the class of matrix generalized inverse Gaussian distributions. The new distributions are shown to arise as conditional distributions of Peirce components of Riesz random matrices. Things are explained in the modern framework of symmetric cones and simple Euclidean Jordan algebra.     

Invariance In the linear model – an argument for χ2 and F in nonnormal situations

For a suitable norm, conservation of the distance between expectation and hypothesis may furnish a basis for data reduction by invariance in the linear, not neces-sarily normal, model. If the norm is Euclidean (i.e. based on some inner product), the maximal invariant is a pair of sums of squares. This provides support for traditional χ² (or F) - methods also in nonnormal cases. If the norm is l^p p≠2, or the supnorm, the maximal invariant is, at the best a air of order statistics.     

Some Matrix Norm Kantorovich Inequalities and Their Applications

In this article, we first present four matrix norm Kantorovich-type inequalities involving non negative definite matrix. Then, based on these inequalities, we propose four new efficiency criteria and present their lower bounds to make efficiency comparisons between the ordinary least squares estimator and the best linear unbiased estimator in a singular linear model.     

A cubic smoothing spline based lack of fit test for nonlinear regression models

The use of a statistic based on cubic spline smoothing is considered for testing nonlinear regression models for lack of fit. The statistic is defined to be the Euclidean squared norm of the smoothed residual vector obtained from fitting the nonlinear model, The asymptotic distribution of the statistic is derived under suitable smooth local alternatives and a numerical example is presented.     

Asymptotic Minimax Risk for the White Noise Model on the Sphere

Estimation of an unknown function on the unit sphere of the Euclidean space is considered. The function is observed in Gaussian continuous time white noise. Uniform norm is chosen as a loss function and exact asymptotic minimax risk is derived extending the result of Korostelev (1993). The exact asymptotic minimax risk is also given for the L ₂-loss, applying the result of Pinsker (1980).     

Multidimensional scaling models with inequality and equality constraints

Recently, a lot of attention has been brought to constrained estimation theory in multidimensional scaling models. So far, only equality constraints have been thoroughly studied. In this paper, the optimization theory is extended to general multidi-mensional scaling models with both inequality and equality constraints. A Newton-Raphson based algorithm is developed to produce the constrained least squares estimate. To illustrate the theory, some classical color data are reanalyzed in the context of the linear Euclidean distance model.     

Matrix trace Wielandt inequalities with statistical applications

The vector correlation coefficient and other measures of association play a very important role in statistics and especially in multivariate analysis. In this paper a new measure of association is proposed and its upper bound is presented by using a matrix trace Wielandt inequality. Also given are relevant results involving Wishart matrices widely used in multivariate analysis, and especially a new alternative for the relative gain of the covariance adjusted estimator of a vector of parameters.     

On Stein Identity,Chernoff Inequality,and Orthogonal Polynomials

In this article, we present a general method for deriving Stein-like identity and Chernoff-like inequality based on orthogonal polynomials. In order to illustrate our method, some applications are given with respect to normal, Gamma, Beta, Poisson, binomial, and negative binomial distribution, not only for random variables but also for random vectors, resulting corresponding Stein-like identity and Chernoff-like inequality are obtained consequently. Within our best knowledge, some of our matrix version results are new in the literature. In addition, forward difference formulae of Charlier polynomials, Krawtchouk polynomials and Meixner polynomials, Stein-like identity, and Chernoff-like inequality with respect to Beta distribution, as well as Rodrigues formula of Meixner polynomials are also prepared in the first time within our limited information. Interestingly, as far as normal, Gamma, Beta, Poisson, binomial, and negative binomial distribution are concerned, we found that their Stein-like identity and corresponding Chernoff-like inequality are related closely, by examining their Rodrigues formula.     

Moment Consistency of the Exchangeably Weighted Bootstrap for Semiparametric M‐estimation

The bootstrap variance estimate is widely used in semiparametric inferences. However, its theoretical validity is a well‐known open problem. In this paper, we provide a first theoretical study on the bootstrap moment estimates in semiparametric models. Specifically, we establish the bootstrap moment consistency of the Euclidean parameter, which immediately implies the consistency of t‐type bootstrap confidence set. It is worth pointing out that the only additional cost to achieve the bootstrap moment consistency in contrast with the distribution consistency is to simply strengthen the L₁ maximal inequality condition required in the latter to the L_p maximal inequality condition for p≥1. The general L_p multiplier inequality developed in this paper is also of independent interest. These general conclusions hold for the bootstrap methods with exchangeable bootstrap weights, for example, non‐parametric bootstrap and Bayesian bootstrap. Our general theory is illustrated in the celebrated Cox regression model.     

一种基于加权F-范数的半正定矩阵的逼近方法

正定性是许多金融预测模型的重要假设前提，然而从实际样本中得到的相关系数矩阵并不能保证其正定性。为此在介绍如何根据样本设定相关系数矩阵以及范数逼近原理的基础上，如何根据该原理找到与之最接近的相关系数矩阵，即最接近的单位对角半正定对称矩阵。通过实证，验证了其方法的有效性。     

On matrix variance inequalities

Olkin and Shepp [2005, A matrix variance inequality. J. Statist. Plann. Inference 130, 351-358] presented a matrix form of Chernoff's inequality for Normal and Gamma (univariate) distributions. We extend and generalize this result, proving Poincaré-type and Bessel-type inequalities, for matrices of arbitrary order and for a large class of distributions.     

Properties of a randomization test for multifactor comparisons of groups

Multivariate hypothesis testing in studies of vegetation is likely to be hindered by unrealistic assumptions when based on conventional statistical methods. This can be overcome by randomization tests. In this paper, the accuracy and power of a MANOVA randomization test are evaluated for one and two factors with interaction with simulated data from three distributions. The randomization test is based on the partitioning of sum of squares computed from Euclidean distances. In one-factor designs, sample size and variance inequality were evaluated. The results showed a high level of accuracy. The power curve was higher with normal distribution, lower with uniform, intermediate with lognormal and was sensitive to variance inequality. In two-factor designs, three methods of permutations and two statistics were compared. The results showed that permutation of the residuals with F _pseudo is accurate and can give good power for testing the interaction and restricted permutation for testing main factors.     

析主成分分析和因子分析——以Frobenius范数下矩阵近似的新视角

本文通过引入数据阵在Frobenius范数下的最优近似等概念来重新探讨主成分和因子分析。我发现，主成分分析中主成分和因子分析中因子得分（通过主成分解因子载荷，然后用最小二乘解因子得分）的估计为数据阵的最优近似（在Frobenius范数下）在不同正交坐标方向矩阵下的坐标。两种方法分别采用了不同的约束条件分解的最优近似（在Frobenius范数下），因为该分解并不唯一。     

Sur un test d'égalité des autocovariances de deux séries chronologiques

We propose a test for the equality of the autocovariance functions of two independent and stationary time series. The test statistic is a quadratic form in the vector of differences of the first J + 1 autocovariances. Its asymptotic distribution is derived under the null hypothesis, and the finite-sample properties of the test, namely the bias and the power, are investigated by Monte Carlo methods. A by-product of this study is a new estimator of the covariance between two sample autocovariances which provides a positive definite covariance matrix. We establish the convergence of this estimator in the L₁ norm.     

Estimating and testing generalized linear models under inequality restrictions

We consider maximum likelihood estimation and likelihood ratio tests under inequality restrictions on the parameters. A special case are order restrictions, which may appear for example in connection with effects of an ordinal qualitative covariate. Our estimation approach is based on the principle of sequential quadratic programming, where the restricted estimate is computed iteratively and a quadratic optimization problem under inequality restrictions is solved in each iteration. Testing for inequality restrictions is based on the likelihood ratio principle. Under certain regularity assumptions the likelihood ratio test statistic is asymptotically distributed like a mixture of χ², where the weights are a function of the restrictions and the information matrix. A major problem in theory is that in general there is no unique least favourable point. We present some empirical findings on finite-sample behaviour of tests and apply the methods to examples from credit scoring and dentistry.     

Exploring the variability of DNA molecules via principal geodesic analysis on the shape space

Most of the linear statistics deal with data lying in a Euclidean space. However, there are many examples, such as DNA molecule topological structures, in which the initial or the transformed data lie in a non-Euclidean space. To get a measure of variability in these situations, the principal component analysis (PCA) is usually performed on a Euclidean tangent space as it cannot be directly implemented on a non-Euclidean space. Instead, principal geodesic analysis (PGA) is a new tool that provides a measure of variability for nonlinear statistics. In this paper, the performance of this new tool is compared with that of the PCA using a real data set representing a DNA molecular structure. It is shown that due to the nonlinearity of space, the PGA explains more variability of the data than the PCA.     

A New Method for Statistical Multidimensional Unfolding

ABSTRACT

Consider the problem of estimating the positions of a set of targets in a multidimensional Euclidean space from distances reported by a number of observers when the observers do not know their own positions in the space. Each observer reports the distance from the observer to each target plus a random error. This statistical problem is the basic model for the various forms of what is called multidimensional unfolding in the psychometric literature. Multidimensional unfolding methodology as developed in the field of cognitive psychology is basically a statistical estimation problem where the data structure is a set of measures that are monotonic functions of Euclidean distances between a number of observers and targets in a multidimensional space. The new method presented in this article deals with estimating the target locations and the observer positions when the observations are functions of the squared distances between observers and targets observed with an additive random error in a two-dimensional space. The method provides robust estimates of the target locations in a multidimensional space for the parametric structure of the data generating model presented in the article. The method also yields estimates of the orientation of the coordinate system and the mean and variances of the observer locations. The mean and the variances are not estimated by standard unfolding methods which yield targets maps that are invariant to a rotation of the coordinate system. The data is transformed so that the nonlinearity due to the squared observer locations is removed. The sampling properties of the estimates are derived from the asymptotic variances of the additive errors of a maximum likelihood factor analysis of the sample covariance matrix of the transformed data augmented with bootstrapping. The robustness of the new method is tested using artificial data. The method is applied to a 2001 survey data set from Turkey to provide a real data example.     

AFFINE‐EQUIVARIANT SPATIAL MEDIAN AND ITS USE IN THE MULTIVARIATE MULTI‐SAMPLE LOCATION PROBLEM

The classical spatial median is not affine‐equivariant, which often turns out to be an unfavourable property. In this paper, the asymptotic properties of an affine‐equivariant modification of the spatial median are investigated. It is shown that under some weak regularity conditions, the modified spatial median computed by means of the sample norming matrix is asymptotically equivalent to the one computed by means of the population norming matrix, which yields its asymptotic normality. A consistent estimate of the asymptotic covariance matrix of the modified spatial median is also presented. These results are implemented in a scheme, where the sample norm is determined by means of the sample Dümbgen scatter matrix. The results are utilized in the construction of affine‐invariant test statistics for testing the multi‐sample hypothesis of equality of location parameters. The performance of the proposed tests is demonstrated through a simulation study.