# NOT RUN {
require(gamair);require(mgcv)
data(mpg)
b <- gam(list(city.mpg~fuel+style+drive+s(weight)+s(hp)+s(make,bs="re"),
hw.mpg~fuel+style+drive++s(weight)+s(hp)+s(make,bs="re")),
family=mvn(d=2),data=mpg)
plot(b,pages=1,scheme=1)
# }
Run the code above in your browser using DataLab