(A, ABj, ABk, BAj, B)
return Vjk - Sj - Sk
Second order estimator following Saltelli 2002
Second order estimator following Saltelli 2002
Second order estimator following Saltelli 2002
def second_order(A, ABj, ABk, BAj, B):
    """Second order estimator following Saltelli 2002"""
    Vjk = np.mean(BAj * ABk - A * B, axis=0) / np.var(np.r_[A, B], axis=0)
    Sj = first_order(A, ABj, B)
    Sk = first_order(A, ABk, B)

    return Vjk - Sj - Sk
def second_order ( A , ABj , ABk , BAj , B ) : Vjk = np . mean ( BAj * ABk - A * B , axis = 0 ) / np . var ( np . r_ [ A , B ] , axis = 0 ) Sj = first_order ( A , ABj , B ) Sk = first_order ( A , ABk , B ) return Vjk - Sj - Sk