Hi,
with small matrices eigen works as expected:
> eigen(cbind(c(1,4),c(4,7)), only.values = TRUE)
$values
[1] 9 -1
$vectors
NULL
> eigen(cbind(c(1,4),c(4,7)))
$values
[1] 9 -1
$vectors
[,1] [,2]
[1,] 0.4472136 -0.8944272
[2,] 0.8944272 0.4472136
> eigen(cbind(c(1,-1),c(1,-1)))
$values
[1] -3.25177e-17+1.570092e-16i -3.25177e-17-1.570092e-16i
$vectors
[,1] [,2]
[1,] 0.7071068+0i 0.7071068+0i
[2,] -0.7071068+0i -0.7071068-0i
With large(?) matrices eigen produces on several of my systems
segfaults:
> data(iris)
> D <- dist(iris[,-5])
> str(D)
Class 'dist' atomic [1:11175] 0.539 0.51 0.648 0.141 0.616 ...
..- attr(*, "Size")= int 150
..- attr(*, "Diag")= logi FALSE
..- attr(*, "Upper")= logi FALSE
..- attr(*, "method")= chr "euclidean"
..- attr(*, "call")= language dist(x = iris[, -5])
> eigen(D)
*** caught segfault ***
address (nil), cause 'unknown'
Traceback:
1: .Call("La_rs", x, only.values, PACKAGE = "base")
2: eigen(D)
All systems are Gentoo systems, i.e. R-2.12.1 and blas-atlas-3.9.23-r4
is installed by compiling the sources. Recompiling R and blas-atlas did
not solve the issue. This issue seems the reason that example(svm)
creates segfaults, too.
Any hint is appreciated.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
