NA / NaN / Inf 调用外部函数 (arg 5)

Question

我有以下问题，我正在生成随机值，当我运行语法时，我收到以下错误“NA / NaN / Inf in a call to an external function (arg 5)”但是当有基础中的文本（至少我读过），在这种情况下它不应该发生。

需要注意的是，Ars 命令指定了变量 x 的范围，即在区间 0

library(ars)

f<-function(x){ 
               125*log(.5+(.25*x))+34*log(x)+38*log(1-x) +72*log(.25)
}


fprima<-function(x){ 
                    (125/(.2+.25*x))+(34/x)-(38/(1-x)) 
                    }


theta<-ars(n=10000,f,fprima,lb=TRUE,ub=TRUE,xlb=0,xub=1)


theta

Answer 1

将message(x) 添加到f 以查看它传递的x 值。在您的原始代码中，您指定了一个起始值x = 0，并且无论指定的边界如何，它都会传递给f。如果您使用默认起始值​​（如代码的当前迭代），则相同。

如果你指定了一个合理的起始值，代码就会运行：

theta<-ars(n=10000,f,fprima,lb=TRUE,ub=TRUE,xlb=0,xub=1, m = 1, x = 0.5)

但是，您随后会收到以下内部错误：

Error in sobroutine sample_...
ifault= 5 

根据文档，这意味着“检测到非凹面”。

您的fprima 函数似乎不正确。根据文档，它应该是d/du log(f(u,...))。基于D(quote(125*log(.5+(.25*x))+34*log(x)+38*log(1-x) +72*log(.25)), "x")，我已经更正了这一点：

fprima<-function(x){ 
  125 * (0.25/(0.5 + (0.25 * x))) + 34 * (1/x) - 38 * (1/(1 - x))
}

这样，代码运行良好。