Hi,

When Aetkin linear model is used, problem holder has to provide weight matrix which is defined as $\Sigma^{-1/2}$. As far as the covariance matrix is always positive-defined the raising to the power $1/2$ is just decomposition like $L^T L$ which is performed by Cholesky technique. In my case the covariance matrix is not known in advance, but there are data to provide an estimator of covariance matrix. One can use the estimator instead of the matrix.

I am trying to implement this in a straightforward manner. I get the estimator of covariance matrix, then I convert it to correlation matrix. A dimension of my problem is close to 100 and input data are quite correlated. This leads to that ratio of maximum eigenvalue to minimal one is quite large. Furthermore, some of eigenvalues are negative due to rounding errors which means absence of positive-definiteness. The Cholesky method can't deal with such bad matrix.

Obviously, Kahan summation may increase accuracy by means of rounding error impact but only a little. As far as I know Page method is not exposed to the such issue because it doesn't require inverting of covariance matrix. If there is a way to deal with correlated input data in painless manner?

upd: I found paper[1] which describes several solutions to my problem. But all of them are `ad-hoc`

in my opinion. There are three naive ways discussed: put off-diagonal elements to zero, rotate the source problem to eigen-coordinates and then crop equations which regard to small eigen-values, put small eigen-values to some value by hands.