[R-br] Plot (x,y) com desvios erro padrão da média como linha vertical

[Daniel Marcelino]

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 em 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 em yahoo.com.br
>         alexandre.santos em cas.ifmt.edu.br
> ======================================================================
>
>
>
> _______________________________________________
> 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.



-- 
"Small steps toward a much better world"

\begin{signature}
Daniel Marcelino
☁ dm.silva em umontreal.ca
☎ (514) 343 6111 #3799
✎ 3200 Jean Brillant,  Office C5071
Montreal, QC; H3T 1N8
Canada
\end{signature}