如何解决R中数组的块对角矩阵
我们如何从 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]))
这本质上简化为单行但使用列表。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。