2014-09-02 122 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