3dLME - REML vs ML

Hi AFNI Team,

A while back Gang Chen graciously added a flag to 3dLME that produces a bucket for log likelihood values, which is necessary if one wants to compare two models in a voxel-wise manner. That said, it appears that 3dLME runs using Restricted Maximum Likelihood (REML) by default. While this is ideal for fitting ‘final’ models, this becomes problematic when attempting to compare intermediate models, which requires fitting all models with ML instead of REML. Would it be possible to modify this to allow the user to specify if the model is fit with REML or ML?