Reproducibility of pymc calculation results

4

It has to do with random seeds. In numpy you also get different results depending on which variables you call first:

import numpy as np

seed = 123
rng = np.random.default_rng(seed)
print(rng.normal(), rng.uniform())

rng = np.random.default_rng(seed)
print(rng.uniform(), rng.normal())

-0.9891213503478509, 0.053821018802222675
0.6823518632481435, -0.3677866514678832

My command was just a trick to make sure the order was the same