R/polysimul1_stats.R

In INESSS-QC/polymed1: Polymédication

#' Statistiques descriptives
#'
#' Statistiques descriptives pour chaque ID
#'
#' @param polysimul1_data dataset créé par la fonction `polysimul1_data()`.
#'
#' @import data.table
#' @export
polysimul1_stats <- function(polysimul1_data){
  x <- copy(polysimul1_data); rm(polysimul1_data)
  x <- polysimul1_statsParID(x)
  return(x)
}

#' Structure data_Rx
#'
#' Statistiques descriptives pour chaque individu.
#'
#' @keywords internal
#' @export
polysimul1_statsParID <- function(x){
  x <- copy(x)
  x <- x[  # statistiques descriptives
    , .(Moyenne = mean(nRx),
        `Écart-type` = sd(nRx),
        Min = min(nRx),
        P5 = quantile(nRx, 0.05),
        P10 = quantile(nRx, 0.10),
        Q1 = quantile(nRx, 0.25),
        `Médiane` = median(nRx),
        Q3 = quantile(nRx, 0.75),
        P90 = quantile(nRx, 0.90),
        P95 = quantile(nRx, 0.95),
        Max = max(nRx)),
    .(ID)
  ]
  x <- cbind(
    x[, .(ID)],
    x[, lapply(.SD, as.numeric), .SDcols = 2:12]
  )
  return(x)
}