Poor Accuracy of ADVI for Linear Regression

If you call approx = pm.ADVI() first, then do approx.fit(n=10000), it should separate the set up with the actual optimization.