search for: ottoboni

Before the next release we really should to sort out the bias issue in sample() reported by Ottoboni and Stark in https://www.stat.berkeley.edu/~stark/Preprints/r-random-issues.pdf and filed aa a bug report by Duncan Murdoch at https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17494. Here are two examples of bad behavior through current R-devel: set.seed(123) m <- (2/5) * 2^32...

...an Wichmann-Hill, it would be great to see the above behavior also for Mersenne-Twister. Thanks for considering. Best regards Kirill On 20.02.19 08:01, Gabriel Becker wrote: > Luke, > > I'm happy to help with this. Its great to see this get tackled (I've cc'ed > Kelli Ottoboni who helped flag this issue). > > I can prepare a patch for the RNGkind related stuff and the doc update. > > As for ???, what are your (and others') thoughts about the possibility of > a) a reproducibility API which takes either an R version (or maybe > alternatively a date) a...

...g. >> >> >> Best regards >> >> Kirill >> >> >> On 20.02.19 08:01, Gabriel Becker wrote: >>> Luke, >>> >>> I'm happy to help with this. Its great to see this get tackled (I've >>> cc'ed >>> Kelli Ottoboni who helped flag this issue). >>> >>> I can prepare a patch for the RNGkind related stuff and the doc update. >>> >>> As for ???, what are your (and others') thoughts about the possibility of >>> a) a reproducibility API which takes either an R versio...

Luke, I'm happy to help with this. Its great to see this get tackled (I've cc'ed Kelli Ottoboni who helped flag this issue). I can prepare a patch for the RNGkind related stuff and the doc update. As for ???, what are your (and others') thoughts about the possibility of a) a reproducibility API which takes either an R version (or maybe alternatively a date) and sets the RNGkind to the d...

...or Mersenne-Twister. > Thanks for considering. > > > Best regards > > Kirill > > > On 20.02.19 08:01, Gabriel Becker wrote: >> Luke, >> >> I'm happy to help with this. Its great to see this get tackled (I've >> cc'ed >> Kelli Ottoboni who helped flag this issue). >> >> I can prepare a patch for the RNGkind related stuff and the doc update. >> >> As for ???, what are your (and others') thoughts about the possibility of >> a) a reproducibility API which takes either an R version (or maybe >>...

I'm getting ready to submit an update of survival, and is my habit I run the checks on all packages that depend/import/suggest? survival.? I am getting some very odd behaviour wrt non-reproducability.? It came to a head when some things failed on one machine and worked on another.?? I found that the difference was that the failure was using the 3/27 release and the success was still on a

...rs from quantization effects that can get quite large when n is just below 2^31). A better method, used in Python, is to generate ceil(log_2(n)) pseudo-random bits, add 1, and discard values bigger than n. Attached is a short document explaining the issue in more detail. Best, Kellie -- Kellie Ottoboni Ph.D. Statistics '19, University of California, Berkeley Fellow at Berkeley Institute for Data Science Mobile: (650) 520-5056 Website: www.stat.berkeley.edu/~kellieotto -------------- next part -------------- A non-text attachment was scrubbed... Name: sample-bug.pdf Type: application/pdf Size...

Could this be related to "SIGNIFICANT USER-VISIBLE CHANGES The default method for generating from a discrete uniform distribution (used in sample(), for instance) has been changed. This addresses the fact, pointed out by Ottoboni and Stark, that the previous method made sample() noticeably non-uniform on large populations. See PR#17494 for a discussion. The previous method can be requested using RNGkind() or RNGversion() if necessary for reproduction of old results. Thanks to Duncan Murdoch for contributing the patch and Ga...

...y default produces files in serialization format version 2 (unless the original is already in format version 3). * The default method for generating from a discrete uniform distribution (used in sample(), for instance) has been changed. This addresses the fact, pointed out by Ottoboni and Stark, that the previous method made sample() noticeably non-uniform on large populations. See PR#17494 for a discussion. The previous method can be requested using RNGkind() or RNGversion() if necessary for reproduction of old results. Thanks to Duncan Murdoch for...

...y default produces files in serialization format version 2 (unless the original is already in format version 3). * The default method for generating from a discrete uniform distribution (used in sample(), for instance) has been changed. This addresses the fact, pointed out by Ottoboni and Stark, that the previous method made sample() noticeably non-uniform on large populations. See PR#17494 for a discussion. The previous method can be requested using RNGkind() or RNGversion() if necessary for reproduction of old results. Thanks to Duncan Murdoch for...

...y default produces files in serialization format version 2 (unless the original is already in format version 3). * The default method for generating from a discrete uniform distribution (used in sample(), for instance) has been changed. This addresses the fact, pointed out by Ottoboni and Stark, that the previous method made sample() noticeably non-uniform on large populations. See PR#17494 for a discussion. The previous method can be requested using RNGkind() or RNGversion() if necessary for reproduction of old results. Thanks to Duncan Murdoch for...