2 results
Search Results
Now showing 1 - 2 of 2
- ItemBayesian inference for spectral projectors of the covariance matrix(Ithaca, NY : Cornell University Library, 2018) Silin, Igor; Spokoiny, VladimirLet X1,…,Xn be an i.i.d. sample in Rp with zero mean and the covariance matrix Σ∗. The classical PCA approach recovers the projector P∗J onto the principal eigenspace of Σ∗ by its empirical counterpart ˆPJ. Recent paper [24] investigated the asymptotic distribution of the Frobenius distance between the projectors ∥ˆPJ−P∗J∥2, while [27] offered a bootstrap procedure to measure uncertainty in recovering this subspace P∗J even in a finite sample setup. The present paper considers this problem from a Bayesian perspective and suggests to use the credible sets of the pseudo-posterior distribution on the space of covariance matrices induced by the conjugated Inverse Wishart prior as sharp confidence sets. This yields a numerically efficient procedure. Moreover, we theoretically justify this method and derive finite sample bounds on the corresponding coverage probability. Contrary to [24, 27], the obtained results are valid for non-Gaussian data: the main assumption that we impose is the concentration of the sample covariance ˆΣ in a vicinity of Σ∗. Numerical simulations illustrate good performance of the proposed procedure even on non-Gaussian data in a rather challenging regime.
- ItemGaussian processes with multidimensional distribution inputs via optimal transport and Hilbertian embedding(Ithaca, NY : Cornell University Library, 2020) Bachoc, François; Suvorikova, Alexandra; Ginsbourger, David; Loubes, Jean-Michel; Spokoiny, VladimirIn this work, we propose a way to construct Gaussian processes indexed by multidimensional distributions. More precisely, we tackle the problem of defining positive definite kernels between multivariate distributions via notions of optimal transport and appealing to Hilbert space embeddings. Besides presenting a characterization of radial positive definite and strictly positive definite kernels on general Hilbert spaces, we investigate the statistical properties of our theoretical and empirical kernels, focusing in particular on consistency as well as the special case of Gaussian distributions. A wide set of applications is presented, both using simulations and implementation with real data.