if (FALSE) ## 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(a, which = "samples")
Run the code above in your browser using DataLab