115
6

Sharp detection of smooth signals in a high-dimensional sparse matrix with indirect observations

Abstract

We consider a matrix-valued Gaussian sequence model, that is, we observe a sequence of high-dimensional M×NM \times N matrices of heterogeneous Gaussian random variables xij,kx_{ij,k} for i{1,...,M}i \in\{1,...,M\}, j{1,...,N}j \in \{1,...,N\} and kZk \in \mathbb{Z}. The standard deviation of our observations is \epks\ep k^s for some \ep>0\ep >0 and s0s \geq 0. We give sharp rates for the detection of a sparse submatrix of size m×nm \times n with active components. A component (i,j)(i,j) is said active if the sequence {xij,k}k\{x_{ij,k}\}_k have mean {θij,k}k\{\theta_{ij,k}\}_k within a Sobolev ellipsoid of smoothness τ>0\tau >0 and total energy kθij,k2\sum_k \theta^2_{ij,k} larger than some r\ep2r^2_\ep. Our rates involve relationships between m,n,Mm,\, n, \, M and NN tending to infinity such that m/Mm/M, n/Nn/N and \ep\ep tend to 0, such that a test procedure that we construct has asymptotic minimax risk tending to 0. We prove corresponding lower bounds under additional assumptions on the relative size of the submatrix in the large matrix of observations. Except for these additional conditions our rates are asymptotically sharp. Lower bounds for hypothesis testing problems mean that no test procedure can distinguish between the null hypothesis (no signal) and the alternative, i.e. the minimax risk for testing tends to 1.

View on arXiv
Comments on this paper