广义相加模型GAM的R语言实现

library(mgcv)
set.seed(2) ## simulate some data... 
dat <- gamSim(1,n=400,dist="normal",scale=2)
b <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat)
summary(b)

这是官方给出的代码,我没搞明白答题dat<-gamSim是什么意思,而且我自己运行的时候就报错。
有大神可以分享一下完整代码及注释吗

不知道你这个问题是否已经解决, 如果还没有解决的话:

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 以帮助更多的人 ^-^