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)

我收到错误:

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]总是只有一列.)