mcp
mcp copied to clipboard
Set seed
There should be a seed-argument for all functions involving randomness cf this question (and comment-answer) on StackOverflow: https://stackoverflow.com/questions/73416610/why-is-set-seed-not-working-for-mcp-r-package
Suggested behavior for mcp()
:
-
mcp(..., seed = 42)
sets a seed before producingchains
"random" numbers, passing them as seeds to rjags chains. -
mcp(..., seed = c(1,2,3))
requireslength(seed) == length(chains)
and passes these seeds directly to rjags chains. -
mcp(..., seed = NULL)
(default) lets rjags set seeds for the chains.
tidy_samples()
, fitted()
, plot()
, etc. should just accept an atomic numeric seed
or NULL.