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)

我收到错误：

if（nrow（M）== 2）{：参数长度为零
时出错

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）{：参数长度为零时出错

1 个答案: