swh:1:snp:0c004a03453a29b80f921a24433f7e780b9ceb53
Raw File
Tip revision: c1069ff942f364c8bbda771405176d5f6173b3c8 authored by vobencha on 30 October 2018, 15:54:39 UTC
bump x.y.z versions to odd y after creation of RELEASE_3_8 branch
Tip revision: c1069ff
calculatedAnnotation.R
#' Create calculated annotation
#'
#' \code{calculatedAnnotation} adds a column calculated by operation
#'
#' @param es ExpressionSet object.
#'
#' @param operation Name of the operation to perform calculation
#'
#' @param columns List of specified columns' indices (optional),
#' indices start from 0#'
#'
#' @param rows List of specified rows' indices (optional), indices start from 0
#'
#' @param name Name of the new annotation
#'
#' @param isColumns Apply fn to columns
#'
#' @return Nothing. Annotated dataset will be assigned to es in environment
#'
#' @import Biobase
#'

calculatedAnnotation <- function (es, operation, rows = c(),
                                  columns = c(), isColumns = FALSE, name = NULL) {
    rows <- getIndicesVector(rows, nrow(exprs(es)))
    columns <- getIndicesVector(columns, ncol(exprs(es)))

    fn <- tolower(operation)
    if (!isColumns) {
        fData(es)[[name]] <- NA
        fData(es)[[name]][rows] <- apply(exprs(es[rows,columns]), 1, fn)
    } else {
        phenoData(es)[[name]] <- NA
        phenoData(es)[[name]][columns] <- apply(exprs(es[rows,columns]), 2, fn)
    }

    assign("es", es, envir = parent.frame())

}
back to top