# (Based on belief that noise/10 is in decibel units)
noise.lm <- lm(noise/10 ~ size * type * side, data = auto.noise)
# Interaction plot of predictions
emmip(noise.lm, type ~ size | side)
# Confidence intervals
plot(emmeans(noise.lm, ~ size | side*type))
Run the code above in your browser using DataLab