[SciPy-User] [SciPy-user] scipy.odr - Correct form of covx, covy matrix
Tue Feb 14 08:55:47 CST 2012
On Tue, Feb 14, 2012 at 9:29 AM, Robert Kern <email@example.com> wrote:
> On Tue, Feb 14, 2012 at 02:38, Elfnor <firstname.lastname@example.org> wrote:
>> I am trying to solve a generalised least squares problem of the form A*x = b
>> where A is the problem design matrix and b is a vector of correlated
>> observations with covariance V. V is not diagonal.
to be precise V is the correlation of the error b-A*x, I assume that's
what you mean.
>> That is find x that minimizes (b - A*x)'*inv(V)*(b - A*x) similar to the
>> matlab function lscov.
> I think you can distribute the inv(V) by doing a Cholesky
> factorization, then doing a Cholesky-solve on both A and b, then doing
> a linear solve on the transformed A and b.
> from scipy import linalg
> cho = linalg.cho_factor(V)
> Ac = linalg.cho_solve(cho, A)
> bc = linalg.cho_solve(cho, b)
> x = linalg.solve(Ac, bc)
That's pretty much what statsmodels GLS does. The transformed
variables have a w prefix, wendog, wexog, wresid
I don't think we have many examples with the full sized (nobs, nobs) V matrix.
> Robert Kern
> "I have come to believe that the whole world is an enigma, a harmless
> enigma that is made terrible by our own mad attempt to interpret it as
> though it had an underlying truth."
> -- Umberto Eco
> SciPy-User mailing list
More information about the SciPy-User