<div dir="ltr">Pessoal, gostaria de fazer um gráfico conforme o anexo, cujo o intervalo de<br>confiança fique com formato de sombra dentro do gráfico. Conforme figura<br>anexada:  <br><div><br></div><div>URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190904/bfde02f0/attachment-0001.png" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190904/bfde02f0/attachment-0001.png</a>>  <br></div><div><br></div><div>Como tentativa usei o  *polygon*, mas não sei se existe um outro recurso<br>melhor. Se alguém souber e puder contribuir, a lista do R agradece. Segue o<br>código:<br><br> x <- seq(1:30)<br>Y <- x + rnorm(30, 0,3)<br>model <-  lm(Y ~ x)<br>pred <- data.frame(x1 = x)<br>pred$y1 <- predict(model, newdata = pred, interval = "confidence")<br>plot(consumo, type = "n",<br>     main = "Consumo médio, kg",<br>     xlab = "Consumo estimado, kg ",<br>     ylab = "Consumo observado, kg" )<br>icfit = pred$y1<br>yy = c(as.numeric(as.vector(icfit[,2])),<br>rev(as.numeric(as.vector(icfit[,3]))))<br>xx = c(x, rev(x))<br>polygon(xx, yy, col = "lightsteelblue2", border = NA)<br>points(consumo)<br>abline(model, col="blue", lwd= 2)<br>points(seq(2:15), seq(2:15), type = "l", col = "green", lwd = 2)<br>legend("topright", legend = c("modelo linear base y = x",<br>                              "modelo ajustado", "Bandas de confiança, 95%"<br>),       lty = c(1,1,1), col = c("green", "blue", "lightsteelblue2"),<br>       lwd = c(2,2,NA), bty = "n", pch = c(NA,NA, 15))  <br></div></div>