Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: Mata not computing eigenvectors/values


From   Matthew J Baker <[email protected]>
To   [email protected]
Subject   Re: st: Mata not computing eigenvectors/values
Date   Sun, 9 Jan 2011 21:22:19 -0500 (EST)

While I am not sure about why the symeigensystem command isn't working, the problem might have something to do with the "A" matrix being rank deficient. In any case, I find that the commands:

eigensystem(A,X=.,L=.)

and

lefteigensystem(A,Y=.,M=.)

both generate values for eigenvalues and eigenvectors for the aforementioned A matrix (collected in X and L, and Y and M, respectively). 

Matt Baker

---- Original message ----
>Date: Sun, 9 Jan 2011 17:45:19 -0500
>From: [email protected] (on behalf of Zachary Neal <[email protected]>)
>Subject: st: Mata not computing eigenvectors/values  
>To: [email protected]
>
>Running the commands below, mata does not yield any eigenvectors or
>eigenvalues.  Computing eigenvectors/values in stata using matrix
>symeigen seems to work without a problem.  Any help would be
>appreciated.
>
>mata
>A = (0,0,0,0,0,1,0,1,0,0,0,0,0 \ ///
>     0,0,0,0,0,0,0,1,0,1,0,0,0 \ ///
>     0,0,0,1,0,1,1,0,1,0,1,1,1 \ ///
>     0,0,1,0,0,0,1,0,0,0,0,0,0 \ ///
>     0,0,0,0,0,0,0,1,0,0,0,0,0 \ ///
>     1,0,1,0,0,0,0,0,0,0,0,1,0 \ ///
>     0,0,1,1,0,0,0,0,0,0,0,0,0 \ ///
>     1,1,0,0,1,0,0,0,0,0,0,0,0 \ ///
>     0,0,1,0,0,0,0,0,0,0,0,0,0 \ ///
>     0,1,0,0,0,0,0,0,0,0,0,0,0 \ ///
>     0,0,1,0,0,0,0,0,0,0,0,0,0 \ ///
>     0,0,1,0,0,1,0,0,0,0,0,0,0 \ ///
>     0,0,1,0,0,0,0,0,0,0,0,0,0)
>eigenvector = 0
>eigenvalue = 0
>symeigensystem(network,eigenvector,eigenvalue)
>eigenvector
>eigenvalue
>end
>
>-- 
>Zachary Neal, PhD
>Department of Sociology
>Michigan State University
>[email protected]
>http://www.msu.edu/~zpneal
>*
>*   For searches and help try:
>*   http://www.stata.com/help.cgi?search
>*   http://www.stata.com/support/statalist/faq
>*   http://www.ats.ucla.edu/stat/stata/
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index