3

R svd() 関数 (LAPACK を使用) が負の特異値を与える 399 x 399 の行列があります。これは起こるべきではありません - 誰もこれを見たことがありますか? LINPACK オプションを使用すると発生しないので、LAPACK svd のバグだと思います。

ganymede: R --vanilla

R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)

R is free software, etc ...

> load('A.dat')
> ls()
[1] "A"
> dim(A)
[1] 399 399
> 
> L1 <- svd(A)
> any( L1$d < 0 )
[1] TRUE
> L1$d[1:4]
[1] 80.18833 68.93905 61.62659 57.62883
> L1$d[396:399]
[1]  3.777844e-15  3.582460e-15  3.175665e-15 -6.512578e+00
> 
> L2 <- svd(A,LINPACK=TRUE)
> any( L2$d < 0 )
[1] FALSE
> L2$d[1:4]
[1] 80.18833 68.93905 61.62659 57.62883
> L2$d[396:399]
[1] 8.565532e-32 3.254162e-32 3.484425e-47 5.411232e-48
>
4

1 に答える 1

0

このバグレポートによると、問題はLAPACK3.4.1に対して再リンクすることで修正できる可能性があります。汎用バージョン3.1.1のアップグレードは簡単ではないようで、手順はOSによって異なります。

于 2012-07-12T18:26:44.193 に答える