[Rd] lda and svd error (PR#630)

Torsten.Hothorn@rzmail.uni-erlangen.de Torsten.Hothorn@rzmail.uni-erlangen.de
Thu, 10 Aug 2000 15:11:56 +0200 (MET DST)


  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.
  Send mail to mime@docserver.cac.washington.edu for more info.

--1133331701-1931357602-965913218=:3057
Content-Type: TEXT/PLAIN; charset=US-ASCII


For a matrix A:

> A
            [,1]       [,2]
 [1,] -2.4491897  0.2345687
 [2,]  0.1530018 -0.5513182
 [3,] -1.3877140  1.0925618
 [4,]  0.5521078  0.7036401
 [5,]  0.4744903 -1.0868091
 [6,] -0.4296282  0.3931673
 [7,] -0.2413387 -0.7461586
 [8,]  1.0567564 -1.1807899
 [9,]  0.7903647  0.8294825
[10,]  1.6456857 -0.2484038
[11,] -1.1795527 -0.6904602
[12,]  2.6470244  0.7050178
[13,]  1.8169865 -0.1057402
[14,]  1.9720294 -0.1091834

and groupings [1] 0 0 0 0 0 0 0 1 1 1 1 1 1 1

> library(MASS)
> lda(A, groupings)
Error in svd(X, nu = 0) : error  2  in dsvdc

is reported. This is hard to reproduce on other machines, even 

> write.table(A, "dummy")
> B <- as.matrix(read.table("dummy"))
> lda(B, groupings)

on the same machine works fine. I add my .RData file renamed 
as DataR (maybe this can help to reproduce the error). The error occured
during a simulation, I added the sources as "svderr.R".

Torsten

--please do not edit the information below--

Version:
 platform = i686-pc-linux-gnu
 arch = i686
 os = linux-gnu
 system = i686, linux-gnu
 status = 
 major = 1
 minor = 1.0
 year = 2000
 month = June
 day = 15
 language = R

Search Path:
 .GlobalEnv, Autoloads, package:base

--1133331701-1931357602-965913218=:3057
Content-Type: APPLICATION/octet-stream; name=DataR
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.LNX.4.21.0008101513380.3057@artemis.imbe.med.uni-erlangen.de>
Content-Description: 
Content-Disposition: attachment; filename=DataR

UkRYMQoAAAAFAAAAAAAAAAMAAAADZGltAAAAAAEAAAABQQAAAAAAAAgAAAAI
Z3JvdXBpbmcAAAAMAAAADC5SYW5kb20uc2VlZAAAAAoAAAAKLlRyYWNlYmFj
awAAAAAAAgAAAAAAAAAAAAAAAQAAAAAAAAAAAAAABf////8AAAACAAAAAAAA
AAD/////AAAAEAAAAAAAAAAAAAAAAQAAACoAAAAqc3RvcChwYXN0ZSgiZXJy
b3IgIiwgeiRpbmZvLCAiIGluIGRzdmRjIikpAAD/////AAAAAgAAAAAAAAAA
/////wAAABAAAAAAAAAAAAAAAAEAAAAOAAAADnN2ZChYLCBudSA9IDApAAD/
////AAAAAgAAAAAAAAAA/////wAAABAAAAAAAAAAAAAAAAEAAAAYAAAAGGxk
YS5kZWZhdWx0KEEsIGdyb3VwaW5nKf////8AAAACAAAAAAAAAAD/////AAAA
EAAAAAAAAAAAAAAAAQAAABEAAAARTmV4dE1ldGhvZCgibGRhIikAAAD/////
AAAAAgAAAAAAAAAA/////wAAABAAAAAAAAAAAAAAAAEAAAAXAAAAF2xkYS5t
YXRyaXgoQSwgZ3JvdXBpbmcpAP////8AAAACAAAAAAAAAAD/////AAAAEAAA
AAAAAAAAAAAAAQAAABAAAAAQbGRhKEEsIGdyb3VwaW5nKf//////////////
////////////////////////////AAAAAgAAAAAAAAAAAAAAAQAAAAAAAAAA
AAAABP////8AAAANAAAAAAAAAAAAAAADAAAAASGIE+ETc44V/////wAAAAIA
AAAAAAAAAAAAAAEAAAAAAAAAAAAAAAP/////AAAADgAAAAAAAAAAAAAADgAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAP/AAAAAAAAA/8AAAAAAAAD/wAAAAAAAAP/AAAAAAAAA/8AAA
AAAAAD/wAAAAAAAAP/AAAAAAAAD/////AAAAAgAAAAAAAAAAAAAAAQAAAAAA
AAAAAAAAAv////8AAAAOAAAAAAAAAAAAAAAcwAOX8MHIaeA/w5WQE34/vb/2
NBOMsXdCP+Gq3gvUFnM/3l4MaikJlL/bfwcqUk/Gv87kL9G6OKo/8Oh5XMER
GD/pSqrXVI3mP/pUuo1+UDG/8t9ys0UCeEAFLRsm5nqQP/0SYHbpeIM//41u
ti/9pT/OBliWQmStv+GkZimVGgA/8XsiHJT4IT/mhDhKsAtrv/FjkfP5/5w/
2SmnFC117r/n4IgDNauUv/Lkg+5c7Wg/6ose7jvUyb/Py7Hhj7m3v+YYP/an
MfY/5o+BfEI5N7+7Ecl7N0bPv7vzcZLmB6UAAAACAAAAAAAAAAAAAAABAAAA
AAAAAAAAAAAB/////wAAAA0AAAAAAAAAAAAAAAIAAAAOAAAAAv//////////
////////////////////////////////
--1133331701-1931357602-965913218=:3057
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="svderr.R"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.LNX.4.21.0008101513381.3057@artemis.imbe.med.uni-erlangen.de>
Content-Description: 
Content-Disposition: attachment; filename="svderr.R"

bGlicmFyeShNQVNTKQ0KDQpTU2xkYSA8LSBmdW5jdGlvbihYLCBZLCBxKQ0K
ew0KCW1YIDwtIHJlcCgxLCBucm93KFgpKSUqJXQoYXBwbHkoWCwgMiwgbWVh
bikpDQoJbVkgPC0gcmVwKDEsIG5yb3coWSkpJSoldChhcHBseShZLCAyLCBt
ZWFuKSkNCglYIDwtIHJiaW5kKFgsWSkNCglteCA8LSByYmluZChtWCwgbVkp
DQoJWCA8LSBYIC0gbXggDQogCVcgPC0gKHQoWCklKiVYKQ0KICAgICAgICBl
ZyA8LSBlaWdlbihzb2x2ZShkaWFnKGRpYWcoVykpKSUqJVcpDQoJcSA8LSBs
ZW5ndGgoZWckdmFsdWVzW2VnJHZhbHVlcyA+IDFdKQ0KICAgICAgICBEIDwt
IGVnJHZlY3RvcnNbLDE6cV0NCgkjIEQgPC0gMS9zcXJ0KGRpYWcoVykpDQoJ
cmV0dXJuKEQpDQp9DQoNCkVycm9yUmF0ZSA8LSBmdW5jdGlvbihkaXNjcmlt
LCBuKQ0Kew0KCWZvbmUgPC0gc3VtKGFzLmludGVnZXIoZGlzY3JpbVsxOm5d
KS0xKS9uDQoJZm51bGwgPC0gKG4gLSBzdW0oYXMuaW50ZWdlcihkaXNjcmlt
WyhuKzEpOigyKm4pXSktMSkpL24NCgkjIGNhdCgiZm51bGw6ICIsIGZudWxs
LCAiIGZvbmUgIiwgZm9uZSwgIlxuIikgDQoJKGZudWxsICsgZm9uZSkvMg0K
fQ0KDQpzaW0xIDwtIGZ1bmN0aW9uKE4sbSxuLCBtdTAsIG11MSwgUykNCnsN
CgllcnJTdGFuZCA8LSBjKCkNCgllcnJTUyA8LSBjKCkNCg0KCWZvciAoaSBp
biAxOk4pDQoJew0KCQlsZWFyblggPC0gbXZybm9ybShtLCBtdTAsIFMpDQoJ
CWxlYXJuWSA8LSBtdnJub3JtKG0sIG11MSwgUykNCgkJZ3JvdXBpbmcgPC0g
YyhyZXAoMCwgbSksIHJlcCgxLCBtKSkNCgkJbGRhU3RhbmQgPC0gbGRhKHJi
aW5kKGxlYXJuWCwgbGVhcm5ZKSwgZ3JvdXBpbmcpDQoJCUQgPC0gU1NsZGEo
bGVhcm5YLCBsZWFyblksIDQpDQoJCUEgPC0gcmJpbmQobGVhcm5YLCBsZWFy
blkpJSolRA0KCQlsZGFTUyA8LSBsZGEoQSwgZ3JvdXBpbmcpDQoJCXRlc3RY
IDwtIG12cm5vcm0obiwgbXUwLCBTKQ0KCQl0ZXN0WSA8LSBtdnJub3JtKG4s
IG11MSwgUykNCgkJZGlzY3JTdGFuZCA8LSBwcmVkaWN0LmxkYShsZGFTdGFu
ZCwgcmJpbmQodGVzdFgsIHRlc3RZKSkkY2xhc3MNCgkJZGlzY3JTUyA8LSBw
cmVkaWN0LmxkYShsZGFTUywgcmJpbmQodGVzdFgsIHRlc3RZKSUqJUQpJGNs
YXNzDQoJCWVyclN0YW5kIDwtIGMoZXJyU3RhbmQsIEVycm9yUmF0ZShkaXNj
clN0YW5kLCBuKSkNCgkJZXJyU1MgPC0gYyhlcnJTUywgRXJyb3JSYXRlKGRp
c2NyU1MsIG4pKQ0KCQljYXQoIlN0YW5kOiAiLCBlcnJTdGFuZFtpXSwgIiBT
UzogIiwgZXJyU1NbaV0sICJcbiIpDQogCX0NCgljYXQoIkVycm9yIFJhdGUg
U3RhbmRhcmQgTERBOiAiLCBtZWFuKGVyclN0YW5kKSwgIlxuIikNCgljYXQo
IkVycm9yIFJhdGUgU1MgTERBOiAiLCBtZWFuKGVyclNTKSwgIlxuIikNCgly
ZXR1cm4obGlzdChlcnJTdGFuZCwgZXJyU1MpKQ0KfQ0KDQojIFMgPC0gbWF0
cml4KGMoMSwwLCAwLjUsIDAsIDAsIDEsIDAsIDAuOSwgMC41LCAwLCAyLCAw
LjMsIDAsIDAuOSwgMC4zLCAxKSwgNCwgNCkNClMgPC0gbWF0cml4KGMoMSww
LDAsMC45LDAsMSwwLDAsMCwwLDEsMCwwLjksMCwwLDEpLCA0LCA0KQ0KbXUw
IDwtIHJlcCgwLCA0KQ0KbXUxIDwtIHJlcCgxLCA0KQ0KbSA8LSBjKDUsNywx
MCwxNSwyMCw1MCwxMDApDQpuIDwtIDEwMDANCk4gPC0gMTAwDQoNCmVyciA8
LSBwbm9ybSgtMC41KnNxcnQodChtdTEgLSBtdTApJSolc29sdmUoUyklKiUo
bXUxIC0gbXUwKSkpDQpjYXQoIkVycm9yIHJhdGU6ICIsIGVyciwgIlxuIikN
Cg0KZUwgPC0gYygpDQplU1MgPC1jKCkNCg0KZm9yIChqIGluIG0pDQp7DQoJ
ZXN0ZXJyIDwtIHNpbTEoTiwgaiwgbiwgbXUwLCBtdTEsIFMpDQoJZUwgPC0g
YyhlTCwgbWVhbihlc3RlcnJbWzFdXSkpDQoJZVNTIDwtIGMoZVNTLCBtZWFu
KGVzdGVycltbMl1dKSkNCn0NCg0K
--1133331701-1931357602-965913218=:3057--

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._