Let\'s make a replicable example:
This is my initial matrix
d<- matrix(1:80,,5) d [,1] [,2] [,3] [,4] [,5] [1,] 1 17 33 49 65 [2,
We can do this directly
do.call(rbind,lapply(1:ncol(d),function(t) matrix(d[,t],ncol = 4)))