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}