When trying to do orthogonal distance regression for a linear model, the scipy.odr package fails whenever I hand it errors (in x and y) in array shape (as opposed to scalar). The Error is “Problem is not full rank at solution”. It works, as long as these parameters are scalars, but that is not my requirement.