Runnable code:
library(cola)
library(scRNAseq)
data = LengESCData()
mat = assays(data)$normcounts
mat = log2(mat + 1)
mat = adjust_matrix(mat)
anno = colData(data)
anno = as.data.frame(anno)
rh = hierarchical_partition(mat, subset = 500, cores = 4, anno = anno)
saveRDS(rh, file = "LengESC_cola_rh.rds")
cola_report(rh, output = "LengESC_cola_rh_report", title = "cola Report for Hierarchical Partitioning - 'LengESC'")
The HTML report is here. (generated by cola version 2.0.0, R 4.1.0.)