Major update

Added options for modeling first-order interactions, as well as
quadratic terms.  Graphs are now optional (useful for large model
frames).  Robust standard errors can also now be specified using
sandwich/lmtest.