Muchas gracias Manuel. Creo que ya he solucionado el problema. No sé si hay un camino más corto; pero como necesitaba utilizar solo la suma de cuadrados residual he añadido un par de líneas. Adjunto mi solución:
AnovaVariableX <- aov(AnteriorX ~ Grupo, data=Anacova)
resultanova<-summary(AnovaVariableX)
> str(resultanova)
List of 1
$ :Classes ‘anova’ and 'data.frame': 2 obs. of 5 variables:
..$ Df : num [1:2] 2 21
..$ Sum Sq : num [1:2] 36.7 56.9
..$ Mean Sq: num [1:2] 18.37 2.71
..$ F value: num [1:2] 6.78 NA
..$ Pr(>F) : num [1:2] 0.00533 NA
- attr(*, "class")= chr [1:2] "summary.aov" "listof"
sumacuadrados<-resultanova[[1]]['Sum Sq']
> sumacuadrados$`Sum Sq`[2]
[1] 56.875
Muchas gracias, de nuevo