A method for discovering gene regulatory models using relational fuzzy models. A method is provided that includes clustering gene expression data into a set of clusters and identifying a representative subset of genes from the set of clusters; and building a relational fuzzy model using the representative subset.