Felipe e Rodrigo, obrigado pela ajuda!
Com a sugestão do Rodrigo apresentava aquele erro mencionado, o que foi resolvido utilizando a sintaxe escrita pelo Felipe substituindo o "~" por "by"

aggregate(produtividade , by = list(cultivar), mean, data = exemplo)


Muito obrigado pessoal.

Em 25 de novembro de 2015 16:27, Felipe <felipe.e.barletta@gmail.com> escreveu:
João,

A variável ano está como factor?

ano <- as.factor(ano)
aggregate(produtividade , by = list(cultivar,ano), mean, data = exemplo)

--
Atenciosamente
Felipe E. Barletta Mendes
Estatístico - Conre3 9766-A
+55 (41)-92077191
+55 (41)-33287216

_______________________________________________
R-br mailing list
R-br@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 fornea cdigo mnimo reproduzvel.



--
João Rodrigo de Castro
Programa de Pós-Graduação em Meteorologia
Bolsista Laboratório de Agrometeorologia - Embrapa Clima Temperado
Universidade Federal de Pelotas