1
I am interested in applying this model described in the variable d for each of the 1024 elements described in my second for and finally store the percentiles of each of these distributions in a vector I called vetorPercent.
I believe that with this code I am only generating a distribution d and not a distribution d[i] for each element i of my sample:
precomed = vector(mode = "numeric", length = 1581L)
for (x in 1:1581) {
precomed[x] = (dados1[x] + dados2[x])/2
}
vetorPercent <- list()
for (i in 1:1024) {
d <- dados3[i]*rnorm(1024, precomed[i],
sqrt(var(dados3)))
vetorPercent[[i]] <- quantile(d, c(.05, .95))
}