Original repository for the paper MR-DoC2, however there is a custom function implemented in the R package umx()