BetaPrime Distribution

[Naba7]

b1 = BetaPrime()
rand(rng, b1) === quantile(b1, rand(rng))

this returns false, while this says otherwise.

[Elrod]

Make sure not to call rand separate times, but I’m assuming you didn’t when actually testing.

More importantly here is the definition being called:

julia> @which rand(MersenneTwister(1), BetaPrime())
rand(rng::AbstractRNG, d::BetaPrime) in Distributions at /home/chriselrod/.julia/packages/Distributions/jEqbk/src/univariate/continuous/betaprime.jl:113

You can also use @less to look at it, or @edit to open the file in a text editor.

[Naba7]

thanks! This helps.