r - Skellam Distribution -
as u guys know skellam bundle removed cran (dont inquire why). after time of net research, couldnt find skellam pmf function, wrote myself.
skellam <- function(k,mu1,mu2){ return(exp(-mu1-mu2)*((mu1/mu2)^(k/2))*besseli(2*sqrt(mu1*mu2),k)) }
r statistics
No comments:
Post a Comment