#standard deviation
circ_sd <- function(theta){
theta <- na.omit(theta)
n<-length(theta)
C1 <- sum(cos(theta))
S1 <- sum(sin(theta))
R.bar <- sqrt(C1^2+S1^2)/n
sd <- sqrt(-2*log(R.bar))
sd
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.