[R] R 1.5.1 on AIX 5.1

[Peter Dalgaard BSA]

"Andy Pierce" <apierce at stny.rr.com> writes:


>      [,1] [,2]
> [1,]    1    1
> [2,]    2    4
> [3,]    3    9
> [4,]    4   16
> [5,]    5   25
> > m1 <- crossprod(x)
> > m1
>      [,1] [,2]
> [1,]   NA  225
> [2,]  225  979
...
> This gets beyond the failing test case (but does nothing to fix the crossprod) and
> then the trees test case fails with:
...
 
> Any tips on what might be wrong or how to go about debugging these
> errors? 

Start with the first one. You're not going to trust a system that
miscomputes a crossproduct! 

You need to figure out exactly how that NA sneaks into m1. That's
about as close as I can get since I don't actually know what debuggers
you have available. With gdb, I'd set a breakpoint in symcrossprod()
inside src/main/array.c.

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._