Revision 66357a4654bfe45c5fe4a9565bdf67d7b5a6f05c authored by clevermx on 18 January 2022, 11:09:39 UTC, committed by Space on 18 January 2022, 11:09:39 UTC
Add meta-files generation for new counts processing.
Add support for 8+ versions of ARCHs4.

Co-authored-by: Alexey Sergushichev <alserg@itmo.ru>
Merge-request: PHANTASUS-MR-21
Merged-by: Maksim Kleverov <klevermx@gmail.com>
1 parent d230978
Raw File
createES.R
#' Create ExpressionSet.
#'
#' \code{createES} function produces an ExpressionSet object from given data,
#'     and exports it to global scope.
#'
#' @param data Gene expression matrix.
#'
#' @param pData Matrix with phenotypical data.
#'
#' @param varLabels Names of phenoData columns.
#'
#' @param fData Matrix with feature data.
#'
#' @param fvarLabels Names of featureData columns.
#'
#' @param eData List with experimentData
#'
#' @return produced ExpressionSet object
#'
#' @import Biobase
#' @importFrom methods new
#'
#' @examples
#' \dontrun{
#' data <- matrix(1:15, 5, 3)
#' pData <- c("A", "B", "C")
#' varLabels <- "cat"
#' fData <- c("p", "r", "s", "t", "u")
#' fvarLabels <- "id"
#' eData <- list(name="", lab="", contact="", title="", url="", other=list(), pubMedIds="")
#' createES(data, pData, varLabels, fData, fvarLabels, eData)
#' }
#'
createES <- function(data, pData, varLabels, fData, fvarLabels, eData) {
    phenoData <- AnnotatedDataFrame(data.frame(pData, stringsAsFactors = FALSE))
    varLabels(phenoData) <- varLabels

    featureData <- AnnotatedDataFrame(data.frame(fData, stringsAsFactors = FALSE))
    varLabels(featureData) <- fvarLabels

    ed <- new ("MIAME",
              name=eData$name,
              lab=eData$lab,
              title=eData$title,
              contact=eData$contact,
              pubMedIds=eData$pubMedIds,
              url=eData$url,
              other=eData$other)

    es <- ExpressionSet(assayData = data,
                        phenoData = phenoData,
                        featureData = featureData,
                        experimentData = ed)
    assign("es", es, envir = parent.frame())
    es
}
back to top