将附加参数传递给 pracma::integral2

Question

我正在使用pracma::integral2 集成一个有时需要额外参数的功能，但integral2 似乎无法识别我提供的参数。 考虑这两个函数：

fun1 <- function(x, y) exp(-x^2 - y^2)
fun2 <- function(x, y, a) (1 / a) * exp(-x^2 - y^2)

如果我再跑

integral2(fun1,
      xmin = 0,
      xmax = 1,
      ymin = 0,
      ymax = 1)$Q

我按预期得到 0.5577463。但如果我跑：

   integral2(fun2,
          xmin = 0,
          xmax = 1,
          ymin = 0,
          ymax = 1,
          a = 1)$Q

我明白了：

Error in fun(x, y, ...) : argument "a" is missing, with no default 

引用显示如下：

4. fun(x, y, ...)
5. FUN(X, Y)
6. .tensor(xmin, xmax, thetaL, thetaR, phiB, phiT, FUN, phiBvar, phiTvar，矢量化 = 矢量化，奇异 = 奇异）
7. integral2(fun2, xmin = 0, xmax = 1, ymin = 0, ymax = 1, a = 1)

我不知道.tensor是什么，但是看起来它的产品FUN(X,Y)已经失去了...。

这里出了什么问题？

Answer 1

a 与integral2 的abstol 参数相匹配。只需将 a 重命名为 fun2 中的其他名称，如下面的解决方案 1 中所示，或者明确指定 abstol，这样它就不会像下面的解决方案 2 中那样混淆。

library(pracma)

# solution 1
fun2aa <- function(x, y, aa) (1 / aa) * exp(-x^2 - y^2)
integral2(fun2aa,
          xmin = 0,
          xmax = 1,
          ymin = 0,
          ymax = 1,
          aa = 1)$Q
## [1] 0.5577463

# solution 2
fun2 <- function(x, y, a) (1 / a) * exp(-x^2 - y^2) # same as in question
integral2(fun2,
          xmin = 0,
          xmax = 1,
          ymin = 0,
          ymax = 1,
          abstol = 0, # added
          a = 1)$Q
## [1] 0.5577463

请注意，如果integral2 在abstolabstol 之前放置了点.这里integral2a 与integral2 相同，只是点点点放在ymax 之后。

# solution 3

fun2 <- function(x, y, a) (1 / a) * exp(-x^2 - y^2) # same as in question

integral2a <- function (fun, xmin, xmax, ymin, ymax, ..., 
  sector = FALSE, reltol = 1e-06, abstol = 0, maxlist = 5000, 
  singular = FALSE, vectorized = TRUE) {}
body(integral2a) <- body(integral2)
environment(integral2a) <- environment(integral2)

integral2a(fun2,
          xmin = 0,
          xmax = 1,
          ymin = 0,
          ymax = 1,
          a = 1)$Q
## [1] 0.5577463