R中truncdist的积分函数的“奇怪”行为
托梅克
从数学上讲，以下是不可能的
library(truncdist)
q = function(x, L, R ) dtrunc(x, "exp", rate=0.1, a=L,b=R)
integrate(q, L=2, R=3, lower =0, upper = 27 )
integrate(q, L=2, R=3, lower =0, upper = 29 )
integrate(q, L=2, R=3, lower =27, upper = 29 )
integrate(q, L=2, R=3, lower =0, upper = 30 )
我们发现第一个积分给出一个正数，第二个积分通过添加将自身积分为零的第三个间隔而得出零。这是一个问题integrate还是truncdist？
我们可以使用以下内容查找更多此类问题
z=numeric()
for(i in 1:50){
z[i]=integrate(q, L=2, R=3, lower =0, upper = i)$value
}
我需要怎么做才能找到正确的积分（从0到积分时全为1 i>=3）？
罗兰
来自help("integrate")：
像所有数值积分例程一样，这些函数会在一组有限的点上评估函数。如果函数在几乎所有范围内都近似恒定（尤其是零），则结果和误差估计可能会严重错误。
您找到了一个示例：
curve(q(x, 2, 3), from = -1, to = 30)
您不应该对分布密度函数进行数字积分。使用累积分布函数：
diff(ptrunc(c(0, 29), "exp", rate = 0.1, a = 2, b = 3))
#[1] 1
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
编辑于
相关文章
Related 相关文章
- 1
R中匹配函数的奇怪行为
- 2
R 3.0.0中by（）函数的奇怪行为？
- 3
R中匹配函数的奇怪行为
- 4
R中的积分函数
- 5
JavaScript中奇怪的函数行为
- 6
JavaScript中奇怪的函数行为
- 7
Javascript函数中的奇怪行为
- 8
递归函数中的奇怪/奇怪行为
- 9
R中奇怪的strptime行为
- 10
R中列表的奇怪行为
- 11
R中的积分函数返回错误
- 12
Javascript新Date函数中的奇怪行为
- 13
dplyr中do（）函数的奇怪行为
- 14
虚拟函数中final的奇怪行为
- 15
关于C ++中构造函数的奇怪行为
- 16
递归函数在PHP中的行为很奇怪
- 17
data.table v1.9.5（R）中shift（）函数的奇怪行为
- 18
在R函数中管道传递给return（）时会出现奇怪的行为？
- 19
从R中的stdin读取时的奇怪行为
- 20
R中gtools :: mixedsort包的奇怪行为
- 21
R中的附加/分离行为非常奇怪
- 22
MATCH函数的奇怪行为
- 23
等例函数的奇怪行为
- 24
递归函数的奇怪行为
- 25
.capitalize函数的奇怪行为
- 26
jQuery函数奇怪的行为
- 27
'constexpr'函数奇怪的行为？
- 28
递归函数的奇怪行为
- 29
C ++中的函数指针转换中的奇怪行为
我来说两句