Skellam Distribution

正如你們知道Skellam軟件包已從CRAN中移除（不問爲什麼）。經過一段時間的互聯網研究，我找不到一個skellam pmf函數，所以我自己寫了。Skellam Distribution

2013-02-19 Windys

當有人說「不問爲什麼」時，我很想問「爲什麼？」。爲什麼？ – Spacedman 2013-02-19 12:46:03

好吧，我不是來自CRARN，你知道;）MB是它記錄不完整或不適用 – Windys 2013-02-19 14:14:01

您可以隨時從存檔中獲取舊版本。然後你可以問維護人員Jerry W. Lewis <[email protected]>。檔案上的版本通過R CMD CHECK給我。並安裝。並運行... – Spacedman 2013-02-19 16:51:58

skellam <- function(k,mu1,mu2){ 
    return(exp(-mu1-mu2)*((mu1/mu2)^(k/2))*besselI(2*sqrt(mu1*mu2),k)) 
}

2013-02-19 11:07:26 Windys

請參閱'VGAM'中的'dskellam'。分佈任務視圖對找到這種功能很有幫助。 http://cran.r-project.org/web/views/Distributions.html – 2013-02-19 11:49:55

該包缺少累積質量函數，以防您想要貢獻。 – 2013-02-19 11:50:43

這很好，但考慮到這個功能的簡單性，有什麼意義呢？任何「應該」使用skellam函數的人都可以直接從wikipedia中複製它。 – 2013-02-19 13:08:00

回答