<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Obrigado Prof. Paulo e Rodrigo,<br>
<br>
Problema resolvido,<br>
<br>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
<br>
<br>
<div class="moz-cite-prefix">Em 12/11/2013 10:28, Paulo Justiniano
escreveu:<br>
</div>
<blockquote
cite="mid:alpine.DEB.2.00.1311121114430.6890@pataxo.est.ufpr.br"
type="cite">Naso era para funcionar mesmo o inicialmente tentando
<br>
<br>
Acrescentando detalhes:
<br>
seus objetos originais nao possuem os nomes de coluna que voce
usou no tapply()
<br>
<br>
procure organizar dados em data-frames e certifique-se que os
nomes de colunas estao definidos
<br>
<br>
d<-data.frame(trat=trat,id=id,y=y)
<br>
names(d)
<br>
<br>
<br>
Outra coisa evite criar tantos objetos de depois juntar no
data.frame()
<br>
para nao criar confusao de nomes e objetos no seu workspace
<br>
<br>
- crie direto no data-frame
<br>
<br>
d <- data.frame(
<br>
y = c(rnorm(100,5), rnorm(80,20), rnorm(100,15)),
<br>
trat = c(rep(1,100), rep(2,80), rep(3,100)),
<br>
id = c(rep(60,100), rep(84,180))
<br>
)
<br>
<br>
revise sua chamada de funcoes, na.rm nao é um argumento de
length()
<br>
<br>
talvez vc queira isto:
<br>
ep.b3<-with(d, tapply(y, list(trat,id), sd, na.rm = T))/with(d,
sqrt(tapply(y, list(trat,id), length)))
<br>
<br>
ou de forma mais direta:
<br>
<br>
ep.b3<-with(d, tapply(y, list(trat,id), function(x) sd(x, na.rm
= T)/length(x)))
<br>
<br>
<br>
<br>
ou
<br>
- apague os objetos nao necessários apos combina-los no data.frame
<br>
<br>
On Tue, 12 Nov 2013, Rodrigo Coster wrote:
<br>
<br>
<blockquote type="cite">Não é mais fácil usar aggregate?
<br>
<br>
<br>
2013/11/12 ASANTOS <a class="moz-txt-link-rfc2396E" href="mailto:alexandresantosbr@yahoo.com.br"><alexandresantosbr@yahoo.com.br></a>
<br>
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
<br>
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
<br>
erro padrão
<br>
ep.b3
<br>
<br>
--
<br>
======================================================================
<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: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<br>
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<br>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
<br>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
<br>
======================================================================
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível.
<br>
<br>
<br>
<br>
<br>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>