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

尼古拉基斯

我有一个对象列表，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函数创建这个新矩阵？

阿克伦

我们可以将 to 转换list为array然后使用applywith MARGINasc(1, 2)并按IQR元素应用

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通过IQR循环应用list

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

本文收集自互联网，转载请注明来源。

如有侵权，请联系 [email protected] 删除。