library(mgcv)
set.seed(2) ## simulate some data...
dat <- gamSim(1,n=400,dist="normal",scale=2)
set.mgcv.options(gam.method="REML",plot.scheme=2)
b <- gam(y~s(x0)+s(x1)+te(x2,x3),data=dat)
b
plot(b,pages=1)
set.mgcv.options() ## reset to 'factory fresh'
b <- gam(y~s(x0)+s(x1)+te(x2,x3),data=dat)
b
plot(b,pages=1)
Run the code above in your browser using DataLab