R/remove_QC_higher_peaks.R

In lucanard/GH_OMIC: The package deconvolutes Q-exactive data in matrix-like data using Yann Guitton scripts and then cleans the dataset eliminating redundant peaks

remove_QC_higher_peaks <- function (clean_peaklist, QC = QCss) {
  target <- grep("mzML.1", colnames(clean_peaklist))
  samples <- clean_peaklist[,target]
  max.QC <- apply(samples[,QC], 1, max)
  max.sam <- apply(samples[,-QC], 1, max)
  clean_peaklist1 <- clean_peaklist[which(max.QC <= max.sam),]
  return(clean_peaklist1)
}