[Rd] A different error in sample()

Martin Maechler m@echler @ending from @t@t@m@th@ethz@ch
Thu Sep 20 17:46:18 CEST 2018


>>>>> Martin Maechler 
>>>>>     on Thu, 20 Sep 2018 09:20:46 +0200 writes:

>>>>> Wolfgang Huber 
>>>>>     on Thu, 20 Sep 2018 08:47:47 +0200 writes:

    >> FWIW, I suspect this is related to the function
    >> R_unif_index that was introduced in src/main/RNG.c around
    >> revision 72356, or the way this function is used in
    >> do_sample in src/main/random.c.

    > Yes, it is just the use of 'dn' instead of 'n'
    > - a one letter thinko I'd say.

    > But *no*, it's much older than revision 72356; e.g., it's already in

    > R version 3.0.0 (2013-04-03) -- "Masked Marvel"

    > but not yet in

    > R version 2.15.3 (2013-03-01) -- "Security Blanket"

    > ----

    > Here, I clearly think we see a regression bug, and hopefully not
    > one that should trigger often in practice...
    > and  -- without any statistics about the consequences out in
    > package space --
    > I do think we should fix this in code and let the documentation
    > become "great again" ;-)

We have agreed that this is simply a regression and should be
fixed without a change to the documenation.

Consequently, ~ 5 minutes ago

$ svn log -v -c75338
------------------------------------------------------------------------
r75338 | maechler | 2018-09-20 17:38:46 +0200 (Thu, 20 Sep 2018) | 1 line
Changed paths:
   M /trunk/doc/NEWS.Rd
   M /trunk/src/main/random.c
   M /trunk/tests/reg-tests-1d.R

revert sample.int(<non-integer>, k, replace=TRUE) to sane pre_R-3.0.0 behaviour
------------------------------------------------------------------------

This is now back to "correct" behaviour  in  "R-devel (>= 75338)"

(and, as Duncan Murdoch also said by choosing this thread's
 Subject, this is really a different issue than the  "Bias in R's....")

Martin



More information about the R-devel mailing list