LRTs now work, both LRT(model0, model1) and and LRT(sim::Simulation).  Plotting now works for an LRT.  Everything looks good for one example with an isotropic fitted with constant, iso, and DTI.