BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Donald Richards

Donald Richards appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.STStatistics Theorymath.PRApplicationsComputation

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2026arXiv

Stein's method for the matrix normal distribution

This work presents the first systematic development of Stein's method for matrix distributions. We establish the basic essential ingredients of Stein's method for matrix normal approximation: we derive a generator-based Stein identity from a matrix Ornstein--Uhlenbeck diffusion with two-sided scales, provide an explicit semigroup representation for the solution of the Stein equation, and obtain regularity estimates for the solution. The new methodology is illustrated with three statistical applications, these being smooth Wasserstein distance bounds to quantify the matrix central limit theorem, a Wasserstein distance bound for the matrix normal approximation of the centered matrix $T$ distribution, and the derivation of Stein's method-of-moments estimators for scale parameters of the matrix normal distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Basic Treatment of the Distance Covariance

The distance covariance of Székely, et al. [23] and Székely and Rizzo [21], a powerful measure of dependence between sets of multivariate random variables, has the crucial feature that it equals zero if and only if the sets are mutually independent. Hence the distance covariance can be applied to multivariate data to detect arbitrary types of non-linear associations between sets of variables. We provide in this article a basic, albeit rigorous, introductory treatment of the distance covariance. Our investigations yield an approach that can be used as the foundation for presentation of this important and timely topic even in advanced undergraduate- or junior graduate-level courses on mathematical statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Continuous-Time Markov Chain Model for the Spread of COVID-19

Since late 2019 the novel coronavirus, also known as COVID-19, has caused a pandemic that persists. This paper shows how a continuous-time Markov chain model for the spread of COVID-19 can be used to explain, and justify to undergraduate students, strategies now being used in attempts to control the virus. The material in the paper is written at the level of students who are taking an introductory course on the theory and applications of stochastic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hoffmann-Jørgensen Inequalities for Random Walks on the Cone of Positive Definite Matrices

We consider random walks on the cone of $m \times m$ positive definite matrices, where the underlying random matrices have orthogonally invariant distributions on the cone and the Riemannian metric is the measure of distance on the cone. By applying results of Khare and Rajaratnam (Ann. Probab., 45 (2017), 4101--4111), we obtain inequalities of Hoffmann-Jørgensen type for such random walks on the cone. In the case of the Wishart distribution $W_m(a,I_m)$, with index parameter $a$ and matrix parameter $I_m$, the identity matrix, we derive explicit and computable bounds for each term appearing in the Hoffmann-Jørgensen inequalities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Product Inequalities for Multivariate Gaussian, Gamma, and Positively Upper Orthant Dependent Distributions

The Gaussian product inequality is an important conjecture concerning the moments of Gaussian random vectors. While all attempts to prove the Gaussian product inequality in full generality have been unsuccessful to date, numerous partial results have been derived in recent decades and we provide here further results on the problem. Most importantly, we establish a strong version of the Gaussian product inequality for multivariate gamma distributions in the case of nonnegative correlations, thereby extending a result recently derived by Genest and Ouimet [5]. Further, we show that the Gaussian product inequality holds with nonnegative exponents for all random vectors with positive components whenever the underlying vector is positively upper orthant dependent. Finally, we show that the Gaussian product inequality with negative exponents follows directly from the Gaussian correlation inequality.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shrinkage Estimation for the Diagonal Multivariate Exponential Families

We study shrinkage estimation of the mean parameters of a class of multivariate distributions for which the diagonal entries of the corresponding covariance matrix are certain quadratic functions of the mean parameter. This class of distributions includes the diagonal multivariate natural exponential families. We propose two classes of semi-parametric shrinkage estimators for the mean and construct unbiased estimators of the corresponding risk. We establish the asymptotic consistency and convergence rates for these shrinkage estimators under squared error loss as both $n$, the sample size, and $p$, the dimension, tend to infinity. Next, we specialize these results to the diagonal multivariate natural exponential families, which have been classified as consisting of the normal, Poisson, gamma, multinomial, negative multinomial, and hybrid classes of distributions. We establish the consistency of our estimators in the normal, gamma, and negative multinomial cases subject to the condition that $p n^{-1/3} (\log{n})^{4/3} \to 0$, and in the Poisson and multinomial cases if $p n^{-1/2} \to 0$, as $n,p \to \infty$. Simulation studies are provided to evaluate the performance of our estimators and we illustrate that, in the gamma and Poisson cases, our estimators achieve lower risk than the maximum likelihood estimator, thereby demonstrating the superiority of our estimators over the maximum likelihood estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Distance Correlation Coefficients for Lancaster Distributions

We consider the problem of calculating distance correlation coefficients between random vectors whose joint distributions belong to the class of Lancaster distributions. We derive under mild convergence conditions a general series representation for the distance covariance for these distributions. To illustrate the general theory, we apply the series representation to derive explicit expressions for the distance covariance and distance correlation coefficients for the bivariate normal distribution and its generalizations of Lancaster type, the multivariate normal distributions, and the bivariate gamma, Poisson, and negative binomial distributions which are of Lancaster type.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exact formulas for the normalizing constants of Wishart distributions for graphical models

Gaussian graphical models have received considerable attention during the past four decades from the statistical and machine learning communities. In Bayesian treatments of this model, the G-Wishart distribution serves as the conjugate prior for inverse covariance matrices satisfying graphical constraints. While it is straightforward to posit the unnormalized densities, the normalizing constants of these distributions have been known only for graphs that are chordal, or decomposable. Up until now, it was unknown whether the normalizing constant for a general graph could be represented explicitly, and a considerable body of computational literature emerged that attempted to avoid this apparent intractability. We close this question by providing an explicit representation of the G-Wishart normalizing constant for general graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Maximum Likelihood Estimation for Linear Gaussian Covariance Models

We study parameter estimation in linear Gaussian covariance models, which are $p$-dimensional Gaussian models with linear constraints on the covariance matrix. Maximum likelihood estimation for this class of models leads to a non-convex optimization problem which typically has many local maxima. Using recent results on the asymptotic distribution of extreme eigenvalues of the Wishart distribution, we provide sufficient conditions for any hill-climbing method to converge to the global maximum. Although we are primarily interested in the case in which $n>\!\!>p$, the proofs of our results utilize large-sample asymptotic theory under the scheme $n/p \to γ> 1$. Remarkably, our numerical simulations indicate that our results remain valid for $p$ as small as $2$. An important consequence of this analysis is that for sample sizes $n \simeq 14 p$, maximum likelihood estimation for linear Gaussian covariance models behaves as if it were a convex optimization problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chi-Square Mixture Representations for the Distribution of the Scalar Schur Complement in a Noncentral Wishart Matrix

We show that the distribution of the scalar Schur complement in a noncentral Wishart matrix is a mixture of central chi-square distributions with different degrees of freedom. For the case of a rank-1 noncentrality matrix, the weights of the mixture representation arise from a noncentral beta mixture of Poisson distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gaussian Random Particles with Flexible Hausdorff Dimension

Gaussian particles provide a flexible framework for modelling and simulating three-dimensional star-shaped random sets. In our framework, the radial function of the particle arises from a kernel smoothing, and is associated with an isotropic random field on the sphere. If the kernel is a von Mises--Fisher density, or uniform on a spherical cap, the correlation function of the associated random field admits a closed form expression. The Hausdorff dimension of the surface of the Gaussian particle reflects the decay of the correlation function at the origin, as quantified by the fractal index. Under power kernels we obtain particles with boundaries of any Hausdorff dimension between 2 and 3.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Generalization of an Integral Arising in the Theory of Distance Correlation

We generalize an integral which arises in several areas in probability and statistics and which is at the core of the field of distance correlation, a concept developed by Székely, Rizzo and Bakirov (2007) to measure dependence between random variables. Let $m$ be a positive integer and let ${\cos_m}(u)$, $u \in \mathbb{R}$, be the truncated Maclaurin expansion of ${\cos}(u)$, where the expansion is truncated at the $m$th summand. For $t, x \in \mathbb{R}^d$, let $\langle t,x\rangle$ and $\|x\|$ denote the standard Euclidean inner product and norm, respectively. We establish the integral formula: For $α\in \mathbb{C}$ and $x \in \mathbb{R}^d$, $\int_{\mathbb{R}^d} [\cos_m(\langle t,x\rangle) - \cos(\langle t,x\rangle)] \,{\rm d}t/{\|t\|^{d+α}} = C(d,α) \, \|x\|^α$, with absolute convergence if and only if $2(m-1) < \Re(α) < 2m$. Moreover, the constant $C(d,α)$ does not depend on $m$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

EM algorithms for estimating the Bernstein copula

A method that uses order statistics to construct multivariate distributions with fixed marginals and which utilizes a representation of the Bernstein copula in terms of a finite mixture distribution is proposed. Expectation-maximization (EM) algorithms to estimate the Bernstein copula are proposed, and a local convergence property is proved. Moreover, asymptotic properties of the proposed semiparametric estimators are provided. Illustrative examples are presented using three real data sets and a 3-dimensional simulated data set. These studies show that the Bernstein copula is able to represent various distributions flexibly and that the proposed EM algorithms work well for such data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The affinely invariant distance correlation

Székely, Rizzo and Bakirov (Ann. Statist. 35 (2007) 2769-2794) and Székely and Rizzo (Ann. Appl. Statist. 3 (2009) 1236-1265), in two seminal papers, introduced the powerful concept of distance correlation as a measure of dependence between sets of random variables. We study in this paper an affinely invariant version of the distance correlation and an empirical version of that distance correlation, and we establish the consistency of the empirical quantity. In the case of subvectors of a multivariate normally distributed random vector, we provide exact expressions for the affinely invariant distance correlation in both finite-dimensional and asymptotic settings, and in the finite-dimensional case we find that the affinely invariant distance correlation is a function of the canonical correlation coefficients. To illustrate our results, we consider time series of wind vectors at the Stateline wind energy center in Oregon and Washington, and we derive the empirical auto and cross distance correlation functions between wind vectors at distinct meteorological stations.

0 citations0 reviewsNo code linkedOpen access