[2018年最新整理]R包mgcv的广义加模型函数gam().docxVIP

R包mgcv的广义加模型函数gam（）library(mgcv)set.seed(0) ## 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)plot(b,pages=1,residuals=TRUE)??## show partial residuals[＃显示部分残差]plot(b,pages=1,seWithMean=TRUE) ## `with intercept CIs[＃拦截“证明书]## same fit in two parts [＃同样适合两部分]G - gam(y~s(x0)+s(x1)+s(x2)+s(x3),fit=FALSE,data=dat)b - gam(G=G)print(b)## change the smoothness selection method to REML[＃REML法改变平滑的选择方法]b0 - gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat,method=REML)plot(b0,pages=1,scheme=1)## Would a smooth interaction of x0 and x1 be better?[＃X0和X1的流畅交互更好？]## Use tensor product smooth of x0 and x1, basis [＃使用张量积X0和X1顺利，基础]## dimension 49 (see ?te for details, also ?t2).[＃尺寸49（见详情TE，也T2）。]bt - gam(y~te(x0,x1,k=7)+s(x2)+s(x3),data=dat, method=REML)plot(bt,pages=1)?plot(bt,pages=1,scheme=2) ## alternative visualization[＃替代的可视化]AIC(b0,bt) ## interaction worse than additive[＃交互比添加剂差]## If it is believed that x0 and x1 are naturally on [＃如果它被认为X0和X1自然]## the same scale, and should be treated isotropically [＃相同的规模，应被视为各向同性]## then could try...[＃然后可以尝试...]bs - gam(y~s(x0,x1,k=50)+s(x2)+s(x3),data=dat, method=REML)plot(bs,pages=1)AIC(b0,bt,bs) ## additive still better. [＃添加剂更好。]## Now do automatic terms selection as well[＃现在做选择自动条款，以及]b1 - gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat, method=REML,select=TRUE)plot(b1,pages=1)## set the smoothing parameter for the first term, estimate rest ...[＃设置的第一任期内，估计截断的平滑参数...]bp - gam(y~s(x0)+s(x1)+s(x2)+s(x3),sp=c(0.01,-1,-1,-1),data=dat)plot(bp,pages=1,scheme=1)## alternatively...[＃或者...]bp - gam(y~s(x0,sp=.01)+s(x1)+s(x2)+s(x3),data=dat)# set lower bounds on smoothing parameters [平滑参数设定的下限]bp-gam(y~s(x0)+s(x1)+s(x2)+s(x3),?min.sp=c(0.001,0.01,0,10),data=dat)?print(b);print(bp)# same with REML[同与REML法]bp-gam(y~s(x0)+s(x1)+s(x2)+s(x3),?min.sp=c(0.1,0.1,0,10),data=dat,method=REML)?print(b0);print(bp)## now a GAM with 3df regression spline term amp; 2 penalized terms[＃现在与的3DF回归样条