Estimating Skellam distribution and regression parameters in Stata

Abstract

The Skellam distribution is a discrete probability distribution related to the difference between two independent Poisson-distributed random variables. It has been used in a variety of contexts, including sports or supply and demand imbalances in shared transportation. Stata does not support the Skellam distribution or Skellam regression. We present a command, skellamreg, to estimate the parameters of a Skellam distribution and Skellam regression model using Mata’s optimize function.

Keywords

st0748 skellamreg skellamreg postestimation Skellam distribution Skellam regression modified Bessel function of the first kind maximum likelihood optimize

Get full access to this article

View all access options for this article.

References

Abramowitz

Stegun

I. A.

1972. Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables. New York: Dover.

Dobbie

J. M.

1961. Letter to the editor—A doubled-ended queuing problem of Kendall. Operations Research 9: 755–757. https://doi.org/10.1287/opre.9.5.755.

Irwin

J. O.

1937. The frequency distribution of the difference between two independent variates following the same Poisson distribution. Journal of the Royal Statistical Society, A ser., 100: 415–416. https://doi.org/10.2307/2980526.

Jann

2005. Making regression tables from stored estimates. Stata Journal 5: 288–308. https://doi.org/10.1177/1536867X0500500302.

Jann

2007. Making regression tables simplified. Stata Journal 7: 227–244. https://doi.org/10.1177/1536867X0700700207.

Karlis

Ntzoufras

2003. Analysis of sports data using bivariate Poisson models. Journal of the Royal Statistical Society, D ser., 52: 381–393. https://doi.org/10.1111/1467-9884.00366.

Kendall

D. G.

1951. Some problems in the theory of queues. Journal of the Royal Statistical Society, B ser., 13: 151–185. https://doi.org/10.1111/j.2517-6161.1951.tb00080.x.

Lewis

J. W.

Brown

P. E.

Tsagris

2017. skellam: Densities and sampling for the Skellam distribution version 0.2.1. R-Forge. https://rdrr.io/rforge/skellam/.

Liu

Pelechrinis

2021a. xinliupitt/skellam_regression. GitHub. https://github.com/xinliupitt/skellam_regression.

10.

Liu

Pelechrinis

2021b. Excess demand prediction for bike sharing systems. PLOS ONE 16: e0252894. https://doi.org/10.1371/journal.pone.0252894.

11.

Moreau

J.-P

2011. Program to calculate the first kind modified Bessel function of integer order N, for any real X, using the function BESSI(N,X).

12.

Skellam

J. G.

1946. The frequency distribution of the difference between two Poisson variates belonging to different populations. Journal of the Royal Statistical Society, A ser., 109: 296. https://doi.org/10.2307/2981372.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB