[R] vectorizing a matrix computation

From: Christoph Lehmann <christoph.lehmann_at_gmx.ch>
Date: Tue 20 Jul 2004 - 19:31:04 EST

Dear R users

I have a 4-dimensional matrix (actually several 3d (x,y, slices) matrices appended over time (volumes))

say, e.g. I want to z-transform the data (subtract the mean and divide by the std-deviation)

for (slice in 1:slices) {

     for (x in 1:x.dim) {
         for (y in 1:y.dim) {
         t <- as.matrix(my.matrix[x,y,slice,1:volumes])
         for (vol in 1:volumes) {
             my.matrix.transformed[x,y,slice,vol] <- 
(my.matrix[x,y,slice,vol] - mean(t))/sqrt(var(t))


how can I vectorize such a function using, one of the *apply functions?

many thanks



R-help@stat.math.ethz.ch mailing list

https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Jul 20 19:44:35 2004

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 02:37:09 EST