r语言模拟不同次数下同时抛掷两枚硬币正面朝上概率并绘图
时间: 2023-08-23 21:06:08 浏览: 66
好的,可以使用以下R代码来模拟不同次数下同时抛掷两枚硬币正面朝上的概率,并绘制图表:
```R
# 设定模拟次数
n <- c(10, 50, 100, 500, 1000, 5000, 10000)
# 创建一个空的向量来存储结果
prob <- numeric(length(n))
# 进行模拟
for(i in 1:length(n)) {
coin1 <- sample(c("H", "T"), size = n[i], replace = TRUE)
coin2 <- sample(c("H", "T"), size = n[i], replace = TRUE)
prob[i] <- mean(coin1 == "H" & coin2 == "H")
}
# 绘制结果图表
plot(n, prob, type = "b", pch = 19, xlim = c(0, 10000), ylim = c(0, 1),
xlab = "Number of Trials", ylab = "Probability of Both Heads")
```
在这段代码中,我们首先设定了不同的模拟次数,并创建一个空的向量来存储每个模拟次数下同时抛掷两枚硬币正面朝上的概率。然后,我们通过一个for循环来进行模拟,每次模拟都把两枚硬币抛掷n[i]次,然后计算同时正面朝上的概率。最后,我们使用plot函数来绘制结果图表。
图表的横轴是模拟次数,纵轴是同时抛掷两枚硬币正面朝上的概率。通过这个图表,我们可以看到,随着模拟次数的增加,我们得到的结果越来越接近于真实概率,即0.25。