require(lsmeans)
noise.lm <- lm(noise ~ size * type * side, data = auto.noise)
# Interaction plot of predictions
lsmip(noise.lm, type ~ size | side)
# Confidence intervals
plot(lsmeans(noise.lm, ~ size | side*type))
Run the code above in your browser using DataLab