r - 斯凯拉姆分布

Question

正如你们所知，Skellam 包已从 CRAN 中删除（不要问为什么）。经过一段时间的互联网研究，我找不到 skellam pmf 函数，所以我自己编写了它。

Answer 1

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