R/removeSig.R

#' removeSig
#' @export

removeSig <-
function(peak.mat,out_path){ # remove top ranked signals
	fs_out <- read.table(file=out_path, sep=",",
                     na.strings = "", stringsAsFactors=F, head=T)
  fs_out <- lapply(as.list(fs_out), function(x) x[!is.na(x)])
  fs_out_pos <- fs_out[[1]]
  idx  <- !(colnames(peak.mat$p) %in% fs_out_pos)
  peak.mat$p <- peak.mat$p[,idx,drop=F]
  fs_out_neg <- fs_out[[2]]
  idx  <- !(colnames(peak.mat$n) %in% fs_out_neg)
  peak.mat$n <- peak.mat$n[,idx,drop=F]
	return(peak.mat)
}
jasenfinch/OrbiFIEmisc documentation built on May 18, 2019, 4:49 p.m.