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 <matthew.baker@hunter.cuny.edu>
To   statalist@hsphsun2.harvard.edu
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: owner-statalist@hsphsun2.harvard.edu (on behalf of Zachary Neal <zpneal@gmail.com>)
>Subject: st: Mata not computing eigenvectors/values  
>To: statalist@hsphsun2.harvard.edu
>
>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
>zpneal@msu.edu
>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