[Rd] svd error (PR#631)

[polzehl@wias-berlin.de]

--=====================_24736660==_
Content-Type: text/plain; charset="iso-8859-1"; format=flowed
Content-Transfer-Encoding: quoted-printable

SVD-Error on

R 1.1.0
Windows 98


I get the following error applying  svd  on a positive definite matrix :

 > sk2
               [,1]         [,2]          [,3]          [,4]          [,5]
[1,]  1.0460139783  0.084356992 -2.810553e-04 -1.303237e-03  2.011809e-03
[2,]  0.0843569924  1.154650879 -5.152560e-04 -2.389213e-03  3.688230e-03
[3,] -0.0002810553 -0.000515256  1.000002e+00  7.960227e-06 -1.228821e-05
[4,] -0.0013032374 -0.002389213  7.960227e-06  1.000037e+00 -5.697973e-05
[5,]  0.0020118087  0.003688230 -1.228821e-05 -5.697973e-05  1.000088e+00
 > svd(sk2)
Error in svd(sk2) : error  4  in dsvdc

The problem seems to be connected to accuracy since

 > svd(sk2-.0000000000000000001)

works correctly.

I've saved the matrix  sk2  in the attached file.

Regards,

J=F6rg Polzehl
--=====================_24736660==_
Content-Type: application/octet-stream; name="sk2"
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="sk2"

UkRYMQoAAAACAAAAAAAAAAMAAAADZGltAAAAAAMAAAADc2syAAAAAAIAAAAAAAAAAAAAAAEAAAAA
AAAAAAAAAAL/////AAAADgAAAAAAAAAAAAAAGT/wvHknO/XaP7WYa3uwMRC/MmtTLMsBmL9VWix1
04ctP2B7EYuQF0Y/tZhre7AxED/yeXMzKfWjv0DiR8RYcTa/Y5KK1dPLfT9uNsdP3h7YvzJrUyzL
AZi/QOJHxFhxNj/wAAHM0oaKPuCxnUzM7+G+6cUur1Qr0r9VWix104ctv2OSitXTy30+4LGdTMzv
4T/wACa0QRFCvw3fsLftHfg/YHsRi5AXRj9uNsdP3h7YvunFLq9UK9K/Dd+wt+0d+D/wAFw7fg5q
AAAAAgAAAAAAAAAAAAAAAQAAAAAAAAAAAAAAAf////8AAAANAAAAAAAAAAAAAAACAAAABQAAAAX/
/////////////////////////w==
--=====================_24736660==_
Content-Type: text/plain; charset="iso-8859-1"; format=flowed
Content-Transfer-Encoding: quoted-printable


=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Dr. J=F6rg Polzehl
WIAS, Mohrenstr. 39
D-10117 Berlin, Germany
email: polzehl@wias-berlin.de
Tel. +49 30/20372 481
Fax: +49 30/2044975
WWW: http://www.wias-berlin.de/~polzehl

--=====================_24736660==_--


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._