def METHOD_NAME(self):
    obs_r, obs_p, obs_perm_stats = hommola_cospeciation(
        self.hdist, self.pdist, self.interact, 0)
    exp_p = np.nan
    exp_r = 0.83170965463247915
    exp_perm_stats = np.array([])
    npt.assert_equal(obs_p, exp_p)
    self.assertAlmostEqual(obs_r, exp_r)
    npt.assert_equal(obs_perm_stats, exp_perm_stats)