基于R中列表对象的矩阵生成

Question

我在 R 中有一个对象列表，如下所示：

set.seed(1234)
data <- matrix(rnorm(3*4,mean=0,sd=1), 3, 4) 
results <- lapply(1:ncol(data), function(i) outer(data[, i], data[, i]))

这将产生以下矩阵

[[1]]
           [,1]        [,2]       [,3]
[1,]  1.4570077 -0.33487534 -1.3089918
[2,] -0.3348753  0.07696698  0.3008557
[3,] -1.3089918  0.30085569  1.1760127

[[2]]
          [,1]       [,2]       [,3]
[1,]  5.502298 -1.0065968 -1.1870541
[2,] -1.006597  0.1841480  0.2171611
[3,] -1.187054  0.2171611  0.2560926

[[3]]
          [,1]      [,2]      [,3]
[1,] 0.3303260 0.3141712 0.3244131
[2,] 0.3141712 0.2988064 0.3085474
[3,] 0.3244131 0.3085474 0.3186061

[[4]]
          [,1]      [,2]      [,3]
[1,] 0.7921673 0.4247196 0.8886017
[2,] 0.4247196 0.2277129 0.4764227
[3,] 0.8886017 0.4764227 0.9967755

如何创建包含相应元素之间的四分位间距 (IQR) 的单个矩阵？这意味着新矩阵将是一个 B=3X3 矩阵，例如 B(1,1)=IQR(c(1.4570077,5.502298,0.3303260,0.7921673))=1.791623、B(1,2)=IQR(c(-0.33487534, -1.0065968, 0.3141712, 0.4247196))= 0.844614。如何使用 apply 函数创建这个新矩阵？

Answer 1

我们可以将 list 转换为 array 然后使用 apply 和 MARGIN 作为 c(1, 2) 并应用 IQR elementwise

apply(simplify2array(results), c(1, 2), IQR)

-输出

      [,1]       [,2]       [,3]
[1,] 1.791623 0.84461397 1.68299882
[2,] 0.844614 0.08813351 0.07058422
[3,] 1.682999 0.07058422 0.73860710

或者另一种选择是 transpose list，然后通过遍历 list

来应用 IQR

library(purrr)
`dim<-`(map_dbl(transpose(results), IQR), dim(results[[1]]))
         [,1]       [,2]       [,3]
[1,] 1.791623 0.84461397 1.68299882
[2,] 0.844614 0.08813351 0.07058422
[3,] 1.682999 0.07058422 0.73860710

基于R中列表对象的矩阵生成

Matrix generation based on list's objects in R

r

list

matrix

apply

matrix-multiplication