# NOT RUN {
code <- nimbleCode({
mu ~ dnorm(0, 1)
x ~ dnorm(mu, 1)
})
output <- MCMCsuite(code,
data = list(x=3),
inits = list(mu=0),
niter = 10000,
monitors = 'mu',
MCMCs = c('nimble', 'nimble_RW'),
summaryStats = c('mean', 'sd', 'max', 'function(x) max(abs(x))'),
makePlot = FALSE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab