使用 R 中的索引向量的子集矩阵

Question

我在尝试理解代码时遇到了一些问题。我有一个代码，它应该加载到矩阵中，并在下面的列中调用顶部的 n 功能。

TopFeatures <- function(matrix, ncols = 1, nfeatures = 20){
   mtx <- matrix[ , ncols, drop = FALSE]
   data <- mtx[order(mtx, decreasing = TRUE), , drop = FALSE]
   features <- head(x = rownames(x = data), n = nfeatures
   return(features)
} 

但是，每当我尝试超过 1 列时，我都会收到如下错误。

> TopFeatures(so, ncols = 1:3, nfeatures = 10)
Error in mtx[order(mtx, decreasing = TRUE), , drop = FALSE] : 
  subscript out of bounds

我相信错误来自这一行：

data <- mtx[order(mtx, decreasing = TRUE), , drop = FALSE]

但是，我真的不明白这条线应该做什么。

我的理解是，这里的目的是对矩阵中的所有单元格进行排序，然后按降序排序，然后挑出最靠前的n 特征。函数order 对您的矩阵进行排序并返回一个索引向量，该向量等于矩阵中单元格的数量，但您不能使用索引向量进行子集化，因为该数字将超过行数。我想这就是为什么我有这个错误。但是，我不知道如何解决这个问题。

有人可以指导我如何解决这个问题吗？

非常感谢。

Answer 1

您应该order 仅基于一列。因此，当您拥有一列时，它可以工作，但当您拥有多于一列时，它就不起作用。

假设第一列是你想order你的数据，你可以试试：

TopFeatures <- function(matrix, ncols = 1, nfeatures = 20){
   mtx <- matrix[ , ncols, drop = FALSE]
   data <- mtx[order(mtx[, 1], decreasing = TRUE), , drop = FALSE]
   features <- head(x = rownames(data), n = nfeatures)
   return(features)
} 

---

在这一行

mtx[order(mtx, decreasing = TRUE), , drop = FALSE]

您正在尝试根据mtx 矩阵中的值排列行。

假设您有一个单列矩阵。

mtx <- matrix(1:10)
data <- mtx[order(mtx, decreasing = TRUE), , drop = FALSE] #gives
data

#      [,1]
# [1,]   10
# [2,]    9
# [3,]    8
# [4,]    7
# [5,]    6
# [6,]    5
# [7,]    4
# [8,]    3
# [9,]    2
#[10,]    1

而使用 2 列矩阵，您会得到：

mtx <- matrix(1:10, ncol = 2)
data <- mtx[order(mtx, decreasing = TRUE), , drop = FALSE]

mtx[order(mtx, reduction = TRUE), , drop = FALSE] 中的错误： 下标越界

即因为

order(mtx, decreasing = TRUE) #returns
#[1] 10  9  8  7  6  5  4  3  2  1

但您的数据中只有 5 行，因此出现错误。

如果你 order 只在第一列，你会得到

order(mtx[, 1], decreasing = TRUE)
#[1] 5 4 3 2 1

这不会给你一个子集错误。