Re: [Rd] pnorm problem (PR#7302)

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Fri 22 Oct 2004 - 07:08:32 EST


On Thu, 21 Oct 2004, Uwe Ligges wrote:

> terra@gnome.org wrote:
> > Full_Name: Morten Welinder
> > Version: 2
> > OS: Solaris(sparc)
> > Submission from: (NULL) (65.213.85.208)
> >
>
> Works for me with R-2.0.0, Solaris 5.7 (UltraSparc), gcc-3.2.3, 32-bit.
> So you have to specify compiler version, OS version, R version etc. much
> more precisely.

And for me under Solaris 8, gcc-3.4.2 or Sun cc Forte 7, 32- or 64-bit, with options as recommended in R-admin.

I suspect this is caused by unsafe compiler options.

> > (gdb) p pnorm(-100000.1, 0.0, 1.0, 0, 1)
> > $42 = NaN(0xfffffffffffff)
> >
> > Expected: -0
> >
> > (gdb) p pnorm(-100000.0, 0.0, 1.0, 0, 1)
> > $43 = -0
> >
> > Good.
> >
> > (I know this is not a typical usage.)

-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Received on Fri Oct 22 18:19:26 2004

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 09:00:53 EST