The 8-parameter Fisher–Bingham distribution on the sphere
Tóm tắt
The Fisher–Bingham distribution (
$$\mathrm {FB}_8$$
) is an eight-parameter family of probability density functions (PDF) on the unit sphere that, under certain conditions, reduce to spherical analogues of bivariate normal PDFs. Due to difficulties in its interpretation and estimation, applications have been mainly restricted to subclasses of
$$\mathrm {FB}_8$$
, such as the Kent (
$$\mathrm {FB}_5$$
) or von Mises–Fisher (vMF) distributions. However, these subclasses often do not adequately describe directional data that are not symmetric along great circles. The normalizing constant of
$$\mathrm {FB}_8$$
can be numerically integrated, and recently Kume and Sei showed that it can be computed using an adjusted holonomic gradient method. Both approaches, however, can be computationally expensive. In this paper, I show that the normalization of
$$\mathrm {FB}_8$$
can be expressed as an infinite sum consisting of hypergeometric functions, similar to that of the
$$\mathrm {FB}_5$$
. This allows the normalization to be computed under summation with adequate stopping conditions. I then fit the
$$\mathrm {FB}_8$$
to two datasets using a maximum-likelihood approach and show its improvements over a fit with the more restrictive
$$\mathrm {FB}_5$$
distribution.
Tài liệu tham khảo
Abramowitz M, Stegun IA (1964) Handbook of mathematical functions with formulas, graphs, and mathematical tables, National Bureau of Standards Applied Mathematics Series, vol 55, 10th edn. Department of Commerce, National Bureau of Standards, Gaithersburg
Bingham C, Mardia KV (1978) A small circle distribution on the sphere. Biometrika 65(2):379–389. https://doi.org/10.1093/biomet/65.2.379
Chirkin D (2013) Event reconstruction in IceCube based on direct event re-simulation. In: Proceedings, 33rd international cosmic ray conference (ICRC2013): Rio de Janeiro, Brazil, July 2–9, 2013, p 0581. http://inspirehep.net/record/1413016/files/icrc2013-0581.pdf. Accessed 14 June 2019
Fisher R (1953) Dispersion on a sphere. Proc R Soc Lond Ser A 217:295–305. https://doi.org/10.1098/rspa.1953.0064
Fraenkel E (2014) From radio pulse to elusive particle. PhD thesis, University of Groningen
Gorski KM, Hivon E, Banday AJ, Wandelt BD, Hansen FK, Reinecke M, Bartelmann M (2005) HEALPix: a framework for high-resolution discretization and fast analysis of data distributed on the sphere. Astrophys J 622(2):759–771. https://doi.org/10.1086/427976
Hoffleit D, Warren Jr W (1995) Vizier online data catalog: Bright star catalogue. VizieR Online Data Catalog 5050
Kent JT (1982) The Fisher–Bingham distribution on the sphere. J R Stat Soc Ser B 44(1):71–80
Kim B, Huckemann S, Schulz J, Jung S (2019) Small-sphere distributions for directional data with application to medical imaging. Scand J Stat. https://doi.org/10.1111/sjos.12381
Kume A, Sei T (2018) On the exact maximum likelihood inference of Fisher–Bingham distributions using an adjusted holonomic gradient method. Stat Comput 28(4):835–847. https://doi.org/10.1007/s11222-017-9765-3
Kume A, Wood ATA (2005) Saddlepoint approximations for the Bingham and Fisher–Bingham normalising constants. Biometrika 92(2):465–476. https://doi.org/10.1093/biomet/92.2.465
Mardia KV (1972) Statistics of directional data, probability and mathematical statistics. Academic Press, London, New York, p 13
Mardia KV (1975) Statistics of directional data. J R Stat Soc Ser B 37(3):349–393
Mardia KV, Jupp PE (2000) Directional statistics. Wiley series in probability and statistics. Wiley, Chichester
Nocedal J, Wright SJ (2006) Numerical optimization, 2nd edn. Springer, New York
Oliphant TE (2006) A guide to NumPy, vol 1. Trelgol Publishing, New York
Rivest LP (1984) On the information matrix for symmetric distributions on the hypersphere. Ann Stat 12(3):1085–1089. https://doi.org/10.1214/aos/1176346724
Van Der Walt S, Colbert SC, Varoquaux G (2011) The NumPy array: a structure for efficient numerical computation. Comput Sci Eng 13(2):22
Virtanen P, Gommers R, Oliphant TE, Haberland M, Reddy T, Cournapeau D, Burovski E, Peterson P, Weckesser W, Bright J, van der Walt SJ, Brett M, Wilson J, Jarrod Millman K, Mayorov N, Nelson ARJ, Jones E, Kern R, Larson E, Carey C, Polat I, Feng Y, Moore EW, VanderPlas J, Laxalde D, Perktold J, Cimrman R, Henriksen I, Quintero EA, Harris CR, Archibald AM, Ribeiro AH, Pedregosa F, van Mulbregt P et al (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nature Methods 17:261–272. https://doi.org/10.1038/s41592-019-0686-2
Zonca A, Singer LP, Lenz D, Reinecke M, Rosset C, Hivon E, Gorski KM (2019) healpy: equal area pixelization and spherical harmonics transforms for data on the sphere in Python. J Open Source Softw 4(35):1298. https://doi.org/10.21105/joss.01298