Как смоделировать загруженные кости в R

Question

Рассмотрим загруженный кубик так, что вероятность получить результат 1 равна 2p / 3, вероятность получения 2, 3, 4 или 5 равна p каждый, а вероятность получения 6 равна 3p / 2.

1. Создайте компьютерную программу для имитации загруженных выше кубиков. Генерировать вывод.
2. Используйте выходные данные для построения графика, чтобы убедиться, что ваша программа работает правильно.

Я обнаружил, что p равно 6/37, но я не знаю, как кодировать загруженные кости.

Я знаю, как симулировать честные кости:

sample(1:6, size=1, replace = TRUE)

Answer 1

Вот один из способов сделать это: вы можете установить вероятности внутри выборки, указав аргумент prob.

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

p <- 6/37
myprobs <- c(2/3 * p, p, p, p, p, 3/2*p)

set.seed(1)
x <- sample(1:6, size = 100, replace = TRUE, prob = myprobs)

barplot(table(x))