Added plots with 95% confidence intervals.
Found cool new function called 'allEffects' which can simultaneously plot partial residuals for all variables, including interactions, with 95% CIs
Plot produced by plot(allEffects(best.mod))