There are a few post in general about dealing with Bad initial energy:
In your case, the observed value is outside of the bound of between1, between2… (in general, these kind of parameters you need to make sure they are ordered constrained). For example, p1 is bounded between lower=np.min(dataset1) and upper=between1, but the default of between1 is actually < np.max(dataset1).
I would suggest you to remove the TruncatedNormal and use Normal instead.