Statistical inference for simultaneous clustering of gene expression data