0

Let $f(x)$ be the probability density function (pdf) of the standard beta distribution on $(0,1)$. And let $f_d(x)$ be the pdf of the generalized beta distribution on $(0,d)$. I know that,

$$f_d(x) = d \cdot f(\frac{x}{d})$$

The cumulative distribution function (cdf) of the generalized beta distribution is,

$$F_d(x) = \int^x_{u=0} f_d(u) \, du.$$

I am tempted to replace $f_d(u)$ by $\left( d \cdot f(\frac{x}{d}) \right)$ in the integral to get,

$$F_d(x) = \int^x_{u=0} df(\frac{u}{d}) \, du = d \cdot \int^x_{u=0} f(\frac{u}{d}) \, du = d \cdot F(\frac{x}{d}).$$

Where $F$ is the cdf of the standard beta distribution. This is wrong because taking $x = d$ gives me $F_d(d) = d$ since $F(1) = 1$.

Question: Where have I gone wrong? Why I can't I replace $f_d(u)$ by $\left( d \cdot f(\frac{x}{d}) \right)$ in the integral?

Legendre
  • 2,845

1 Answers1

2

OP wrote:


Let $f(x)$ be the probability density function (pdf) of the standard beta distribution on $(0,1)$. And let $f_d(x)$ be the pdf of the generalized beta distribution on $(0,d)$. I know that,

$$f_d(x) = d \cdot f(\frac{x}{d})$$


Scaled Beta

If $X$~$Beta(a,b)$ (a standard Beta on (0,1)), then the pdf of the scaled Beta ... $Y = u(X) = c*X$ will be:

$$f_c(y) = \frac{1}{c} \cdot f(\frac{y}{c}) ..... \text{not .... } c \cdot f(\frac{y}{c})$$

This is because, by the method of transformations, the pdf of $Y=u(X)$ is:

$$ |J| f( u^{-1}(y)) $$

where $x=u^{-1}(y)$ is the inverse of the transformation equation $Y = c X$, and $J=\frac{d u^{-1}(y)}{dy}$ denotes the Jacobian of the transformation . You have calculated the Jacobian of the transformation incorrectly.

P.S. What is described here is more correctly termed a scaled Beta, not a generalised Beta (which has generalised upper and lower bounds).

Addendum

In response to the comments below:

Let $X$ ~ standard Beta with pdf $f(x)$:

        f = (x^(a-1) (1-x)^(b-1))/Beta[a, b];
 domain[f] = {x, 0, 1} && {a > 0, b > 0}; 

... then, using Mathematica here, you wish to confirm:

Integrate[(1/c) f /. x -> y/c, {y, 0, c}, Assumptions -> {a > 0, b > 0, c > 0}] 

1

wolfies
  • 5,174