Hi No problem with eigen (also in the asymmetric case, unless you set the option symmetric=T for an asymmetric matrix) or svd. I think your problem is here:
Pi2<-Pi %*% t(Pi) #that means Pi2 is symmetric and the eigenvalues should
#be the squared eigenvalues of Pi This not true generally but only for symmetric Pi (i.e. Pi=t(Pi)) Yours Marcel -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._