Pois bem Walmes, o que eu quero saber na verdade é o seguinte:

Dá pra fazer:

m0 <- lm(y~ 1a+2b+3c+4ab+5ac+6bc) 
anova(m0) 

?

Eu acho que realmenteo não consigo explicar o problema. É que eu ja tenho os coeficientes. Acho que deve ser até um erro absurdo esta minha pergunta. Por que até onde eu sei, o correto seria pedir algo do tipo:

m0 <- lm(y~ a+b+c+ab+ac+bc) 
anova(m0)  

Sem os coeficientes, dái eu obteria os coeficientes e depois a análise certo? parece que estou querendo fazer o contrário.