Displaying 1 result from an estimated 1 matches for "pgbeta2".
Did you mean:
pbeta2
2005 Dec 11
1
Quantile function for the generalized beta distribution of the 2nd kind
I have succeded in defining the cdf of the generalized
beta of the second kind, eg.
pgbeta2 <- function(quint,b,a,p1,p2) {
integrate(function(x)
{exp(log(a)+(a*p1-1)*log(x)-(a*p1)*log(b)-log(beta(p1,p2))-(p1+p2)*log(1+(x/b)^a))},0,quint)$value
}
but I'm facing problems with the quantile function. I
tried something like
qgbeta2 <- function(proba,b,a,p1,p2) {
optimize(functio...