Deprecated/deprecatedClasses.R

In osenan/cliqueMS: Annotation of Isotopes, Adducts and Fragmentation Adducts for in-Source LC/MS Metabolomics Data

#' @export
#' @title 'anClique-class' for annotating isotopes and adducts
#'
#' @aliases anClique-class
#' @description
#' S3 Class \code{anClique-class} for annotating isotopes and adducts
#' in processed m/z data. Features are first
#' grouped based on a similarity network algorithm and then
#' annotation of isotopes and adducts is performed in each group.
#' @details anClique-class contains the following elements:
#' 
#' 'peaklist'
#' Is a data.frame with m/z, retention time
#' and intensity information for each feature. It also contains
#' adduct and isotope information if annotation has been performed.
#'
#' 'network'
#' Is an igraph undirected network of similarity 
#' used to compute groups of features before annotation.
#' 
#' 'cliques'
#' Is a list that contains the groups of features.
#' Each id corresponds to a row in the peaklist.
#'
#' 'isotopes'
#' Is a data.frame with the column 'feature' for feature id, 
#' column 'charge' for the charge, column 'grade' that starts with 0
#' and it is 1 for the first isotope, 2 for the second and so on and
#' column 'cluster' which labels each group of features that are
#' isotopes.
#'
#' 'cliquesFound' is
#' TRUE if clique groups have been computed,
#' 'isoFound' is
#' TRUE if isotopes have been annotated,
#' 'anFound' is
#' TRUE if annotation of adducts have been computed.
#' @examples
#' mzfile <- system.file("standards.mzXML", package = "cliqueMS")
#' library(xcms)
#' mzraw <- readMSData(files = mzfile, mode = "onDisk")
#' cpw <- CentWaveParam(ppm = 15, peakwidth = c(5,20), snthresh = 10)
#' mzData <- findChromPeaks(object = mzraw, param = cpw)
#' ex.anClique <- createanClique(mzData = mzData)
#' summary(ex.anClique)
#' @seealso \code{\link{createanClique}}
anClique <- structure(list("peaklist" = data.frame(),
    "network" = igraph::empty_graph(),
    "cliques" = list(),
    "cliquesFound" = FALSE,
    "isotopes" = data.frame,
    "isoFound" = FALSE,
    "anFound" = FALSE),
    class = "anClique")


#' @export
#' @title 'createanClique' generic function to create an object
#' of class 'anClique'.
#'
#' @description
#' \code{createanClique} creates an 'anClique' object from processed m/z data.e
#' @param mzData An object with processed m/z data. See methods for
#' valid class types.
#' @return
#' An 'anClique' S3 object with all elements to perform clique grouping,
#' isotope annotation and adduct annotation.
#' @examples
#' mzfile <- system.file("standards.mzXML", package = "cliqueMS")
#' library(xcms)
#' mzraw <- readMSData(files = mzfile, mode = "onDisk")
#' cpw <- CentWaveParam(ppm = 15, peakwidth = c(5,20), snthresh = 10)
#' mzData <- findChromPeaks(object = mzraw, param = cpw)
#' ex.anClique <- createanClique(mzData = mzData)
#' summary(ex.anClique)
#' @seealso \code{\link{anClique-class}}
createanClique <- function(mzData) UseMethod("createanClique")

#' @export
#' @title 'createanClique.xcmsSet' produces an object of class 'anClique'.
#'
#' @description
#' \code{createanClique.xcmsSet} creates an 'anClique' object from 'xcmsSet'
#' processed m/z data.
#' @param mzData A 'xcmsSet' object with processed m/z data.
#' @return
#' An 'anClique' S3 object with all elements to perform clique grouping,
#' isotope annotation and adduct annotation.
#' @details CAMERA package has to be installed to use this method.
#' @examples
#' mzfile <- system.file("standards.mzXML", package = "cliqueMS")
#' msSet <- xcms::xcmsSet(files = mzfile, method = "centWave",
#' ppm = 15, peakwidth = c(5,20), snthresh = 10)
#' ex.anClique <- createanClique(msSet)
#' summary(ex.anClique)
#' @seealso \code{\link{anClique-class}}
createanClique.xcmsSet <- function(mzData) {
    if (!requireNamespace("CAMERA", quietly = TRUE)) {
        stop("Package CAMERA needed for 'xcmsSet' processed data. Please use
'XCMSnExp' objects or install package CAMERA.",
    call. = FALSE)
    }
    if(is(mzData,"xcmsSet") == FALSE) {
        stop("mzData should be of class xcmsSet") }
    peaklist = as.data.frame(mzData@peaks)
    cliques = list()
    isotopes = matrix()
    return(structure(list("peaklist" = peaklist,
    "network" = igraph::empty_graph(),
    "cliques" = cliques,
    "cliquesFound" = FALSE,
    "isotopes" = list(),
    "isoFound" = FALSE,
    "anFound" = FALSE),
    class = "anClique"))
}

#' @export
#' @title 'createanClique.XCMSnExp' produces an object of class 'anClique'.
#'
#' @description
#' \code{createanClique.XCMSnExp} creates an 'anClique' object from
#' 'XCMSnExp' processed m/z data.
#' @param mzData A 'XCMSnExp' object with processed m/z data.
#' @return
#' An 'anClique' S3 object with all elements to perform clique grouping,
#' isotope annotation and adduct annotation.
#' @examples
#' require(xcms)
#' mzfile <- system.file("standards.mzXML", package = "cliqueMS")
#' rawMS <- readMSData(files = mzfile, mode = "onDisk")
#' cpw <- CentWaveParam(ppm = 15, peakwidth = c(5,20), snthresh = 10)
#' mzData <- findChromPeaks(rawMS, cpw)
#' ex.anClique <- createanClique(mzData)
#' summary(ex.anClique)
#' @seealso \code{\link{anClique-class}}
createanClique.XCMSnExp <- function(mzData) {
    if(is(mzData,"XCMSnExp") == FALSE) {
        stop("mzData should be of class XCMSnExp")
    }
    peaklist = as.data.frame(xcms::chromPeaks(mzData))
    cliques = list()
    isotopes = matrix()
    return(structure(list("peaklist" = peaklist,
    "network" = igraph::empty_graph(),
    "cliques" = cliques,
    "cliquesFound" = FALSE,
    "isotopes" = list(),
    "isoFound" = FALSE,
    "anFound" = FALSE),
    class = "anClique"))
}

#' @export
summary.anClique <- function(object, ...)
{
    message(paste("anClique object with",nrow(object$peaklist),
    "features"), sep = " ")
    if(object$cliquesFound) {
        message(paste("Features have been splitted into",
        length(object$cliques), "cliques", sep = " "))
    } else {
        message("No computed clique groups")
    }
    if(object$isoFound) {
        if( sum(is.na(unlist(object$isotopes))) ==
            length(unlist(object$isotopes)) ) {
            message("0 Features are isotopes")
        } else {
        message(paste(nrow(object$isotopes), "Features are isotopes", sep = " "))
        }
    } else {
        message("No isotope annotation")
    }
    if(object$anFound) {
        pos1 = which(!is.na(object$peaklist$an1))
        pos2 = which(!is.na(object$peaklist$an2))
        pos3 = which(!is.na(object$peaklist$an3))
        pos4 = which(!is.na(object$peaklist$an4))
        pos5 = which(!is.na(object$peaklist$an5))
        anFeatures = unique(c(pos1,pos2,pos3,pos4,pos5))
        message(paste(length(anFeatures), "features annotated", sep = " "))
    } else {
        message("No adduct annotation")
    }
}