1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68
|
library("matrixStats")
rowDiffs_R <- function(x, lag = 1L, differences = 1L, ..., useNames = TRUE) {
ncol2 <- ncol(x) - lag * differences
if (ncol2 <= 0) {
y <- matrix(x[integer(0L)], nrow = nrow(x), ncol = 0L)
# Preserve names attribute
if (useNames && !is.null(rownames(x))) rownames(y) <- rownames(x)
return(y)
}
suppressWarnings({
y <- apply(x, MARGIN = 1L, FUN = diff, lag = lag, differences = differences)
})
y <- t(y)
# Preserve dimnames attribute
dim(y) <- c(nrow(x), ncol2)
if (useNames && !is.null(dimnames(x))) {
colnames <- colnames(x)
if (!is.null(colnames)) {
len <- length(colnames)
colnames <- colnames[(len - ncol2 + 1):len]
}
dimnames(y) <- list(rownames(x), colnames)
}
else dimnames(y) <- NULL
y
}
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Subsetted tests
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
source("utils/validateIndicesFramework.R")
x <- matrix(runif(6 * 6, min = -6, max = 6), nrow = 6, ncol = 6)
storage.mode(x) <- "integer"
# To check dimnames attribute
dimnames <- list(letters[1:6], LETTERS[1:6])
# Test with and without dimnames on x
for (setDimnames in c(TRUE, FALSE)) {
if (setDimnames) dimnames(x) <- dimnames
else dimnames(x) <- NULL
count <- 0L
for (rows in index_cases) {
for (cols in index_cases) {
count <- count + 1L
for (lag in 1:2) {
for (differences in 1:3) {
# Check dimnames attribute
useNames <- c(TRUE, FALSE)
useNames <- useNames[count %% length(useNames) + 1]
validateIndicesTestMatrix(x, rows, cols,
ftest = rowDiffs, fsure = rowDiffs_R,
lag = lag, differences = differences, useNames = useNames)
validateIndicesTestMatrix(x, rows, cols,
ftest = function(x, rows, cols, ..., useNames) {
t(colDiffs(t(x), rows = cols, cols = rows, ..., useNames = useNames))
}, fsure = rowDiffs_R,
lag = lag, differences = differences, useNames = useNames)
}
}
}
}
}
|