Thanks for the reply @jonsedar. Here is a stripped down version of my code. It works fine with the pm.sample_prior_predictive line commented out. df[“OOB Config”] is my variable with missing values. I have tried specifying “shape = 1” but it just gives a slightly different error.
test = pm.Model()
with test:
OOB = pm.Normal("OOB", mu=df["OOB Config"].mean(), sigma=10.0, observed=df["OOB Config"])
noise = pm.Normal("eps", mu=0, sigma=1.0)
y = pm.Deterministic("y", OOB + noise)
trace = pm.sample(10)
prior = pm.sample_prior_predictive()
Best,
-Brendan