Note to self:
From R 3.6.0, base R’s
sample function works differently. This is to fix a bug in the previous version of the function (which caused some integers to be generated more often than others). A side effect of this is that simulations run with
sample() before v3.6.0 will give different results from those run with >= v3.6.0.
For backward compatibility, you can add the following line to the top of the script:
sample‘s behavior to the v3.6.0 default, use the following line of code:
More details from this Revolutions blog post.