Search

This paper provides the first explicit formula for the expectation of the product of two disjoint principal minors of a Wishart random matrix, solving a part of a broader problem put forth by Samuel S. Wilks in 1934 in the Annals of Mathematics. The proof makes crucial use of hypergeometric functions of matrix argument and their Laplace transforms. Additionally, a Wishart generalization of the Gaussian product inequality conjecture is formulated and a stronger quantitative version is proved to hold in the case of two minors., Comment: 14 pages, 0 figures

This paper introduces a local linear smoother for regression surfaces on the simplex. The estimator solves a least-squares regression problem weighted by a locally adaptive Dirichlet kernel, ensuring excellent boundary properties. Asymptotic results for the bias, variance, mean squared error, and mean integrated squared error are derived, generalizing the univariate results of Chen (2002). A simulation study shows that the proposed local linear estimator with Dirichlet kernel outperforms its only direct competitor in the literature, the Nadaraya-Watson estimator with Dirichlet kernel due to Bouzebda, Nezzal, and Elhattab (2024)., Comment: 20 pages, 2 tables, 3 figures

In 1934, the American statistician Samuel S. Wilks derived remarkable formulas for the joint moments of embedded principal minors of sample covariance matrices in multivariate Gaussian populations, and he used them to compute the moments of sample statistics in various applications related to multivariate linear regression. These important but little-known moment results were extended in 1963 by the Australian statistician A. Graham Constantine using Bartlett's decomposition. In this note, a new proof of Wilks' results is derived using the concept of iterated Schur complements, thereby bypassing Bartlett's decomposition. Furthermore, Wilks' open problem of evaluating joint moments of disjoint principal minors of Wishart random matrices is related to the Gaussian product inequality conjecture., Comment: 6 pages, 0 figures

The large-sample behavior of non-degenerate multivariate $U$-statistics of arbitrary degree is investigated under the assumption that their kernel depends on parameters that can be estimated consistently. Mild regularity conditions are given which guarantee that once properly normalized, such statistics are asymptotically multivariate Gaussian both under the null hypothesis and sequences of local alternatives. The work of Randles (1982, \emph{Ann. Statist.}) is extended in three ways: the data and the kernel values can be multivariate rather than univariate, the limiting behavior under local alternatives is studied for the first time, and the effect of knowing some of the nuisance parameters is quantified. These results can be applied to a broad range of goodness-of-fit testing contexts, as shown in two specific examples., Comment: 22 pages, 3 figures

This paper extends various results related to the Gaussian product inequality (GPI) conjecture to the setting of disjoint principal minors of Wishart random matrices. This includes product-type inequalities for matrix-variate analogs of completely monotone functions and Bernstein functions of Wishart disjoint principal minors, respectively. In particular, the product-type inequalities apply to inverse determinant powers. Quantitative versions of the inequalities are also obtained when there is a mix of positive and negative exponents. Furthermore, an extended form of the GPI is shown to hold for the eigenvalues of Wishart random matrices by virtue of their law being multivariate totally positive of order~2 ($\mathrm{MTP}_2$). A new, unexplored avenue of research is presented to study the GPI from the point of view of elliptical distributions., Comment: 27 pages, 0 figures

This note presents a refined local approximation for the logarithm of the ratio between the negative multinomial probability mass function and a multivariate normal density, both having the same mean-covariance structure. This approximation, which is derived using Stirling's formula and a meticulous treatment of Taylor expansions, yields an upper bound on the Hellinger distance between the jittered negative multinomial distribution and the corresponding multivariate normal distribution. Upper bounds on the Le Cam distance between negative multinomial and multivariate normal experiments ensue., Comment: 12 pages, 0 figures

This paper introduces a novel density estimator supported on $d$-dimensional half-spaces. It stands out as the first asymmetric kernel smoother for half-spaces in the literature. Using the multivariate inverse Gaussian (MIG) density from Minami (2003) as the kernel and incorporating locally adaptive parameters, the estimator achieves desirable boundary properties. To analyze its mean integrated squared error (MISE) and asymptotic normality, a local limit theorem and probability metric bounds are established between the MIG and the corresponding multivariate Gaussian distribution with the same mean vector and covariance matrix, which may also be of independent interest. Additionally, a new algorithm for generating MIG random vectors is developed, proving to be faster and more accurate than Minami's algorithm based on a Brownian hitting location representation. This algorithm is then used to discuss and compare optimal MISE and likelihood cross-validation bandwidths for the estimator in a simulation study under various target distributions. As an illustration, the MIG asymmetric kernel is used to smooth the posterior distribution of a generalized Pareto model fitted to large electromagnetic storms., Comment: 45 pages, 6 figures, 2 tables

Bayes spaces were initially designed to provide a geometric framework for the modeling and analysis of distributional data. It has recently come to light that this methodology can be exploited to provide an orthogonal decomposition of bivariate probability distributions into an independent and an interaction part. In this paper, new insights into these results are provided by reformulating them using Hilbert space theory and a multivariate extension is developed using a distributional analog of the Hoeffding-Sobol identity. A connection between the resulting decomposition of a multivariate density and its copula-based representation is also provided.

This note reports partial results related to the Gaussian product inequality (GPI) conjecture for the joint distribution of traces of Wishart matrices. In particular, several GPI-related results from Wei (2014) and Liu et al. (2015) are extended in two ways: by replacing the power functions with more general classes of functions, and by replacing the usual Gaussian and multivariate gamma distributional assumptions by the more general trace-Wishart distribution assumption. These findings suggest that a Kronecker product form of the GPI holds for diagonal blocks of any Wishart distribution., Comment: 10 pages, 0 figures

A combinatorial proof of the Gaussian product inequality (GPI) is given under the assumption that each component of a centered Gaussian random vector $\boldsymbol{X} = (X_1, \ldots, X_d)$ of arbitrary length can be written as a linear combination, with coefficients of identical sign, of the components of a standard Gaussian random vector. This condition on $\boldsymbol{X}$ is shown to be strictly weaker than the assumption that the density of the random vector $(|X_1|, \ldots, |X_d|)$ is multivariate totally positive of order $2$, abbreviated MTP${}_2$, for which the GPI is already known to hold. Under this condition, the paper highlights a new link between the GPI and the monotonicity of a certain ratio of gamma functions., Comment: 9 pages, 1 figure

This paper considers the asymptotic behavior in $\beta$-H\"older spaces, and under $L^p$ losses, of a Dirichlet kernel density estimator proposed by Aitchison and Lauder (1985) for the analysis of compositional data. In recent work, Ouimet and Tolosana-Delgado (2022) established the uniform strong consistency and asymptotic normality of this estimator. As a complement, it is shown here that the Aitchison-Lauder estimator can achieve the minimax rate asymptotically for a suitable choice of bandwidth whenever $(p,\beta) \in [1, 3) \times (0, 2]$ or $(p, \beta) \in \mathcal{A}_d$, where $\mathcal{A}_d$ is a specific subset of $[3, 4) \times (0, 2]$ that depends on the dimension $d$ of the Dirichlet kernel. It is also shown that this estimator cannot be minimax when either $p \in [4, \infty)$ or $\beta \in (2, \infty)$. These results extend to the multivariate case, and also rectify in a minor way, earlier findings of Bertin and Klutchnikoff (2011) concerning the minimax properties of Beta kernel estimators., Comment: 18 pages, 1 figure

We answer an open question posed by the second author at the Salzburg workshop on Dependence Models and Copulas in 2016 concerning the size of the family $\mathcal{A}^\rho_v$ ($\mathcal{A}^\tau_v$) of all Pickands dependence functions $A$ whose corresponding Extreme-Value-Copulas have Spearman $\rho$ (Kendall $\tau$) equal to some arbitrary, fixed value $v \in [0,1]$. After determining compact sets $\Omega^\rho_v, \Omega^\tau_v \subseteq [0,1] \times [\frac{1}{2},1]$ containing the graphs of all Pickands dependence functions from the classes $\mathcal{A}^\rho_v$ and $\mathcal{A}^\tau_v$ respectively, we then show that both sets are best possible.

Continuation refers to the operation by which the cumulative distribution function of a discontinuous random vector is made continuous through multilinear interpolation. The copula that results from the application of this technique to the classical empirical copula is either called the multilinear or the checkerboard copula. As shown by Genest and Ne\v{s}lehov\'{a} (Astin Bull. 37 (2007) 475-515) and Ne\v{s}lehov\'{a} (J. Multivariate Anal. 98 (2007) 544-567), this copula plays a central role in characterizing dependence concepts in discrete random vectors. In this paper, the authors establish the asymptotic behavior of the empirical process associated with the multilinear copula based on $d$-variate count data. This empirical process does not generally converge in law on the space $\mathcal {C}([0,1]^d)$ of continuous functions on $[0,1]^d$, equipped with the uniform norm. However, the authors show that the process converges in $\mathcal{C}(K)$ for any compact $K\subset\mathcal{O}$, where $\mathcal{O}$ is a dense open subset of $[0,1]^d$, whose complement is the Cartesian product of the ranges of the marginal distribution functions. This result is sufficient to deduce the weak limit of many functionals of the process, including classical statistics for monotone trend. It also leads to a powerful and consistent test of independence which is applicable even to sparse contingency tables whose dimension is sample size dependent., Comment: Published in at http://dx.doi.org/10.3150/13-BEJ524 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

It is often reasonable to assume that the dependence structure of a bivariate continuous distribution belongs to the class of extreme-value copulas. The latter are characterized by their Pickands dependence function. In this paper, a procedure is proposed for testing whether this function belongs to a given parametric family. The test is based on a Cram\'{e}r--von Mises statistic measuring the distance between an estimate of the parametric Pickands dependence function and either one of two nonparametric estimators thereof studied by Genest and Segers [Ann. Statist. 37 (2009) 2990--3022]. As the limiting distribution of the test statistic depends on unknown parameters, it must be estimated via a parametric bootstrap procedure, the validity of which is established. Monte Carlo simulations are used to assess the power of the test and an extension to dependence structures that are left-tail decreasing in both variables is considered., Comment: Published in at http://dx.doi.org/10.3150/10-BEJ279 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

Martin Bradbury Wilk was born on December 18, 1922, in Montr\'{e}al, Qu\'{e}bec, Canada. He completed a B.Eng. degree in Chemical Engineering in 1945 at McGill University and worked as a Research Engineer on the Atomic Energy Project for the National Research Council of Canada from 1945 to 1950. He then went to Iowa State College, where he completed a M.Sc. and a Ph.D. degree in Statistics in 1953 and 1955, respectively. After a one-year post-doc with John Tukey, he became Assistant Director of the Statistical Techniques Research Group at Princeton University in 1956--1957, and then served as Professor and Director of Research in Statistics at Rutgers University from 1959 to 1963. In parallel, he also had a 14-year career at Bell Laboratories, Murray Hill, New Jersey. From 1956 to 1969, he was in turn Member of Technical Staff, Head of the Statistical Models and Methods Research Department, and Statistical Director in Management Sciences Research. He wrote a number of influential papers in statistical methodology during that period, notably testing procedures for normality (the Shapiro--Wilk statistic) and probability plotting techniques for multivariate data. In 1970, Martin moved into higher management levels of the American Telephone and Telegraph (AT&T) Company. He occupied various positions culminating as Assistant Vice-President and Director of Corporate Planning. In 1980, he returned to Canada and became the first professional statistician to serve as Chief Statistician. His accomplishments at Statistics Canada were numerous and contributed to a resurgence of the institution's international standing. He played a crucial role in the reinstatement of the Cabinet-cancelled 1986 Census., Comment: Published in at http://dx.doi.org/10.1214/08-STS272 the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Conditions are given under which the empirical copula process associated with a random sample from a bivariate continuous distribution has a smaller asymptotic covariance function than the standard empirical process based on observations from the copula. Illustrations are provided and consequences for inference are outlined., Comment: 14 pages, 2 figures

November 27, 2004, marked the 250th anniversary of the death of Abraham De Moivre, best known in statistical circles for his famous large-sample approximation to the binomial distribution, whose generalization is now referred to as the Central Limit Theorem. De Moivre was one of the great pioneers of classical probability theory. He also made seminal contributions in analytic geometry, complex analysis and the theory of annuities. The first biography of De Moivre, on which almost all subsequent ones have since relied, was written in French by Matthew Maty. It was published in 1755 in the Journal britannique. The authors provide here, for the first time, a complete translation into English of Maty's biography of De Moivre. New material, much of it taken from modern sources, is given in footnotes, along with numerous annotations designed to provide additional clarity to Maty's biography for contemporary readers., Comment: Published at http://dx.doi.org/10.1214/088342306000000268 in the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Deheuvels [J. Multivariate Anal. 11 (1981) 102--113] and Genest and R\'{e}millard [Test 13 (2004) 335--369] have shown that powerful rank tests of multivariate independence can be based on combinations of asymptotically independent Cram\'{e}r--von Mises statistics derived from a M\"{o}bius decomposition of the empirical copula process. A result on the large-sample behavior of this process under contiguous sequences of alternatives is used here to give a representation of the limiting distribution of such test statistics and to compute their relative local asymptotic efficiency. Local power curves and asymptotic relative efficiencies are compared under familiar classes of copula alternatives., Comment: Published at http://dx.doi.org/10.1214/009053606000000984 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Consider a continuous random pair $(X,Y)$ whose dependence is characterized by an extreme-value copula with Pickands dependence function $A$. When the marginal distributions of $X$ and $Y$ are known, several consistent estimators of $A$ are available. Most of them are variants of the estimators due to Pickands [Bull. Inst. Internat. Statist. 49 (1981) 859--878] and Cap\'{e}ra\`{a}, Foug\`{e}res and Genest [Biometrika 84 (1997) 567--577]. In this paper, rank-based versions of these estimators are proposed for the more common case where the margins of $X$ and $Y$ are unknown. Results on the limit behavior of a class of weighted bivariate empirical processes are used to show the consistency and asymptotic normality of these rank-based estimators. Their finite- and large-sample performance is then compared to that of their known-margin analogues, as well as with endpoint-corrected versions thereof. Explicit formulas and consistent estimates for their asymptotic variances are also given., Comment: Published in at http://dx.doi.org/10.1214/08-AOS672 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)