R/MAD.R

In Ortem-a/mypackage: Calculate main statistic values

# Считает среднее абсолютное отклонение входящей таблицы
MAD <- function(df) {
  #' Calculate mean absolute deviation
  #'
  #' @description This function calculate the mean absolute deviation
  #' of whole input table.
  #' @param df Numeric 2 dimension table
  #' @return Good way: one number - mean absolute deviation of whole input table
  #'
  #' Bad way: -1 if input data are invalid
  #' @details The inputs can be also numeric vector.
  #' @examples
  #' MAD(data.frame)
  #' MAD(matrix)
  #' MAD(vector)
  #' MAD(List)
  #' @export
  
  # посчитать кол-во эл-тов
  if (nrow(df) == 0){
    n <- ncol(df)
  } else if (ncol(df) == 0) {
    n <- nrow(df)
  } else {
    n <- nrow(df) * ncol(df)
  }
  if (n == 0){
    return(-1)
  }
  # найти сумму квадратов разности элемента и средней арифметической
  # найти среднее арифметическое
  sum <- 0 # сумма всех элементов
  for (i in 1:nrow(df)){
    for (j in 1:ncol(df)){
      sum <- sum + df[i,j] #df[[i]][j]
    }
  }
  m <- sum / n # средняя арифметическая
  sum <- 0 # сумма
  for (i in 1:nrow(df)){
    for (j in 1:ncol(df)){
      sum <- sum + abs(df[i,j] - m)
    }
  }
  result <- sum / n
  return(result)
}