If I wanted to have a von Mises distribution that is bounded with a lower bound of 350 degrees and an upper bound of 30 degrees (so encompassing 40 degrees on a 0-359 continuous scale), how would I implement that? It looks like the documentation says that von Mises go from a scale of -np.pi to np.pi, so what happens if I enter a value outside of that? Does it automatically get converted into the -np.pi to np.pi scale?
For example would the following code work correctly?
a = pm.Uniform('a', lower = np.deg2rad(-10.), upper = np.deg2rad(30.)) b = pm.VonMises('b', mu=a, kappa = 30.)