I’m still a pymc3 newbie, so I cannot explain exactly why this fixed it, but yes that’s what I’m saying. When I had data_y.values.shape = (300000,), then this line would crash with a memory error:
Y_obs = pm.Normal(‘Y_obs’, mu=mu, sd=sigma, observed=data_y.values)
and when I changed the shape to (300000,1), it stopped crashing.