r - 使用 R 3.2.2 从 mgcv 包运行 GAM 函数

Question

我最近遇到了从以前工作的代码运行 GAM 模型的问题。我相信它与更新的 R 版本和 mgcv 包的更新版本有关。所以很高兴知道是否有人有同样的问题或有解决方案。

我目前正在运行：R 版本 3.2.2 (2015-08-14) -- Windows 上的“消防安全”。我正在使用 mgcv 包 1.8-7。

下面是在我的计算机上运行时重新生成错误消息的示例代码。

###Load package
library(mgcv)

这是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())

这行得通。

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

给出以下错误消息：

.C(C_gdi1, X = as.double(x[good, ]), E = as.double(Sr), Eb = as.double(Eb), 中的错误：参数数量不正确 (48)，预期为 47 'gdi1'

谢谢你的帮助！

Answer 1

我已经重新安装了 R 和 mgcv 包，似乎这已经解决了这个问题。