BetaPrime Distribution

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

this returns false, while this says otherwise.

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.

1 Like

thanks! This helps.

1 Like