R - if (nrow(M) == 2) { 错误：参数长度为零

Question

在以下代码中：

det<-function(M)
{
    if (nrow(M) == 2)
    {
        return(M[1,1]*M[2,2]-M[1,2]*M[2,1])
    }
    else
    {
      d<-0
      for (i in 1:ncol(M))
      {
        d<-d+(-1)^(i+1)*M[1][i]*det(M[-1,i])
      }
      return(d)
   }
}

执行我确定的功能

 x<-matrix(c(1:9),nrow=3)

执行函数时

 det(x)

我收到错误：

Error in if (nrow(M) == 2) { : argument is of length zero

Answer 1

您正在使用递归函数，问题出在以下行中：

d<-d+(-1)^(i+1)*M[1][i]*det(M[-1,i])

此处您使用 M[-1, i] 对矩阵进行子集化，R 将矩阵转换为数值向量（R 在进行子集化时始终转换为最简单的数据类型）。您可以设置参数 drop=FALSE 来避免这种情况。有关详细信息，请参阅 ?"[" 和 ?drop。

x <- matrix(1:9, nrow=3)
x[-1, 1]
# [1] 2 3
x[-1, 1, drop=FALSE]
#      [,1]
# [1,]    2
# [2,]    3

（请注意，即使您修复了此问题，您的函数也会失败。在您的 if 语句中，您访问了第二列，但 M[-1, i] 将始终只有一列。）

R - if (nrow(M) == 2) { 错误：参数长度为零

R - Error in if (nrow(M) == 2) { : argument is of length zero

arguments

if-statement

r