Eu estou trabalhando com a linguagem de programação R. Eu estou tentando construir um laço que executa o seguinte :
Passo 1: Manter a geração de dois números aleatórios "a" e "b", até que tanto "a" e "b" são maiores do que 12
Passo 2: o controle de quantos números aleatórios tinha que ser gerado até que demorou para o Passo 1 para ser concluída
Passo 3: Repita o Passo 1 e o Passo 2 100 vezes
Desde que eu não sei como manter a geração de números aleatórios até que uma condição seja satisfeita, eu tentei gerar uma grande quantidade de números aleatórios na esperança de que a condição for satisfeita (não é provavelmente a melhor maneira de escrever isso):
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)
Problema: Quando eu olhar para os resultados, notei que o loop é incorretamente considerando-se a condição a ser satisfeita, por exemplo:
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
Por exemplo, em cada uma destas linhas - as "a" e "b" são menores do que 12.
Alguém sabe porque isso está acontecendo e alguém por favor pode me mostrar como corrigir este problema?
Obrigado!