ClassEx/mysd.R

In agracy2246/MATH4753grac0009: Adam Gracy MATH4753

#' Calculate the standard deviation
#'
#' @param x Numeric Vector
#'
#' @return Numeric standard deviation
#'
#' @examples
#' L <- 1:20
#' mysd(L)
#'
#' @export
mysd <- function(x){
  N        <- length(x)
  variance <- sum((x - mean(x))^2) / (N - 1)
  std      <- sqrt(variance)

  cat("Variance: ", variance, "\nDeviation: ",std, sep = "")

}