## Not run: ------------------------------------
# ## Simulate and run model with MCMC.
# set.seed(123)
# d <- GAMart()
# b <- bamlss(num ~ s(x1) + s(x2) + s(x3) + te(lon,lat), data = d)
#
# ## Continue sampling.
# a <- continue(b)
#
# ## Plot all samples.
# plot(samples(a, combine = FALSE))
## ---------------------------------------------
Run the code above in your browser using DataLab