Bom dia.

Eu tenho um experimento em látice 8x8 com duas repetições. Fiz a análise de variância utilizando a função "PBIB" - Analysis of the Partially Balanced Incomplete Block Design. 

Script:

exp01.gy<-PBIB.test(bloco,gen,rep,GY,k=8,group=TRUE, console=TRUE, method="REML")

Minha dúvida é como fazer para estimar os resíduos desse modelo, uma vez que não estou conseguindo com a função "residuals". O que você me sugerem fazer?

Desde já agradeço.

Rodrigo.