R Under development (unstable) (2023-11-06 r85483 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library("matrixStats")
> 
> ## Always allow testing of the 'center' argument (as long as it's not defunct)
> options(matrixStats.center.onUse = "ignore")
> 
> rowMads_R <- function(x, na.rm = FALSE, ..., useNames = NA) {
+   suppressWarnings({
+     res <- apply(x, MARGIN = 1L, FUN = mad, na.rm = na.rm)
+   })
+   if (is.na(useNames) || !useNames) names(res) <- NULL
+   res
+ }
> 
> colMads_R <- function(x, na.rm = FALSE, ..., useNames = NA) {
+   suppressWarnings({
+     res <- apply(x, MARGIN = 2L, FUN = mad, na.rm = na.rm)
+   })
+   if (is.na(useNames) || !useNames) names(res) <- NULL
+   res
+ }
> 
> rowMads_center <- function(x, rows = NULL, cols = NULL, na.rm = FALSE, ..., useNames = NA) {
+   center <- rowMedians(x, cols = cols, na.rm = na.rm, useNames = FALSE)
+   rowMads(x, rows = rows, cols = cols, center = center, na.rm = na.rm, useNames = useNames)
+ }
> 
> colMads_center <- function(x, rows = NULL, cols = NULL, na.rm = FALSE, ..., useNames = NA) {
+   center <- colMedians(x, rows = rows, na.rm = na.rm, useNames = FALSE)
+   colMads(x, rows = rows, cols = cols, center = center, na.rm = na.rm, useNames = useNames)
+ }
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # A 3x3 matrix (no ties)
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> for (mode in c("integer", "double")) {
+   x <- matrix(c(1, 2, 3, 2, 3, 4, 3, 4, 5) + 0.1, nrow = 3, ncol = 3)
+   cat("mode: ", mode, "\n", sep = "")
+   storage.mode(x) <- mode
+   str(x)
+   
+   # To check name attribute
+   dimnames <- list(letters[1:3], LETTERS[1:3])
+   
+   # Test with and without dimnames on x
+   for (setDimnames in c(TRUE, FALSE)) {
+     if (setDimnames) dimnames(x) <- dimnames
+     else dimnames(x) <- NULL
+     tx <- t(x)
+     # Check names attribute
+     for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("rowMads():\n")
+       center <- rowMedians(x, na.rm = TRUE, useNames = useNames)
+       r0 <- rowMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- rowMads(x, na.rm = TRUE, useNames = useNames)
+       r1b <- rowMads_center(x, na.rm = TRUE, useNames = useNames)
+       r1c <- rowMads(x, center = center, na.rm = TRUE, useNames = useNames)
+       r2 <- colMads(tx, na.rm = TRUE, useNames = useNames)
+       r2b <- colMads_center(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r1b, r1))
+       stopifnot(all.equal(r1c, r1))
+       stopifnot(all.equal(r2, r0))
+       stopifnot(all.equal(r2b, r2))
+       
+       cat("colMads():\n")
+       center <- colMedians(x, na.rm = TRUE, useNames = useNames)
+       r0 <- colMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- colMads(x, na.rm = TRUE, useNames = useNames)
+       r1b <- colMads_center(x, na.rm = TRUE, useNames = useNames)
+       r1c <- colMads(x, center = center, na.rm = TRUE, useNames = useNames)
+       r2 <- rowMads(tx, na.rm = TRUE, useNames = useNames)
+       r2b <- rowMads_center(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r1b, r1))
+       stopifnot(all.equal(r1c, r1))
+       stopifnot(all.equal(r2, r0))
+       stopifnot(all.equal(r2b, r2))
+     }
+   }
+ }
mode: integer
 int [1:3, 1:3] 1 2 3 2 3 4 3 4 5
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
mode: double
 num [1:3, 1:3] 1.1 2.1 3.1 2.1 3.1 4.1 3.1 4.1 5.1
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Ties: a 4x4 matrix
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> for (mode in c("integer", "double")) {
+   x <- matrix(c(1:4, 2:5, 3:6, 4:7) + 0.1, nrow = 4, ncol = 4)
+   cat("mode: ", mode, "\n", sep = "")
+   storage.mode(x) <- mode
+   str(x)
+   tx <- t(x)
+   
+   # To check name attribute
+   dimnames <- list(letters[1:4], LETTERS[1:4])
+   
+   # Test with and without dimnames on x
+   for (setDimnames in c(TRUE, FALSE)) {
+     if (setDimnames) dimnames(x) <- dimnames
+     else dimnames(x) <- NULL
+     tx <- t(x)
+     # Check names attribute
+     for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("rowMads():\n")
+       r0 <- rowMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- rowMads(x, na.rm = TRUE, useNames = useNames)
+       r2 <- colMads(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+       
+       cat("colMads():\n")
+       r0 <- colMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- colMads(x, na.rm = TRUE, useNames = useNames)
+       r2 <- rowMads(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+     }
+   }
+ }
mode: integer
 int [1:4, 1:4] 1 2 3 4 2 3 4 5 3 4 ...
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
mode: double
 num [1:4, 1:4] 1.1 2.1 3.1 4.1 2.1 3.1 4.1 5.1 3.1 4.1 ...
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
> tx <- NULL
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # No ties: a 3x3 matrix with an NA value
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> for (mode in c("integer", "double")) {
+   x <- matrix(c(1, 2, 3, 2, 3, 4, 3, 4, 5) + 0.1, nrow = 3, ncol = 3)
+   x[2, 2] <- NA_real_
+ 
+   cat("mode: ", mode, "\n", sep = "")
+   storage.mode(x) <- mode
+   str(x)
+   tx <- t(x)
+   
+   # To check name attribute
+   dimnames <- list(letters[1:3], LETTERS[1:3])
+   
+   # Test with and without dimnames on x
+   for (setDimnames in c(TRUE, FALSE)) {
+     if (setDimnames) dimnames(x) <- dimnames
+     else dimnames(x) <- NULL
+     tx <- t(x)
+     # Check names attribute
+     for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("rowMads():\n")
+       r0 <- rowMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- rowMads(x, na.rm = TRUE, useNames = useNames)
+       r2 <- colMads(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+       
+       cat("colMads():\n")
+       r0 <- colMads_R(x, na.rm = TRUE, useNames = useNames)
+       r1 <- colMads(x, na.rm = TRUE, useNames = useNames)
+       r2 <- rowMads(tx, na.rm = TRUE, useNames = useNames)
+       stopifnot(all.equal(r1, r2))
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+     }
+   }
+ }
mode: integer
 int [1:3, 1:3] 1 2 3 2 NA 4 3 4 5
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
mode: double
 num [1:3, 1:3] 1.1 2.1 3.1 2.1 NA 4.1 3.1 4.1 5.1
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
rowMads():
colMads():
> tx <- NULL
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # With and without some NAs
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> for (add_na in c(FALSE, TRUE)) {
+   cat("add_na = ", add_na, "\n", sep = "")
+ 
+   x <- matrix(1:50, nrow = 10L, ncol = 5L)
+   if (add_na) {
+     x[3:7, c(2, 4)] <- NA_real_
+   }
+   
+   # To check name attribute
+   dimnames <- list(letters[1:10], LETTERS[1:5])
+   
+   # Test with and without dimnames on x
+   for (setDimnames in c(TRUE, FALSE)) {
+     if (setDimnames) dimnames(x) <- dimnames
+     else dimnames(x) <- NULL
+     tx <- t(x)
+     # Row/column ranges
+     for (na.rm in c(FALSE, TRUE)) {
+       # Check names attribute
+       for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("na.rm = ", na.rm, "\n", sep = "")
+         cat("rowMads():\n")
+         center <- rowMedians(x, na.rm = na.rm, useNames = useNames)
+         r0 <- rowMads_R(x, na.rm = na.rm, useNames = useNames)
+         r1 <- rowMads(x, na.rm = na.rm, useNames = useNames)
+         r2 <- rowMads(x, center = center, na.rm = na.rm, useNames = useNames)
+         r3 <- colMads(tx, na.rm = na.rm, useNames = useNames)
+         stopifnot(all.equal(r1, r0))
+         stopifnot(all.equal(r2, r0))
+         stopifnot(all.equal(r2, r1))
+         stopifnot(all.equal(r3, r0))
+         stopifnot(all.equal(r3, r1))
+         stopifnot(all.equal(r3, r2))
+         
+         cat("colMads():\n")
+         center <- colMedians(x, na.rm = na.rm, useNames = useNames)
+         r0 <- colMads_R(x, na.rm = na.rm, useNames = useNames)
+         r1 <- colMads(x, na.rm = na.rm, useNames = useNames)
+         r2 <- colMads(x, center = center, na.rm = na.rm, useNames = useNames)
+         r3 <- rowMads(tx, na.rm = na.rm, useNames = useNames)
+         stopifnot(all.equal(r1, r0))
+         stopifnot(all.equal(r2, r0))
+         stopifnot(all.equal(r2, r1))
+         stopifnot(all.equal(r3, r0))
+         stopifnot(all.equal(r3, r1))
+         stopifnot(all.equal(r3, r2))
+       }
+     }
+   }
+   tx <- NULL
+ } # for (add_na ...)
add_na = FALSE
na.rm = FALSE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
add_na = TRUE
na.rm = FALSE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = FALSE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
na.rm = TRUE
rowMads():
colMads():
> 
> 
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # All NAs
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> x <- matrix(NA_real_, nrow = 10L, ncol = 5L)
> # Test with and without dimnames on x
> for (setDimnames in c(TRUE, FALSE)) {
+   if (setDimnames) dimnames(x) <- dimnames
+   else dimnames(x) <- NULL
+   tx <- t(x)
+   # Row/column ranges
+   for (na.rm in c(FALSE, TRUE)) {
+     # Check names attribute
+     for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("na.rm = ", na.rm, "\n", sep = "")
+       r0 <- rowMads_R(x, na.rm = na.rm, useNames = useNames)
+       if (na.rm) r0[is.na(r0)] <- NaN
+       r1 <- rowMads(x, na.rm = na.rm, useNames = useNames)
+       r2 <- colMads(tx, na.rm = na.rm, useNames = useNames)
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+       stopifnot(all.equal(r1, r2))
+     }
+   }
+ }
na.rm = FALSE
na.rm = FALSE
na.rm = TRUE
na.rm = TRUE
na.rm = FALSE
na.rm = FALSE
na.rm = TRUE
na.rm = TRUE
> tx <- NULL
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # A 1x1 matrix
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> x <- matrix(0, nrow = 1L, ncol = 1L)
> dimnames <- list("a", "A")
> # Test with and without dimnames on x
> for (setDimnames in c(TRUE, FALSE)) {
+   if (setDimnames) dimnames(x) <- dimnames
+   else dimnames(x) <- NULL
+   tx <- t(x)
+   # Row/column ranges
+   for (na.rm in c(FALSE, TRUE)) {
+     # Check names attribute
+     for (useNames in c(if (!matrixStats:::isUseNamesNADefunct()) NA, TRUE, FALSE)) {
+       cat("na.rm = ", na.rm, "\n", sep = "")
+       r0 <- rowMads_R(x, na.rm = na.rm, useNames = useNames)
+       r1 <- rowMads(x, na.rm = na.rm, useNames = useNames)
+       r2 <- colMads(tx, na.rm = na.rm, useNames = useNames)
+       stopifnot(all.equal(r1, r0))
+       stopifnot(all.equal(r2, r0))
+       stopifnot(all.equal(r1, r2))
+     }
+   }
+ }
na.rm = FALSE
na.rm = FALSE
na.rm = TRUE
na.rm = TRUE
na.rm = FALSE
na.rm = FALSE
na.rm = TRUE
na.rm = TRUE
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # A 0x0 matrix
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> x <- matrix(double(0), nrow = 0, ncol = 0)
> tx <- t(x)
> for (na.rm in c(FALSE, TRUE)) {
+   cat("na.rm = ", na.rm, "\n", sep = "")
+   r0 <- rowMads_R(x, na.rm = na.rm)
+   r1 <- rowMads(x, na.rm = na.rm)
+   r2 <- colMads(tx, na.rm = na.rm)
+   stopifnot(all.equal(r1, r2))
+   stopifnot(all.equal(r1, r0))
+   stopifnot(all.equal(r2, r0))
+ }
na.rm = FALSE
na.rm = TRUE
> 
> proc.time()
   user  system elapsed 
   0.31    0.07    0.39