Displaying 1 result from an estimated 1 matches for "bugmatrix".
Did you mean:
bigmatrix
2005 Jul 31
1
eigen gives NaN in $vectors (PR#8041)
...2005-06-20)
OS: GNU/Linux 2.6.8-2-386, Debian testing
Submission from: (NULL) (65.24.56.41)
Relevant Bugs Ids : 7987, 7989
H is a 100x100 singular but symmetric matrix (a matrix defining the
neighbourhood structure for a spatial data) available from..
http://www.stat.ohio-state.edu/~pkapat/miscl/bugMatrix.RData
> load(file="bugMatrix.RData")
> ls()
"H"
> I <- diag(rep(1,100))
> g <- 0.1318467
> M <- I - g*H
> y <- eigen(M)
> y$vec[,91]
## all NaNs
Interestingly, as suggested in bug id 7989, this is very sensitive to the
scaling factor g. The er...