FindDEACC: Calculate precision, recall, and TNR of identified DE genes

In xnnba1984/DoubletCollection: Find doublets in scRNA-seq data

Calculate precision, recall, and TNR of identified DE genes

Description

Calculate precision, recall, and TNR of identified DE genes across different doublet-detection methods and DE methods.

Usage

FindDEACC(DE.list, DE.truth, gene.all)

Arguments

DE.list	A list of identified DE genes.
DE.truth	A list of ground-truth DE genes.
gene.all	A vector of all gene names.

Value

A list of precision, recall, and TNR across different doublet-detection methods and DE methods.

Examples

table.DE.all <- data.frame()
for(DE.method in c('MAST', 'wilcox', 'bimod')){
 DE.list <- FindDE(data.removal.list = data.removal.list, method.de = DE.method)
 DE.acc.list <- FindDEACC(DE.list=DE.list, DE.truth=data.de$gene.de, gene.all=rownames(data.de$count))
 table.DE <- ListToDataframe(l = DE.acc.list, type = 'barplot')
 table.DE[['DE_method']] <- DE.method
 table.DE.all <- rbind(table.DE.all, table.DE)
}

xnnba1984/DoubletCollection documentation built on Dec. 10, 2022, 11:13 a.m.