compare function in indicator species analysis works but can be made
more efficient using np.where()