/usr/share/doc/r-cran-matrixstats/tests/rowSds.R is in r-cran-matrixstats 0.10.3-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | library("matrixStats")
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# With and without some NAs
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
for (addNA in c(FALSE, TRUE)) {
cat("addNA=", addNA, "\n", sep="")
x <- matrix(1:100, nrow=20, ncol=5)
if (addNA) {
x[13:17,c(2,4)] <- as.double(NA)
}
# Row/column ranges
for (na.rm in c(FALSE, TRUE)) {
cat("na.rm=", na.rm, "\n", sep="")
r0 <- apply(x, MARGIN=1, FUN=sd, na.rm=na.rm)
r1 <- rowSds(x, na.rm=na.rm)
r2 <- colSds(t(x), na.rm=na.rm)
stopifnot(all.equal(r1, r2))
stopifnot(all.equal(r1, r0))
stopifnot(all.equal(r2, r0))
}
} # for (addNA ...)
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# All NAs
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
x <- matrix(as.double(NA), nrow=20, ncol=5)
for (na.rm in c(FALSE, TRUE)) {
cat("na.rm=", na.rm, "\n", sep="")
suppressWarnings({
r0 <- apply(x, MARGIN=1, FUN=sd, na.rm=na.rm)
})
r1 <- rowSds(x, na.rm=na.rm)
r2 <- colSds(t(x), na.rm=na.rm)
stopifnot(all.equal(r1, r2))
stopifnot(all.equal(r1, r0))
stopifnot(all.equal(r2, r0))
}
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# A 1x1 matrix
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
x <- matrix(0, nrow=1, ncol=1)
for (na.rm in c(FALSE, TRUE)) {
cat("na.rm=", na.rm, "\n", sep="")
suppressWarnings({
r0 <- apply(x, MARGIN=1, FUN=sd, na.rm=na.rm)
})
r1 <- rowSds(x, na.rm=na.rm)
r2 <- colSds(t(x), na.rm=na.rm)
stopifnot(all.equal(r1, r2))
stopifnot(all.equal(r1, r0))
stopifnot(all.equal(r2, r0))
}
|