R: Condizioni Logiche Non Essere Rispettati

0

Domanda

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!

data-manipulation loops r while-loop
2021-11-24 02:45:09
1

Migliore risposta

4

Come circa questo modo? Come tag while-loopHo provato a utilizzarlo.

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
2021-11-24 02:52:21

In altre lingue

Questa pagina è in altre lingue

Русский
..................................................................................................................
Polski
..................................................................................................................
Română
..................................................................................................................
한국어
..................................................................................................................
हिन्दी
..................................................................................................................
Français
..................................................................................................................
Türk
..................................................................................................................
Česk
..................................................................................................................
Português
..................................................................................................................
ไทย
..................................................................................................................
中文
..................................................................................................................
Español
..................................................................................................................
Slovenský
..................................................................................................................