[R-br] fatorial 5x3 + testemunha

Preciso rodar uma análise fatorial 5x3 + testemunha e estou usando i seguinte script segue em anexo o conjunto de dados para tornar o script reproduzível mamona=read.csv2(file.choose(),h=T) mamona.f1=TukeysHSD(aov(IVG~palha:test,data=mamona)) mamona.f2=TukeysHSD(aov(IVG~posicao:test,data=mamona)) mamona.f1.f2=TukeysHSD(aov(IVG~palha*posicao:test,data=mamona)) Vcs poderiam me informar se o erro está no script ou se meu conjunto de dados está organizado errado? Muito onrigado.   Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil