Coskewness Explained

In probability theory and statistics, coskewness is a measure of how much three random variables change together. Coskewness is the third standardized cross central moment, related to skewness as covariance is related to variance. In 1976, Krauss and Litzenberger used it to examine risk in stock market investments.^[1] The application to risk was extended by Harvey and Siddique in 2000.^[2]

If three random variables exhibit positive coskewness they will tend to undergo extreme deviations at the same time, an odd number of which are in the positive direction (so all three random variables undergoing extreme positive deviations, or one undergoing an extreme positive deviation while the other two undergo extreme negative deviations). Similarly, if three random variables exhibit negative coskewness they will tend to undergo extreme deviations at the same time, an even number of which are in the positive direction (so all three random variables undergoing extreme negative deviations, or one undergoing an extreme negative deviation while the other two undergo extreme positive deviations).

Types

There are two different measures for the degree of coskewness in data.

Coskewness

For three random variables X, Y and Z, the non-trivial coskewness statistic is defined as:^[3]

S(X,Y,Z)=

	\operatorname{E
	\left[(X-\operatorname{E}[X])(Y-\operatorname{E}[Y])(Z-\operatorname{E}[Z])\right]}

{\sigma_X\sigma_Y\sigma_Z}

where E[''X''] is the expected value of X, also known as the mean of X, and

\sigma_X

is the standard deviation of X.

Standardized rank coskewness

Bernard, Chen, Rüschendorf and Vanduffel defined the standardized rank coskewness of three random variables X, Y and Z as:^[4]

RS(X,Y,Z)=32\operatorname{E}

\left[\left(F

X(X)-	1
	2

\right)\left(F

Y(Y)-	1
	2

\right)\left(F

Z(Z)-	1
	2

\right)\right]

where F_X (X), F_Y (Y) and F_Z (Z) are the cumulative distribution functions of X, Y and Z, respectively.

Properties

Skewness is a special case of the coskewness when the three random variables are identical:

S(X,X,X)=

	\operatorname{E
	\left[(X

	3}
\operatorname{E}[X])
	X

={\operatorname{skewness}[X]},

For two random variables, X and Y, the skewness of the sum, X + Y, is

S_X+Y={1\over

	3}{\left[
\sigma
	X+Y

	3S
\sigma
	X

	2\sigma
3\sigma
	YS(X,X,Y)

+3\sigma_X\sigma

	2S(X,Y,Y)

	Y

	3S
\sigma
	Y

\right]},

where S_X is the skewness of X and

\sigma_X

is the standard deviation of X. It follows that the sum of two random variables can be skewed (S_X+Y ≠ 0) even if both random variables have zero skew in isolation (S_X = 0 and S_Y = 0).

The standardized rank coskewness RS(X, Y, Z) satisfies the following properties:^[4]

(1) −1 ≤ RS(X, Y, Z) ≤ 1.

(2) The upper bound of 1 is obtained by the copula given in (3.3) in Bernard, Chen, Rüschendorf and Vanduffel (2023). The lower bound of −1 is obtained by the copula (3.5) in the same paper.

(3) It is invariant under strictly increasing transformations, i.e., when f_i, i = 1, 2, 3, are arbitrary strictly increasing functions, RS(X, Y, Z) = RS(f₁ (X), f₂ (Y), f₃ (Z)).

(4) RS(X, Y, Z) = 0 if X, Y and Z are independent.

Example

Let X be standard normally distributed and Y be the distribution obtained by setting X=Y whenever X<0 and drawing Y independently from a standard half-normal distribution whenever X>0. In other words, X and Y are both standard normally distributed with the property that they are completely correlated for negative values and uncorrelated apart from sign for positive values. The joint probability density function is

f_X,Y(x,y)=

	-x^2/2
e

\sqrt{2\pi

} \left(H(-x)\delta(x-y) + 2H(x)H(y) \frac\right)where H(x) is the Heaviside step function and δ(x) is the Dirac delta function. The third moments are easily calculated by integrating with respect to this density:

S(X,X,Y)=S(X,Y,Y)=-

	1
	\sqrt{2\pi

} \approx -0.399

Note that although X and Y are individually standard normally distributed, the distribution of the sum X+Y is significantly skewed. From integration with respect to density, we find that the covariance of X and Y is

\operatorname{cov}(X,Y)=

	1
	2

	1
	\pi

from which it follows that the standard deviation of their sum is

\sigma_X+Y=\sqrt{3+

	2
	\pi

}Using the skewness sum formula above, we have

S_X+Y=-

	3\sqrt{2
	\pi}{(2+3\pi)

^3/2

} \approx -0.345This can also be computed directly from the probability density function of the sum:

f_X+Y(u)=

	-u^2/8
e

2\sqrt{2\pi

} H(-u) + \frac \operatorname\left(\frac\right) H(u)

Bernard, Chen, Rüschendorf and Vanduffel (2023) found risk bounds on coskewness for some popular marginal distributions as shown in the following table.^[4]

Marginal distributions

Minimum coskewness

Maximum coskewness

\mu_i

	2
\sigma
	i

)

-	2\sqrt{2\pi

}

	2\sqrt{2\pi

}

Student(

\nu

\nu>3

4(\nu-2)\sqrt{(\nu-2

)\pi\Gamma(	\nu+1
	2

2)\pi\Gamma(	\nu
	2

)}{(3-4\nu+\nu

)}

4(\nu-2)\sqrt{(\nu-2

)\pi\Gamma(	\nu+1
	2

2)\pi\Gamma(	\nu
	2

)}{(3-4\nu+\nu

)}

Laplace(

\mu_i

b_i

)

-	3\sqrt{2

}

	3\sqrt{2

}

a_i

b_i

)

-	3\sqrt{3

}

	3\sqrt{3

}

where

\Gamma(x)

is the gamma function.

Notes and References

Friend . Irwin . Randolf Westerfield . Co-Skewness and Capital Asset Pricing . The Journal of Finance . 1980 . 35 . 4 . 897–913 . 10.1111/j.1540-6261.1980.tb03508.x.
Book: Jondeau . Eric . Ser-Huang Poon. Michael Rockinger. Financial Modeling Under Non-Gaussian Distributions. 2007 . Springer . 978-1-84628-696-4. 31–32.
Book: Miller. Michael B.. Mathematics and Statistics for Financial Risk Management . Chapter 3. Basic Statistics . 2014 . John Wiley & Sons, Inc.. Hoboken, New Jersey. 978-1-118-75029-2. 53–56. 2nd. http://www.wiley.com/WileyCDA/WileyTitle/productCd-1118750292.html.
Bernard . Carole . Jinghui . Chen . Rüschendorf . Ludger . Vanduffel . Steven . Coskewness under dependence uncertainty . Statistics and Probability Letters . 5 May 2023 . 199 . 8.