2014-09-02 133 views
1

我只是在学习如何使用R.我练了一些东西的统计,如正态分布,泊松等逼近结果

当我试图计算概率,答案是一个数字非常接近为零(0),程序显示为结果0,所以我看不到完整的答案,我需要完整的答案。总有一个概率,即使是一个小的!

我的问题是:我可以关闭自我逼近或可以使用哪个代码来获得完整答案?

实施例:

1-pbinom(q =10, size = 10,prob = 0.8) 

结果:

0 

回答

3

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=)设置在控制台中进行舍入值,但这不是这里发生的情况。

0

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