set.seed(42)
n <- 500
individus <- data.frame(
satisfaction = pmin(10, pmax(0, round(rnorm(n, 5.8, 2.1)))),
milieu = sample(c("urbain","rural"), n, TRUE),
quintile = sample(paste0("Q", 1:5), n, TRUE),
poids = runif(n, 0.8, 1.3)
)
satisfaction_vie(individus, "satisfaction", poids="poids",
sous_groupes=c("milieu","quintile"))
Run the code above in your browser using DataLab