时间戳

首页 > 解决方案 > 使用矩阵索引列矩阵时如何理解错误“subscipts out of bound”

r - 使用矩阵索引列矩阵时如何理解错误“subscipts out of bound”

问题描述

考虑以下代码示例：

目标是用 mypoints 中的相应值替换 mymatrix 中的索引。

# example 1--------------------------------------------- 
mymatrix <- c(3,4,1,2, 2,1,4,3, 4,3,2,1)
dim(mymatrix) <- c(4,3)
mypoints <- matrix(c(10, 20, 30, 40), nrow(mymatrix), 1)
# loop
mysummary <- mypoints[mymatrix[,1]]
for (i in 2:ncol(mymatrix) ) {
  mysummary <- cbind(mysummary, mypoints[mymatrix[,i]])
}
# ------------------------------------------------------

结果是：

     mysummary      
[1,]        30 20 40
[2,]        40 10 30
[3,]        10 40 20
[4,]        20 30 10

可以用更简洁的方式替换循环：

> mysummary      <- mypoints[mymatrix]
> dim(mysummary) <- dim(mymatrix)

结果：

     [,1] [,2] [,3]
[1,]   30   20   40
[2,]   40   10   30
[3,]   10   40   20
[4,]   20   30   10

但是，这在以下较小的示例中不起作用：

# example 3---------------------------------------------
> mymatrix <- c(2,1, 2,1)
> dim(mymatrix) <- c(2,2)
> mypoints <- matrix(c(10, 20), nrow(mymatrix), 1)
> mysummary      <- mypoints[mymatrix]
Error in mypoints[mymatrix] : subscript out of bounds

我不知道哪个下标超出范围以及在哪个数组中。以及如何防止这个错误。

标签： rmatrixvectorizationsubscript

解决方案

如果索引矩阵mymatrix有两列，则 R 假定您想使用第一列作为行索引，第二列作为列索引。mymatrix所以它失败了，因为你的第二个例子中的第一行是2,2，它超出了mypoints.

在您的第一个示例中，mymatrix有超过 2 列，因此 R 将其视为向量，这是在这种情况下唯一有意义的方法。

推荐阅读