f1<-c(rep("A",20),rep("B",20))

f2<-rep(c(gl(4,5)),2)

res1<-c(0,0, 20,  0, 20, 20,  0,  0, 20, 40, 30, 40, 40, 20, 40, 30, 20, 30, 30, 50, 40, 40, 40, 70, 60, 40, 40, 40, 60, 20, 30, 30, 50, 40, 40, 50, 30, 40, 30, 30)

res2<-c(0,0,1,0,1,1,0,0,1,1.75,1.333333,1,1.5,1,1.75,1.333333,1.5,1.666667,1.333333,1.666667,1,1,1.25,1,1.166667,1.25,1.25,1.25,1.333333,1,2.333333,1,1.2,1,1,2.2,1.333333,1.5,

f1<-factor(f1)

f2<-factor(f2)

y<-cbind(res1,res2)

require(mvoutlier)

aq.plot(y, alpha=0.05) # O que devo fazer com os outlines multivariados?

mav<-manova(y~f1*f2)

summary(mav)

# Como testar os contrastes multivariados de interesse no R?

Ex: Nível A de f1 vs Nível B de f1 dentro do nível um de f2

# Com fazer a análise Stepdown no R?