R/pop_sd.r

In bmotif: Motif Analyses of Bipartite Networks

pop_sd <- function(v) {
  # compute population sd instead of sample sd
  # built in sd uses denominator n - 1
  n <- length(v)
  if(n == 1){
    return(0)
  } else {
    return(stats::sd(v) * sqrt((n - 1) / n))
  }
}