R/circ_sd.R

Defines functions circ_sd

  #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
    
  }
  
Gi-Seop/ODA documentation built on Jan. 6, 2020, 12:49 p.m.