A major goal of biomedical research is to identify molecular features associated with a biological or clinical class of interest. Differential expression analysis has long been used for this purpose, however, conventional methods perform poorly when applied to data with high within class heterogeneity.