gráficos com bandas de confiança no R - com sombra
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/attac...
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))
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_enAR816A... 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@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/attac...
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@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.
Considere esses exemplos lattice: http://leg.ufpr.br/~walmes/cursoR/mgest/1medidas-influen.html ggplot2: http://leg.ufpr.br/~walmes/analises/CECarducci/cafe_pedotrans.html Usando a graphics com *código reproduzível* m0 <- lm(dist ~ poly(speed, degree = 2), data = cars) grid <- with(cars, data.frame(speed = seq(min(speed), max(speed), length.out = 31))) mat <- predict(m0, newdata = grid, interval = "confidence") grid <- cbind(grid, as.data.frame(mat)) str(grid) plot(dist ~ speed, data = cars, type = "n") with(grid, polygon(x = c(speed, rev(speed)), y = c(lwr, rev(upr)), col = "gray", border = NA)) lines(fit ~ speed, data = grid) points(dist ~ speed, data = cars) À disposição. Walmes.
Tenta Hmisc::rm.boot Pedro Brasil Em sex, 6 de set de 2019 às 09:57, Walmes Zeviani por (R-br) < r-br@listas.c3sl.ufpr.br> escreveu:
Considere esses exemplos
lattice: http://leg.ufpr.br/~walmes/cursoR/mgest/1medidas-influen.html ggplot2: http://leg.ufpr.br/~walmes/analises/CECarducci/cafe_pedotrans.html
Usando a graphics com *código reproduzível*
m0 <- lm(dist ~ poly(speed, degree = 2), data = cars)
grid <- with(cars, data.frame(speed = seq(min(speed), max(speed), length.out = 31)))
mat <- predict(m0, newdata = grid, interval = "confidence")
grid <- cbind(grid, as.data.frame(mat)) str(grid)
plot(dist ~ speed, data = cars, type = "n") with(grid, polygon(x = c(speed, rev(speed)), y = c(lwr, rev(upr)), col = "gray", border = NA)) lines(fit ~ speed, data = grid) points(dist ~ speed, data = cars)
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@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.
participantes (4)
-
Nicolay Cunha -
Nilson Guiotoku -
Pedro Emmanuel Alvarenga Americano do Brasil -
Walmes Zeviani