Monte Carlo 连续 3 张或更多面孔

Question

我编写此代码是为了在 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 的更正提出建议。

Answer 1

添加到我的评论中，您可以使用函数 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
  }
}