diff --git a/src/qibo/states.py b/src/qibo/states.py
index 6cb0493f8c..2f8827e555 100644
--- a/src/qibo/states.py
+++ b/src/qibo/states.py
@@ -299,7 +299,11 @@ def __array__(self):
         return self.state(numpy=True)
 
     def probabilities(self, qubits=None):
-        """Calculates measurement probabilities by tracing out qubits.
+        """Calculates measurement probabilities by tracing out qubits. When
+        a noise model is applied to the circuit this method will return the
+        probabilities from the last execution shot, which may or may not
+        correspond to the noise-free state depending on whether the noise was
+        applied or not in that specific execution.
 
         Args:
             qubits (list, set): Set of qubits that are measured.