Closed TomMaullin closed 1 year ago
REML is now implemented following PR #46 for the PSFS algorithm. It would be nice to add this functionality to the other algorithms but this is a low priority. It would also be nice to give the user the option of whether to use ML or REML but at present, I am unsure whether this functionality is operational. Until this is resolved this issue will remain open.
REML should be possible as it allows unbiased variance estimation. However, as the 3D code is designed for high n, this is a low priority.