powered by
This is an implementation of the algorithm given by Owen via repeated integration by parts.
OwensQOwen(nu, t, delta, a=0, b)
degree of Owen<U+2019>s Q
parameter t
parameter delta
lower integration limit. Only a=0 implemented, other values give an error.
a=0
upper integration limit
Numeric value of Owen<U+2019>s Q function.
Owen DB. A special case of a bivariate non-central t-distribution. Biometrika. 1965;52(3/4):437--46. 10.2307/2333696
OwensQ, OwensT
OwensQ
OwensT
# NOT RUN { # comparison of the results of both implementations # both should give [1] 0.0731726 OwensQ(2, 2.92, 4.2135, 0, 2.0407) OwensQOwen(2, 2.92, 4.2135, 0, 2.0407) # }
Run the code above in your browser using DataLab