我只是在学习如何使用R.我练了一些东西的统计,如正态分布,泊松等逼近结果
当我试图计算概率,答案是一个数字非常接近为零(0),程序显示为结果0,所以我看不到完整的答案,我需要完整的答案。总有一个概率,即使是一个小的!
我的问题是:我可以关闭自我逼近或可以使用哪个代码来获得完整答案?
实施例:
1-pbinom(q =10, size = 10,prob = 0.8)
结果:
0
我只是在学习如何使用R.我练了一些东西的统计,如正态分布,泊松等逼近结果
当我试图计算概率,答案是一个数字非常接近为零(0),程序显示为结果0,所以我看不到完整的答案,我需要完整的答案。总有一个概率,即使是一个小的!
我的问题是:我可以关闭自我逼近或可以使用哪个代码来获得完整答案?
实施例:
1-pbinom(q =10, size = 10,prob = 0.8)
结果:
0
的pbinom
函数给出的累积密度函数。这是一个值小于或等于特定值的概率。所以,一个像10二项式分布离散分布绘制
pbinom(10, 10, .8)
# [1] 1
告诉你,有一个100%的变化,你会看到10次或更少的成功。
也许你会想的概率密度函数(或概率密度函数,因为这是一个离散分布)dbinom
dbinom(10, 10, .8)
# [1] 0.1073742
意味着有一个约11%的机会,所有的画会成功。这也是事实,
sum(dbinom(0:10, 10, .8))
# [1] 1
,通过获得0的概率之和正好是1
这些情况下,你所得到的确切答案左右。 R根据options(digits=)
设置在控制台中进行舍入值,但这不是这里发生的情况。
pbinom
是二项分布的分布函数,它是离散的,因此可以精确为1(如你的例子)。你可能一直在想连续分布像正常或伽玛分布。在这种情况下,舍入可能会导致结果被截断,例如
> 1 - pnorm(10, 0, 1)
[1] 0
然而,p[dist]
功能有旨在解决这一问题的争论lower.tail=FALSE
:
> pnorm(10, 0, 1, lower.tail=FALSE)
[1] 7.619853e-24