Giters

compbiomed / celda

Bayesian Hierarchical Modeling for Clustering Single Cell Genomic Data

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

reorder.celda_* Error when all cluster labels are the same

definitelysean opened this issue · comments

commented

In working with simulateCells.celda_CG, I realized if alpha is very small, all of the cells will often be pushed into the same cluster. This causes some downstream errors in reorder.celda_CG(), in which hclust() is run, but there aren't enough classes!

Should add an explicit check to the reorder function to ensure that not only is K > 2, but that the number of unique z is > 1. Same with L/y, and should be in all models.