我正在试图整合以下函数:h:
功能:
integralpos <- function(h){
(h)^(m-1)*exp(-x[4]*h-(r[v]-h-x[5]*x[3] + 0.5*(x[6]^2)*x[3])^2/
(2*x[6]^2*x[3]))}
CALL:
integrate(integralneg, -Inf, 0, abs.tol = 1e-08)$val
如果我有所有变量的示例值(exept h),这可以正常工作。
当我尝试将积分函数嵌套在变量变化的循环中时,会出现问题。我最大限度地使用bobyqa。根据x的似然估计。 bobyqa函数改变了x [。]的值,但由于integralpos未被定义为x的函数,因此它不会输入&#34;那里。
有人可以帮助我吗?有没有办法运行与变化变量的集成?还有其他功能吗?
非常感谢
答案 0 :(得分:0)
通过传递变量解决了这个问题:
integrate(integralneg, -Inf, 0, abs.tol = 1e-08,x=x,r=r,m=m)$val
请求帮助