R/polysimul2_stats.R

In INESSS-QC/polymed1: Polymédication

#' Statistiques descriptives
#'
#' Statistiques descriptives des statistiques descriptives de `polysimul1()`.
#'
#' @param data_polysimul1 dataset créé par la fonction `polysimul1_stats()`.
#'
#' @import data.table
#' @export
polysimul2_stats <- function(data_polysimul1){
  x <- copy(data_polysimul1); rm(data_polysimul1)
  stats_desc <- lapply(names(x)[2:ncol(x)], function(col){  # statistiques descriptives pour chaque colonne (sauf ID)
    x[  # statistiques descriptives
      , .(Moyenne = mean(x[[col]]),
          `Écart-type` = sd(x[[col]]),
          Min = min(x[[col]]),
          P5 = quantile(x[[col]], 0.05),
          P10 = quantile(x[[col]], 0.10),
          Q1 = quantile(x[[col]], 0.25),
          `Médiane` = median(x[[col]]),
          Q3 = quantile(x[[col]], 0.75),
          P90 = quantile(x[[col]], 0.90),
          P95 = quantile(x[[col]], 0.95),
          Max = max(x[[col]]),
          n = uniqueN(x$ID))
      ]
  })
  names(stats_desc) <- names(x)[2:ncol(x)]  # nommer les éléments de la liste créée
  return(stats_desc)
}