Fits models of gene expression evolution to expression data from coregulated groups of genes, assuming inverse gamma distributed rate variation.