Using functions as variables in pymc model

You shouldn’t use eval, that’s just for debugging. Regarding your error it means you probably have something with shape=(1,) that should be () on the second argument of set_subtensor