789676424d9767a82de32b02bb0db474ab06e972,skbio/stats/evolve/_hommola.py,,hommola_cospeciation,#Any#Any#Any#Any#,18

Before Change

    sort = (x + y).argsort()

    // calculate the observed correlation coefficient for this host/symbionts
    r = pearsonr(x[sort], y[sort])[0]

    // now do permutatitons. Initialize index lists of the appropriate size.
    mp = np.arange(num_pars)

After Change

            x_p = _get_dist(hosts_k_labels, hosts_t_labels, host_dist.data, mh)

            // calculate shuffled correlation coefficient
            perm_stats[i] = pearsonr(x_p, y_p)[0]

        p_value = ((perm_stats >= corr_coeff).sum() + 1) / (permutations + 1)

    return corr_coeff, p_value, perm_stats

Instances