Computes |mean(exp(i * (phi - theta)))|, a scalar summary of circular
alignment between two angular variables.
See also
Other directional dependence indicators:
conditional_variation_index(),
local_circ_moment()
Examples
dat <- simulate_toroidal(n = 80, scenario = "diagonal", seed = 1)
alignment_index(dat$theta, dat$phi)
#> [1] 0.9545508