Thank you I suspected it might be something of that nature but lacked experience with PyMC.
I will have a go with those two proposed alternatives and update.
I’ll try using pm.Data first but naively it doesn’t make sense to me that a single discrete value provided via data should be computationally more efficient than a single discrete value provided as a constant.