#circular-circular correlation
cc_cor <- function(theta1, theta2){
s1 <- sin(theta1-circ_mean(theta1))
s2 <- sin(theta2-circ_mean(theta2))
pcc <-sum(s1*s2)/sqrt(sum(s1^2)*sum(s2^2))
pcc
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.