test for negative RTs

Author: poldrack

rtanalysis/rtanalysis.py

@@ -56,6 +56,10 @@ def fit(self, rt, accuracy, verbose=True):
         rt = rt.mask(~accuracy)
         self.meanrt_ = rt.mean()
 
+        try:
+            assert rt.min() >  0
+        except:
+            raise ValueError( "negative response times found")
         if verbose:
             print(f"mean RT: {self.meanrt_}")
             print(f"mean accuracy: {self.meanacc_}")
@@ -75,3 +79,5 @@ def _ensure_series_type(var):
         if type(var) is not pd.core.series.Series:
             var = pd.Series(var)
         return var
+
+# %%