使用gam s()包装的ExtractVars中的模型公式无效
我想做的很简单。尝试拟合gam模型,在该模型中,我可以确定smoomthing样条曲线的df(公式部分中的“s”函数)。调用gam模型应该会返回拟合,但它会给出错误。我做错了什么使用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)+
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工作得更好,但是避免任何问题更简单,在任何一个会话中只加载一个包。是的,我注意到了..谢谢