如何在 R 中模拟加载的骰子

How to simulate a loaded dice in R

考虑一个已加载的骰子,结果为 1 的概率为 2p/3,获得 2、3、4 或 5 的概率各为 p,获得 6 的概率为 3p/2。

  1. 创建一个计算机程序来模拟上面加载的骰子。生成输出。
  2. 使用输出绘制图表以检查您的程序是否正常运行。

我发现 p 是 6/37,但我不知道如何为已加载的骰子编码。

我知道如何模拟公平的骰子:

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

这是一种方法:您可以通过指定 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))