使用R 3.2.2从mgcv软件包运行GAM功能

Question

I have recently run into a problem running a GAM model from a previously working code. I believe it is related to an updated R-Version and an updated Version of the mgcv package. So it would be great to know if anyone has the same problem or has a solution to it.

I am currently running: R version 3.2.2 (2015-08-14) -- "Fire Safety" on Windows. And I am using the mgcv Package 1.8-7.

Below is an example code that re-produces the error message, when run on my computer.

###Load package
library(mgcv)

This is mgcv 1.8-7.

###Simulate some example data
set.seed(2) ## simulate some data... 
dat <- gamSim(1,n=400,dist="normal",scale=2)

###Run normal model
b <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat, family=gaussian())

This works.

###change the smoothness selection method to REML
b0 <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat,method="REML")

Gives the following error message:

Error in .C(C_gdi1, X = as.double(x[good, ]), E = as.double(Sr), Eb = as.double(Eb), : Incorrect number of arguments (48), expecting 47 for 'gdi1'

Thanks for your help!

Answer 1

我已经重新安装了R和mgcv软件包，看来这已经解决了问题。