um exemplo reproduzivel... set.seed(1) da = data.frame(Trat=sample(6, 1000, rep=T), x=rnorm(1000), y=rexp(1000, 1)) aggregate(da[,-1], by=list(Trat=da[['Trat']]), function(x) c(media=mean(x, na.rm=T), stderr=sd(x, na.rm=T)))
b _______________________________________________ R-br mailing list [email protected] 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.
