Prezados,

Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados:

#ANOVA:
> mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod)
Warning message:
In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) :
  modelo Error() é singular

Error: PS[, 1]:PS[, 3]
          Df Sum Sq Mean Sq F value Pr(>F)  
PS[, 1]    1    8.1   8.100    11.3 0.0152 *
Residuals  6    4.3   0.717                 
---
Signif. codes:  
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: Within
                Df Sum Sq Mean Sq F value   Pr(>F)    
PS[, 2]          4  108.9   27.23   24.94 3.02e-08 ***
PS[, 1]:PS[, 2]  4  496.9  124.23  113.79 3.09e-15 ***
Residuals       24   26.2    1.09                     
---
Signif. codes:  
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

#Resíduos:

> residuals(mod)
NULL
#Teste de Shapiro-Wilk:
> res<-residuals(mod)
> shapiro.test(res)
Error: is.numeric(x) is not TRUE

Obrigada.
Att.,
Angélica.