Last chance! 50% off unlimited learning
Sale ends in
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.
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
# 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