Combine matrices to block diagonal structure
blockdiag(x, ..., pad = 0)
x
...
pad
A <- diag(3)+1 blockdiag(A,A,A,pad=NA)
Klaus K. Holst
Useful links