# 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 DataCamp Workspace