Added new example data from Shipley 2013
Added model control arguments

All functions now have control = list() where control arguments (e.g.,
lmeControl(), lmerControl() etc.) can be used to help with model
convergence of d-sep models, partial models, etc.