作业：使用 R 模拟抛硬币直到连续正面朝上

我是 R 新手，因此在这里询问（尚未找到非常有用的详细模拟教程。）

问题陈述是这样的

1. 模拟抛硬币20次并记录正面的数量 和最长的头。
2. 模拟抛硬币并记录所需的抛掷次数，直到依次（连续）出现 2、3、4 次正面（负二项式？）
3. 使用不同的种子进行 100 次运行，以找出记录的项目的分布。

如何用 R 编程语言解决这个问题？

对于 1，我执行了以下操作：

n=20                #no of trials
y=NULL              #initializing a vector of NULL values   
for (i in 1:n) {
    x=runif(1)      #random uniform
    if (x<0.5) {    #if condition for assigning heads / tails
        y[i]='H'
    } else {
       y[i]='T'
    }
}
y                   #print the vector storing the heads and tails.

对于 2，我理解这是负二项式的情况，因为它说“直到 2,3,4 头依次出现”。但我不确定如何继续在 R 中编写代码或提出逻辑。

对于 3，我猜测必须计算 100 次运行中获得的正面数量，即正面的运行长度。但这意味着什么不同的种子？我对此不确定。

note是的，这是一个家庭作业问题（我无法如此标记它）。我不需要代码，只需一些指针/更正/澄清/建议就会有帮助。

edit为了模拟 2，我在 MATLAB 上尝试了以下代码 -

head_count = 0;
flip_count = 0;

while (head_count < 3)                  #counting for 3 multiple heads
    flip = rand(1);
    flip_count = flip_count+1;
    if (flip > 0.5)
        head_count = head_count + 1;
    else 
         head_count = 0;       #invalidate head counts if not successive.
    end
end
disp(flip_count)

请告诉我这个逻辑是否正确。我可以使用 R 重新创建它。

由于这是家庭作业，我想为您提供一些关于#2 和#3 的指导。

对于#2，创建一个循环，不断抛硬币并计算连续出现正面的数量。每次抛掷时，抛掷次数加 1，当达到要求的正面数时，只需返回抛掷次数。

为了帮助您开始，这将执行 nbTosses tossesL

coin <- c('h','t')

ComputeNbTosses <- function(targetTosses) {
  nbtosses = 0;
  while (nbtosses < targetTosses) {
    nbtosses = nbtosses + 1
  }
  return(nbtosses);
} 

您只需修改它以接受参数 targetHeads 并仅在到达 tagetHeads 时退出 while 循环。

对于 #3，将 #2 放入函数中，然后使用set.seed(integer)，设置种子。计算机上的随机数并不是真正随机的。大多数随机生成器允许您设置种子。如果使用相同的种子，您将获得两次相同的伪随机数序列，这允许调试。

编辑：我将#2从递归更改为普通循环，为他编写代码会更容易。

编辑（提示）： 这将计算您连续有多少个正面，现在您必须在达到正面目标值时停止循环：

coin <- c('h','t')

ComputeNbTosses <- function(targetTosses) {
  nbTosses = 0;
  nbHeadsInRow = 0;
  allTosses = c();
  while (nbTosses < targetTosses) {
    toss = sample(coin,1,T)
    allTosses = c(allTosses, toss);
    if (toss == 'h') {
      nbHeadsInRow = nbHeadsInRow + 1;
    } else {
       nbHeadsInRow = 0;
    }
    nbTosses = nbTosses + 1
  }
  ret = list();
  ret$nbTosses = nbTosses;
  ret$nbHeadsInRowAtEnd = nbHeadsInRow;
  ret$allTosses = allTosses;
  return(ret);
} 

ComputeNbTosses(5);
ComputeNbTosses(8);

这将打印当它到达 nbTosses 时，您在行中有多少个头。这不是你想要的。你想要的是：

• 将参数从targetTosses to targetHeadsInRow
• 改变条件中的while以便它停止时nbHeadsInRow达到targetHeadsInRow.
• 你需要返回的是nbTosses（达到你的条件的抛掷次数）

如果有任何不清楚的地方请告诉我。