nanoraw_stats.py 文件源码

def calc_fishers_method(pos_pvals, offset):
    pvals_np = np.empty(pos_pvals[-1][1] + 1)
    pvals_np[:] = np.NAN
    pvals_np[[list(zip(*pos_pvals)[1])]] = np.maximum(
        zip(*pos_pvals)[0], nh.SMALLEST_PVAL)

    fishers_pvals = [
        _calc_fm_pval(pvals_np[pos - offset:pos + offset + 1])
        if pos - offset >= 0 and
        pos + offset + 1 <= pvals_np.shape[0] and
        not np.any(np.isnan(pvals_np[pos - offset:pos + offset + 1])
                   ) else 1.0
        for _, pos, _, _ in pos_pvals]

    return fishers_pvals