qbeta.R(0.6, 2, 3) # 0.4445
qbeta.R(0.6, 2, 3) - qbeta(0.6, 2,3) # almost 0
qbetaRV <- Vectorize(qbeta.R, "alpha") # now can use
curve(qbetaRV(x, 1.5, 2.5))
curve(qbeta (x, 1.5, 2.5), add=TRUE, lwd = 3, col = adjustcolor("red", 1/2))
## an example of disagreement (and doubt, as borderline, close to underflow):
qbeta.R(0.5078, .01, 5) # -> 2.77558e-15 # but
qbeta (0.5078, .01, 5) # now gives 4.651188e-31 -- correctly!
qbeta (0.5078, .01, 5, ncp=0)# ditto
## which is because qbeta() now works in log-x scale here:
curve(pbeta(x, .01, 5), 1e-40, 1, n=10001, log="x", xaxt="n")
sfsmisc::eaxis(1); abline(h=.5078, lty=3); abline(v=4.651188e-31,col=2)
Run the code above in your browser using DataLab