## Not run:
# ## generate some data
# set.seed(111)
# n <- 500
#
# ## regressors
# dat <- data.frame(x = runif(n, -3, 3), z = runif(n, -3, 3),
# w = runif(n, 0, 6), fac = factor(rep(1:10, n/10)))
#
# ## response
# dat$y <- with(dat, 1.5 + sin(x) + cos(z) * sin(w) +
# c(2.67, 5, 6, 3, 4, 2, 6, 7, 9, 7.5)[fac] + rnorm(n, sd = 0.6))
#
# ## estimate model
# b <- bayesx(y ~ sx(x) + sx(z, w, bs = "te") + fac,
# data = dat, method = "MCMC")
#
# ## now show summary statistics
# summary(b)
# ## End(Not run)
Run the code above in your browser using DataLab