Ошибка в GAM Пуассона с нулевым завышением - неопределенная штрафная вероятность - PullRequest
Новая
       43

Ошибка в GAM Пуассона с нулевым завышением - неопределенная штрафная вероятность

0 голосов
/ 29 мая 2020

Я построил отрицательный биномиальный GAM в mgcv. 

> library(mgcv); library(mgcViz)
> gam14 <- gam(n ~ State + method + State*method + s(month, bs = "cc", k = 12, by = method) + s(time, by = method), data = control_l_gp_pad, family = nb(), method = "REML", knots = knots, select = TRUE)

Эта модель работает нормально, и я получаю следующие остаточные графики. 

> k <- getViz(gam14)
> check(k, a.qq = list(method = "auto"), a.hist = list(bins = 20))

Residual plots

Мои остаточные графики показывают, что моя модель плохо предсказывает количество 0, а также в некоторой степени подсчет 1, 2 и 3.

Затем я попытался последовать за ним Сообщение в блоге Гэвина Симпсона для соответствия GAM с нулевым надуванием Пуассона. Однако я получаю следующую ошибку. 

> gam14 <- gam(list(n ~ State + method + State*method + s(month, bs = "cc", k = 12, by = method) + s(time, by = method), ~ State + method), data = control_l_gp_pad, family = ziplss, method = "REML", knots = knots, select = TRUE)
Error in gam.fit5(x, y, sp, Sl = Sl, weights = weights, offset = offset,  : 
  indefinite penalized likelihood in gam.fit5

Я искал в Интернете решение этой ошибки, но пока ничего не добился.

У меня два вопроса:

1) Что означает эта ошибка? (Я могу догадаться, но я предпочитаю слышать людей, которые знают, о чем они говорят)

2) Есть ли у вас какие-либо предлагаемые исправления для этой ошибки?

Заранее спасибо.

...