The Owen’s T function is
presented in four new ways, one of them as a series similar to the Euler’s
arctangent series divided by 2π, which is its majorant series. All
possibilities enable numerically stable and fast convergent computation of the bivariate normal integral with simple
recursion. When tested? computation on a random sample of one million
parameter triplets with
uniformly distributed components and using double precision arithmetic, the
maximum absolute error was 3.45×10-16.
In additional testing, focusing on cases with correlation coefficients close to
one in absolute value, when the computation may be very sensitive to small
rounding errors, the accuracy was retained.
In rare potentially critical cases, a simple adjustment to the computation
procedure was performed—one potentially critical computation was
replaced with two equivalent non-critical ones. All new series are suitable for
vector and high-precision computation, assuming they are supplemented with appropriate efficient and accurate
computation of the arctangent and standard normal cumulative
distribution functions. They are implemented by the R package Phi2rho
References
[1]
Abramowitz, M. and Stegun, I.A. (1972) Handbook of Mathematical Functions. 9th Edition, Dover Publications, Inc., New York.
[2]
Gupta, S.S. (1963) Bibliography on the Multivariate Normal Integrals and Related Topics. The Annals of Mathematical Statistics, 34, 829-838. https://doi.org/10.1214/aoms/1177704005
[3]
Donnelly, T.G. (1973) Algorithm 462: Bivariate Normal Distribution. Communications of the ACM, 16, 638. https://doi.org/10.1145/362375.362414
[4]
Owen, D.B. (1956) Tables for Computing Bivariate Normal Probabilities. The Annals of Mathematical Statistics, 27, 1075-1090. https://doi.org/10.1214/aoms/1177728074
[5]
Divgi, D.R. (1979) Calculation of Univariate and Bivariate Normal Probability Functions. The Annals of Statistics, 7, 903-910. https://doi.org/10.1214/aos/1176344739
[6]
Drezner, Z. and Wesolowsky, G.O. (1990) On the Computation of the Bivariate Normal Integral. Journal of Statistical Computation and Simulation, 35, 101-107. https://doi.org/10.1080/00949659008811236
[7]
Patefield, M. and Tandy, D. (2000) Fast and Accurate Calculation of Owen’s T Function. Journal of Statistical Software, 5, 1-25. https://doi.org/10.18637/jss.v005.i05
[8]
Genz, A. (2004) Numerical Computation of Rectangular Bivariate and Trivariate Normal and t Probabilities. Statistics and Computing, 14, 251-260. https://doi.org/10.1023/B:STCO.0000035304.20635.31
[9]
Meyer, C. (2013) Recursive Numerical Evaluation of the Cumulative Bivariate Normal Distribution. Journal of Statistical Software, 52, 1-14. https://doi.org/10.18637/jss.v052.i10
[10]
Fayed, H.A. and Atiya, A.F. (2014) A Novel Series Expansion for the Multivariate Normal Probability Integrals Based on Fourier Series. Mathematics of Computation, 83, 2385-2402. https://doi.org/10.1090/S0025-5718-2014-02844-5
[11]
Terza, J.V. and Welland, U. (1991) A Comparison of Bivariate Normal Algorithms. Journal of Statistical Computation and Simulation, 39, 115-127. https://doi.org/10.1080/00949659108811343
[12]
Gai, J. (2001) A Computational Study of the Bivariate Normal Probability Function. Master’s Thesis, Queen’s University, Kingston.
[13]
Ağca, Ş.A. and Chance, D.M. (2003) Speed and Accuracy Comparison of Bivariate Normal Distribution Approximations for Option Pricing. The Journal of Computational Finance, 6, 61-96. https://doi.org/10.21314/JCF.2003.102
[14]
Owen, D.B. (1980) A Table of Normal Integrals. Communications in Statistics—Simulation and Computation, 9, 389-419. https://doi.org/10.1080/03610918008812164
[15]
Yee, T.W. (2023) VGAM: Vector Generalized Linear and Additive Models, R Package Version 1.1-8. https://CRAN.R-project.org/package=VGAM
[16]
Laurent, S. (2023) OwenQ: Owen Q-Function, R Package Version 1.0.7. https://CRAN.R-project.org/package=OwenQ
[17]
Genz, A. and Bretz, F. (2009) Computation of Multivariate Normal and t Probabilities. Lecture Notes in Statistics, Springer-Verlag, Heidelberg. https://doi.org/10.1007/978-3-642-01689-9
[18]
Maechler, M. (2023) Rmpfr: R MPFR—Multiple Precision Floating-Point Reliable, R Package Version 0.9-3. https://CRAN.R-project.org/package=Rmpfr
[19]
Kendall, M.G. (1941) Proof of Relations Connected with the Tetrachoric Series and Its Generalization. Biometrika, 32, 196-198. https://doi.org/10.1093/biomet/32.2.196
[20]
Vasicek, O.A. (1998) A Series Expansion for the Bivariate Normal Integral. The Journal of Computational Finance, 1, 5-10. https://doi.org/10.21314/JCF.1998.015
[21]
Plackett, R.L. (1954) A Reduction Formula for Normal Multivariate Integrals. Biometrika, 41, 351-360. https://doi.org/10.1093/biomet/41.3-4.351
[22]
Gradshteyn, I.S. and Ryzhik, I.M. (2015) Table of Integrals, Series, and Products. 8th Edition, Academic Press, Amsterdam.
[23]
Alzer, H. (1997) On Some Inequalities for the Incomplete Gamma Function. Mathematics of Computation, 66, 771-778. https://doi.org/10.1090/S0025-5718-97-00814-4
[24]
Gautschi, W. (1998) The Incomplete Gamma Functions since Tricomi. In: Tricomi’s Ideas and Contemporary Applied Mathematics, Atti Convegni Lincei, Vol. 147, Accademia Nazionale dei Lincei, Rome, 203-237.
Alzer, H. (1998) Inequalities for the Chi Square Distribution Function. Journal of Mathematical Analysis and Applications, 223, 151-157. https://doi.org/10.1006/jmaa.1998.5965
[27]
Olver, F.W.J., Lozier, D.W., Boisvert, R.F. and Clark C.W. (2010) NIST Handbook of Mathematical Functions. National Institute of Standards and Technology & Cambridge University Press, Cambridge.
[28]
Marsaglia, G. (2004) Evaluating the Normal Distribution. Journal of Statistical Software, 11, 1-11. https://doi.org/10.18637/jss.v011.i04
[29]
R Core Team (2023) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna. https://www.R-project.org/
[30]
Komelj, J. (2023) Phi2rho: Owen’s T Function and Bivariate Normal Integral, R Package Version 1.0. https://CRAN.R-project.org/package=Phi2rho
[31]
Robitzsch, A. (2020) pbv: Probabilities for Bivariate Normal Distribution, R Package Version 0.4-22. https://CRAN.R-project.org/package=pbv
[32]
Fayed, H.A. and Atiya, A.F. (2014) An Evaluation of the Integral of the Product of the Error Function and the Normal Probability Density with Application to the Bivariate Normal Integral. Mathematics of Computation, 83, 235-250. https://doi.org/10.1090/S0025-5718-2013-02720-2
[33]
Fousse, L., Hanrot, G., Lefèvre, V., Pélissier, P. and Zimmermann, P. (2007) MPFR: A Multiple-Precision Binary Floating-Point Library with Correct Rounding. ACM Transactions on Mathematical Software, 33, 13-es. https://doi.org/10.1145/1236463.1236468