如何在 R 中模拟加载的骰子
How to simulate a loaded dice in R
考虑一个已加载的骰子,结果为 1 的概率为 2p/3,获得 2、3、4 或 5 的概率各为 p,获得 6 的概率为 3p/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))
考虑一个已加载的骰子,结果为 1 的概率为 2p/3,获得 2、3、4 或 5 的概率各为 p,获得 6 的概率为 3p/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))