R中包含乘积项的函数的积分

Posted 2023-03-24

【中文标题】R中包含乘积项的函数的积分

【英文标题】：Integration of a function in R containing product term

【发布时间】：2020-07-15 11:05:54

【问题描述】：

我需要在 R 的形式中集成一个函数

f<-function(x,a)prod((2-log(a))^x)
a<-c(0.5,0.7,0.9,1.2)
integrate(f,0,1,a=a)

运行上述代码后，出现以下错误。

Error in integrate(f, 0, 1, a = a) :
evaluation of function gave a result of wrong length
In addition: Warning message:
In (2 - log(a))^x :
longer object length is not a multiple of shorter object length

请提供问题的解决方案。

【参考方案1】：

问题是integrate 为向量x 评估f(x,a)。

你可以这样做：

a <- c(0.5,0.7,0.9,1.2)

g <- function(x, a)
  sapply(x, function(xi)
    prod((2-log(a))^xi)
  )


integrate(g, 0, 1, a = a)

【讨论】：

非常感谢@Stephen。它解决了我的问题。