In light of #219, we should test out https://github.com/arviz-devs/arviz-stats/blob/main/src/arviz_stats/base/stats_utils.py#L164 to see if there are any pain points with this numpy-only version compared to the scipy implementation.