<div dir="ltr">Não é mais fácil usar aggregate? <div><br></div><div><div>d<-data.frame(trat,id,y) # data.frame() é melhor do que cbind()</div><div>aggregate(data.frame(dados=d$y), list(trat=d$trat, id=d$id), function(x) sd(x)/sqrt(length(x)))</div>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/11/12 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caros membros,<br>
<br>
Tenho variáveis respostas com diferentes tamanhos de vetores e na hora de calcular o erro padrão da média usando o with() e tapply(), em função dos diferentes comprimentos de vetores , não estou conseguindo calcular, alguém poderia me ajudar, segue CRM:<br>
<br>
y1<-rnorm(100,5)### Variável resposta<br>
y2<-rnorm(80,20)<br>
y3<-rnorm(100,15)<br>
y<-c(y1,y2,y3)<br>
<br>
<br>
trat1<-rep(1,100)### Variável explicativa tratamento<br>
trat2<-rep(2,80)<br>
trat3<-rep(3,100)<br>
trat<-c(trat1,trat2,trat3)<br>
<br>
id1<-rep(60,100)### Variável explicativa idade<br>
id2<-rep(84,180)<br>
id<-c(id1,id2)<br>
<br>
d<-cbind(trat,id,y)<br>
<br>
ep.b3<-with(d, tapply(y, list(trat,id), sd, na.rm = T)) /sqrt(with(d, tapply(y, list(trat,id), length, na.rm = T))) ### Calculando o erro padrão<br>
ep.b3<br>
<br>
-- <br>
==============================<u></u>==============================<u></u>==========<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@<u></u>yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.<u></u>br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/<u></u>1360403201088680</a><br>
==============================<u></u>==============================<u></u>==========<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>