蒙特卡洛连续3张或以上的面孔

如何解决蒙特卡洛连续3张或以上的面孔

我编写了这段代码，以100000次迭代的模拟检查3个或更多连续的面，并进行了5次模切。我认为这是正确的方向，但我缺少一些东西。我不断收到缺失值错误：

nrep = 500000
count = 0
for (i in 1:nrep) {
  roll = sample(6,5)
  print(roll)
  if (roll[i] == roll[i+1] & roll[i+1] == roll[i+2]) count = count + 1
}
print(count)

请仅建议使用基数R进行更正。

解决方法

在我的评论中，您可以使用函数rle()计算向量中等值游程的长度和值。您可以执行以下操作

nrep = 500000
count = 0
for (i in 1:nrep) {
  roll = sample(6,5,replace = TRUE)
  roll_rle = rle(roll)
  if (any(roll_rle$lengths >= 3)) {
    print(roll)
    count = count + 1
  }
}