search for: xbval

...<- sort(y) # Ordene todas as médias geradas LI <- y[((1-nc)/2)*rep] # Limite inferior, LS <- y[rep-(((1-nc)/2)*rep)] # Limite superior ### PARTE GRÁFICA ### x <- mean(y) xvals <- seq(-LI, LI, length.out=5000) dvals <- dnorm(xvals,mean(y), sd(y))[1:5000] xbvals <- seq(LS, LS*2, length.out=5000) dbvals <- dnorm(xbvals,mean(y), sd(y))[1:5000] ahist <- hist(y, freq=FALSE, col="lightblue", main="Intervalo de confiança") polygon(c(xvals,LI,LI), c(dvals,dnorm(LI,mean(y), sd(y)),min(dvals)), col="orange", lty=0) polygon(...