Boa tarde prezados, estou tentando fazer um gráfico de barras considerando 3 fatores, como p gráfico que o comando "bargraph.CI" permite fazer com dois fatores. Mas, não estou conseguindo. Alguém já fez este gráfico ou poderia me ajudar? Obrigada -- Natália Martins
Envie um código mínimo reproduzível para podermos nos orientar e lhe ajudar Em 27-01-2014 17:22, Natalia Martins escreveu:
Boa tarde prezados,
estou tentando fazer um gráfico de barras considerando 3 fatores,ie
como p gráfico que o comando "bargraph.CI" permite fazer com dois fatores.
Mas, não estou conseguindo. Alguém já fez este gráfico ou poderia me ajudar?
Obrigada
--
Natália Martins
_______________________________________________ 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.
Olá, Natalia! Veja se é algo nessa linha... ### <code r> require(psych) dd <- data.frame(ID=as.factor(rep(1:5,4)),A=sample(10:50,20),B=sample(10:50,20),C=sample(10:50,20)); dd error.bars.by(dd[2:4],dd[1], alpha=.05, bar=TRUE) ?error.bars.by ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 27 de janeiro de 2014 17:42, Fernando <nandodesouza@gmail.com> escreveu:
Envie um código mínimo reproduzível para podermos nos orientar e lhe ajudar Em 27-01-2014 17:22, Natalia Martins escreveu:
Boa tarde prezados,
estou tentando fazer um gráfico de barras considerando 3 fatores,ie
como p gráfico que o comando "bargraph.CI" permite fazer com dois fatores.
Mas, não estou conseguindo. Alguém já fez este gráfico ou poderia me ajudar?
Obrigada
--
Natália Martins
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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 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.
Um outro exemplo, dessa vez usando sciplot::bargraph.CI... ### <code r> dURL <- 'http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5d.zip' download.file(dURL, basename(dURL), mode='wb') unzip(basename(dURL),'lentils.csv') dat <- read.table('lentils.csv', head=T, sep=',') attach(dat) require(sciplot) bargraph.CI(x.factor=FARM, group=VARIETY, response=YIELD, legend=T, xlab="Location", ylab="Yield (kg/ha)", ylim=c(0,800), cex.names=1, cex.lab=1, x.leg=6, y.leg=700, cex.leg=1, leg.lab=c("Var A","Var B","Var C"), col=grey.colors(3), uc=T, lc=T, err.width=0.1, err.col="black", err.lty=1) ### </code> Fonte/Source: http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5D.pdf Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
O código que usei foi este: require(sciplot) bargraph.CI(pH, SHAM, group = time, data = sham, xlab = "Tratamento", ylab = "Sham", cex.lab = 1.2, x.leg = 1, col = c("gray18", "gray37", "gray60"), density = c(0,0.8), legend.text = NULL) legend("topright", c("2","4","8"), pch = c(15,15,15), col=c("gray18","gray37","gray60"), title="Time (h)") sendo igual ao que você enviou Éder. Porém eu ainda gostaria de acrescentar os niveis de mais um fator no eixo x. 2014-01-27 Éder Comunello <comunello.eder@gmail.com>
Um outro exemplo, dessa vez usando sciplot::bargraph.CI...
### <code r> dURL <- 'http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5d.zip' download.file(dURL, basename(dURL), mode='wb') unzip(basename(dURL),'lentils.csv') dat <- read.table('lentils.csv', head=T, sep=',') attach(dat)
require(sciplot) bargraph.CI(x.factor=FARM, group=VARIETY, response=YIELD, legend=T, xlab="Location", ylab="Yield (kg/ha)", ylim=c(0,800), cex.names=1, cex.lab=1, x.leg=6, y.leg=700, cex.leg=1, leg.lab=c("Var A","Var B","Var C"), col=grey.colors(3), uc=T, lc=T, err.width=0.1, err.col="black", err.lty=1) ### </code>
Fonte/Source: http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5D.pdf
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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.
-- Natália Martins Contato: (35) 91812482
Natalia, bom dia! Os códigos que postei são para postar três variáveis ou mais, justapostas, adicionando as barras de erros em cada uma delas. Poderia indicar o que você precisa diferenciar em relação aos gráficos gerados pela execução do código? Adicionalmente, sugiro encaminhar uma amostra do seu dataset pra rodar o código. Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Segue: time pH sinal SHAM 1 2 5.8 1 0.002 2 2 5.8 1 0.003 3 2 5.8 1 0.008 4 2 4.5 1 0.004 5 2 4.5 1 0.002 6 2 4.5 1 0.006 7 2 4.0 1 0.043 8 2 4.0 1 0.024 9 2 4.0 1 0.039 10 2 5.0 1 0.013 11 2 5.0 1 0.008 12 2 5.0 1 0.004 13 4 5.8 1 0.007 14 4 5.8 1 0.001 15 4 5.8 1 0.007 16 4 4.5 1 0.017 17 4 4.5 1 0.005 . . . 63 8 5.8 2 0.019 64 8 4.5 2 0.184 65 8 4.5 2 0.188 66 8 4.5 2 0.198 67 8 4.0 2 0.119 68 8 4.0 2 0.134 69 8 4.0 2 0.231 70 8 5.0 2 0.072 71 8 5.0 2 0.135 72 8 5.0 2 0.123 o objetivo é fazer um grafico de barras em que no eixo y esteja a var. sham, no x pH e sinal e "para a legenda" time. Em 28 de janeiro de 2014 10:33, Éder Comunello <comunello.eder@gmail.com>escreveu:
Natalia, bom dia!
Os códigos que postei são para postar três variáveis ou mais, justapostas, adicionando as barras de erros em cada uma delas. Poderia indicar o que você precisa diferenciar em relação aos gráficos gerados pela execução do código?
Adicionalmente, sugiro encaminhar uma amostra do seu dataset pra rodar o código.
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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.
-- Natália Martins Contato: (35) 91812482
Natalia, Uma tentativa... ### <code r> dat <- structure(list(time = c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L), pH = c(5.8, 5.8, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5, 5.8, 5.8, 5.8, 4.5, 4.5, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5), sinal = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), SHAM = c(0.002, 0.003, 0.008, 0.004, 0.002, 0.006, 0.043, 0.024, 0.039, 0.013, 0.008, 0.004, 0.007, 0.001, 0.007, 0.017, 0.005, 0.019, 0.184, 0.188, 0.198, 0.119, 0.134, 0.231, 0.072, 0.135, 0.123)), .Names = c("time", "pH", "sinal", "SHAM"), class = "data.frame", row.names = c(NA, -27L)) dat$new <- with(dat, paste0(format(dat$pH, dig=2), " (", sinal, ")")); dat require(sciplot) bargraph.CI(x.factor=new, response=SHAM, group = time, data = dat, xlab = "Tratamento", ylab = "Sham", cex.lab = 1.2, x.leg = 1, col = c("gray18", "gray37", "gray60"), density = c(0,0.8), legend.text = NULL) legend("topright", c("2","4","8"), pch = c(15,15,15), col=c("gray18","gray37","gray60"), title="Time (h)") ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 28 de janeiro de 2014 09:42, Natalia Martins <nsmbarreto@gmail.com>escreveu:
Segue:
time pH sinal SHAM 1 2 5.8 1 0.002 2 2 5.8 1 0.003 3 2 5.8 1 0.008 4 2 4.5 1 0.004 5 2 4.5 1 0.002 6 2 4.5 1 0.006 7 2 4.0 1 0.043 8 2 4.0 1 0.024 9 2 4.0 1 0.039 10 2 5.0 1 0.013 11 2 5.0 1 0.008 12 2 5.0 1 0.004 13 4 5.8 1 0.007 14 4 5.8 1 0.001 15 4 5.8 1 0.007 16 4 4.5 1 0.017 17 4 4.5 1 0.005 . . . 63 8 5.8 2 0.019 64 8 4.5 2 0.184 65 8 4.5 2 0.188 66 8 4.5 2 0.198 67 8 4.0 2 0.119 68 8 4.0 2 0.134 69 8 4.0 2 0.231 70 8 5.0 2 0.072 71 8 5.0 2 0.135 72 8 5.0 2 0.123
o objetivo é fazer um grafico de barras em que no eixo y esteja a var. sham, no x pH e sinal e "para a legenda" time.
Em 28 de janeiro de 2014 10:33, Éder Comunello <comunello.eder@gmail.com>escreveu:
Natalia, bom dia!
Os códigos que postei são para postar três variáveis ou mais, justapostas, adicionando as barras de erros em cada uma delas. Poderia indicar o que você precisa diferenciar em relação aos gráficos gerados pela execução do código?
Adicionalmente, sugiro encaminhar uma amostra do seu dataset pra rodar o código.
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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.
--
Natália Martins Contato: (35) 91812482
_______________________________________________ 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.
Perfeito, deu certinho. Muito obrigada!!! Em 28 de janeiro de 2014 11:27, Éder Comunello <comunello.eder@gmail.com>escreveu:
Natalia,
Uma tentativa...
### <code r> dat <- structure(list(time = c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L), pH = c(5.8, 5.8, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5, 5.8, 5.8, 5.8, 4.5, 4.5, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5), sinal = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), SHAM = c(0.002, 0.003, 0.008, 0.004, 0.002, 0.006, 0.043, 0.024, 0.039, 0.013, 0.008, 0.004, 0.007, 0.001, 0.007, 0.017, 0.005, 0.019, 0.184, 0.188, 0.198, 0.119, 0.134, 0.231, 0.072, 0.135, 0.123)), .Names = c("time", "pH", "sinal", "SHAM"), class = "data.frame", row.names = c(NA, -27L))
dat$new <- with(dat, paste0(format(dat$pH, dig=2), " (", sinal, ")")); dat
require(sciplot)
bargraph.CI(x.factor=new, response=SHAM, group = time, data = dat, xlab = "Tratamento", ylab = "Sham", cex.lab = 1.2, x.leg = 1, col = c("gray18", "gray37", "gray60"), density = c(0,0.8), legend.text = NULL)
legend("topright", c("2","4","8"), pch = c(15,15,15), col=c("gray18","gray37","gray60"), title="Time (h)")
### </code>
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Em 28 de janeiro de 2014 09:42, Natalia Martins <nsmbarreto@gmail.com>escreveu:
Segue:
time pH sinal SHAM 1 2 5.8 1 0.002 2 2 5.8 1 0.003 3 2 5.8 1 0.008 4 2 4.5 1 0.004 5 2 4.5 1 0.002 6 2 4.5 1 0.006 7 2 4.0 1 0.043 8 2 4.0 1 0.024 9 2 4.0 1 0.039 10 2 5.0 1 0.013 11 2 5.0 1 0.008 12 2 5.0 1 0.004 13 4 5.8 1 0.007 14 4 5.8 1 0.001 15 4 5.8 1 0.007 16 4 4.5 1 0.017 17 4 4.5 1 0.005 . . . 63 8 5.8 2 0.019 64 8 4.5 2 0.184 65 8 4.5 2 0.188 66 8 4.5 2 0.198 67 8 4.0 2 0.119 68 8 4.0 2 0.134 69 8 4.0 2 0.231 70 8 5.0 2 0.072 71 8 5.0 2 0.135 72 8 5.0 2 0.123
o objetivo é fazer um grafico de barras em que no eixo y esteja a var. sham, no x pH e sinal e "para a legenda" time.
Em 28 de janeiro de 2014 10:33, Éder Comunello <comunello.eder@gmail.com>escreveu:
Natalia, bom dia!
Os códigos que postei são para postar três variáveis ou mais, justapostas, adicionando as barras de erros em cada uma delas. Poderia indicar o que você precisa diferenciar em relação aos gráficos gerados pela execução do código?
Adicionalmente, sugiro encaminhar uma amostra do seu dataset pra rodar o código.
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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.
--
Natália Martins Contato: (35) 91812482
_______________________________________________ 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.
_______________________________________________ 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.
-- Natália Martins Contato: (35) 91812482
participantes (3)
-
Fernando -
Natalia Martins -
Éder Comunello