Optimal Estimation of Bacterial Growth Rates Based on Permuted Monotone Matrix
Rong Ma, Tony Cai, and Hongzhe Li
Abstract:
Motivated by the problem of estimating the bacterial growth rates for genome assemblies from shotgun metagenomic data, we consider the permuted monotone matrix model Y = Θ Π + Z, where Y∈ ℝn×p is observed, Θ ∈ ℝn×p is an unknown approximately rank-one signal matrix with monotone rows, Π ∈ ℝn×p is an unknown permutation matrix, and Z∈ ℝn×p is the noise matrix. This paper studies the estimation of the extreme values associated to the signal matrix Θ, including its first and last columns, as well as their difference (the range vector). Treating these estimation problems as compound decision problems, minimax rate-optimal estimators are constructed using the spectral column sorting method. Numerical experiments through simulated and synthetic microbiome metagenomic data are presented, showing the superiority of the proposed methods over the alternatives. The methods are illustrated by comparing the growth rates of gut bacteria between inflammatory bowel disease patients and normal controls.