Oi Nilson, Infelizmente seu CMR ainda não é totalmente reproduzível. Note que você chama plot(consumo...), mas não há nenhum objeto "consumo" até então. Substituindo "consumo" por "x" e "Y" o código funciona.
Sugiro que você tente fazer o que quer no ggplot e use geom_ribbon para isso. Há muitos exemplos disponíveis, creio que algum servirá para o seu propósito. https://www.google.com/search?q=geom_ribbon&newwindow=1&rlz=1C5CHFA_enAR816AR817&sxsrf=ACYBGNRFJCCuwK4wOGpQ_UvWpF0fLQ6XWw:1567719453513&source=lnms&tbm=isch&sa=X&ved=0ahUKEwjhhomJ0rrkAhX2CrkGHQL2B9gQ_AUIEigB&biw=1432&bih=784&dpr=2 Um abraço Nicolay ----------------------------------------------------------------------- Grupo de Ecología de la Polinización (https://sites.google.com/view/ecopol/home) INIBIOMA, CONICET-Universidad Nacional del Comahue Quintral 1250 8400 San Carlos de Bariloche Rio Negro, Argentina ----------------------------------------------------------------------- On Thu, Sep 5, 2019 at 6:12 PM Nilson Guiotoku por (R-br) < [email protected]> wrote: > 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: > > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190904/bfde02f0/attachment-0001.png > > > > 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)) > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
