Tuesday, May 30 |
07:30 - 09:00 |
Breakfast (Restaurant at your assigned hotel) |
09:00 - 09:45 |
Vladimir Koltchinskii: Efficient Estimation of Smooth Functionals of Covariance Operators ↓ A problem of efficient estimation of a smooth functional of unknown
covariance operator Σ of a mean zero Gaussian random variable in a Hilbert
space based on a sample of its i.i.d. observations will be discussed.
The goal is to find an estimator whose distribution is approximately
normal with a minimax optimal variance in a setting when either the dimension of the space, or so called effective
rank of the covariance operator are allowed to be large (although much smaller than
the sample size). This problem has been recently solved in our joint paper with Loeffler and Nickl in the case
of estimation of a linear functional of unknown eigenvector of Σ
corresponding to its largest eigenvalue (the top principal component).
The efficient estimator developed in this paper does not coincide
with the naive estimator based on the top principal component of
sample covariance which is not efficient due to its large bias.
An approach to a more general problem of efficient estimation of a functional
⟨f(Σ),B⟩ for a given sufficiently smooth function f:R↦R
and given operator B will be also discussed. (Conference Room San Felipe) |
09:45 - 10:30 |
Florence Merlevède: On strong approximations for some classes of random iterates ↓ his talk is devoted to strong approximations in the dependent setting. The famous results of Koml\'os, Major and Tusn\'ady (1975-1976) state that it
is possible to approximate almost surely the partial sums of size n of i.i.d. centered random
variables in Lp (p>2) by a Wiener process with an
error term of order o(n1/p). In the case of functions of random iterates generated by an iid sequence, we
we shall give new dependent conditions, expressed in terms of a natural coupling (in L∞ or in L1), under which the strong approximation result holds with rate
o(n1/p). The proof is an adaptation of the recent construction given in Berkes, Liu and Wu (2014).
As we shall see our conditions are well adapted to a large variety of models, including left random
walks on GLd(R), contracting iterated random functions, autoregressive Lipschitz processes, and some ergodic Markov chains.
We shall also provide some examples showing that our L1-coupling condition is in some sense optimal. This talk is based on a joint work with J. Dedecker and C. Cuny. (Conference Room San Felipe) |
10:30 - 11:00 |
Coffee Break (Conference Room San Felipe) |
11:00 - 11:45 |
Radoslaw Adamczak: Uncertainty relations for high dimensional unitary matrices ↓ I will present various types of uncertainty relations satisfied by Haar distributed random unitary matrices in high dimensions. If time permits I will also discuss their applications to quantum information theory (in particular to the problems of information locking and data hiding) and to special cases of the Dvoretzky theorem. (Conference Room San Felipe) |
11:45 - 12:30 |
Stanislav Minsker: Random Matrices with Heavy-Tailed Entries: Tight Mean Estimators and Applications to Statistics ↓ Estimation of the covariance matrix has attracted significant attention of the statistical research community over the years, partially due to important applications such as Principal Component Analysis. However, frequently used empirical covariance estimator (and its modifications) is very sensitive to outliers, or ``atypical’’ points in the sample.
As P. Huber wrote in 1964, “...This raises a question which could have been asked already by Gauss, but which was, as far as I know, only raised a few years ago (notably by Tukey): what happens if the true distribution deviates slightly from the assumed normal one? As is now well known, the sample mean then may have a catastrophically bad performance…”
Motivated by Tukey's question, we develop a new estimator of the (element-wise) mean of a random matrix, which includes covariance estimation problem as a special case. Assuming that the entries of a matrix possess only finite second moment, this new estimator admits sub-Gaussian or sub-exponential concentration around the unknown mean in the operator norm.
Our arguments rely on generic chaining techniques applied to operator-valued stochastic processes, as well as bounds on the trace moment-generating function.
We will discuss extensions of our approach to matrix-valued U-statistics and examples such as matrix completion problem.
Part of the talk will be based on a joint work with Xiaohan Wei. (Conference Room San Felipe) |
13:30 - 15:00 |
Lunch (Restaurant Hotel Hacienda Los Laureles) |
15:00 - 15:45 |
Victor Perez Abreu: Some extensions of the chain: From Hermitian Brownian motion to Dyson-Brownian process to free Brownian motion (Conference Room San Felipe) |
15:45 - 16:30 |
Adelaide Olivier: Estimation of the division rate in growth-fragmentation models ↓ This talk is concerned with growth-fragmentation models, implemented for investigating the growth of a population of cells. From a stochastic point of view, we are dealing with the evolution of a system of particles. The evolution of the system is then driven by two phenomenons. First, particles evolve on a deterministic basis: they age, they grow. Secondly, particles split randomly: a particle of age a or size x splits into two particles (of age 0, of size at birth x/2), at a rate B depending on the age or on the size of the splitting particle. The division rate B is an unknown function, on which we focus our attention.
The main goal of this statistical work is to get a nonparametric estimate of the division rate. To do so, various observation schemes may be considered: {\bf (i)} Observation {\it up to a given generation} in the genealogical tree of the population. {\bf (ii)} Continuous time observation {\it up to time T}. This scheme, which differs radically from those previously used, entails specific difficulties -- mainly a bias selection.
As different as these two schemes may be, I will try to highlight their common features. Mainly, a competition occurs between the growth of the tree (measured by the Malthus parameter) and the convergence to equilibrium of the branching process (measured by some ρB say). We still have to deal with open questions in these toy models: adaptativity with respect to the smothness of B - which requires some deviation inequalities (availaible to treat the observation scheme (i) but not (ii)), but also adaptivity with respect to the parameter ρB measuring the convergence to equilibrium. (Conference Room San Felipe) |
16:30 - 17:00 |
Coffee Break (Conference Room San Felipe) |
17:00 - 17:45 |
Rafal Meller: Two-sided moment estimates for random chaoses. ↓ Let X1,…,Xn be random variables such that there exists a constant C>1 satisfying ‖ for every p \geq 1.
We define random chaos S=\sum a_{i_1,...,i_d} X_{i_1}\cdots X_{i_d}. We will show two-sided deterministic bounds on ||S||_p, with constant depending only on C and d in two cases:
1) X_1,\ldots,X_n are nonnegative and a_{i_1,...,i_d}\geq 0.
2) X_1,\ldots ,X_n are symmetric, d=2. (Conference Room San Felipe) |
19:00 - 21:00 |
Dinner (Restaurant Hotel Hacienda Los Laureles) |