问题描述
我们如何从 R 中的 3 维数组构造块对角矩阵？从矩阵列表（例如，Reduce(magic::adiag,list_of_matrices)）或单个矩阵（例如，magic::adiag(matrix1,matrix2)）开始时，有多种可能性。但是，当我们从数组开始时，我找不到任何东西：
matrices <- array(NA,c(3,3,2))
matrices[,1] <- diag(1,3)
matrices[,2] <- matrix(rnorm(9),3)
是否有任何有效的解决方案来构建相应的 9x9 块矩阵，或者只是转换为列表并使用 magic::adiag 是否更好？后者似乎效率相对较低，尤其是在矩阵数量较多的情况下。
解决方法
我想转换为列表并使用 magic::adiag 是最快的方法。尝试以下几行代码，它很短而且我经常使用：
library(magic)
arr <- array(1:8,c(2,2,3))
do.call("adiag",lapply(seq(dim(arr)[3]),function(x) arr[,x]))
这本质上简化为单行但使用列表。