[R] R-0.64.0 installation on SGI IRIX 6.5 machine

Jayshree Sarma jsarma at cs.gmu.edu
Fri Apr 30 16:37:42 CEST 1999


On 29 Apr 1999, Peter Dalgaard BSA wrote:

> Jayshree Sarma <jsarma at cs.gmu.edu> writes:
> 
> > gmake[1]: *** [simple-tests.Rout] Error 1
> > gmake[1]: Leaving directory `/usr/local/R-0.64.0/tests'
> > gmake: [test-All] Error 2 (ignored)
> > ***********************************
> > 
> > The test that fails from the simple-test.Rout file is
> > 
> > > Meps <- .Machine $ double.eps
> > > abs(1- .Machine$double.xmax * 
> >        10^(-.Machine$double.max.exp*log10(2)))/Meps < 1e3
> 
> The interesting bit of information missing from your post is what the 
> actual value of the left hand side of that comparison is. On
> Linux/i386 I get 560. If you get 1001 (say), I wouldn't worry too
> much, but if it's some wild number, then you may be in deeper trouble.
> (e.g. if log10(2) is not calculating the right value)
> -- 

  Given below is the output from R.  I do not get 560 but a rather large
number.  I ran the individual pieces too so you can compare the output
with the output from your machine and maybe tell me which is one that is
wrong. Is the code within abs()  correct?

***************output from R*************
Type    "q()" to quit R.

> log10(2)
[1] 0.30103
> -.Machine$double.max.exp
[1] -1024
> 10^(-.Machine$double.max.exp*log10(2))
[1] 0
> .Machine$double.xmax
[1] 1.797693e+308
> .Machine$double.xmax * 10^(-.Machine$double.max.exp*log10(2))
[1] 0
> abs(1- .Machine$double.xmax * 10^(-.Machine$double.max.exp*log10(2)))
[1] 1
> Meps <- .Machine $ double.eps
> .Machine $ double.eps
[1] 2.220446e-16
> abs(1- .Machine$double.xmax * 10^(-.Machine$double.max.exp*log10(2)))/Meps 
[1] 4.5036e+15
> abs(1- .Machine$double.xmax * 10^(-.Machine$double.max.exp*log10(2)))/Meps < 1e3
[1] FALSE

******************************end **************

Thanks,

Jayshree Sarma

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



More information about the R-help mailing list