A model-based clustering via mixture of hierarchical models with covariate adjustment for detecting differentially expressed genes from paired design Article in BMC Bioinformatics (November 2023)