pgeom accuracy (PR#6792) - R-devel

Mon, Apr 19, 2004 10:13 AM #

Full_Name: Morten Welinder
Version: snapshot
OS: 
Submission from: (NULL) (65.213.85.218)


This should fix the remaining two 1-p cancellation issues.

  double l_rt = log1p (-p) * (x + 1);
  if (log_p)
    return R_DT_Clog (l_rt);
  else
    return lower_tail ? -expm1 (l_rt) : exp (l_rt);

Peter Dalgaard

Mon, Apr 19, 2004 10:59 AM #

terra@gnome.org writes:

Never mind the actual content, but snapshot of what? You are being
lazy in a way that potentially causes more work for other people.  

Please state the version exactly, including the date (of the snapshot,
not the submission!) and whether it is of R-devel or R-patched.

O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907

terra@diku.dk

Mon, Apr 19, 2004 11:05 AM #

Apologies.  Snapshot from

    rsync --stats -z -rC --delete rsync.r-project.org::r-devel R

yesterday which contains the fragment pow(1-p,...) which is problematic
for small p.

(If the snapshot I get with the above command has an official name I'll
be happy to use that in future reports.)

Morten

Peter Dalgaard

Mon, Apr 19, 2004 11:19 AM #

terra@gnome.org (Morten Welinder) writes:

The sources should be time-stamped. Take a look at

?version
?bug.report

O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907

Martin Maechler

Mon, Apr 19, 2004 12:49 PM #

PD> terra@gnome.org (Morten Welinder) writes:
    >> Apologies.  Snapshot from
    >> 
    >> rsync --stats -z -rC --delete
    >> rsync.r-project.org::r-devel R
    >> 
    >> yesterday which contains the fragment pow(1-p,...) which
    >> is problematic for small p.
    >> 
    >> (If the snapshot I get with the above command has an
    >> official name I'll be happy to use that in future
    >> reports.)
    >> 
    >> Morten

    PD> The sources should be time-stamped. Take a look at

    PD>  ?bug.report
    PD>  ?version

if Morten has compiled the snapshot version,...

Otherwise, use  
	   <source>/VERSION
and	   <source>/date-stamp