Merci pour vos réponses, j'ai finalement utilisé l'estimation de la fonction de répartition et la méthode d'inversion de la fonction de répartition: y <- c(rlnorm(100), rpois(100, 5), rnorm(100), rnorm(100, 3), rbinom(100, 1, 0.5), rbeta(100, 0.5, 5)) Ecdf(y) p=Ecdf(y) f=function(u){ tempo=p$x[which...