Merci beaucoup pour la suggestion, mais cela ne semble pas produire l'effet attendu. Voici le nouveau code: tablestat <- function(formula, data= Data, FUN= alist(min, mean, median, max, var, sd, coefvar, sem, confint), prob= 0.95, ...) { TABLE <- aggregate(formula, data, sub...