# NOT RUN {
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) # -> 1.776357e-15 now gives 4.651188e-31 !!!
qbeta (0.5078, .01, 5, ncp=0) # also gives 4.651188e-31
# }
Run the code above in your browser using DataCamp Workspace