Model-based biclustering for expression data