Re: [Rd] Randomness not due to seed

2011-07-25 Thread Paul Johnson
On Tue, Jul 19, 2011 at 8:13 AM, jeroen00ms wrote: > I am working on a reproducible computing platform for which I would like to > be able to _exactly_ reproduce an R object. However, I am experiencing > unexpected randomness in some calculations. I have a hard time finding out > exactly how it oc

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Dirk Eddelbuettel
On 20 July 2011 at 18:02, peter dalgaard wrote: | | On Jul 20, 2011, at 15:38 , Dirk Eddelbuettel wrote: | | > | > On 20 July 2011 at 14:03, Jeroen Ooms wrote: | > | >> I think Bill Dunlap's answer addressed it: the claim appears to be false. | > | | > | Here is another example where there i

Re: [Rd] Randomness not due to seed

2011-07-20 Thread peter dalgaard
On Jul 20, 2011, at 15:38 , Dirk Eddelbuettel wrote: > > On 20 July 2011 at 14:03, Jeroen Ooms wrote: > | >> I think Bill Dunlap's answer addressed it: the claim appears to be > false. > | > | Here is another example where there is randomness that is not due to > | the seed. On the same machi

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Paul Gilbert
It does not look like your calculation is using the random number generator, so the other responses are probably more to the point. However, beware that setting the seed is not enough to guarantee the same random numbers. You need to also make sure you are using the same uniform RNG and any o

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Duncan Murdoch
to:r-devel-boun...@r-project.org] On Behalf Of Duncan Murdoch Sent: 20 July 2011 14:47 To: Jeroen Ooms Cc: r-devel@r-project.org Subject: Re: [Rd] Randomness not due to seed I would guess the error below is because of Java messing around in the hardware. It's pretty common on Windows for DLL

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Martyn Byng
-project.org [mailto:r-devel-boun...@r-project.org] On Behalf Of Duncan Murdoch Sent: 20 July 2011 14:47 To: Jeroen Ooms Cc: r-devel@r-project.org Subject: Re: [Rd] Randomness not due to seed I would guess the error below is because of Java messing around in the hardware. It's pretty comm

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Duncan Murdoch
I would guess the error below is because of Java messing around in the hardware. It's pretty common on Windows for DLLs to attempt to change the precision setting on the floating point processor; I hadn't seen that before on Linux, but that would be my guess as to the cause. It's also possibl

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Dirk Eddelbuettel
On 20 July 2011 at 14:03, Jeroen Ooms wrote: | >> I think Bill Dunlap's answer addressed it:  the claim appears to be false. | | Here is another example where there is randomness that is not due to | the seed. On the same machine, the same R binary, but through another | interface. First directly

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Jeroen Ooms
>> I think Bill Dunlap's answer addressed it:  the claim appears to be false. Here is another example where there is randomness that is not due to the seed. On the same machine, the same R binary, but through another interface. First directly in the shell: > sessionInfo() R version 2.13.1 (2011-0

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Duncan Murdoch
On 11-07-19 8:01 PM, Mike Marchywka wrote: Date: Tue, 19 Jul 2011 06:13:01 -0700 From: jeroen.o...@stat.ucla.edu To: r-devel@r-project.org Subject: [Rd] Randomness not due to seed I am working on a reproducible computing platform for which I would

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Mike Marchywka
> Date: Tue, 19 Jul 2011 06:13:01 -0700 > From: jeroen.o...@stat.ucla.edu > To: r-devel@r-project.org > Subject: [Rd] Randomness not due to seed > > I am working on a reproducible computing platform for which I would like to > be able to _exactly_ rep

Re: [Rd] Randomness not due to seed

2011-07-19 Thread William Dunlap
Did you actually see some rounding differences? The lm objects made in the calls to maklm will differ in the environments attached to the formula (because you made the formula in the function). If I change both copies of that .Environment attribute to .GlobalEnv (or any other environment), then i