I think I figured it out, the model wasn’t converging because the offsets were occasionally changing by one or two and the lambdas had to adjust when that happened. If I model once to get the approximate offsets, then model again with the offsets fixed, I get good convergence.