Так что, похоже, вы действительно хотите усеченное нормальное распределение. Я бы предложил truncdist
. Затем вы можете сделать
library(truncdist)
set.seed(123)
x <- rtrunc(1000, "norm", a = 5, mean = 10, sd = 4)
, чтобы получить нужные вам образцы:
summary(x)
# Min. 1st Qu. Median Mean 3rd Qu. Max.
# 5.009 8.268 10.440 10.780 13.002 23.091
plot(density(x, from = 5))