[R-br] gráficos com bandas de confiança no R - com sombra

[Nicolay Cunha]

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) <
r-br em listas.c3sl.ufpr.br> 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
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190905/43f26340/attachment.html>