m <- 2 n <- 3 matlist <- matrix(vector("list"), m, n) for (i in 1:m) for (j in 1:n) matlist[[i,j]] <- matrix(sample(12), 3, 4) matlist[[1,2]] ## access with double[[ subscripts