Sorry @Evgenii for making you try that, I had not realised that my formulation implied a different model! Thanks @nhuurre and @andrjohns for explaining that so well, perhaps it will be valuable to others anyway. I’m going to edit my post to mention that already, so that people won’t use my code wrongly if they don’t read to the end.