From 6ced409bb42d127ebac85e563ee5ff01071b583f Mon Sep 17 00:00:00 2001
From: Chris Fonnesbeck <fonnesbeck@gmail.com>
Date: Tue, 18 Oct 2016 15:49:00 -0500
Subject: [PATCH] Replaced multinomial probability bounds test in logp

---
 pymc3/distributions/multivariate.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/pymc3/distributions/multivariate.py b/pymc3/distributions/multivariate.py
index 09d7bf3876..5f5d389cef 100755
--- a/pymc3/distributions/multivariate.py
+++ b/pymc3/distributions/multivariate.py
@@ -305,13 +305,13 @@ def logp(self, x):
         else:
             x_sum = x
             n_sum = n
-
         return bound(
             factln(n_sum) + tt.sum(x_sum * tt.log(p) - factln(x_sum)),
             tt.all(x >= 0),
             tt.all(x <= n),
             tt.eq(tt.sum(x_sum), n_sum),
-            tt.isclose(p.sum(), 1),
+            tt.all(p <= 1),
+            tt.eq(p.sum(), 1),
             n >= 0)