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