Predict
m1 <- lm(y ~ A*B, data = mydata)
p  <- predict(m1, type="terms", terms="A:B")