Ok, thanks. I’ve modified the code to do that and now it works fine.
data = [0.49131995252993826, 0.2774397121516236, 0.5381935236905224, 0.19753121107715765]
with pm.Model() as basic_model:
theta = pm.Uniform('theta', lower=max(data), upper=1)
y = pm.Uniform('y', lower=0, upper=theta,
observed=data)
print(basic_model.check_test_point())
trace = pm.sample(1000, progressbar=True)