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))
}