## Not run:
# full_model_fits <- fitModel(HSMM[sample(nrow(fData(HSMM_filtered)), 100),], modelFormulaStr="expression~sm.ns(Pseudotime)")
# expression_curve_matrix <- responseMatrix(full_model_fits)
# clusters <- clusterGenes(expression_curve_matrix, k=4)
# plot_clusters(HSMM_filtered[ordering_genes,], clusters)
# ## End(Not run)
Run the code above in your browser using DataLab