set.seed(42)
n <- 1000
ages_base <- sample(10:89, n, TRUE)
# Attraction sur les ages ronds
ages <- ifelse(runif(n) < 0.20,
ages_base - (ages_base %% 10),
ages_base)
individus <- data.frame(age = pmax(10, pmin(89, ages)))
myers(individus, "age")
Run the code above in your browser using DataLab