In [63]:
%matplotlib inline
import matplotlib
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import vonmises
import numpy as np
In [65]:
x = np.linspace(-1, 1)
In [66]:
#1D Case
k1 = vonmises.pdf(x, kappa=1)
k10 = vonmises.pdf(x, kappa=10)
k100 = vonmises.pdf(x, kappa=100)
X, Y = np.meshgrid(k1, k10)
fig1 = plt.figure()
ax = fig1.add_subplot(111)
ax.plot(k1,k10)
Out[66]:
[<matplotlib.lines.Line2D at 0x7ff8569d2f10>]
In [67]:
## 2D Case
## The plot of k1 vs k10 should be the same as above
K = vonmises.pdf([x,x], kappa=[[1],[10]])
fig1 = plt.figure()
ax = fig1.add_subplot(111)
ax.plot(K[0],K[1])
Out[67]:
[<matplotlib.lines.Line2D at 0x7ff85681d590>]
In [ ]: