R：应用于 nXm 矩阵的每一列，返回 n*n*m 数组答案

【问题标题】：R: Apply over every column of nXm matrix, return n*n*m arrayR：应用于 nXm 矩阵的每一列，返回 n*n*m 数组
【发布时间】：2018-10-15 17:45:50
【问题描述】：

在矩阵/nXm 数组的每一列上应用返回 n*n 矩阵的函数并返回 nXnXm 数组的适当函数是什么？我可以让它与列表一起使用：

dm <- structure(c(0.205292439079469, 0.222642769676055, 0.222786614980172, 
              0.210917067409095, 0.123955673451174, 0.10591537361648, 0.0887665882561804, 
              0.0725263089185602, 0.174088645938512, 0.171936521036433, 0.159619467190219, 
              0.143492610130578, 0.23303504847129, 0.269389294469438, 0.279103519495381, 
              0.240824583373471, 0.192637398624569, 0.201551950015116, 0.195748345926477, 
              0.181448686311511), .Dim = 4:5)
sapply(1:ncol(dm), function(x){dm[, x] %*% t(dm[, x])^-1}, simplify = F)

谢谢！

【问题讨论】：

标签： r apply

【解决方案1】：

使用simplify = "array" 而不是simplify = F。

sapply(1:ncol(dm), function(x){dm[, x] %*% t(dm[, x])^-1}, simplify = "array")

【讨论】：