[SciPy-User] [SciPy-user] scipy.odr - Correct form of covx, covy matrix

Elfnor elfnor@gmail....
Mon Feb 13 20:38:03 CST 2012




Robert Kern-2 wrote:
> 
> 
> covx and covy aren't covariance matrices across observations, i.e. not
> the 10x10 matrices that you have here. Rather, what you need to
> provide is a length-10 array of covariance matrices, one for each
> observation. For the case where each of your observations are
> 1-dimensional, you don't construct covariance matrices at all, just
> provide the standard deviations:
> 
>   p_dat = RealData(p_x, p_y, sx=p_sx, sy=p_sy)
> 
> https://github.com/scipy/scipy/blob/master/scipy/odr/odrpack.py#L27
> 
> -- 
> 


OK. I gave an example with diagonal covariance matrices to make the problem
easy to explain. 
covy in ODR I see would be 10 2 by 2 matrices for the pearson's example,
each giving the covaraiance between the parameters at that observation.

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.

That is find x that minimizes (b - A*x)'*inv(V)*(b - A*x) similar to the
matlab function lscov.

I'm now looking at the GLS fit from scikits.statsmodels.

Thanks Eleanor
-- 
View this message in context: http://old.nabble.com/scipy.odr---Correct-form-of-covx%2C-covy-matrix-tp33312093p33319485.html
Sent from the Scipy-User mailing list archive at Nabble.com.



More information about the SciPy-User mailing list