Update: I have tried refitting the model due to finding this past question Loo moment matching not working - #10 by sergej but now R aborts as soon as I try loo2 <- loo(fit, moment_match = TRUE). This particular model is extremely quick to fit (4 chains of 2000 run in less than a minute), so I don’t see why this would break?