[R] Is this a bug of pweibull()? (Follow up)

ripley@stats.ox.ac.uk ripley at stats.ox.ac.uk
Sun May 12 11:44:17 CEST 2002

On Sun, 12 May 2002, Lu Chi-Hsien Joseph wrote:

> Please allow me to add just a little more about this:
> nothing wrong with pweibull(), namely, the two cases I reported:
> pweibull(3:10, 2) and pweibull(3:10, 2.1),
> in rw1041 and earlier version.
> I wonder this might just due to the change from rw1041 to rw1050,
> however, I can't find anything relevant (seems to me) in the News
> or Readme.

It's there in NEWS: search for pweibull:

    o	New function expm1() both in R and in C API, for accurate exp(x)-1;
	precision improvement in pexp() and pweibull() in some cases.

I have posted at least twice (and copied to you) that the bug is the
the C code for expm1 used on Windows, new in rw1050.  It's fixed in

Let's not flog this dead horse.

> Thanks Sundar for the suggestion of using 1 - pweibull(..., lower.tail=TRUE),
> but I think I rather stay with rw1041 and wait for correction comming since
> it would be a not-so-worth effort for me to modify many functions I wrote
> at this moment.

There will be no correction until rw1051 unless you build R-patched

Brian D. Ripley,                  ripley at 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 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

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