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. 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. Best regards, Joseph Lu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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. (PR#1334-5) 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 R-patched. 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 yourself. -- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Dear Prof. Ripley,
Extremely sorry for my ignorance on the new function expm1()
that affects pweibull() in some cases!
Since I am not using C for a long time that also makes me to
ignore the news/changes that involving C codes.
Thanks very much for your response and explanation.
BTW, thanks Prof. Fox for his testing on Win2000.
Best regards,
Joseph Lu
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._