Am Montag, den 21.02.2011, 11:11 +0100 schrieb Juergen
Rose:> 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.
Still some additional information:
root at moose:/root(5)# R
R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> sessionInfo()
R version 2.12.1 (2010-12-16)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods
base