aggregate(), não seria o caso?
sapply(OCTOPUS_S_2003[, -(1:3)],
tapply(OCTOPUS_S_2003[,2:3], mean, na.rm=TRUE))
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@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br