R/utils.R

In vjcitn/BiocOncoTK: Bioconductor components for general cancer genomics

#' @importFrom SummarizedExperiment colData
#
.load_mock = function(stub) {
 get(load(system.file(sprintf("mocks/%s.rda", stub), package="BiocOncoTK")))
}
#' utilities for mock data (not involving internet access for vignette)
#' @rdname utils
#' @note These functions are provided only for avoiding reliance on internet connectivity for document production.
#' @return a list of DRProfSet instances
#' @examples
#' load_ccleNRAS()
#' dim(load_nrasdf())
#' @export
load_ccleNRAS = function() {
 .load_mock("ccleNRAS")
}
#' @rdname utils
#' @return a data.frame with fields 'Cell_line_primary_name', 'RMA_normalized_expression', 'HGNC_gene_symbol'
#' @export
load_NRAS_AHR = function() {
 .load_mock("NRAS_AHR")
}
#' @rdname utils
#' @return a data.frame
#' @export
load_nrasdf = function() {
 .load_mock("nrasdf")
}

#' bind MSI data to a SummarizedExperiment
#' @param se SummarizedExperiment instance
#' @param useDing logical(1) if TRUE, use MSIsensor outputs from Ding et al. Cell 2018, otherwise use firehose labelings msi-h,msi-l
#' @param onlyHL logical(1) if TRUE, retain only msi-h, msi-l records;
#' ignored if useDing is TRUE
#' @return SummarizedExperiment instance with expanded colData,
#' samples limited to those with microsatellite instability values.
#' The additional variable is called 'msiTest' and is numerical
#' if useDing is TRUE and is character (msi-{h,l,s}) otherwise.
#' @note This function adds the column \code{msiTest} to
#' \code{colData(se)}.  The contents of the column are given by
#' \code{\link{fireMSI}}.  Samples in \code{se}
#' that do not correspond to a row of \code{\link{fireMSI}}
#' are dropped.  If there is already a column named \code{msiTest}
#' in \code{colData(se)}, it is replaced and samples are filtered
#' as described, and a message is given.  If none of the
#' samples in \code{se} have rows in \code{\link{fireMSI}},
#' an error is thrown.  *OF NOTE:* The MSIsensor data from
#' Ding's cell paper (see help(dingMSI) for URL) provides the
#' participant barcode.  The participant barcode is a substring of
#' the sample barcode.  Be sure to filter the input SummarizedExperiment
#' to include only tumor samples, using the substr(colnames(se),14,15)
#' (values "10"..."14" correspond to normal, non-tumor samples.)
#' Additionally, bindMSI will only work if the colnames of the (filtered)
#' SummarizedExperiment have been truncated to the participant barcode,
#' that is, the first 12 characters of the sample barcode.
#' @examples
#' bindMSI
#' @export
bindMSI = function(se, useDing=TRUE, onlyHL=TRUE) {
 inpat = colnames(se)
 if ("msiTest" %in% names(colData(se)))
   message("'msiTest' found among colData columns, will replace.")
 if (useDing) {
    msi = BiocOncoTK::dingMSI
    msivbl = "MSIsensor.score"
    msipat = msi$participant_barcode
    rownames(msi) = msipat
    }
 else {
    msi = BiocOncoTK::fireMSI
    msivbl = "msiTest"
    msipat = msi$patient_barcode
    rownames(msi) = msipat
    }
 ok = intersect(inpat, msipat)
 if (length(ok)==0) stop("no patients in SE found in MSI resource")
 se = se[, ok]
 se$msiTest = msi[colnames(se),msivbl]
 if (!useDing) se$msiTest = as.character(se$msiTest)
 if (!useDing && onlyHL) se = se[, which(se$msiTest %in% c("msi-h", "msi-l"))]
 se
}