(A, AB, B)
return 0.5 * np.mean((A - AB) ** 2, axis=0) / np.var(np.r_[A, B], axis=0)
Total order estimator following Saltelli et al. 2010 CPC, normalized by
    sample variance
Total order estimator following Saltelli et al. 2010 CPC, normalized by
    sample variance
Total order estimator following Saltelli et al . 2010 CPC normalized by sample variance
def total_order(A, AB, B):
    """
    Total order estimator following Saltelli et al. 2010 CPC, normalized by
    sample variance
    """
    return 0.5 * np.mean((A - AB) ** 2, axis=0) / np.var(np.r_[A, B], axis=0)
def total_order ( A , AB , B ) : return 0.5 * np . mean ( ( A - AB ) ** 2 , axis = 0 ) / np . var ( np . r_ [ A , B ] , axis = 0 )