Base R

lm(response ~ A*B*C*D, data)
plot(xdata, ydata)

Package doetools

Installing

install.packages("devtools")
devtools::install_github("jensenlab/doetools")

Commands

get_effects(model)
show_effects(model, ordered="abs")
show_model(model)
farplot(data, factors, response)
add_dispersion(data, factors, response)
effect_dotplot(model)

Package daewr

Installing

install.packages("daewr")
library("daewr")

Half-normal plots

halfnorm(get_effects(model))

If you can an error about if (logC), you probably have NA effects in your model. Remove them with na.omit.

halfnorm(na.omit(get_effects(model)))