tests/rowTabulates.R

In matrixStats: Methods that apply to rows and columns of a matrix

library("matrixStats")

N <- 6L
K <- 5L
J <- 5L
data <- matrix(1:J, nrow=N, ncol=K)


modes <- c("integer", "raw")
for (mode in modes) {
  cat(sprintf("Mode: %s...\n", mode))

  x <- data
  storage.mode(x) <- mode
  print(x)

  y <- rowTabulates(x)
  print(y)
  stopifnot(identical(dim(y), c(N,J)))
  
  y <- colTabulates(x)
  print(y)
  stopifnot(identical(dim(y), c(K,J)))
  
  # Count only certain values
  y <- rowTabulates(x, values=1:3)
  print(y)
  stopifnot(identical(dim(y), c(N,3L)))
  
  y <- colTabulates(x, values=1:3)
  print(y)
  stopifnot(identical(dim(y), c(K,3L)))

  cat(sprintf("Mode: %s...done\n", mode))
} # for (mode ...)