Circular-linear correlation

Implementation of the circular-linear correlation measure introduced by Mardia (1976) and Johnson and Wehrly (1977) as cited in Jammalamadaka & Sengupta (2001).

Usage

circ_lin_corr(circ_x, lin_x, na.rm = FALSE)

Arguments

circ_x

circular variable

lin_x

linear variable

na.rm

a logical value indicating whether NA values should be removed before the computation proceeds

Value

circular-linear correlation measure

Details

This measure is computed as r^2 = (r_xc^2+r_xs^2-2 r_xc r_xsr_cs)/(1-r_cs^2) where r_xc = corr(x, cos()), r_xs = corr(x, sin()), r_cs = corr(cos(), sin()), and and x are the circular and linear variables, respectively.

References

Jammalamadaka, S. R., & SenGupta, A. (2001). Topics in Circular Statistics. WORLD SCIENTIFIC. doi:10.1142/4031

Examples

x <- rnorm(50)
a <- as.vector(circular::rvonmises(50, 0, 5))
#> Warning: an object is coerced to the class 'circular' using default value for the following components:
#>   type: 'angles'
#>   units: 'radians'
#>   template: 'none'
#>   modulo: 'asis'
#>   zero: 0
#>   rotation: 'counter'
#> conversion.circularmuradians0counter
circ_lin_corr(x + a, x)
#> [1] 0.8024563