Cokurtosis Explained

In probability theory and statistics, cokurtosis is a measure of how much two random variables change together. Cokurtosis is the fourth standardized cross central moment.^[1] If two random variables exhibit a high level of cokurtosis they will tend to undergo extreme positive and negative deviations at the same time.

Definition

For two random variables X and Y there are three non-trivial cokurtosis statistics^[2]

K(X,X,X,Y)={\operatorname{E}{[(X-\operatorname{E}[X])^3(Y-\operatorname{E}[Y])]}\over

	3
\sigma
	X

\sigma_Y},

K(X,X,Y,Y)={\operatorname{E}{[(X-\operatorname{E}[X])^2(Y-\operatorname{E}[Y])^2]}\over

	2
\sigma
	X

	2},
\sigma
	Y

and

K(X,Y,Y,Y)={\operatorname{E}{[(X-\operatorname{E}[X])(Y-\operatorname{E}[Y])^3]}\over\sigma_X

	3},
\sigma
	Y

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.

Properties

Kurtosis is a special case of the cokurtosis when the two random variables are identical:

K(X,X,X,X)={\operatorname{E}{[(X-\operatorname{E}[X])^4]}\over

	4}
\sigma
	X

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

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

\begin{align} K_X+Y={1\over

	4}
\sigma
	X+Y

	4K
\sigma
	X

	3\sigma
4\sigma
	YK(X,X,X,Y)

	2K(X,X,Y,Y)
6\sigma
	Y

\\ &{}+4\sigma_X\sigma

	3K(X,Y,Y,Y)

	Y

	4K
\sigma
	Y

], \end{align}

where

K_X

is the kurtosis of X and

\sigma_X

is the standard deviation of X.

It follows that the sum of two random variables can have kurtosis different from 3 (

K_X+Y ≠ 3

) even if both random variables have kurtosis of 3 in isolation (

K_X=3

and

K_Y=3

The cokurtosis between variables X and Y does not depend on the scale on which the variables are expressed. If we are analyzing the relationship between X and Y, the cokurtosis between X and Y will be the same as the cokurtosis between a + bX and c + dY, where a, b, c and d are constants.

Examples

Bivariate normal distribution

Let X and Y each be normally distributed with correlation coefficient ρ. The cokurtosis terms are

K(X,X,Y,Y)=1+2\rho²

K(X,X,X,Y)=K(X,Y,Y,Y)=3\rho

Since the cokurtosis depends only on ρ, which is already completely determined by the lower-degree covariance matrix, the cokurtosis of the bivariate normal distribution contains no new information about the distribution. It is a convenient reference, however, for comparing to other distributions.

Nonlinearly correlated normal distributions

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 fourth moments are easily calculated by integrating with respect to this density:

K(X,X,Y,Y)=2

K(X,X,X,Y)=K(X,Y,Y,Y)=

	3	+
	2

	2
	\pi

≈ 2.137

It is useful to compare this result to what would have been obtained for an ordinary bivariate normal distribution with the usual linear correlation. From integration with respect to density, we find that the linear correlation coefficient of X and Y is

\rho=

	1
	2

	1
	\pi

≈ 0.818

A bivariate normal distribution with this value of ρ would have

K(X,X,Y,Y) ≈ 2.455

and

K(X,X,X,Y) ≈ 2.339

. Therefore, all of the cokurtosis terms of this distribution with this nonlinear correlation are smaller than what would have been expected from a bivariate normal distribution with ρ=0.818.

Note that although X and Y are individually standard normally distributed, the distribution of the sum X+Y is platykurtic. The standard deviation of the sum is

\sigma_X+Y=\sqrt{3+

	2
	\pi

}Inserting that and the individual cokurtosis values into the kurtosis sum formula above, we have

K_X+Y=

	2\pi(8+15\pi)
	(2+3\pi)²

≈ 2.654

This 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)

Notes and References

Book: Miller. Michael B.. Mathematics and Statistics for Financial Risk Management. 2014. John Wiley & Sons, Inc.. Hoboken, New Jersey. 978-1-118-75029-2. 53–56. 2nd.
Book: Meucci. Attilio. Risk and Asset Allocation. 2005. Springer-Verlag. Berlin. 978-3642009648. 58–59.