set.seed(42)
n <- 500
individus <- data.frame(
bonheur = sample(1:4, n, TRUE, prob=c(0.20,0.45,0.25,0.10)),
milieu = sample(c("urbain","rural"), n, TRUE),
poids = runif(n, 0.8, 1.3)
)
bonheur_declare(individus, "bonheur", codes_heureux=c(1,2),
poids="poids", sous_groupes="milieu")
Run the code above in your browser using DataLab