Transformations to symmetry are frequently used in statistics to improve the accuracy of statistical models that assume the normal distribution or at least symmetrically distributed errors. Typically the transformation is chosen after the data has been observed but it is of theoretical interest to see what is the best symmetrizing transformation for various possible distributions.

We show the probability density function for the transformed random variable,

,

,

where

and

is a random variable with support on

. This table shows the value of

for each of the distributions that makes skewness of

close to zero for the distributions used.

For the first three distributions, the usual skewness coefficient is defined for

, so the best

to symmetrize the distribution was found using

*Mathematica*'s built-in function

FindMinimum to minimize the absolute skewness coefficient. The skewness coefficient may be computed using

*Mathematica*'s built-in function

NIntegrate to obtain the first three moments. In the case of the inverse gamma distribution the mean and variance are undefined so a different definition of skewness is needed. The simplest definition for this problem is to take the difference in areas with respect to the mode. Using this definition,

was obtained. The values given in the table agree visually with the display.

In general, the value of

to symmetrize the distribution may depend on the shape parameter. Of course, for some distributions, such as bimodal ones, there may not exist a symmetrizing transformation of the type discussed in this Demonstration.

It is of interest to note that increasing

pushes the right tail of the distribution out toward

and pulls the left tail in towards 0. Decreasing

reverses this effect so the right tail is pulled in and the left tail is pushed out toward

. This effect is sometimes referred to as the ladder of transformations with increasing/decreasing

corresponding to going up/down the ladder.