Julia, ja tentou aprimorar algo dessa forma:

modelo<-aov(Ciclo ~ Rep + Die*Id + Blo/Rep + Error(Blo/Die), data=dados)

Eu tentei aqui e saiu algum resultado parecido, mas falta alguma coisa....

Eu conheço essa tese do Iemma, estou fazendo algo como o contrário do que ele fez.

Augusto