inst/benchmarks/scale.r

In wrathematics/pcapack: Principal Components Analysis Tools and Utilities

library(pcapack)
library(rbenchmark)

set.seed(1234)

m <- 2e4
n <- 500
mean <- 0
sd <- 1
x <- matrix(rnorm(m*n, sd=sd, mean=mean), m, n)

reps <- 10

cols <- c("test", "replications", "elapsed", "relative")

print("center=TRUE, scale=FALSE")
benchmark(R=A <- scale(x, TRUE, FALSE), pcapack=B <- scale2(x, TRUE, FALSE), replications=reps, columns=cols) 
all.equal(A, B, check.attributes=FALSE)

print("center=FALSE, scale=TRUE")
benchmark(R=A <- scale(x, FALSE, TRUE), pcapack=B <- scale2(x, FALSE, TRUE), replications=reps, columns=cols)
all.equal(A, B, check.attributes=FALSE)

print("center=TRUE, scale=TRUE")
benchmark(R=A <- scale(x, TRUE, TRUE), pcapack=B <- scale2(x, TRUE, TRUE), replications=reps, columns=cols)
all.equal(A, B, check.attributes=FALSE)