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);
terra@gnome.org writes:> Full_Name: Morten Welinder > Version: snapshot > OS: > Submission from: (NULL) (65.213.85.218)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
>>>>> "PD" == Peter Dalgaard <p.dalgaard@biostat.ku.dk> >>>>> on 19 Apr 2004 18:15:14 +0200 writes: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