modified RNA-seq.R in R (3349f8ff) · Commits · Chaos / chaos.tools

R/RNA-seq.R

+9 −2

Original line number	Diff line number	Diff line
		@@ -354,6 +354,10 @@ c.volcano <- function(
		special_gene <- NULL
		}

		if(length(compare_label) != n_compare){
		compare_label <- selected_compare
		}

		dt <- de.dt(
		.data,
		fold_change,
		@@ -380,10 +384,13 @@ c.volcano <- function(
		][
		,label_special := fcase(gene_name %in% special_gene,gene_name)
		][
		data.table(compare,compare_label),
		on = "compare"
		data.table(compare = selected_compare,compare_label),
		on = .(compare)
		]




		stat <- merge(
		dt[,.(main = .N),.(compare,regulation)],
		dt[,.(sub = .N),.(compare,regulation = sub_regulation)],