Estou tentando realizar um cálculo e necessito da ajuda de vocês. Primeiro eu faço uma sequência com distribuição normal com 420 números, depois eu faço uma sequência de números que chamo de mci
com 42 números, depois o próximo cálculo tem uma condição if
.
Nesta parte eu preciso que na minha condição sempre compare a primeira lista de vetor do mci
com o primeiro número dos valores que gerei da minha distribuição normal, depois todo o processo se repita com o segundo numero até completar os 420.
Por exemplo: Para rp = (28, 27, 26...), bwg = (conta com mci[i]) > rp
Quero que esse rp
seja sempre o 28 até acabar a primeira sequencia de 42 numeros, e depois comparar com 27 da segunda sequencia do meu vetor.
Observe que a condição sempre compara sempre o mesmo valor de rp
, e acho que o R não está fazendo isso, comparando por exemplo o valor 26,912 com 0,435 onde ele deveria comparar com o 0,429
max <- 30
K <- 0.0118
Xm <- 21
SD <- 0.851636356306513
mean.b <- 28
rp <- rnorm(420, mean = mean.b, sd = SD)
BW <- 0.0223*rp^0.8944
mci <- seq(200, 270, 1.7)
bwg <- ifelse((max*(1-exp(-K*(MCi-(BW*Xm))))) > rp, yes = rp, no = (max*(1-exp(-K*(MCi-(BW*Xm))))))