Thanks. I think the plot you generated helps me a lot. Now the only warnings I have is a future warning and a The number of effective samples is smaller than 25% for some parameters. I’m not sure of the relative importance of this last warning. Could it be related to the fact you mentioned that n=10 (scalar)? I wrote obs = [10,12,19,9,15,7] and launched it and the warning remained. Perhaps with a bigger vector?