0

我在 R 中遇到了“集成”命令的问题,我不知道如何解决它,因此非常感谢如何处理这个问题的任何帮助!

我有一些数据,我使用 nls 命令拟合了一个函数。我想在一定范围内(例如 2-20)集成这个函数,在我的例子中,它说明了水柱中 2-20 米的荧光——可以是任何函数。

我尝试这样做的方式:

depth <- seq(0,100, by=.1)  
#### First making a vector that simulates depth (Have the real once in my own data set).

flu <- 0.216 + 0.140*depth + (-0.01538*(depth^2)) + 0.0004134*(depth^3)  
### Fitted function.

integrand <- function(depth) {flu}

integrate(integrand, lower= 2, upper= 20)

当我这样做时,R 说: 积分错误(整数,下 = 2,上 = 20):函数的评估给出了错误长度的结果

我也尝试过对流感函数进行矢量化,然后再次集成,但这没有帮助。也许它是带有浮点的东西?但我不知道如何处理。

希望你能帮助我 - 感谢您的时间和提前帮助!!!索伦

4

1 回答 1

3
 depth <- seq(0,100, by=.1)  
 flu <- function(depth) {0.216 + 0.140*depth + (-0.01538*(depth^2)) + 0.0004134*(depth^3) }
 integrate(flu, lower= 2, upper= 20)
#7.170026 with absolute error < 8e-14
于 2012-06-04T16:02:00.587 回答