Incorrect sample statement

[storopoli]

All of the sample statements are incorrect for Turing Models:

This implies that Turing will sample 2_000 iterations after discarding the nadapts default in NUTS() (which is 1_000). So this gives you 3_000 iterations with 1_000 as warmup.

sample(model, NUTS(), MCMCThreads(), 2_000, 4)

The correct approach would be:

sample(model, NUTS(1_000, 0.8), MCMCThreads(), 1_000, 4)