如何在 R 中定義概率質量 function

Question

我想在處理另一個計算時將概率質量 function 定義（實現）到 R 環境中。 由於沒有太多的 R 經驗，我不知道如何將帶有這些參數的方程式翻譯成代碼（已經嘗試到處搜索）。 感謝幫助。 [![此方程給出的概率質量 function][1]][1]

Answer 1

這可以通過Beta 分布表示

你的表達式 = BetaPDF(r, x) t (1-t) / x

P <- function(x, theta) {
    shape1 <- 2
    shape2 <- x
    dbeta(theta, shape1, shape2)*theta*(1.0-theta)/x
}

請檢查數學，另請參見此處。

更新

或者表示為更高次冪的Beta PDF

你的表達式 = BetaPDF(r+1, x+1) r /((r+x+1)*(r+x))

Answer 2

它已經是 R 中的內置 pmf - 當您安裝 R 時它就在那里。 如果您已經定義了r和theta ，只需調用dnbinom

dnbinom(x, size=r, prob=theta)

查看幫助?dnbinom

當有代碼已經調試過並且每天由大量知識淵博的用戶（例如許多研究統計學家使用它）進行測試時，不要浪費您的時間調試自己的代碼——他們當然應該發現除了非常罕見的問題之外的任何東西. 它也可能比你寫的在數值上更穩定，更有效地實現。