Rosiana,

O script abaixo é da análise de um experimento DQL com 4 amostras por parcela. Adapte pro seu caso.

http://www.leg.ufpr.br/~walmes/ensino/af722-2014-01/af722-2014-01-05dqlrept.R

Se eu entendi bem, o seu modelo seria

y_{ijk} = trat_i+parcela_{i/j}+peixe_{i/j/k}

em que trat é teu fator experimental, parcela é o efeito aleatório da j parcela dentro do i trat e k peixe dentro da j parcela dentro do i trat.

À disposição.
Walmes.​