Displaying 1 result from an estimated 1 matches for "601x601".
2006 Mar 03
1
NA in eigen()
...about
my eigen decomposition of this matrix -- and if so, is it my matrix or the
program?
I have no idea what's causing it, and I can't get a reproducible example,
other than with my large matrix. My original matrix has no NAs in it. Here
is code, but of course it requires my original, 601x601 symmetric matrix
called mat
> any(is.na(mat))
[1] FALSE
> any(is.na(d))
[1] FALSE
> dim(mat)
[1] 601 601
> length(which(d==0))
[1] 5
> d<-rowSums(mat)
> temp1<-eigen(diag(d)-mat,symmetric=T)
> temp2<-eigen(diag(d)-mat,symmetric=T,EISPACK=T)
> any(is.na(tem...