我正在尝试在R中运行一些相当深的递归代码,但它一直给我这个错误:
错误:C堆栈使用率太接近限制
我的输出CStack_info()
是:
Cstack_info()
size current direction eval_depth
67108864 8120 1 2
我的机器上有足够的内存,我只是想弄清楚如何增加R的CStack。
编辑:有人要求一个可复制的例子。这是导致问题的一些基本示例代码。多次运行f(1,1),您会收到错误消息。请注意,我已经设置了--max-ppsize = 500000和options(expressions = 500000),因此,如果您未设置这些参数,则可能会收到关于这两件事之一的错误。如您所见,递归在这里可以进行得很深入,我不知道如何使它一致地工作。谢谢。
f <- function(root=1,lambda=1) {
x <- c(0,1);
prob <- c(1/(lambda+1),lambda/(lambda+1));
repeat {
if(root == 0) {
break;
}
else {
child <- sample(x,2,replace=TRUE,prob);
if(child[1] == 0 && child[2] == 0) {
break;
}
if(child[1] == 1) {
child[1] <- f(root=child[1],lambda);
}
if(child[2] == 1 && child[1] == 0) {
child[2] <- f(root=child[2],lambda);
}
}
if(child[1] == 0 && child[2] == 0) {
break;
}
if(child[1] == 1 || child[2] == 1) {
root <- sample(x,1,replace=TRUE,prob);
}
}
return(root)
}
options(expressions = somethinglarge)