Let $\widetilde X_{M\times N}$ be a rectangular data matrix with entries with
independent real valued entries $[\widetilde x_{ij}]$ satisfying $\mathbb E
\,\widetilde{x}_{ij} = 0$ and $\mathbb E \,\widetilde{x}^2_{ij} = {1 \over M}$,
$N, M\to \infty$ and these entries have a sub-exponential decay at the tails.
We will be working in the regime $N/M = d_N \in (0,\infty), \lim_{N\to
\infty}d_N \neq 1$.
We prove the universality of covariance matrices of the form $H_{N \times N}
= {1 \over N} \tp{X}X$ where $[X]_{M \times N}$ is a rectangular matrix with
independent real valued entries $[x_{ij}]$ satisfying $\E \,x_{ij} = 0$ and $\E
\,x^2_{ij} = {1 \over M}$, $N, M\to \infty$. Furthermore it is assumed that
these entries have sub-exponential tails. We will study the asymptotics in the
regime $N/M = d_N \in (0,\infty), \lim_{N\to \infty}d_N \neq 1$.
We propose a class of scale mixture of uniform distributions to generate
shrinkage priors for the covariance matrix. This new class of priors enjoys a
number of advantages over the traditional scale mixture of normal priors,
including its simplicity in characterizing the prior density based on its
first-order derivative and computationally efficiency based on a Gibbs sampler.
We first discuss the theory and computational details of this new approach for
the covariance matrix estimation.
We study random walk based algorithms for posterior simulation in a large of
class of Bayesian nonparametric problems with Gaussian random field or Gaussian
process priors. Our emphasis is both on developing practical guidelines for the
design and implementation of efficient algorithms for these naturally high
dimensional problems, and on the development of rigorous underpinning theory.
We illustrate via an example that, in designing algorithms for nonparametric
problems, it is important to take advantage of the infinite dimensional
structure inherent in both the prior and likelihood.
The Metropolis-adjusted Langevin (MALA) algorithm is a sampling algorithm
which makes local moves by incorporating information about the gradient of the
target density. In this paper we study the efficiency of MALA on a natural
class of target measures supported on an infinite dimensional Hilbert space.
These natural measures have density with respect to a Gaussian random field
measure and arise in many applications such as Bayesian nonparametric
statistics and the theory of conditioned diffusions.
Approximate Bayesian computation (ABC), also known as likelihood-free
methods, have become a favourite tool for the analysis of complex stochastic
models, primarily in population genetics but also in financial analyses. We
advocated in Grelaud et al.
Diffusion limits of MCMC methods in high dimensions provide a useful
theoretical tool for studying computational complexity. In particular they lead
directly to precise estimates of the number of steps required to explore the
target measure, in stationarity, as a function of the dimension of the state
space. However, to date such results have only been proved for target measures
with a product structure, severely limiting their applicability.
We investigate the properties of the Hybrid Monte-Carlo algorithm (HMC) in
high dimensions. HMC develops a Markov chain reversible w.r.t. a given target
distribution $\Pi$ by using separable Hamiltonian dynamics with potential
$-\log\Pi$. The additional momentum variables are chosen at random from the
Boltzmann distribution and the continuous-time Hamiltonian dynamics are then
discretised using the leapfrog scheme. The induced bias is removed via a
Metropolis-Hastings accept/reject rule.
We consider a family of stochastic processes $\{X_t^\epsilon, t \in T\}$ on a
metric space $T$, with a parameter $\epsilon \downarrow 0$. We study the
conditions under which \lim_{\e \to 0} \P \Big(\sup_{t \in T} |X_t^\e| < \delta
\Big) =1 when one has the \textit{a priori} estimate on the modulus of
continuity and the value at one point.
We demonstrate that stochastic differential equations (SDEs) driven by
fractional Brownian motion with Hurst parameter H > 1/2 have similar ergodic
properties as SDEs driven by standard Brownian motion. The focus in this
article is on hypoelliptic systems satisfying H\"ormander's condition. We show
that such systems satisfy a suitable version of the strong Feller property and
we conclude that they admit a unique stationary solution that is physical in
the sense that it does not "look into the future".