switch to multinomial composition for mixture sampling

ljwolf · ljwolf · commit cfc280d6586b · 2016-10-19T21:53:42.000-07:00
diff --git a/sklearn/mixture/base.py b/sklearn/mixture/base.py
@@ -385,7 +385,7 @@ def sample(self, n_samples=1):
 
         _, n_features = self.means_.shape
         rng = check_random_state(self.random_state)
-        n_samples_comp = np.round(self.weights_ * n_samples).astype(int)
+        n_samples_comp = rng.multinomial(n_samples, self.weights_).astype(int)
 
         if self.covariance_type == 'full':
             X = np.vstack([
