Last chance! 50% off unlimited learning
Sale ends in
eigen()
seemed
to be difficult, an unscaled version of this has been posted to the
web, accompanying an E-mail to R-help
(data(CAex)
dgCMatrix
.eigen(CAex)
fell into an
infinite loop (whereas eigen(CAex, EISPACK=TRUE)
has
been okay).data(CAex)
str(CAex) # of class "dgCMatrix"
image(CAex)# -> it's a simple band matrix with 5 bands
## and the eigen values are basically 1 (42 times) and 0 (30 x):
zapsmall(ev <- eigen(CAex, EISPACK=TRUE, only.values=TRUE)$values)
## i.e., the matrix is symmetric, hence
sCA <- as(CAex, "symmetricMatrix")
## and
stopifnot(class(sCA) == "dsCMatrix",
as(sCA, "matrix") == as(CAex, "matrix"))
Run the code above in your browser using DataLab