The beta prime distribution is the distribution of the odds ratio associated with a random variable with the beta distribution. Since variables with beta distributions are often used to model random probabilities and proportions, the corresponding odds ratios occur naturally as well.
Suppose that \( U \) has the beta distribution with shape parameters \( a, \, b \in (0, \infty) \). Random variable \( X = U \big/ (1 - U) \) has the beta prime distribution with shape parameters \( a \) and \( b \).
The special case \( a = b = 1 \) is known as the standard beta prime distribution and corresponds to the odds ratio of a variable with the standard uniform distribution. Since \( U \) has a continuous distribution on the interval \( (0, 1) \), random variable \( X \) has a continuous distribution on the interval \( (0, \infty) \).
Suppose that \( X \) has the beta prime distribution with shape parameters \( a, \, b \in (0, \infty) \), and as usual, let \( B \) denote the beta function.
\( X \) has probability density function \( f \) given by \[ f(x) = \frac{1}{B(a, b)} \frac{x^{a - 1}}{(1 + x)^{a + b}}, \quad x \in (0, \infty) \]
First, recall that the beta PDF \( g \) with parameters \( a \) and \( b \) is \[ g(u) = u^{a-1} (1 - u)^{b-1}, \quad u \in (0, 1) \] The transformation \( x = u \big/ (1 - u) \) maps \( (0, 1) \) onto \( (0, \infty) \) and is increasing. The inverse transformation is \( u = x \big/ (x + 1) \), and \( 1 - u = 1 \big/ (x + 1) \) and \( du / dx = 1 \big/ (x + 1)^2 \). Thus, by the change of variables formula, \[ f(x) = g(u) \frac{du}{dx} = \frac{1}{B(a, b)} \left(\frac{x}{x+1}\right)^{a-1} \left(\frac{1}{x + 1}\right)^{b-1} \frac{1}{(x + 1)^2} = \frac{1}{B(a, b)} \frac{x^{a-1}}{(x + 1)^{a + b}}, \quad x \in (0, \infty) \]
If \( a \ge 1 \), the probability density function is defined at \( x = 0 \), so in this case, it's customary add this endpoint to the domain. In particular, for the standard beta prime distribution, \[ f(x) = \frac{1}{(1 + x)^2}, \quad x \in [0, \infty) \] Qualitatively, the first order properties of the probability density function \( f \) depend only on \( a \), and in particular on how \( a \) compares in order to 1.
The probability density function \( f \) satisfies the following properties:
These properties follow from standard calculus. The first derivative of \( f \) is \[ f^\prime(x) = \frac{1}{B(a, b)} \frac{x^{a-2}}{(1 + x)^{a + b + 1}} [(a - 1) - x(b + 1)], \quad x \in (0, \infty) \]
Qualitatively, the second order properties of \( f \) also depend only on \( a \), with transitions at \( a = 1 \) and \( a = 2 \).
For \( a \gt 1 \), define \begin{align} x_1 & = \frac{(a - 1)(b + 2) - \sqrt{(a - 1)(b + 2)(a + b)}}{(b + 1)(b + 2)} \\ x_2 & = \frac{(a - 1)(b + 2) + \sqrt{(a - 1)(b + 2)(a + b)}}{(b + 1)(b + 2)} \end{align} The probability density function \( f \) satisfies the following properties:
These results follow from standard calculus. The second derivative of \( f \) is \[ f^{\prime\prime}(x) = \frac{1}{B(a, b)} \frac{x^{a-3}}{(1 + x)^{a + b + 2}}\left[(a - 1)(a - 2) - 2 (a - 1)(b + 2) x + (b + 1)(b + 2)x^2\right], \quad x \in (0, \infty) \]
Open the special distribution simulator and select the beta prime distribution. Vary the parameters and note the shape of the probability density function. For selected values of the parameters, run the simulation 1000 times and compare the empirical density function to the probability density function.
Because of the definition of the beta prime variable, the distribution function of \( X \) has a simple expression in terms of the beta distribution function with the same parameters, which in turn is the regularized incomplete beta function. So let \( G \) denote the distribution function of the beta distribution with parameters \( a, \, b \in (0, \infty) \), and recall that \[ G(x) = \frac{B(x; a, b)}{B(a, b)}, \quad x \in (0, 1) \]
\( X \) has distribution function \( F \) given by \[ F(x) = G\left(\frac{x}{x + 1}\right), \quad x \in [0, \infty) \]
Similarly, the quantile function of \( X \) has a simple expression in terms of the beta quantile function \( G^{-1} \) with the same parameters.
\( X \) has quantile function \( F^{-1} \) given by \[ F^{-1}(p) = \frac{G^{-1}(p)}{1 - G^{-1}(p)}, \quad p \in [0, 1) \]
Open the quantile app and choose the beta prime distribution. Vary the parameters and note the shape of the distribution function. For selected values of the parameters, find the median and the first and third quartiles.
For certain values of the parameters, the distribution and quantile functions have simple, closed form expressions.
If \( a \in (0, \infty) \) and \( b = 1 \) then
For \( a \gt 0 \) and \( b = 1 \), \( G(u) = u^a \) for \( u \in [0, 1] \) and \( G^{-1}(p) = p^{1/a} \) for \( p \in [0, 1] \)
If \( a = 1 \) and \( b \in (0, \infty) \) then
For \( a = 1 \) and \( b \gt 0 \), \( G(u) = 1 - (1 - u)^b \) for \( u \in [0, 1] \) and \( G^{-1}(p) = 1 - (1 - p)^{1/b} \) for \( p \in [0, 1] \).
If \( a = b = \frac{1}{2} \) then
For \( a = b = \frac{1}{2} \), \( G(u) = \frac{2}{\pi} \arcsin\left(\sqrt{u}\right) \) for \( u \in (0, 1) \) and \( G^{-1}(p) = \sin^2\left(\frac{\pi}{2} p\right) \) for \( p \in [0, 1] \)
When \( a = b = \frac{1}{2} \), \( X \) is the odds ratio for a variable with the standard arcsine distribution.
As before, \( X \) denotes a random variable with the beta prime distribution, with parameters \( a, \, b \in (0, \infty) \). The moments of \( X \) have a simple expression in terms of the beta function.
If \( t \in (-a, b) \) then \[ \E\left(X^t\right) = \frac{B(a + t, b - t)}{B(a, b)} \] If \(t \in (-\infty, -a] \cup [b, \infty)\) then \( \E(X^t) = \infty \).
Once again, let \( g \) denote the beta PDF with parameters \( a \) and \( b \). With the transformation \( x = u \big/ (1 - u) \), as in the proof of , we have \( f(x) dx = g(u) du \). Hence \[ \int_0^\infty x^t f(x) dx = \int_0^1 \left(\frac{u}{1 - u}\right)^t g(u) du = \frac{1}{B(a, b)} \int_0^1 u^{a + t - 1} (1 - u)^{b - t - 1} du \] If \( t \le -a \) the improper integral diverges to \( \infty \) at 0. If \( t \ge b \) the improper integral diverges to \( \infty \) at 1. If \( -a \lt t \lt b \) the integral is \( B(a + t, b - t) \) by definition of the beta function.
Of course, we are usually most interested in the integer moments of \( X \). Recall that for \( x \in \R \) and \( n \in \N \), the rising power of \( x \) of order \( n \) is \( x^{[n]} = x (x + 1) \cdots (x + n - 1) \).
Suppose that \( n \in \N \). If \( n \lt b \) Then \[ \E\left(X^n\right) = \prod_{k=1}^n \frac{a + k - 1}{b - k} \] If \( n \ge b \) then \( \E\left(X^n\right) = \infty \).
From , \[ \E(X^n) = \frac{B(a + n, b - n)}{B(a, b)} = \frac{\Gamma(a + n) \Gamma(a - n)}{\Gamma(a + b)} \frac{\Gamma(a + b)}{\Gamma(a) \Gamma(b)} = \frac{\Gamma(a + n)}{\Gamma(a)} \frac{\Gamma(b - n)}{\Gamma(b)} = \frac{a^{[n]}}{(b - n)^{[n]}} \] by a basic property of the gamma function.
As a corollary, we have the mean and variance.
If \( b \gt 1 \) then \[ \E(X) = \frac{a}{b - 1} \]
If \( b \gt 2 \) then \[ \var(X) = \frac{a (a + b - 1)}{(b - 1)^2 (b - 2)} \]
Open the Special Distribution Simulator and select the beta prime distribution. Vary the parameters and note the size and location of the mean\(\pm\)standard deviation bar. For selected values of the parameters, run the simulation 1000 times and compare the empirical mean and standard deviation to the distribution mean and standard deviation.
Finally, the general moment result in leads to the skewness and kurtosis of \( X \).
If \( b \gt 3 \) then \[ \skw(X) = \frac{2 (2 a + b - 1)}{b - 3} \sqrt{\frac{b - 2}{a (a + b - 1)}}\]
In particular, the distibution is positively skewed for all \( a \gt 0 \) and \( b \gt 3 \).
If \( b \gt 4 \) then \[ \kur(X) = \frac{3 a^3 b^2 + 69 a^3 b - 30 a^3 + 6 a^2 b^3 +12 a^2 b^2 -78 a^2 b +60 a^2 + 3 a b^4 + 9 a b^3 - 69 a b^2 + 99 a b - 42 a + 6 b^4 - 30 b^3 + 54 b^2- 42 b + 12}{(a + b - 1)(b - 3)(b - 4)} \]
The most important connection is the one between the beta prime distribution and the beta distribution given in definition . We repeat this for emphasis.
Suppose that \( a, \, b \in (0, \infty) \).
The beta prime family is closed under the reciprocal transformation.
If \( X \) has the beta prime distribution with parameters \( a, \, b \in (0, \infty) \) then \( 1 / X \) has the beta prime distribution with parameters \( b \) and \( a \).
A direct proof using the change of variables formula is possible, of course, but a better proof uses a corresponding property of the beta distribution. By definition, we can take \( X = U \big/ (1 - U) \) where \( U \) has the beta distribution with parameters \( a \) and \( b \). But then \( 1 / X = (1 - U) \big/ U \), and \( 1 - U \) has the beta distribution with parameters \( b \) and \( a \). By another application of the definition, \( 1 / X \) has the beta prime distribution with parameters \( b \) and \( a \).
The beta prime distribution is closely related to the \( F \) distribution by a simple scale transformation.
Connections with the \( F \) distributions.
Let \( f \) denote the PDF of \( X \) and \( g \) the PDF of \( Y \).
The beta prime is the distribution of the ratio of independent variables with standard gamma distributions. (Recall that standard here means that the scale parameter is 1.)
Suppose that \( Y \) and \( Z \) are independent and have standard gamma distributions with shape parameters \( a \in (0, \infty) \) and \( b \in (0, \infty) \), respectively. Then \( X = Y / Z \) has the beta prime distribution with parameters \( a \) and \( b \).
Of course, a direct proof can be constructed, but a better approach is to use the previous result. Thus suppose that \( Y \) and \( Z \) are as stated in the theorem. Then \( 2 Y \) and \( 2 Z \) are independent chi-square variables with \( 2 a \) and \( 2 b \) degrees of freedom, respectively. Hence \[ W = \frac{Y / 2a}{Z / 2b} \] has the \( F \) distribution with \( 2 a \) degrees of freedom in the numerator and \( 2 b \) degrees of freedom in the denominator. By the previous result, \[ X = \frac{2 a}{2 b} W = \frac{Y}{Z} \] has the beta prime distribution with parameters \( a \) and \( b \).
The standard beta prime distribution is the same as the standard log-logistic distribution.
The PDF of the standard beta prime distribution is \( f(x) = 1 \big/ (1 + x)^2 \) for \( x \in [0, \infty) \), which is the same as the PDF of the standard log-logistic distribution.
Finally, the beta prime distribution is a member of the general exponential family of distributions.
Suppose that \( X \) has the beta prime distribution with parameters \( a, \, b \in (0, \infty) \). Then \( X \) has a two-parameter general exponential distribution with natural parameters \( a - 1 \) and \( -(a + b) \) and natural statistics \( \ln(X) \) and \( \ln(1 + X) \).
This follows from the definition of the general exponential family, since the PDF can be written in the form \[ f(x) = \frac{1}{B(a, b)} \exp[(a - 1) \ln(x) - (a + b) \ln(1 + x)], \quad x \in (0, \infty) \]