Plot (x,y) com desvios erro padrão da média como linha vertical
Boa noite pessoal, Gostaria de plotar um gráfico com pontos que representam as médias e que o erro padrão da média previamente fornecido aparecesse sobre os pontos, sendo meu CRM: ##Gráfico casal10<-c(16.03636,16.14,16.62,23.05)##Média 10 casais sd.casal10<-(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') # Eu pergunto: tem como sd.casal10 e sd.casal1 aparecerem como uma barra vertical ("arrows") sobre os pontos que representam as média? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
Você precisa arruma o seu exemplo, pois ele não reproduz direito. Se eu entendi corretamente, você quer algo do tipo: plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90) Agora, porque você não faz tudo direto ao invés de re-entrar as médias e desvios? Digamos que o seus dados brutos são iris data(iris) mu <- tapply(Petal.Length,Species,mean) ci95 <- function (x) { return(1.96*sqrt(var(x)/length(x))) } se <- tapply(Petal.Length,Species, ci95) plot(1:3,mu,xlim=c(0,4),ylim=c(1,7),pch=20, las=1, bty="L", xlab="Species", ylab="Petal Length (cm)") arrows( x0=1:3, x1=1:3, y0=mu+se, y1=mu-se, code=3, length=0.1, angle=90) 2013/1/10 ASANTOS <alexandresantosbr@yahoo.com.br>:
Boa noite pessoal,
Gostaria de plotar um gráfico com pontos que representam as médias e que o erro padrão da média previamente fornecido aparecesse sobre os pontos, sendo meu CRM:
##Gráfico casal10<-c(16.03636,16.14,16.62,23.05)##Média 10 casais sd.casal10<-(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') # Eu pergunto: tem como sd.casal10 e sd.casal1 aparecerem como uma barra vertical ("arrows") sobre os pontos que representam as média?
Obrigado
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ 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.
-- "Small steps toward a much better world" \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (514) 343 6111 #3799 ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature}
Bom Dia Daniel, Ñão usei os dados brutos pois fiz uma transformação com os dados para análise. O problema com meu CRM foi a falta de concatenar (c) os erros padrões ficando: casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, + code=3, length=0.1, angle=90) Mesmo assim os arrows não ficaram corretamente sobre os pontos, Obrigado, Alexandre Em 10/01/2013 23:49, Daniel Marcelino escreveu:
plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
No arrows() o teu x0 e x1 nao deveria ser ger ao invés de 1:12? casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal *sd.casal10*<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal (eu troquei pra casal10 para o codigo funcionar) ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(*x0=ger, x1=ger*, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90) 2013/1/11 ASANTOS <alexandresantosbr@yahoo.com.br>
Bom Dia Daniel, Ñão usei os dados brutos pois fiz uma transformação com os dados para análise. O problema com meu CRM foi a falta de concatenar (c) os erros padrões ficando: casal10<-c(0.8402676,0.**69461420,0.5656461,1.479279)#**SD para 10 casais casal1<-c(14.14000,14.21,17.**37,17.26667)## Média 1casal sd.casal1<-c(0.1166190,0.**07520343,1.1059486,0.768739)#**SD para 1 casal
ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),**pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, + code=3, length=0.1, angle=90)
Mesmo assim os arrows não ficaram corretamente sobre os pontos,
Obrigado,
Alexandre
Em 10/01/2013 23:49, Daniel Marcelino escreveu:
plot(ger,casal10,ylim=c(0,25),**pch=16)
lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
-- ==============================**==============================**========== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@**yahoo.com.br<e-mails%3Aalexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.**br<alexandre.santos@cas.ifmt.edu.br> ==============================**==============================**==========
______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.
Sim, pode ser. Você conseguiu plotar afinal? Daniel 2013/1/11 Rodrigo Coster <rcoster@gmail.com>:
No arrows() o teu x0 e x1 nao deveria ser ger ao invés de 1:12?
casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal10<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal (eu troquei pra casal10 para o codigo funcionar)
ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=ger, x1=ger, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
2013/1/11 ASANTOS <alexandresantosbr@yahoo.com.br>
Bom Dia Daniel, Ñão usei os dados brutos pois fiz uma transformação com os dados para análise. O problema com meu CRM foi a falta de concatenar (c) os erros padrões ficando: casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal
ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, + code=3, length=0.1, angle=90)
Mesmo assim os arrows não ficaram corretamente sobre os pontos,
Obrigado,
Alexandre
Em 10/01/2013 23:49, Daniel Marcelino escreveu:
plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ 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.
-- "Small steps toward a much better world" \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (514) 343 6111 #3799 ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature}
Obrigado Daniel e Rodrigo, problema resolvido, segue CRM abaixo: casal10<-c(16.03636,16.14,16.62,23.05)##Média 10 casais sd.casal10<-(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=ger, x1=ger, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90) arrows(x0=ger, x1=ger, y0=casal1+sd.casal1, y1=casal1-sd.casal1, code=3, length=0.1, angle=90, col='gray') # Alexandre Em 11/01/2013 20:12, Daniel Marcelino escreveu:
Sim, pode ser. Você conseguiu plotar afinal?
Daniel
2013/1/11 Rodrigo Coster <rcoster@gmail.com>:
No arrows() o teu x0 e x1 nao deveria ser ger ao invés de 1:12?
casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal10<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal (eu troquei pra casal10 para o codigo funcionar)
ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=ger, x1=ger, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
2013/1/11 ASANTOS <alexandresantosbr@yahoo.com.br>
Bom Dia Daniel, Ñão usei os dados brutos pois fiz uma transformação com os dados para análise. O problema com meu CRM foi a falta de concatenar (c) os erros padrões ficando: casal10<-c(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-c(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal
ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, + code=3, length=0.1, angle=90)
Mesmo assim os arrows não ficaram corretamente sobre os pontos,
Obrigado,
Alexandre
Em 10/01/2013 23:49, Daniel Marcelino escreveu:
plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') arrows(x0=1:12, x1=1:12, y0=casal10+sd.casal10, y1=casal10-sd.casal10, code=3, length=0.1, angle=90)
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ 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.
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
Embora chegando um pouco tarde nesta 'trédi', posso perguntar porque você não usa as funções do pacote gplots plotCI() que faz, a meu ver exatamente o que você quer, e a plotmeans() que eu acho melhor porque a partir dos vetores de dados ela calcula as médias, erros padrão e plota? -- Cesar Rabak 2013/1/10 ASANTOS <alexandresantosbr@yahoo.com.br>:
Boa noite pessoal,
Gostaria de plotar um gráfico com pontos que representam as médias e que o erro padrão da média previamente fornecido aparecesse sobre os pontos, sendo meu CRM:
##Gráfico casal10<-c(16.03636,16.14,16.62,23.05)##Média 10 casais sd.casal10<-(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') # Eu pergunto: tem como sd.casal10 e sd.casal1 aparecerem como uma barra vertical ("arrows") sobre os pontos que representam as média?
Obrigado
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ 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)
-
ASANTOS -
Cesar Rabak -
Daniel Marcelino -
Rodrigo Coster