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

Question

Я построил отрицательный биномиальный 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))

Мои графики остатков предполагают, что моя модель плохо предсказывает счет, равный 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) Есть ли у вас какие-либо предлагаемые способы исправления этой ошибки?

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