Estou com dificuldades em desdobrar a interação em uma regressão em fatorial.

Segue um exemplo abaixo. Como obter um modelo de regressão de r em função de b para a1 e a2??????

a<-as.factor(c(rep(1,9),rep(2,9)))

b<-c(rep(0,3),rep(5,3),rep(10,3), rep(0,3),rep(5,3),rep(10,3))

r<-c(1,1.2,1.3,2,2.3,2.25,4,3.9,3.8,5,4.8,4.7,3.5,3,3.4,1.7,1.5,1.2)

b<-as.factor(b)

m<-lm(r~a*b)

anova(m)

b<-as.numeric(b)

regressão<-lm(r~b)

regressão

anova(regressão)

Como desdobrar b dentro de a1 e a2?????

Obrigado.

Marcelo.