mgcv ziplss predict contrasts 只能应用于具有 2 个或更多级别的因子

Question

我使用 mgcv package 中的 gam() 拟合了一个零膨胀泊松 model

 ziplss.fit.mixed <- gam(
    list(
      n.ind ~ as.factor(data.type) + s(year) +s(enz,bs="re") + s(clc.3,bs="re") + te(lon, lat), 
      ~ as.factor(data.type) + s(year) + s(enz,bs="re") + s(clc.3,bs="re") + te(lon, lat)),
    family=ziplss(),
    data = dat,
    control = gam.control(keepData = TRUE)
  )

每次我尝试进行预测时，我都会收到以下消息

Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels

我什至尝试使用我用于拟合的数据进行预测。 它给出了同样的错误

predict.gam(ziplss.fit.mixed, newdata = ziplss.fit.mixed$data, newdata.guaranteed = TRUE)

如果我改回具有相同公式的常规泊松游戏，那么它就可以工作。

知道为什么 ziplss 选项会触发错误吗？

谢谢你的帮助

Answer 1

由于某种未知的原因，在从 ziplss 家庭拟合进行预测时，您似乎必须格式化变量。

例如，如果您在 model 中写了“as.factor(my.var)”，那么您的新数据必须是 class 因子。