R/make_RPKM.R

In doseR: doseR

#' @title make_RPKM Make RPKM.
#' @description make_RPKM populates RPKM slot of SummarizedExperiment
#' S4 object.
#' @usage make_RPKM(se)
#' @param se A SummarizedExperiment object.
#' @return RPKM populated object
#' @examples
#' data(hmel.se)
#' SummarizedExperiment::assays(se)$rpkm <- make_RPKM(se)
#' @author AJ Vaestermark, JR Walters.
#' @references The "doseR" package, 2018 (in press).

make_RPKM <- function(se) {
return(
sweep(sweep(assays(se)$counts, 2, as.vector(colData(se)$Libsizes/1000000),
FUN = '/'), 1, rowData(se)$seglens/1000, FUN = '/')
)
}