Some minor fixes:

• update parameter for scipy.optimize.brentq as float and not array

• update indexing of pandas series using .iloc

• Add pairwise testing: test each pair of conditions individually
• Remove numpy version restriction
• Refactor main module

• Replace Fortran implementation of gausq2 function with python

• Replace deprecated numpy.distutils

• Remove statannotations package. Plotting p-values is now implemented directly