%0 Journal Article
%T The Bivariate Normal Integral via Owen¡¯s <i>T</i> Function as a Modified Euler¡¯s Arctangent Series
%A Janez Komelj
%J American Journal of Computational Mathematics
%P 476-504
%@ 2161-1211
%D 2023
%I Scientific Research Publishing
%R 10.4236/ajcm.2023.134026
%X 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