递归函数返回错误信息

Question

我正在尝试创建一个递归函数来进行矩阵乘法n 次。

我的代码如下：

R <- function(P, n){

  R(P, n-1) %*% P

}

我希望当这个函数被 n = 3 调用来执行时

(P %*% P) %*% P.

举个例子：

> P
     [,1] [,2] [,3]
[1,]  0.6  0.1  0.3
[2,]  0.2  0.7  0.1
[3,]  0.3  0.3  0.4

但是，当我调用该函数时，我收到一条错误消息。

错误：求值嵌套太深：无限递归/选项（表达式=）？
总结时出错：评估嵌套太深：无限递归/选项（表达式=）？

您能解释一下为什么该功能不起作用以及应该如何修改代码吗？

Answer 1

每个递归函数都应该有一个阈值，在该点递归调用将停止。您定义的函数R 没有这种情况，因为触发的错误也提到了

求值嵌套太深：无限递归

所以，

R <- function(P, n){
    if (n==1) return(P)
    R(P, n-1) %*% P 
}

应该可以解决您的问题。鉴于你的例子：

r1 <- (P %*% P) %*% P
r2 <- R(P,3)
all(r1==r2)
#[1] TRUE