Estimating the number of block boundaries from diagonal blockwise matrices without penalization

Vincent Brault (1) , Maud Delattre (1) , Tristan Mary-Huard (1) , Céline Leduc (1, 2)

1 MIA-Paris - Mathématiques et Informatique Appliquées
2 GQE-Le Moulon - Génétique Quantitative et Evolution - Le Moulon (Génétique Végétale)

In computational biology, numerous recent studies have been dedicated to the analysis of the chromatin structure within the cell by two-dimensional segmentation methods. Motivated by this application, we consider the problem of retrieving the diagonal blocks in a matrix of observations. The theoretical properties of the least squares estimators of both the boundaries and the number of blocks are investigated. More precisely, the contribution of the paper is to establish the consistency of these estimators. A surprising consequence of our results is that, contrary to the one-dimensional case, a penalty is not needed for retrieving the true number of diagonal blocks. Finally, the results are illustrated on synthetic data.

Mots clés: Segmentation; Hi-C data; Model selection

Domaines: Sciences du Vivant [q-bio]; Théorie [stat.TH]

Vincent Brault, Maud Delattre, Tristan Mary-Huard, Céline Leduc. Estimating the number of block boundaries from diagonal blockwise matrices without penalization. Scandinavian Journal of Statistics, 2017, 44 (2), pp.563-580. ⟨10.1111/sjos.12266⟩. ⟨hal-01533843⟩