[R] order by which the eigenvalues are presented
Duncan Murdoch
murdoch.duncan at gmail.com
Wed May 6 23:04:18 CEST 2015
On 06/05/2015 12:58 PM, Luis Borda de Agua wrote:
> Thank you, Bill, for your reply. However, I'm afraid I didn't explain myself properly.
>
>
>
> Imagine you have a 2x2 matrix
>
> Then the eigenvalues lambda_1 and lambda_2 are analytically calculated from
>
>
>
> lambda_1 = (-b+sqrt(delta))/2a
>
> lambda_2 = (-b-sqrt(delta))/2a
>
>
>
> where delta = b^2-4ac
>
>
>
> If delta>0 then lambda_1 > lambda_2 always. Otherwise their Real parts are equal.
The first condition is incomplete. You need a>0 as well.
>
>
>
> If we have a 3x3 matrix the three eigenvalues will have very complicated expressions:
>
>
>
> lambda_1 = f_1
>
> lambda_2 = f_2
>
> lambda_3 = f_3
>
>
>
> where f_1,f_2 and f_3 are functions of the elements of the matrix a11,a12...,a33, which are sampled from a given distribution (e.g. normal(0,1)).
>
>
>
> What I would like to know is from which expression (f_1,f_2 or f_3) comes the largest Re part of the eigenvalues. For example, does it always come from f_1 independently of the sampled values of a11,a12...,a33?
I don't know the answer, but I would expect it to depend on the entries
in the matrix in quite a complicated way.
Duncan Murdoch
More information about the R-help
mailing list