## Not run:
# full_model_fits <- fitModel(HSMM_filtered[sample(nrow(fData(HSMM_filtered)), 100),], modelFormulaStr="expression~VGAM::bs(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