r - 关于在 R 中使用 GAM 模型

Question

目前，我正在使用 GAM 模型复制 Wood（2006 年）关于芝加哥空气污染与死亡率之间关系的练习。

因此，我遵循了他在书中使用的代码。该库称为“gamair”，数据框为“chicago”。

我把这个（就像这本书一样）：

gam(death~s(time,bs="cr",k=200)+pm10median+so2median+o3median+tmpd,data=chicago,family=poisson)

但我得到一个错误：

Error in s(time, bs = "cr", k = 200) : 
  unused argument(s) (bs = "cr", k = 200)

有人可以解释一下吗？

非常感谢。

Answer 1

您是否也加载了mgcv包裹gamair？

library(gamair)
library(mgcv)
data(chicago)
gam(death~s(time,bs="cr",k=200)+pm10median+so2median+o3median+tmpd,data=chicago,family=poisson)

对我来说很好。

如果我将包gam不mgcv用于功能，我可以重现您的错误gam。