Hi @contsili, and also @CharFraza because I realise this could be related to the issues you are describing.
The reason that the BLR plots look different from before is that the new version of the example has the BLR model configured with:
fixed_effect_slope=True, # We also model a fixed effect in the slope of the mean for each individual batch effect
fixed_effect_var_slope=True,
If these are both set to false, the plots are similar as before.
With both true:

With both false:
