Gene set analysis is a powerful tool to study the coordinative change of time-course data. However, most existing methods only model the overall change of a gene set, yet completely overlooked heterogeneous time-dependent changes within sub-sets of genes.