diff --git a/plot.png b/plot.png
index c7f1f09..71bb909 100644
Binary files a/plot.png and b/plot.png differ
diff --git a/spherical_kde/__init__.py b/spherical_kde/__init__.py
index f0af471..86036aa 100644
--- a/spherical_kde/__init__.py
+++ b/spherical_kde/__init__.py
@@ -68,7 +68,7 @@ def __init__(self, phi_samples, theta_samples,
                              "shape as weights ({}!={})".format(
                                  len(self.phi), len(self.weights)))
 
-        sigmahat = VonMises_std(self.theta, self.phi)
+        sigmahat = VonMises_std(self.phi, self.theta)
         self.suggested_bandwidth = 1.06*sigmahat*len(weights)**-0.2
 
     def __call__(self, phi, theta):