關於在R中使用GAM模型

Question

目前，我正在使用GAM模型複製Wood（2006）關於芝加哥空氣污染與死亡率之間關係的練習。

所以，我遵循他在他的書中使用的代碼。該圖書館被稱爲「gamair」，數據框架是「芝加哥」。

我把這個（酷似書）：

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，我可以再現您的錯誤。