\(\newcommand{\R}{\mathbb{R}}\)
\(\newcommand{\N}{\mathbb{N}}\)
\(\newcommand{\E}{\mathbb{E}}\)
\(\newcommand{\P}{\mathbb{P}}\)
\(\newcommand{\var}{\text{var}}\)
\(\newcommand{\sd}{\text{sd}}\)
\(\newcommand{\cov}{\text{cov}}\)
\(\newcommand{\cor}{\text{cor}}\)
\(\newcommand{\skw}{\text{skew}}\)
\(\newcommand{\kur}{\text{kurt}}\)

- Random
- 4. Special Distributions
- The U-Power Distribution

The U-power distribution is a U-shaped family of distributions based on a simple family of power functions.

Random variable \(Z\) has the standard U-power distribution with shape parameter \( k \in \N_+ \) if \(Z\) has a continuous distribution on \( [-1, 1] \) with probability density function \(g\) given by \[g(x) = \frac{2 k + 1}{2} x^{2 k}, \quad x \in [-1, 1] \]

The algebraic form of the probability density function explains the name of the distribution. The most common of the standard U-power distributions is the U-quadratic distribution, which corresponds to \( k = 1 \).

The function \(g\) is a valid probability density function and satisfies the following properties:

- \(g\) is symmetric about \( x = 0 \).
- \(g\) decreases and then increases with minimum value at \( x = 0 \).
- The modes are \(x = \pm 1 \).
- \( g \) is concave upward.

These results follow from simple calculus. In particular, note that \[ \int_{-1}^1 x^{2 k} dx = \frac{2}{2 k + 1} \]

Open the Special Distribution Simulator and select the U-power distribution. Vary the shape parameter but keep the default values for the other parameters. Note the graph of the probability density function. For selected values of the shape parameter, run the simulation 1000 times and compare the emprical density function to the probability density function.

\(Z\) has distribution function \(G\) given by \[ G(x) = \frac{1}{2} \left(1 + x^{2 k + 1}\right), \quad x \in [-1, 1] \]

This follows from the PDF given above and simple calculus.

\(Z\) has quantile function \(G^{-1}\) given by \( G^{-1}(p) = (2 p - 1)^{1/(2 k + 1)} \) for \( p \in [0, 1] \).

- \( G^{-1}(1 - p) = -G^{-1}(p) \)
- The first quartile is \(q_1 = -\frac{1}{2^{1/(2 k + 1)}} \)
- The median is 0.
- The third quartile is \(q_3 = \frac{1}{2^{1/(2 k + 1)}} \)

The formula for the quantile function follows immediately from the CDF given above by solving \(p = G(x)\) for \(x\) in terms of \(p \in [0, 1]\). Property (a) follows from the symmetry of the distribution about 0.

Open the Special Distribution Calculator and select the U-power distribution. Vary the shape parameter but keep the default values for the other parameters. Note the shape of the distribution function. For various values of the shape parameter, compute a few quantiles.

Suppose again that \( Z \) has the standard U-power distribution with parameter \( k \in \N_+ \)

The odd-order moments of \( Z \) (about 0) are 0. The even-order moments of \( Z \) are given by \[ \E\left(Z^{2 n}\right) = \frac{2 k + 1}{2 (n + k) + 1}, \quad n \in \N \]

This result follows from simple calculus. The fact that the even order moments are 0 also follows from the symmetry of the distribution about 0.

Since the mean is 0, the moments about 0 are also the central moments.

Note that \( \var(Z) \to 1 \) as \( k \to \infty \).

Open the Special Distribution Simulator and select the U-power distribution. Vary the shape parameter but keep the default values for the other paramters. Note the position and size of the mean \(\pm \) standard deviation bar. For selected values of the shape parameter, run the simulation 1000 times and compare the empirical mean and stadard deviation to the distribution mean and standard deviation.

The skewness and kurtosis of \(Z\) are

- \(\skw(Z) = 0\)
- \(\kur(Z) = \frac{(2 k + 3)^2}{(2 k + 5)(2 k + 1)}\)

The skewness is 0 by the symmetry of the distribution. Since the mean is 0, the kurtosis is \( \E(Z^4) / [\E(Z^2)]^2 \) and so the result follows from the general moment result above.

Note that \( \kur(Z) \to 1 \) as \( k \to \infty \). The excess kurtosis is \( \kur(Z) - 3 = \frac{(2 k + 3)^2}{(2 k + 5)(2 k + 1)} - 3 \).

The U-power probability density function \( g \) actually makes sense for \( k = 0 \) as well, and in this case the distribution reduces to the uniform distribution on the interval \( [-1, 1] \). But of course, this distribution is not U-shaped, except in a degenerate sense. There are other connections to the uniform distribution. The first is a standard result since the quntile function has a simple, closed representation:

Suppose that \( k \in \N_+ \).

- If \( U \) has the standard uniform distribution then \( Z = (2 U - 1)^{1/(2 k + 1)} \) has the standard U-power distribution with parameter \( k \).
- If \( Z \) has the standard U-power distribution with parameter \( k \) then \( U = \frac{1}{2} \left(1 + Z^{2 k + 1} \right) \) has the standard uniform distribution.

Part (a) of course leads to the random quantile method of simulation.

Open the random quantile simulator and select the U-power distribution. Vary the shape parameter but keep the default values for the other parameters. Note the shape of the distribution and density functions. For selected values of the parameter, run the simulation 1000 times and note the random quantiles. Compare the empirical density function to the probability density function.

The standard U-power distribution with shape parameter \( k \) converges to the discrete uniform distribution on \( \{-1, 1\} \) as \( k \to \infty \).

This follows from the definition of convergence in distribution. The U-power distribution function \( G \) is 0 on \( (-\infty, -1] \), is 1 on \( [1, \infty) \), and is given by the formula in (3) on \( [-1, 1] \). As \( k \to \infty \), \( G(x) \to 0 \) for \( x \in (-\infty, -1) \), \( G(x) \to \frac{1}{2} \) for \( x \in (-1, 1) \), and \( G(x) \to 1 \) for \( x \in (1, \infty) \). This agrees with the distribution function of the discrete uniform distribution on \( \{-1, 1\} \) except at the points of discontinuity \( \pm 1 \).

Like so many standard distributions, the standard U-power distribution is generalized by adding location and scale parameters.

If \(Z\) has the standard U-power distribution with shape parameter \( k \in \N_+ \), and if \(\mu \in \R\) and \(c \in (0, \infty)\), then \(X = \mu + c Z\) has the U-power distribution with shape parameter \( k \), location parameter \( \mu \) and scale parameter \(c\).

Note that \( X \) takes values in the interval \( [a, b] \) where \( a = \mu - c \) and \( b = \mu + c \), so the distribution can also be parameterized by the the shape parameter \( k \) and the endpoints \( a \) and \( b \). With this parametrization, the location parameter is \( \mu = \frac{a + b}{2} \) and the scale parameter is \( c = \frac{b - a}{2} \).

\(X\) has probability density function \(f\) given by \[f(x) = \frac{2 k + 1}{2 c} \left(\frac{x - \mu}{c}\right)^{2 k}, \quad x \in [\mu - c, \mu + c] \]

- \(f\) is symmetric about \( \mu \).
- \(f\) decreases and then increases with minimum value at \( x = \mu \).
- The modes are at \( x = \mu \pm c \).
- \( f \) is concave upward.

Recall that \(f(x) = \frac{1}{c} g\left(\frac{x - \mu}{c}\right)\) where \(g\) is the PDF of the standard U-power distribution with shape parameter \( k \), given above.

Open the Special Distribution Simulator and select the U-power distribution. Vary the parameters and note the shape and location of the probability density function. For various values of the parameters, run the simulation 1000 times and compare the emprical density function to the probability density function.

\(X\) has distribution function \(F\) given by \[F(x) = \frac{1}{2}\left[1 + \left(\frac{x - \mu}{c}\right)^{2 k + 1}\right], \quad x \in [\mu - c, \mu + c] \]

Recall that \(F(x) = G\left(\frac{x - \mu}{c}\right)\) where \(G\) is the CDF of the standard U-power distribution with shape parameter \( k \), given above.

\(X\) has quantile function \(F^{-1}\) given by \(F^{-1}(p) = \mu + c (2 p - 1)^{1/(2 k + 1)} \) for \(p \in [0, 1] \).

- \( F^{-1}(1 - p) = \mu - c F^{-1}(p) \)
- The first quartile is \(q_1 = \mu - c \frac{1}{2^{1/(2 k + 1)}} \)
- The median is \( \mu \).
- The third quartile is \(q_3 = \mu + c \frac{1}{2^{1/(2 k + 1)}} \)

Recall that \(F^{-1}(p) = \mu + c G^{-1}(p)\) where \(G^{-1}\) is the quantile function of the standard U-power distribution with shape parameter \( k \), given above.

Open the Special Distribution Calculator and select the U-power distribution. Vary the parameters and note the graph of the distribution function. For various values of the parameters, compute selected values of the distribution function and the quantile function.

Suppose again that \( X \) has the U-power distribution with shape parameter \( k \in \N_+ \), location parameter \( \mu \in \R \), and scale parameter \( c \in (0, \infty) \).

The mean and variance of \( X \) are

- \(\E(X) = \mu \)
- \(\var(X) = c^2 \frac{2 k + 1}{2 k + 3}\)

These results follow from the representation \( X = \mu + c Z \) where \( Z \) has the standard U-power distribution with shape parameter \( k \), and from the mean and variance of \( Z \) given above.

Note that \( \var(Z) \to c^2 \) as \( k \to \infty \)

Open the Special Distribution Simulator and select the U-power distribution. Vary the parameters and note the size and location of the mean \(\pm\) standard deviation bar. For various values of the parameters, run the simulation 1000 times and compare the empirical mean and stadard deviation to the distribution mean and standard deviation.

The odd-order central moments of \( X \) are 0. The even-order central moments of \( X \) are given by \[ \E\left[(x - \mu)^{2 n}\right] = c^{2 n} \frac{2 k + 1}{2 (n + k) + 1}, \quad n \in \N \]

This follows from the representation \( X = \mu + c Z \) where \( Z \) has the standard U-power distribution with shape parameter \( k \), and the central moments of \( Z \) given above.

The skewness and kurtosis of \(X\) are

- \(\skw(X) = 0\)
- \(\kur(X) = \frac{(2 k + 3)^2}{(2 k + 5)(2 k + 1)}\)

Recall that the skewness and kurtosis are defined in terms of the standard score of \(X\) and hence are invariant under a location-scale transformation. Thus, the results are the same as for the standard distribution given above.

Again, \( \kur(X) \to 1 \) as \( k \to \infty \) and the excess kurtosis is \( \kur(X) - 3 = \frac{(2 k + 3)^2}{(2 k + 5)(2 k + 1)} - 3\)

Since the U-power distribution with a given shape parameter is a location-scale family, it is trivially closed under location-scale transformations.

Suppose that \( X \) has the U-power distribution with shape parameter \( k \in \N_+ \), location parameter \( \mu \in \R \), and scale parameter \( c \in (0, \infty) \). If \( \alpha \in \R \) and \( \beta \in (0, \infty) \), then \( Y = \alpha + \beta X \) has the U-power distribution with shape parameter \( k \), location parameter \( \alpha + \beta \mu \), and scale parameter \( \beta c \).

As before, since the distribution and quantile function have simple forms, we have the usual connections with the standard uniform distribution.

Suppose that \( k \in \N_+ \), \( \mu \in \R \) and \( c \in (0, \infty) \).

- If \( U \) has the standard uniform distribution then \( X = \mu + c (2 U - 1)^{1/(2 k + 1)} \) has the U-power distribution with shape parameter \( k \), location parameter \( \mu \), and scale parameter \( c \).
- If \( X \) has the U-power distribution with shape parameter \( k \), location parameter \( \mu \), and scale parameter \( c \), then \( U = \frac{1}{2} \left[1 + \left(\frac{X - \mu}{c}\right)^{2 k + 1} \right] \) has the standard uniform distribution.

Again, part (a) of course leads to the random quantile method of simulation.

Open the random quantile simulator and select the U-power distribution. Vary the parameters and note the shape of the distribution and density functions. For selected values of the parameters, run the simulation 1000 times and note the random quantiles. Compare the empirical density function to the probability density function.

The U-power distribution with given location and scale parameters converges to the discrete uniform distribution at the endpoints as the shape parameter increases.

The U-power distribution with shape parameter \( k \in \N_+ \), location parameter \( \mu \in \R \), and scale parameter \( c \in (0, \infty) \) converges to the discrete uniform distribution on \( \{\mu - c, \mu + c\} \) as \( k \to \infty \).

This follows from the convergence result for the standard U-power distribution given above and basic properties of convergence in distribution.

The U-power distribution is a general exponential family in the shape parameter, if the location and scale parameters are fixed.

Suppose that \( X \) has the U-power distribution with unspecified shape parameter \( k \in \N_+ \), but with specified location parameter \( \mu \in \R \) and scale parameter \( c \in (0, \infty) \). Then \( X \) has a one-parameter exponential distribution with natural parameter \( 2 k \) and natural statistics \( \ln\left(\frac{X - \mu}{c}\right) \).

This follows from the definition of the general exponential family, since the PDF of the U-power distribution can be written as \[ f(x) = \frac{2 k + 1}{2 c} \exp\left[2 k \ln\left(\frac{x - \mu}{c}\right)\right], \quad x \in [\mu - c, \mu + c] \]

Since the U-power distribution has a bounded probability density function on a bounded support interval, it can also be simulated via the rejection method.

Open the rejection method experiment and select the U-power distribution. Vary the parameters and note the shape of the probability density function. For selected values of the parameters, run the experiment 1000 times and watch the scatterplot. Compare the empirical density function, mean, and standard deviation to their distributional counterparts.