R/cl_cor.R

In Gi-Seop/ODA: Ocean Data Analysis package in R

  #circular-linear correlation
  cl_cor <- function(theta,lin){
    
    rad <- theta/(180/pi)
    rxc <- cor(cos(rad), lin)
    rxs <- cor(sin(rad), lin)
    rcs <- cor(cos(rad), sin(rad))
    R.sq <- (rxc^2+rxs^2-2*rxc*rxs*rcs)/(1-rcs^2)
    sqrt(R.sq)
    
  }