Revision 44656149198007189ebc600ddfbd19539e86be1f authored by Vladislav Kamenev on 28 January 2019, 10:05:53 UTC, committed by Vladislav Kamenev on 28 January 2019, 10:05:53 UTC
1 parent 155b528
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())
}
Computing file changes ...