Sto lavorando con il linguaggio di programmazione R. Sto cercando di creare un ciclo che esegue le seguenti operazioni :
Fase 1: Tenere la generazione di due numeri casuali "a" e "b" fino a quando entrambe le "a" e "b" sono maggiori di 12
Passo 2: tenere Traccia del numero di numeri casuali doveva essere generato fino a quando ha preso per il Passaggio 1 per essere completato
Passo 3: Ripetere i passaggi 1 e 2 fino a 100 volte
Dato che non so come mantenere la generazione di numeri casuali fino a quando una condizione è soddisfatta, ho provato a generare una grande quantità di numeri casuali sperando che la condizione è soddisfatta (probabilmente c'è un modo migliore per scrivere questo):
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 guardo i risultati, ho notato che il ciclo è errato considerare la condizione da soddisfare, per esempio:
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
Per esempio, in ognuna di queste righe - sia "a" e "b" sono inferiori a 12.
Qualcuno sa perché questo sta accadendo e per favore qualcuno può mostrarmi come risolvere questo problema?
Grazie!