The simplest way to create an array is to use the array() function from R and then specify the elements and the dimension as shown in the following code snippet:

        arrayA <- array(1:16, dim=c(2,2,4))

The preceding function will produce an array of four matrices of order 2 x 2. By default, it will create the columns of the matrices and place them next to each other. The output of an array looks like this:

        > arrayA
        , , 1

             [,1] [,2]
        [1,]    1    3
        [2,]    2    4

        , , 2

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

        , , 3

             [,1] [,2]
        [1,]    9   11
        [2,]   10   12

        , , 4

             [,1] [,2]
        [1,]   13   15
        [2,]   14   16