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

Question

应用在矩阵/ 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)

谢谢！

Answer 1

使用simplify = "array"代替simplify = F。

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