如何获取多维R数组中单个单元格的DIMNAME？

Question

我有一个大型多维数组（约1900万个元素），其中包含多个不同属性的联合概率。

数组非常稀疏，我只对概率非零的单元格感兴趣。

但是，在筛选数组中的非零元素时，我无法检索筛选值的维度名称（对应于各种属性值）。

以下是一个玩具示例：

array_dim <- c(2,5,5,4)

array_fill <- runif(prod(array_dim))

array_dimnames <- list(
                    c('strawberry', 'blackberry'),
                    c('cranberry', 'banana', 'pineapple', 'apple', 'tangerine'), 
                    c('orange', 'blueberry', 'kiwi', 'grapes', 'guava'),
                    c('plum', 'fig', 'grapefruit', 'lemon')
                    )


fruits <- array(array_fill, dim=array_dim, dimnames=array_dimnames)

我可以获得与某个标准匹配的单元格的索引值（此处，

> which(fruits %in% fruits[fruits>0.9], arr.ind = TRUE)
 [1]   8  23  25  32  33  35  37  76  77  85  90 101 117 121 123 135 154 197

但我无法使用上述索引值并找出它们是什么水果组合，因为在查找特定单元格值时，DIMNAME会被删除：

> fruits[8]
[1] 0.9590207
> fruits[8, drop=FALSE]
[1] 0.9590207
> dimnames(fruits[8])
NULL
> names(fruits[8])
NULL

我尝试将数组转换为data.frame，并使用drop=FALSE参数：

> fruits.df <- as.data.frame(fruits)
> 
> fruits.df[1,2,drop=FALSE]

           banana.orange.plum
strawberry          0.4003854

但是添加条件过滤器失败，因为fruits.df[fruits.df

作为最后的手段，我可以构造 array_index -

我还研究了listaries包。

提前感谢




Answer 1

我也没有找到直接从数组中获取DIMNAME的简单方法。转换数据结构的简单方法是使用dplyr中的作为_tbl_cube，并将其转换为data.frame（或data.table）以查看dimnames：

set.seed(3)
array_dim <- c(2,5,5,4)
array_fill <- runif(prod(array_dim))
array_dimnames <- list(
    dim1=c('strawberry', 'blackberry'),
    dim2=c('cranberry', 'banana', 'pineapple', 'apple', 'tangerine'), 
    dim3=c('orange', 'blueberry', 'kiwi', 'grapes', 'guava'),
    dim4=c('plum', 'fig', 'grapefruit', 'lemon')
)
fruits <- array(array_fill, dim=array_dim, dimnames=array_dimnames)
which(fruits %in% fruits[fruits>0.9], arr.ind = TRUE)
#>  [1]  28  54  56  73  74  85  90 115 161 198
fruits[198]
#> [1] 0.9065314

library(dplyr)
arr.cube <- as.tbl_cube(fruits)
tail(as.data.frame(arr.cube))
#>           dim1      dim2  dim3  dim4    fruits
#> 195 strawberry pineapple guava lemon 0.7057146
#> 196 blackberry pineapple guava lemon 0.3907374
#> 197 strawberry     apple guava lemon 0.8242374
#> 198 blackberry     apple guava lemon 0.9065314
#> 199 strawberry tangerine guava lemon 0.4171170
#> 200 blackberry tangerine guava lemon 0.2791320

在这个例子中，水果[198]的名字是黑莓、苹果、番石榴和柠檬。