r - 一个 R 程序在模拟一个公平的骰子后计算一个的数量

Question

我想模拟一个实验，你扔一个公平的骰子 100 次并计算数量。我想重复这个实验 10^5 次并保存结果。

这是我掷骰子的n代码

dice <- function(n) {
   sample(c(1:6),n,replace = TRUE)
}

x <- dice(100)

接下来，我想计算 1 的个数并模拟 10^5 次（这部分是错误的）：

x <- numeric(10^5)

for(n in 1:10^5){
  x[n] <- sum(dice(100)) 
} 

hist(x,
  main="100 Fair Rolls",
  xlab="Rolls",
  ylab="Probability",
  xlim=c(0,100),
  breaks=-1:1000+1,
  prob=TRUE)

标签： rstatisticsprobability-distribution

Answer 1

我认为你非常接近。如果您按如下方式更改 for 循环，它应该可以工作。

for(n in 1:10^5){

  x[n]<-sum(dice(100)==1) 

}