Adding lqrt initial commit.

Author: adam2392

.gitignore

@@ -1,5 +1,6 @@
 # Byte-compiled / optimized / DLL files
 __pycache__/
+.idea/*
 *.py[cod]
 *$py.class
 

dabest/_classes.py

@@ -470,6 +470,7 @@ def __init__(self, control, test, effect_size,
         from numpy.random import choice, seed
 
         import scipy.stats as spstats
+        import lqrt
 
         # import statsmodels.stats.power as power
 
@@ -658,6 +659,16 @@ def __init__(self, control, test, effect_size,
                 # in terms of rank (eg. all zeros.)
                 pass
 
+            # Likelihood Q-Ratio test:
+            # Assumes a gross-error model of distributions
+            try:
+                lqrt_ratio = lqrt.lqrtest_rel(control, test)
+                self.__pvalue_lqrt = lqrt_ratio.pvalue
+                self.__statistic_lqrt = lqrt_ratio.statistic
+            except ImportError:
+                # warnings.warn("To get")
+                pass
+
             standardized_es = es.cohens_d(control, test, is_paired=False)
             # self.__power = power.tt_ind_solve_power(standardized_es,
             #                                         len(control),
@@ -874,6 +885,21 @@ def statistic_wilcoxon(self):
         except AttributeError:
             return npnan
 
+    @property
+    def pvalue_lqrt(self):
+        from numpy import nan as npnan
+        try:
+            return self.__pvalue_lqrt
+        except AttributeError:
+            return npnan
+
+    @property
+    def statistic_lqrt(self):
+        from numpy import nan as npnan
+        try:
+            return self.__statistic_lqrt
+        except AttributeError:
+            return npnan
 
 
     @property

setup.py

@@ -42,9 +42,9 @@
             'numpy~=1.15',
             'scipy~=1.2',
             'pandas==0.24',
-            
             'matplotlib~=3.0',
-            'seaborn~=0.9'
+            'seaborn~=0.9',
+            'lqrt>=0.3.2'
         ],
         extras_require={'dev': ['pytest==5.2', 'pytest-mpl==0.10']},
         python_requires='>=3.5',