modified RNA-seq.R in R (41bd55ff) · Commits · Chaos / chaos.tools

R/RNA-seq.R

+2 −6

Original line number	Diff line number	Diff line
		@@ -86,7 +86,7 @@ rnaseq123 <- function(

		de.dt <- function(
		.data,
		fold_change = 2,
		fold_change = 1.5,
		p_value = 0.05,
		hyperbola = FALSE
		){
		@@ -104,11 +104,7 @@ de.dt <- function(
		as.data.table(
		.data
		)[
		,hyperbola_log10P := fifelse(
		log2FC < 0,
		-1/(log2FC + log2(fold_change)) - log10(p_value),
		1/(log2FC - log2(fold_change)) - log10(p_value)
		)
		,hyperbola_log10P := 1/((abs(log2FC)) - log2(fold_change)) - log10(p_value)
		][
		,regulation := fcase(
		log2FC > log2(fold_change) & log10P > hyperbola_log10P, "up",