if (requireNamespace("mgcv", quietly = TRUE)) {
x = 1:150
z = c(1:50, rep(1, 100))
y = as.numeric(BJsales)
m = mgcv::gam(y ~ s(x, bs = "b_bart", k=10) + s(z, bs = "b_bart", k=20))
summary(m)
plot(x, y)
lines(x, fitted(m), type="s", col="blue")
}
Run the code above in your browser using DataLab