A characterization of multivariate normality through univariate projections

This paper introduces a new characterization of multivariate normality of a random vector based on univariate normality of linear combinations of its components.     

A multivariate version of Hoeffding’s Phi-Square

A multivariate measure of association is proposed, which extends the bivariate copula-based measure Phi-Square introduced by Hoeffding [22]. We discuss its analytical properties and calculate its explicit value for some copulas of simple form; a simulation procedure to approximate its value is provided otherwise. A nonparametric estimator for multivariate Phi-Square is derived and its asymptotic behavior is established based on the weak convergence of the empirical copula process both in the case of independent observations and dependent observations from strictly stationary strong mixing sequences. The asymptotic variance of the estimator can be estimated by means of nonparametric bootstrap methods. For illustration, the theoretical results are applied to financial asset return data.     

Concordance measures for multivariate non-continuous random vectors

A notion of multivariate concordance suitable for non-continuous random variables is defined and many of its properties are established. This allows the definition of multivariate, non-continuous versions of Kendall’s tau, Spearman’s rho and Spearman’s footrule, which are concordance measures. Since the maximum values of these association measures are not +1 in general, a special attention is given to the computation of upper bounds. The latter turn out to be multivariate generalizations of earlier findings made by Nešlehová (2007) [9] and Denuit and Lambert (2005) [2]. They are easy to compute and can be estimated from a data set of (possibly) discontinuous random vectors. Corrected versions are considered as well.     

On a model selection problem from high-dimensional sample covariance matrices

Modern random matrix theory indicates that when the population size p is not negligible with respect to the sample size n, the sample covariance matrices demonstrate significant deviations from the population covariance matrices. In order to recover the characteristics of the population covariance matrices from the observed sample covariance matrices, several recent solutions are proposed when the order of the underlying population spectral distribution is known. In this paper, we deal with the underlying order selection problem and propose a solution based on the cross-validation principle. We prove the consistency of the proposed procedure.     

Asymptotic distributions of maxima of complete and incomplete samples from multivariate stationary Gaussian sequences

Let be a sequence of d-dimensional stationary Gaussian vectors, and let denote the partial maxima of . Suppose that there are missing data in each component of and let denote the partial maxima of the observed variables. In this note, we study two kinds of asymptotic distributions of the random vector where the correlation and cross-correlation satisfy some dependence conditions.     

On some testing problems for sparse contingency tables

A class of test statistics are proposed for sparse tables with ordered categories. It is shown that for different testing situations these test statistics asymptotically more powerful tests than Pearson's chi-square.     

Consistent multiple testing for change points

The change point problem for independent normal means is considered as a multiple testing problem. Two stepwise methods are considered. Namely, the binary segmentation method of Vostrikova (1981) [7] and the maximum residual down method of Cohen et al. (2009) [5]. Both of these methods are shown to be consistent. Consistent here means that as sample sizes tend to infinity, the probability of making an error (false rejection or false acceptance) tends to zero.     

Depth notions for orthogonal regression

Global depth, tangent depth and simplicial depths for classical and orthogonal regression are compared in examples, and properties that are useful for calculations are derived. The robustness of the maximum simplicial depth estimates is shown in examples. Algorithms for the calculation of depths for orthogonal regression are proposed, and tests for multiple regression are transferred to orthogonal regression. These tests are distribution free in the case of bivariate observations. For a particular test problem, the powers of tests that are based on simplicial depth and tangent depth are compared by simulations.     

On the limiting spectral distribution of the covariance matrices of time-lagged processes

We consider two continuous-time Gaussian processes, one being partially correlated to a time-lagged version of the other. We first give the limiting spectral distribution for the covariance matrices of the increments of the processes when the span between two observations tends to zero. Then, we derive the limiting distribution of the eigenvalues of the sample covariance matrices. This result is obtained when the number of paths of the processes is asymptotically proportional to the number of observations for each single path. As an application, we use the second moment of this distribution together with auxiliary volatility and correlation estimates to construct an adaptive estimator of the time lag between the two processes. Finally, we provide an asymptotic theory for our estimation procedure.     

Stochastic comparisons of multivariate mixtures

In this paper we establish multivariate hazard rate, multivariate reverse hazard rate, and multivariate likelihood ratio stochastic orderings among multivariate random mapping (mixture) distributions. The new results streamline and simplify the proofs of some partial results that have recently appeared in the literature. Some applications in reliability theory and risk management are described.     

On the layered nearest neighbour estimate, the bagged nearest neighbour estimate and the random forest method in regression and classification

Let be identically distributed random vectors in R^d, independently drawn according to some probability density. An observation is said to be a layered nearest neighbour (LNN) of a point if the hyperrectangle defined by and contains no other data points. We first establish consistency results on , the number of LNN of . Then, given a sample of independent identically distributed random vectors from R^d×R, one may estimate the regression function by the LNN estimate , defined as an average over the Y_i’s corresponding to those which are LNN of . Under mild conditions on r, we establish the consistency of towards 0 as n→∞, for almost all and all p≥1, and discuss the links between r_n and the random forest estimates of Breiman (2001) [8]. We finally show the universal consistency of the bagged (bootstrap-aggregated) nearest neighbour method for regression and classification.     

On signature-based expressions of system reliability

The concept of signature was introduced by Samaniego for systems whose components have i.i.d. lifetimes. This concept proved to be useful in the analysis of theoretical behaviors of systems. In particular, it provides an interesting signature-based representation of the system reliability in terms of reliabilities of k-out-of-n systems. In the non-i.i.d. case, we show that, at any time, this representation still holds true for every coherent system if and only if the component states are exchangeable. We also discuss conditions for obtaining an alternative representation of the system reliability in which the signature is replaced by its non-i.i.d. extension. Finally, we discuss conditions for the system reliability to have both representations.     

Subsampling tests for variance changes in the presence of autoregressive parameter shifts

In this paper, we consider the problem of testing for variance changes in the linear autoregressive processes including AR(p) processes when there are autoregressive parameter shifts. In performing a test, we employ the conventional residual CUSUM of squares test (RCUSQ) statistic. The RCUSQ test is based on the subsampling method introduced by Jach and Kokoszka (2004) [16] to eliminate the influence caused by autoregressive parameter shifts. It is shown that under regularity conditions, the test statistic behaves asymptotically the function of a standard Brownian bridge. We establish the asymptotic validity of this method and assess its performance both theoretically and numerically.     

Limit theorems in the Fourier transform method for the estimation of multivariate volatility

In this paper, we prove some limit theorems for the Fourier estimator of multivariate volatility proposed by Malliavin and Mancino (2002, 2009) [14] and [15]. In a general framework of discrete time observations we establish the convergence of the estimator and some associated central limit theorems with explicit asymptotic variance. In particular, our results show that this estimator is consistent for synchronous data, but possibly biased for non-synchronous observations. Moreover, from our general central limit theorem, we deduce that the estimator can be efficient in the case of a synchronous regular sampling. In the non-synchronous sampling case, the expression of the asymptotic variance is in general less tractable. We study this case more precisely through the example of an alternate sampling.     

The multiple hybrid bootstrap — Resampling multivariate linear processes

The paper reconsiders the autoregressive aided periodogram bootstrap (AAPB) which has been suggested in Kreiss and Paparoditis (2003) [18]. Their idea was to combine a time domain parametric and a frequency domain nonparametric bootstrap to mimic not only a part but as much as possible the complete covariance structure of the underlying time series. We extend the AAPB in two directions. Our procedure explicitly leads to bootstrap observations in the time domain and it is applicable to multivariate linear processes, but agrees exactly with the AAPB in the univariate case, when applied to functionals of the periodogram. The asymptotic theory developed shows validity of the multiple hybrid bootstrap procedure for the sample mean, kernel spectral density estimates and, with less generality, for autocovariances.     

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.     

Bayesian inference for dependent elliptical measurement error models

In this article we provide a Bayesian analysis for dependent elliptical measurement error models considering nondifferential and differential errors. In both cases we compute posterior distributions for structural parameters by using squared radial prior distributions for the precision parameters. The main result is that the posterior distribution of location parameters, for specific priors, is invariant with respect to changes in the generator function, in agreement with previous results obtained in the literature under different assumptions. Finally, although the results obtained are valid for any elliptical distribution for the error term, we illustrate those results by using the student-t distribution and a real data set.     

On the weighted multivariate Wilcoxon rank regression estimate

Zhou (2010) introduced a multivariate Wilcoxon regression estimate which possesses some nice properties: computational ease, asymptotic normality and high efficiency. However, it is sensitive to the leverage points. To circumvent this problem, we propose a weighted multivariate Wilcoxon regression estimate. Under some regularity conditions, the asymptotic normality is established. We further study the robustness of the proposed estimate through the influence function. By properly choosing the weight functions, our results show that the corresponding estimate can have bounded influence function on both response and covariates.     

Marginal parameterizations of discrete models defined by a set of conditional independencies

It is well-known that a conditional independence statement for discrete variables is equivalent to constraining to zero a suitable set of log–linear interactions. In this paper we show that this is also equivalent to zero constraints on suitable sets of marginal log–linear interactions, that can be formulated within a class of smooth marginal log–linear models. This result allows much more flexibility than known until now in combining several conditional independencies into a smooth marginal model. This result is the basis for a procedure that can search for such a marginal parameterization, so that, if one exists, the model is smooth.     

On weighting of bivariate margins in pairwise likelihood

Composite and pairwise likelihood methods have recently been increasingly used. For clustered data with varying cluster sizes, we study asymptotic relative efficiencies for various weighted pairwise likelihoods, with weight being a function of cluster size. For longitudinal data, we also study weighted pairwise likelihoods with weights that can depend on lag. Good choice of weights are needed to avoid the undesirable behavior of estimators with low efficiency. Some analytic results are obtained using the multivariate normal distribution. For clustered data, a practically good choice of weight is obtained after study of relative efficiencies for an exchangeable multivariate normal model; they are different from weights that had previously been suggested. For longitudinal data, there are advantages to only include bivariate margins of adjacent or nearly adjacent pairs in the weighted pairwise likelihood.