Source code for limix.qc.regress

# TODO: finish this or get rid of it
def regress_out(Y, X, return_b=False):
    r"""Regresses out X from Y"""
    from scipy.linalg import pinv

    Xd = pinv(X)
    b = Xd.dot(Y)
    Y_out = Y - X.dot(b)
    if return_b:
        return Y_out, b
    else:
        return Y_out