代码之家 › 专栏 › 技术社区 › DaveR

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

mgcv gam r

0

DaveR · 技术社区 · 6 年前

我想做的很简单。试图拟合一个gam模型,在那里我可以决定smoomthing样条(公式部分中的s函数)的df。调用gam模型应该会返回fit,但它会给出错误。我在这里做错什么了?

    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)

1 回复 | 直到 6 年前

1

0

DaveR 6 年前

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