[Rd] Re: [R] Strange means of numbers drawn from rpois (PR#730)
Peter Dalgaard BSA
08 Nov 2000 18:51:59 +0100
Peter Dalgaard BSA <firstname.lastname@example.org> writes:
> email@example.com writes:
> > On 8 Nov 2000, Peter Dalgaard BSA wrote:
> > >Done...
> > Great!
> However, be careful not to reply to the version of the mail without a
> PR# if there's a Cc: to R-bugs. Now there are *two* reports in the
> I've poked around a bit and it seems that the problem has to do with
> the code that tries to save time by reusing calculations if mu ==
> > range(sapply(1:2000, function(n) mean(rpois(10000, c(15,15+1e-8)))))
>  14.8692 15.1200
AHA! Spotted, I think. It is possible to return from rpois in step N,
in which case the initialisations in step P are not performed in that
invokation and not the next times either because it thinks that it has
been done previously.
I'm not quite sure what is the best way to fix it, though. Either kill
muprev when exiting in step N, or make the initialisations currently
in step P before entering step N. Or maybe save two versions of muprev
.. yes, that might be it - will try.
O__ ---- Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (firstname.lastname@example.org) FAX: (+45) 35327907
r-devel 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: email@example.com