Boa tarde, Cezar Rabak.

Entendi melhor, porém ainda não por completo.

Mas, o meu problema agora é que não consigo mais obter o mesmo resultado. Estou importando exatamente os mesmos dados e utilizando exatamente os mesmos comandos, porém não consigo mais o mesmo resultado que enviei no fórum dia 04/09/16. Quando realizo a ANOVA o resultado gerado está sendo diferente (e está errado) daquele que apresentei anteriormente (este está correto). Gostaria de saber como posso obter o "resíduo" para inserir no comando para o teste de Shapiro Willk: "shapiro.test(resíduo)".

 Seguem abaixo os comandos utilizados (acho que é o 'CMR' a que você se referiu) e o dados, alguém poderia me ajudar a encontrar o meu erro e mostrar uma solução? Já procurei muito, mas não encontro.
_____________________________________________________________
> PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T);PS
   PARC SubP REP var
1     C   30   1  12
2     C   30   2  13
3     C   30   3  14
4     C   30   4  12
5     C   45   1  18
6     C   45   2  17
7     C   45   3  16
8     C   45   4  17
9     C   60   1  22
10    C   60   2  21
11    C   60   3  24
12    C   60   4  23
13    C   75   1  24
14    C   75   2  23
15    C   75   3  22
16    C   75   4  24
17    C   90   1  21
18    C   90   2  20
19    C   90   3  21
20    C   90   4  19
21    S   30   1  22
22    S   30   2  21
23    S   30   3  24
24    S   30   4  23
25    S   45   1  21
26    S   45   2  22
27    S   45   3  23
28    S   45   4  21
29    S   60   1  20
30    S   60   2  18
31    S   60   3  19
32    S   60   4  20
33    S   75   1  17
34    S   75   2  15
35    S   75   3  16
36    S   75   4  15
37    S   90   1  12
38    S   90   2  13
39    S   90   3  12
40    S   90   4  11
> attach(PS)

> SubP<-as.factor(SubP);is.factor(SubP)
[1] TRUE

> mod<-aov(var~PARC*SubP+Error(PARC:REP));summary(mod)

Error: PARC:REP
          Df Sum Sq Mean Sq F value Pr(>F) 
PARC       1  6.453   6.453   80.67 0.0706 .
Residuals  1  0.080   0.080                
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: Within
          Df Sum Sq Mean Sq F value   Pr(>F)   
PARC       1    1.7    1.67   1.535    0.226   
SubP       4  108.9   27.23  25.076 6.64e-09 ***
PARC:SubP  4  496.9  124.22 114.418  < 2e-16 ***
Residuals 28   30.4    1.09                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> resíduo<-residuals(mod);resíduo
NULL
___________________________________________________________________
Aguardo retorno.

Obrigada.

Att.,
Angélica.