One of the big changes between 2.3 and 3.0-DEV is that fast=false
fits use a two-stage optimization processes in 2.3, but a single stage in 3.0-DEV. There are a few examples where this leads to very different fits…
You can install the dev version with ]add MixedModels#master