我正在嘗試創建一個函式,它與 R 中的 pbinom 相同,即二項式分布的密度。
我的代碼如下:
MyPBin <- function(X,n,p) {
choose(n,X)*p^X*(1-p)^(n-X)
}
這應該與 R 中 pbinom 中使用的函式匹配: p(x) = choose(n, x) p^x (1-p)^(nx)
但是,我沒有得到正確的結果。我想我已經找到了原因——如果我對所有單獨的選擇函式(n=1、n=2、n=3 等)求和,我得到的結果與 pbinom 相同。
所以我的問題是:當我運行它時,我如何總結該函式使用的所有單獨的選擇函式?
uj5u.com熱心網友回復:
使用總和。這些都給出了相同的結果。(請注意,MyPBin 給出的結果與 dbinom 相同。)
sum(MyPBin(0:3, 10, .3))
## [1] 0.6496107
sum(dbinom(0:3, 10, .3))
## [1] 0.6496107
pbinom(3, 10, .3)
## [1] 0.6496107
轉載請註明出處,本文鏈接:https://www.uj5u.com/net/401153.html
