powered by
This is an implementation of the algorithm given by Owen via repeated integration by parts.
OwensQOwen(nu, t, delta, a=0, b)
Numeric value of Owen’s Q function.
degree of Owen’s Q
parameter t
parameter delta
lower integration limit. Only a=0 implemented, other values give an error.
a=0
upper integration limit
D. Labes
Owen DB. A special case of a bivariate non-central t-distribution. Biometrika. 1965;52(3/4):437--46. tools:::Rd_expr_doi("10.2307/2333696")
OwensQ, OwensT
OwensQ
OwensT
# 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