[R-br] erro bar
Tiago Souza Marçal
tiagosouzamarcal em hotmail.com
Segunda Junho 17 20:42:07 BRT 2013
Geovane 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 em yahoo.com.br
To: r-br em 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 em live.com>
Para: R-br Lista <r-br em listas.c3sl.ufpr.br>; geovane barbosa <geovanecb em 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.
Daniel
2013/6/17 geovane barbosa <geovanecb em yahoo.com.br>
Olá pessoal tudo bem
como faço para construir um erro bar por cada idade tendo uma variável y sendo por exemplo peso.
abraços
De: Marcos Silva <marcosfs2006 em gmail.com>
Para: r-br em listas.c3sl.ufpr.br
Enviadas: Segunda-feira, 17 de Junho de 2013 10:39
Assunto: Re: [R-br]
Dúvida - gdata (notação científica)
Veja se algo como options(scipen=30) funciona.
Abs.
Em 17 de junho de 2013 09:15, Sérgio Henrique almeida da silva ju <sergio.edfisica em gmail.com> escreveu:
Amigos
Estou lendo um banco do excel usando o gdata
require(gdata)
dados <- read.xls(...)
Tenho uma variável com dados numéricos de tamanho maior que 30. Ex:
00000000000000000000000000000083990244081296600001
Quando leio o banco ela fica em notação científica (8.40e+19), porém preciso trabalhar com ela no formato original, qual comando eu posso usar?
--
Sérgio Henrique Almeida da Silva Junior
Doutorando em Epidemiologia em Saúde Pública
Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ
http://lattes.cnpq.br/1611345552843383
Tel: (21) 68463637
http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro
_______________________________________________
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.
--
Marcos F. Silva
http://sites.google.com/site/marcosfs2006
_______________________________________________
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.
_______________________________________________
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.
--
\begin{signature}Daniel Marcelino
☁ dm.silva em umontreal.ca
☎ (514) 343 6111 #3799Skype: d.marcelino
✎ 3200 Jean Brillant, Office C5071Montreal, QC; H3T 1N8
Canada\end{signature}
"Small steps toward a much better world"
_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130618/04ac3ad6/attachment-0001.html>
Mais detalhes sobre a lista de discussão R-br