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

回答