图形的蒙特卡罗模拟

Question

我正在写一个蒙特卡罗模拟来检查有多少次y不是紧挨着另一个y。我变出了一个40 x和10 y的向量，放置在向量中的随机位置。我的目标是计算向量中没有任何相邻y的概率。以下是我尝试过的：

nrep = 100000
count = 0
for (i in 1:nrep) {
  x = sample(c(rep('x', 40), c(rep('y', 10))))
  if (x[i]!=x[i+1] && x[i+1]!=x[i+2]) count = count + 1
}
print(count/nrep)

结果是一个非常小的数字，这对我来说似乎没有意义。

Answer 1

if部分不正确。我们可以使用head/end来检查连续的元素，并查看在一次迭代中是否有任何两个连续的'y'。

nrep = 100000
count = 0
set.seed(2020)

for (i in 1:nrep) {
  x = sample(rep(c('x', 'y'), c(40, 10)))
  if(any(head(x, -1) == 'y' & tail(x, -1) == 'y')) count = count + 1
}

count/nrep
#[1] 0.891