---
title: Projected Rayleigh Statistic (PRS)
url: https://www.emergentmind.com/topics/projected-rayleigh-statistic-prs
type: topic
---

# Projected Rayleigh Statistic (PRS)

The projected Rayleigh statistic (PRS) is a modification of the classic Rayleigh statistic for testing non-uniform relative orientation between two pseudo-vector fields, particularly when the relevant alternatives are preferential parallel or preferential perpendicular alignment. It was introduced by Jow et al. as a statistic tuned to detect a mean direction of zero after mapping axial relative-orientation data to circular data, and it has since been used in analyses of magnetic-field and gas-structure alignment in molecular clouds [1708.04063; 2303.07410].

## 1. Definition and inferential target

Suppose two overlapping pseudo-vector fields on the sky are denoted \(U\) and \(V\). At each position \(i\), one forms the relative angle
\[
\phi_i = \arctan\!\left(\frac{|U_i \times V_i|}{U_i \cdot V_i}\right),
\]
with \(\phi_i \in [-\pi/2,\pi/2]\). The inferential question is whether the \(\phi_i\) are uniformly distributed, corresponding to no preferred alignment, or whether they cluster around \(0^\circ\) or \(90^\circ\), corresponding respectively to parallel or perpendicular relative orientation [1708.04063].

The PRS, denoted \(Z_x\), is specifically tuned to test preferential alignment at \(\phi=0\). A strongly positive \(Z_x\) indicates statistically significant parallel alignment, while a strongly negative \(Z_x\) indicates perpendicular alignment. In this sense, the PRS is not merely a generic non-uniformity test; it is a directional test adapted to the geometry of pseudo-vectors, or spin-2 quantities, for which \(\phi\) and \(\phi+\pi\) are equivalent [1708.04063].

In molecular-cloud applications, the same construction is used with \(\nabla I\), the gradient of a gas-tracer Moment 0 map or intensity map, and \(\hat E\), the unit vector of the polarization electric field. In that setting, \(\phi=0^\circ\) means magnetic field and gas structures are locally parallel, while \(\phi=90^\circ\) means perpendicular [2303.07410].

## 2. Derivation from the classic Rayleigh test

The starting point is the classic Rayleigh statistic for circular data \(\{\theta_i\}\):
\[
Z \equiv \frac{\left(\sum_i \cos \theta_i\right)^2 + \left(\sum_i \sin \theta_i\right)^2}{n}.
\]
Under the null hypothesis of uniform \(\theta_i\), \(Z\to 0\) on average, whereas clustering away from uniformity drives \(Z\) large [1708.04063].

For relative-orientation data between pseudo-vectors, the angles are axial rather than circular. Jow et al. therefore use angle doubling:
\[
\theta_i \equiv 2\phi_i \in [-\pi,\pi].
\]
Under this mapping, parallel alignment at \(\phi=0\) maps to \(\theta=0\), while perpendicular alignment at \(\phi=\pm \pi/2\) maps to \(\theta=\pm \pi\). Rather than using the full two-dimensional Rayleigh vector sum, the PRS projects that sum onto the \(x\)-axis only, because the target alternative is a mean direction of zero. The resulting statistic is
\[
Z_x = \frac{\sum_i \cos \theta_i}{\sqrt{n/2}}
     = \frac{\sum_i \cos 2\phi_i}{\sqrt{n/2}}.
\]

Under the null hypothesis of uniform \(\phi_i\), one has \(E[\cos\theta_i]=0\) and \(\mathrm{Var}[\cos\theta_i]=1/2\). By the central-limit theorem for large \(n\), \(Z_x \mapsto N(0,1)\). Jow et al. also give the measurement-error-driven variance of a single \(Z_x\) as
\[
\sigma^2_{Z_x} = \frac{2\sum_i (\cos\theta_i)^2 - (Z_x)^2}{n}.
\]
This normalization is the basis for interpreting \(Z_x\) directly in units of Gaussian significance [1708.04063].

## 3. Computation, weighting, and significance

In practical use, the computation is straightforward. One first measures two pseudo-vectors at each pixel, such as a column-density gradient direction and a polarization pseudo-vector. One then computes the relative angle \(\phi_i\), doubles it to obtain \(\theta_i = 2\phi_i\), and evaluates the cosine sum. If no weighting is used, the denominator is \(\sqrt{n/2}\), where \(n\) is the number of samples [1708.04063].

Jow et al. also define a weighted form,
\[
Z_x^{*} = \frac{\sum_i w_i \cos \theta_i}{\sqrt{\sum_i w_i^2/2}},
\]
where the per-angle weights \(w_i\) may represent, for example, \(S/N^2\) of polarization. Under uniform \(\phi_i\), \(Z_x^{*}\) also tends to \(N(0,1)\). The weighted form is therefore a direct extension of the same null calibration while allowing heterogeneous uncertainties or data quality across the map [1708.04063].

In the formulation summarized by Mazzei et al., the computation begins from synthetic or observed Stokes \(I,Q,U\) maps, from which one computes the polarization fraction \(p\), the polarization angle \(\chi\), and hence \(\hat E\). The local gradient vector \(\nabla I\) is computed from the same \(I\) map or from a Moment 0 line map, and the relative angle is written as
\[
\phi = \arctan\bigl(|\nabla I \times \hat E|,\;\nabla I \cdot \hat E\bigr).
\]
The PRS is then evaluated using the number of independent samples \(n_{\rm ind}\), with an oversampling correction implemented via the white-noise map protocol. In that application, \(Z_x>0\) is interpreted as a global tendency toward parallel alignment and \(Z_x<0\) as a tendency toward perpendicular alignment; \(|Z_x|\gtrsim 3\) corresponds to a \(\gtrsim 3\sigma\) detection of a preferred orientation [2303.07410].

## 4. Statistical power relative to histogram binning

A central motivation for the PRS is improved statistical efficiency relative to histogram-binning methods. The relevant comparator in Jow et al. is the histogram of relative orientations (HRO) shape statistic
\[
\zeta = \frac{A_c - A_e}{A_c + A_e},
\]
where \(A_c\) counts central angles with \(|\phi|<22.5^\circ\) and \(A_e\) counts edge angles with \(67.5^\circ<|\phi|<90^\circ\). Under the null hypothesis, \(\zeta\) is approximately Gaussian \(N(0,\sigma_\zeta^2(n))\), but it ignores \(50\%\) of the angular range [1708.04063].

Jow et al. compare the two procedures using von Mises alternatives \(\phi \sim \mathrm{VM}(\mu=0,\kappa)\). For \(\kappa\) from \(0\to 0.5\) and detection thresholds of \(3\sigma\), \(4\sigma\), and \(5\sigma\), the PRS power curves lie systematically above those of \(\zeta\). At intermediate \(\kappa\), PRS power can exceed HRO power by \(\simeq 24\%\). In a second experiment with fixed \(\kappa=0.05\) and \(n\) varying from \(5{,}000\to 50{,}000\), the PRS reaches \(99\%\) power at \(n\approx 33{,}000\) for a \(4\sigma\) threshold, whereas HRO requires \(n\approx 39{,}000\) [1708.04063].

The stated conclusion is that the PRS uses every angle continuously and is the Neyman-Pearson optimal test for von Mises alternatives with known mean direction, whereas HRO loses information in binning. This suggests that the primary gain is not a change in the scientific question being asked, but a more efficient use of the same relative-orientation data [1708.04063].

## 5. Empirical performance in the Vela C molecular cloud

Jow et al. apply the PRS to the Vela C molecular cloud using BLASTPol polarization at \(250\), \(350\), and \(500\,\mu\mathrm{m}\) for \(V_i\) and Herschel column-density gradients for \(U_i\). The four sub-regions South-Nest, South-Ridge, Centre-Nest, and Centre-Ridge contain \(n \approx 21{,}000\)–\(50{,}000\) independent \(\phi_i\) [1708.04063].

The explicit values reported for several region-band combinations are:

| Region / Band | HRO \(|\zeta|\) (\(\sigma\)) | PRS \(|Z_x|\) (\(\sigma\)) |
|---|---:|---:|
| South-Nest, \(250\,\mu\mathrm{m}\) | \(5.9 \pm 0.2\) | \(7.7 \pm 0.1\) |
| South-Nest, \(350\,\mu\mathrm{m}\) | \(4.8 \pm 0.2\) | \(6.1 \pm 0.1\) |
| South-Nest, \(500\,\mu\mathrm{m}\) | \(3.9 \pm 0.2\) | \(4.9 \pm 0.1\) |
| South-Ridge, \(250\,\mu\mathrm{m}\) | \(6.1 \pm 0.2\) | \(7.5 \pm 0.1\) |
| South-Ridge, \(350\,\mu\mathrm{m}\) | \(5.0 \pm 0.2\) | \(6.1 \pm 0.1\) |

In every band and region, \(|Z_x| > |\zeta|\), with the largest factor being \(\simeq 1.3\) in South-Nest at \(250\,\mu\mathrm{m}\). Jow et al. state that the PRS reveals the same transition from parallel to perpendicular orientation with increasing \(N_H\), but with tighter uncertainties and higher significance. The Vela C application therefore serves as an observational demonstration of the same efficiency gain already seen in the power calculations [1708.04063].

## 6. Broader uses, physical interpretation, and limitations

Jow et al. describe the PRS as applicable wherever two orientation fields are to be compared, especially when one or both are spin-2 quantities. The examples they give include the relative alignment of galaxy major axes within filaments or sheets in large-scale structure, the orientation of protostellar disks or binary orbital planes versus host-cluster shear, and small-scale polarization versus large-scale shear patterns in galaxy clusters. They also note that one may form a two-dimensional cross-correlation \(Z_x(\Delta x,\Delta y)\) by computing \(Z_x\) between one field and a spatially shifted version of the other, in order to measure correlation length scales [1708.04063].

Mazzei et al. provide a detailed later application in synthetic observations of molecular clouds. In their trans-Alfvénic simulation, which they describe as more strongly magnetized, there is a transition to perpendicular alignment at densities above \(\sim 4\times 10^3\,\mathrm{cm}^{-3}\). This transition is recovered in most synthetic observations of optically thin molecular tracers, whereas for \(^{12}\mathrm{CO}\) the PRS remains in parallel alignment across the whole observer-space. They attribute this difference to optical depth, reporting that the \(\tau=1\) surface for \(^{12}\mathrm{CO}\) is largely in front of the cloud midplane, suggesting that \(^{12}\mathrm{CO}\) mainly probes low volume density gas. In their super-Alfvénic simulation, the magnetic field becomes significantly more tangled and all observed tracers tend toward no preference for perpendicular or parallel alignment [2303.07410].

That study also articulates several limitations. Beam convolution reduces \(|Z_x|\) by washing out small-scale structure, although larger beams can partially counteract this by merging low-density sightlines into dense ones. There is a viewing-angle degeneracy: if the mean field is inclined toward the line of sight, even a strong field can yield \(Z_x \approx 0\), mimicking a weak-field case. Pixel-selection and sensitivity cuts matter because using only the top fraction of pixels can bias the inferred sign of \(Z_x\). The radiative-transfer method, such as LTE versus LVG, changes molecular optical depth and can shift \(Z_x\) accordingly. Mazzei et al. also note that assuming uniform tracer abundance is a limitation, since real chemistry may intensify the perpendicular signature by raising high-density abundances relative to low density [2303.07410].

Across these applications, the PRS functions as a single-number diagnostic of preferential alignment whose sign distinguishes parallel from perpendicular organization and whose magnitude is interpretable on an approximately Gaussian significance scale. Within the scope described by Jow et al. and Mazzei et al., its main advantages are continuous use of the full angular information, compatibility with weighting, and superior power relative to binning-based alternatives when the target alternative is known a priori [1708.04063].

Source: https://www.emergentmind.com/topics/projected-rayleigh-statistic-prs