### ============================================================================
### rowWeightedMads
###
### ----------------------------------------------------------------------------
### Non-exported methods
###
.DelayedMatrix_block_rowWeightedMads <- function(x, w = NULL, rows = NULL,
cols = NULL, na.rm = FALSE,
constant = 1.4826,
center = NULL, ..., useNames = TRUE) {
# Check input type
stopifnot(is(x, "DelayedMatrix"))
DelayedArray:::.get_ans_type(x, must.be.numeric = FALSE)
# Check and subset 'w' (must be either NULL, or a numeric vector of
# length 1 or 'ncol(x)')
if (!is.null(w)) {
stopifnot(is.numeric(w))
if (length(w) != 1L) {
stopifnot(length(w) == ncol(x))
if (!is.null(cols))
w <- w[cols]
}
}
# Check, normalize, and subset 'center'
center <- normarg_center_and_subset(center, nrow(x), "nrow(x)", rows)
# Subset 'x'
x <- ..subset(x, rows, cols)
# Compute result
val <- rowblock_APPLY(x = x,
FUN = .rowWeightedMads_internal,
w = w,
na.rm = na.rm,
constant = constant,
center = center,
...,
useNames = useNames)
if (length(val) == 0L) {
return(numeric(nrow(x)))
}
unlist(val, recursive = FALSE, use.names = useNames)
}
#' @importFrom DelayedArray currentViewport makeNindexFromArrayViewport
.rowWeightedMads_internal <- function(x, center, ..., useNames = TRUE) {
if (!is.null(center) && length(center) != 1L) {
block.env <- parent.frame(2)
vp <- currentViewport(block.env)
subset <- makeNindexFromArrayViewport(vp)[[1]]
if (!is.null(subset)) {
center <- center[as.integer(subset)]
}
}
rowWeightedMads(x, center = center, ..., useNames = useNames)
}
### ----------------------------------------------------------------------------
### Exported methods
###
# ------------------------------------------------------------------------------
# General method
#
#' @inherit MatrixGenerics::rowWeightedMads
#' @importMethodsFrom DelayedArray seed
#' @rdname colWeightedMads
#' @export
#' @examples
#'
#' rowWeightedMads(dm_matrix, w = 3:1)
setMethod("rowWeightedMads", "DelayedMatrix",
function(x, w = NULL, rows = NULL, cols = NULL, na.rm = FALSE,
constant = 1.4826, center = NULL,
force_block_processing = FALSE, ..., useNames = TRUE) {
.smart_seed_dispatcher(x, generic = MatrixGenerics::rowWeightedMads,
blockfun = .DelayedMatrix_block_rowWeightedMads,
force_block_processing = force_block_processing,
w = w,
rows = rows,
cols = cols,
na.rm = na.rm,
constant = constant,
center = center,
...,
useNames = useNames)
}
)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.