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

R/RNA-seq.R

+3 −3

Original line number	Diff line number	Diff line
		@@ -870,13 +870,13 @@ c.subDE <- function(
		compare,
		breaf = TRUE
		){
		necessary_cols <- c("gene_id","gene_name","type","adj.P.Val","log10P",compare)
		if(isTRUE(breaf)){
		common_cols <-c("gene_id","gene_name","type","adj.P.Val")
		final_col <- necessary_cols
		}else{
		common_cols <- c("gene_id","gene_name","type","chr","start","end","strand","AveExpr","F","P.Value","adj.P.Val","log10P")
		common_cols <- c(necessary_cols,"chr","start","end","strand","AveExpr","F","P.Value")
		}

		final_col <- c(common_cols,compare)
		.data[,..final_col] \|>
		setnames(compare,"log2FC") %>%
		.[]