使用gam s()包装的ExtractVars中的模型公式无效

使用gam s()包装的ExtractVars中的模型公式无效,r,gam,mgcv,R,Gam,Mgcv,我想做的很简单。尝试拟合gam模型,在该模型中,我可以确定smoomthing样条曲线的df(公式部分中的“s”函数)。调用gam模型应该会返回拟合,但它会给出错误。我做错了什么 library(mgcv) library(gam) set.seed(2) ## simulate some data... dat <- gamSim(1,n=400,dist="normal",scale=2) b <- gam(y~s(x0,df = 5)+

我想做的很简单。尝试拟合gam模型,在该模型中,我可以确定smoomthing样条曲线的df(公式部分中的“s”函数)。调用gam模型应该会返回拟合,但它会给出错误。我做错了什么

    library(mgcv)
    library(gam)
    set.seed(2) ## simulate some data... 
    dat <- gamSim(1,n=400,dist="normal",scale=2)
    b <- gam(y~s(x0,df = 5)+s(x1)+s(x2)+s(x3),data=dat)
库(mgcv)
图书馆(gam)
种子(2)##模拟一些数据。。。

dat显然与mcgv包和“gam”版本有关。使用gam“1.14”而不使用mcgv是可行的。例如:

    library(ISLR)
    attach(Wage)
    library(gam)
    fit = gam(wage ~ s(year, 4) + s(age, 5) +   education, data = Wage)
    plot(fit)
    summary(fit)

不要同时加载gam和mgcv包。如果要使用
mgcv::gamSim()
但使用
gam::gam()
估计模型,则调用
mgcv::gamSim()
,但不加载包(无
库(“mgcv”)
调用)。我知道如果加载mgcv,开发人员已经做了一些事情使gam工作得更好,但是避免任何问题更简单,在任何一个会话中只加载一个包。是的,我注意到了..谢谢