R-beta: R-0.62.1 under Digital Unix

[Peter Dalgaard BSA]

Albrecht Gebhardt <albrecht.gebhardt at uni-klu.ac.at> writes:

> 
> Now I tried the same compilation procedure as you with 0.62.1 on my alpha
> (with DU 4.0D and DEC Fortran (f77) 5.0) and my resulting binary produces
> correct results for your tests, but gives other curios results:
> 
> > 1:3
>  [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15
> > seq(1,5,by=0.2)
>  [1] 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8
> > seq(1,5)
>  [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15
> > seq(1,20)
>  [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15
> > seq(1,5,length=4)
>  [1]  1.000000  2.333333  3.666667  5.000000  6.333333  7.666667  9.000000
>  [8] 10.333333 11.666667 13.000000 14.333333 15.666667 17.000000 18.333333
> [15] 19.666667 21.000000  5.000000
> 
> 
> I guess this behavior causes man problems, because many functions generate
> internally some index vectors with seq(). This was also my first guess
> after my troubles with R 0.61.2 and R 0.61.3 and the alpha.   

Have you tried pointing a debugger at the seq() function in
src/main/seq.c ? Or just inserting test printouts at appropriate
places? Printing n1,n2 and s1 and s2 just after the calls to asReal could
be revealing. (Use printValue(s1) and ..(s2)).

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