R:逻辑条件未得到尊重

R: Logical Conditions Not Being Respected

提问人:stats_noob 提问时间:11/24/2021 更新时间:11/24/2021 访问量:46

问:

我正在使用 R 编程语言。我正在尝试构建一个执行以下操作的循环:

  • 第 1 步:继续生成两个随机数“a”和“b”,直到“a”和“b”都大于 12

  • 第 2 步:跟踪在完成第 1 步之前必须生成多少个随机数

  • 第 3 步:重复第 1 步和第 2 步 100 次

由于我不知道如何继续生成随机数直到满足条件,因此我尝试生成大量随机数,希望满足该条件(可能有更好的写法):

results <- list()


for (i in 1:100){
  
  # do until break
  repeat {
    
    # repeat many random numbers
    a = rnorm(10000,10,1)
    b = rnorm(10000,10,1)
    
    # does any pair meet the requirement
    if (any(a > 12 & b > 12)) {
      
      # put it in a data.frame
      d_i = data.frame(a,b)
      
      # end repeat
      break
    }
  }
  
  # select all rows until the first time the requirement is met
  # it must be met, otherwise the loop would not have ended
  d_i <- d_i[1:which(d_i$a > 10 & d_i$b > 10)[1], ]
  
  # prep other variables and only keep last row (i.e. the row where the condition was met)
  d_i$index = seq_len(nrow(d_i))
  d_i$iteration = as.factor(i)
e_i = d_i[nrow(d_i),]
  
  results[[i]] <- e_i
  
}

results_df <- do.call(rbind.data.frame, results)

问题:当我查看结果时,我注意到循环错误地考虑了要满足的条件,例如:

head(results_df)

          a        b index iteration
4  10.29053 10.56263     4         1
5  10.95308 10.32236     5         2
3  10.74808 10.50135     3         3
13 11.87705 10.75067    13         4
1  10.17850 10.58678     1         5
14 10.14741 11.07238     1         6

例如,在每一行中,“a”和“b”都小于 12。

有谁知道为什么会发生这种情况,有人可以告诉我如何解决这个问题吗?

谢谢!

R 循环 while 循环 数据操作

评论

0赞 Colombo 11/24/2021
你的第一个条件是,但在第二个条件上,你测试是有原因的吗?>12>10
0赞 stats_noob 11/24/2021
@科伦坡:谢谢你的回复!我想你可能已经发现了我问题中的错误!谢谢!

答:

4赞 Park 11/24/2021 #1

这样怎么样?当你标记时,我尝试使用它。while-loop

res <- matrix(0, nrow = 0, ncol = 3)    

for (j in 1:100){
  a <- rnorm(1, 10, 1)
  b <- rnorm(1, 10, 1)
  i <- 1
  while(a < 12 | b < 12) {
    a <- rnorm(1, 10, 1)
    b <- rnorm(1, 10, 1)
    i <- i + 1
  }
  x <- c(a,b,i)
  res <- rbind(res, x)
}

head(res)
      [,1]     [,2] [,3]
x 12.14232 12.08977  399
x 12.27158 12.01319 1695
x 12.57345 12.42135  302
x 12.07494 12.64841  600
x 12.03210 12.07949   82
x 12.34006 12.00365  782

dim(res)
[1] 100   3