Об использовании моделей GAM в R

Question

В настоящее время я воспроизвожу упражнение Wood (2006) о взаимосвязи между уровнем загрязнения воздуха и смертностью в Чикаго, используя модели GAM.

Итак, я следовал коду, который использовал в своей книге. Библиотека называется «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.