# Regression
x<-1:30
y<-1:30+rnorm(30,0,2)
modele<-lm(y~x)
val.resid(modele,reg=TRUE)
# Analyse de variance
f<-factor(rep(LETTERS[1:3],each=10))
modele2<-lm(y~f)
val.resid(modele2,reg=FALSE,fact="f")
Run the code above in your browser using DataLab