带有 gam s() 包装器的 ExtractVars 中的模型公式无效

invalid model formula in ExtractVars with gam s() wrapper

我想做的很简单。尝试拟合一个 gam 模型,我可以在其中确定平滑样条曲线的 df(公式部分中的 's' 函数)。调用 gam 模型应该 return 适合,但它给出了错误。我在这里做错了什么?

    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)

显然与 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)