Runnable code:
library(cola)
library(scRNAseq)
data = ReprocessedTh2Data()
mat = assays(data)$rsem_tpm
mat = log2(mat + 1)
mat = adjust_matrix(mat)
anno = colData(data)[, c("NREADS")]
anno = as.data.frame(anno)
rh = hierarchical_partition(mat, subset = 500, cores = 4, anno = anno)
saveRDS(rh, file = "ReprocessedTh2_cola_rh.rds")
cola_report(rh, output = "ReprocessedTh2_cola_rh_report", title = "cola Report for Hierarchical Partitioning - 'ReprocessedTh2'")
The HTML report is here. (generated by cola version 2.0.0, R 4.1.0.)