简单的整合

Question

 constant <- function(x){1}

 integrate(constant,lower=1,upper=10)

这可能不仅仅是荒谬简单,而是为什么这不起作用的任何理由？

我应该用不同的方式编写常量函数吗？

Answer 1

您可以使用该Vectorize函数将非向量化函数转换为integrate需要的函数类型:

> constant <- function(x){1}
> Vconst <- Vectorize(constant)
>  integrate(Vconst,lower=1,upper=10)
9 with absolute error < 1e-13