掷硬币蒙特卡罗模拟精确正面

Question

我正在尝试编写一个用于不公平抛硬币的函数，它有四个输入：正面朝上的重量/概率、模拟次数、每个模拟的翻转次数以及我想要找到正面的确切数量获得的概率。

下面的函数返回我不理解的翻转部分的错误：

如果 (percent_heads == nheads) 计数器 1 并且只使用第一个元素 50

counter = 0
percent_heads = NULL
flip_function <- function(p, nheads, nflips, nsim){
  for(i in 1:nsim){
    flips <- sample(c("H","T"), size = nflips, replace = TRUE, prob = c(p, 1-p))
    percent_heads[i] <- length(which(flips == "H")) / nflips
    if(percent_heads == nheads) counter <- counter + 1
    }
  return(counter/nsim)
}

Answer 1

我觉得我是个白痴。经过一番深思熟虑，我相信我得到了它的工作。请比我更聪明的人验证这是否按预期工作：

counter = 0
total_heads = NULL
flip_function <- function(p, nheads, nflips, nsim){
for(i in 1:nsim){
flips <- sample(c("H","T"), size = nflips, replace = TRUE, prob = c(p, 1-p))
total_heads[i] <- length(which(flips == "H"))
if(total_heads[i] == nheads) counter <- counter + 1
}
return(counter/nsim)
}