[R-br] Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo)

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.