Central limit theorem for directional statistics

In probability theory, the central limit theorem states conditions under which the average of a sufficiently large number of independent random variables, each with finite mean and variance, will be approximately normally distributed. Directional statistics is the subdiscipline of statistics that deals with directions (unit vectors in Rⁿ), axes (lines through the origin in Rⁿ) or rotations in Rⁿ. The means and variances of directional quantities are all finite, so that the central limit theorem may be applied to the particular case of directional statistics.

This article will deal only with unit vectors in 2-dimensional space (R²) but the method described can be extended to the general case.

The central limit theorem

A sample of angles [math]\displaystyle{ \theta_i }[/math] are measured, and since they are indefinite to within a factor of [math]\displaystyle{ 2\pi }[/math], the complex definite quantity [math]\displaystyle{ z_i=e^{i\theta_i}=\cos(\theta_i)+i\sin(\theta_i) }[/math] is used as the random variate. The probability distribution from which the sample is drawn may be characterized by its moments, which may be expressed in Cartesian and polar form:

[math]\displaystyle{ m_n=E(z^n)= C_n +i S_n = R_n e^{i \theta_n}\, }[/math]

It follows that:

[math]\displaystyle{ C_n=E(\cos (n\theta))\, }[/math]

[math]\displaystyle{ S_n=E(\sin (n\theta))\, }[/math]

[math]\displaystyle{ R_n=|E(z^n)|=\sqrt{C_n^2+S_n^2}\, }[/math]

[math]\displaystyle{ \theta_n=\arg(E(z^n))\, }[/math]

Sample moments for N trials are:

[math]\displaystyle{ \overline{m_n}=\frac{1}{N}\sum_{i=1}^N z_i^n =\overline{C_n} +i \overline{S_n} = \overline{R_n} e^{i \overline{\theta_n}} }[/math]

where

[math]\displaystyle{ \overline{C_n}=\frac{1}{N}\sum_{i=1}^N\cos(n\theta_i) }[/math]

[math]\displaystyle{ \overline{S_n}=\frac{1}{N}\sum_{i=1}^N\sin(n\theta_i) }[/math]

[math]\displaystyle{ \overline{R_n}=\frac{1}{N}\sum_{i=1}^N |z_i^n| }[/math]

[math]\displaystyle{ \overline{\theta_n}=\frac{1}{N}\sum_{i=1}^N \arg(z_i^n) }[/math]

The vector [[math]\displaystyle{ \overline{ C_1 },\overline{ S_1 } }[/math]] may be used as a representation of the sample mean [math]\displaystyle{ (\overline{m_1}) }[/math] and may be taken as a 2-dimensional random variate. The bivariate central limit theorem states that the joint probability distribution for [math]\displaystyle{ \overline{ C_1 } }[/math] and [math]\displaystyle{ \overline{ S_1 } }[/math] in the limit of a large number of samples is given by:

[math]\displaystyle{ [\overline{C_1},\overline{S_1}] \xrightarrow{d} \mathcal{N}([C_1,S_1],\Sigma/N) }[/math]

where [math]\displaystyle{ \mathcal{N}() }[/math] is the bivariate normal distribution and [math]\displaystyle{ \Sigma }[/math] is the covariance matrix for the circular distribution:

[math]\displaystyle{ \Sigma = \begin{bmatrix} \sigma_{CC} & \sigma_{CS} \\ \sigma_{SC} & \sigma_{SS} \end{bmatrix} \quad }[/math]

[math]\displaystyle{ \sigma_{CC}=E(\cos^2\theta)-E(\cos\theta)^2\, }[/math]

[math]\displaystyle{ \sigma_{CS}=\sigma_{SC}=E(\cos\theta\sin\theta)-E(\cos\theta)E(\sin\theta)\, }[/math]

[math]\displaystyle{ \sigma_{SS}=E(\sin^2\theta)-E(\sin\theta)^2\, }[/math]

Note that the bivariate normal distribution is defined over the entire plane, while the mean is confined to be in the unit ball (on or inside the unit circle). This means that the integral of the limiting (bivariate normal) distribution over the unit ball will not be equal to unity, but rather approach unity as N approaches infinity.

It is desired to state the limiting bivariate distribution in terms of the moments of the distribution.

Covariance matrix in terms of moments

Using multiple angle trigonometric identities

[math]\displaystyle{ C_2= E(\cos(2\theta)) = E(\cos^2\theta-1)=E(1-\sin^2\theta)\, }[/math]

[math]\displaystyle{ S_2= E(\sin(2\theta)) = E(2\cos\theta\sin\theta)\, }[/math]

It follows that:

[math]\displaystyle{ \sigma_{CC}=E(\cos^2\theta)-E(\cos\theta)^2 =\frac{1}{2}\left(1 + C_2 - 2C_1^2\right) }[/math]

[math]\displaystyle{ \sigma_{CS}=E(\cos\theta\sin\theta)-E(\cos\theta)E(\sin\theta)=\frac{1}{2}\left(S_2 - 2 C_1 S_1 \right) }[/math]

[math]\displaystyle{ \sigma_{SS}=E(\sin^2\theta)-E(\sin\theta)^2 =\frac{1}{2}\left(1 - C_2 - 2S_1^2\right) }[/math]

The covariance matrix is now expressed in terms of the moments of the circular distribution.

The central limit theorem may also be expressed in terms of the polar components of the mean. If [math]\displaystyle{ P(\overline{C_1},\overline{S_1})d\overline{C_1}d\overline{S_1} }[/math] is the probability of finding the mean in area element [math]\displaystyle{ d\overline{C_1}d\overline{S_1} }[/math], then that probability may also be written [math]\displaystyle{ P(\overline{R_1}\cos(\overline{\theta_1}),\overline{R_1}\sin(\overline{\theta_1}))\overline{R_1}d\overline{R_1}d\overline{\theta_1} }[/math].