How to bind columns in high dimensional arrays following a specific order?

I have two high dimensional arrays, both with 1000 rows, 3 columns, 10 slices and 4 groups of slices.

I want to bind them but matching the columns between them. Example below:

#dumb data
array1 <- array(1:4, dim = c(2,1,2,2),dimnames = list(NULL,"Ca",1:2,1:2)); array1
array2 <- array(1:4, dim = c(2,1,2,2),dimnames = list(NULL,"Cb",1:2,1:2)); array2

#desired result 
cbind(array1[,,1,1],array2[,,1,1],array1[,,2,1],array2[,,2,1],
array1[,,1,2],array2[,,1,2],array1[,,2,2],array2[,,2,2])

the outcome is:

     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]    1    5    3    7    1    5    3    7
[2,]    2    6    4    8    2    6    4    8

Although column names do not show up, I have ordered manually the columns to be: First column of the first array, first column of the second array, second column of the first array, second column of the second array, and so on.

Is there a other way to do this than manually using cbind?

Thanks in advance!

1 answer

  • answered 2017-08-16 19:25 manotheshark

    Here is one way to accomplish this for the example, but might require some modification for the actual data set:

    m1 <- cbind(array1, array2)
    l1 <- lapply(seq_len(nrow(m1) / 2), function(x) m1[seq(2 * x - 1, 2 * x), ])
    m2 <- do.call(cbind, l1)
    
    > m2
         array1 array2 array1 array2 array1 array2 array1 array2
    [1,]      1      5      3      7      1      5      3      7
    [2,]      2      6      4      8      2      6      4      8