## Not run:
# data("GAMart")
#
# ## normal response
# b <- bayesx(num ~ fac + sx(x1) + sx(x2) + sx(x3) +
# sx(long, lat, bs = "te") + sx(id, bs = "re"),
# data = GAMart)
# summary(b)
# plot(b)
#
# ## binomial response
# b <- bayesx(bin ~ fac + sx(x1) + sx(x2) + sx(x3) +
# sx(long, lat, bs = "te") + sx(id, bs = "re"),
# data = GAMart, family = "binomial", method = "REML")
# summary(b)
# plot(b)
#
# ## categorical response
# b <- bayesx(cat ~ fac + sx(x1) + sx(x2) + sx(x3) +
# sx(long, lat, bs = "te") + sx(id, bs = "re"),
# data = GAMart, family = "cumprobit", method = "REML")
# summary(b)
# plot(b)
# ## End(Not run)
Run the code above in your browser using DataLab