Needed to re-run models with intercept set to 0 to avoid negative predictions of mass. Cleaned up some of the rest of the code.