if (FALSE) {
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)
}
Run the code above in your browser using DataLab