refactor run_model function, so that the analysis is run in the CXGN::MixedModel...