R Как сделать выборку из прерванной кривой перевернутого колокола

Question

Я задал связанный вопрос , перед которым успешно получил ответ. Теперь я хочу сэмплировать значения из кривой перевернутого колокольчика, но исключить диапазон значений, которые находятся посередине, как показано на рисунке ниже:

У меня есть это код, работающий в данный момент:

min <- 1
max <- 20
q <- min + (max-min)*rbeta(10000, 0.5, 0.5)

Как я могу адаптировать его для достижения желаемого результата?

Answer 1

Скажем, вы хотите получить образец 10000 из своего дистрибутива, но не хотите, чтобы в вашем образце были цифры от 5 до 15. Почему бы просто не сделать:

q <- min + (max-min)*rbeta(50000, 0.5, 0.5); 
q <- q[!(q > 5 & q < 15)][1:10000]

Что дает вам это:

hist(q)

Но все равно имеет правильный размер:

length(q)
#> [1] 10000