G
eovane o ultimo valor do seu arquivo (63) só tem um representante, portanto, não é possível calcular o erro padrão ep=sd/sqrt(n). Tente utilizar o algoritmo abaixo:
dados=structure(list(idade = c(35L, 35L, 35L, 36L, 36L, 36L, 36L, 37L,
37L, 40L, 40L), peso = c(59L, 64L, 87L, 52L, 90L, 45L, 87L, 65L,
49L, 26L, 58L)), .Names = c("idade", "peso"), class = "data.frame", row.names = c(NA,
-11L))
attach(dados)
med=tapply(peso,idade,mean)
desv=tapply(peso,idade,sd)
nrep=tapply(idade,idade,length)
erro.padrao=numeric(0)
for(i in 1:length(summary(factor(idade)))){
erro.padrao<-c(erro.padrao,desv[[i]]/sqrt(nrep[[i]]))
}
n=barplot(med,ylim=c(0,100),ylab="Peso",xlab="Idade",col=3)
arrows(n,med-erro.padrao,n,med+erro.padrao,angle=90,code = 3, col = 9, length = 0.1)
Att.
Tiago.
#################################################################
Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES
Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas
#################################################################
Date: Mon, 17 Jun 2013 12:06:50 -0700
From: geovanecb@yahoo.com.br
To: r-br@listas.c3sl.ufpr.br
Subject: [R-br] erro bar
OLá Daniel tudo bem
eu tenho seguinte conjunto de dados, porém tenho mais de 10.000 dados, como poderia construir o erro bar ao longo dessas idades, pegando a média de peso por idade?
idade |
peso |
35 |
59 |
35 |
64 |
35 |
87 |
36 |
52 |
36 |
90 |
36 |
45 |
36 |
87 |
37 |
65 |
37 |
49 |
40 |
26 |
40 |
58 |
41 |
63 |
|
|
####################################
####################################
Prof. Geovane Carlos Barbosa
UCL - Faculdade do Centro Leste
###################################
###################################
De: Daniel Marcelino <dmarcelino@live.com>
Para: R-br Lista <r-br@listas.c3sl.ufpr.br>; geovane barbosa <geovanecb@yahoo.com.br>
Enviadas: Segunda-feira, 17 de Junho de 2013 15:20
Assunto: Re: [R-br] erro bar
Mande um exemplo minimamente reproduzível dos dados.
_______________________________________________
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.