Outlier Removal (5SD to NA)

outs5SD.R

outs5SD <- function(vec, n_outsONLY = FALSE){
  m = mean(vec, na.rm = T)
  s = sd(vec, na.rm = T)
  # if less than 5 sd below or above the mean replace with na's
  print(paste0(length(vec[vec < m - s*5 | vec > m + s*5]), " number of outliers above or below 5SD"))
  vec[vec < m - s*5 | vec > m + s*5] <- rep(NA, length(vec[vec < m - s*5 | vec > m + s*5]))
  if(n_outsONLY == FALSE) {
    return(vec)
  } else {} # no return of vector
}