R GAMLSS错误：从新数据集预测拟合值

我已经建立了一个GAMLSS模型，以调查基于现有数据集的鱼类种群驱动因素，但在尝试预测新数据集的响应值时出现了一个问题

要创建虚拟数据集，请执行以下操作：

Site <- c("Angle Crossing","Angle Crossing Pool","Casuarina Sands","Kambah Pool","Kissops Flat","Lawler Rd","Point Hut Crossing","Retallacks Hole","Scottsdale","Tharwa Sandwash")
Year <- round((rnorm(n=100, mean=2013, sd=2.530846)), digits=0)
LogTurbidity <- (rnorm(n=100, mean=2.026, sd=1.417185))
Datej <- rnorm(n=100, mean=105.0, sd=41.66927)
Catch <- round((rnorm(n=100, mean=1.596, sd=1.895757)), digits=0)
Catch <- ifelse(Catch < 0, 0, Catch)
mc.dummy <- data.frame(Site=Site, Catch=Catch, Year=Year, 
Turbidity=LogTurbidity, Datej=Datej)

---

我使用了上面所有相同的代码，没有

+random（Site）

，一切都很好，所以我的问题是如何让

predict（）

函数在模型中以随机效果工作？

没有数据，因此解释将是高度推测性的。你是在让我们猜测。谢谢你的澄清。现在包含了复制错误消息的数据。它现在是可复制的。我得到了同样的错误，即使有一些小的修改和包含一个数据参数。我也做了谷歌搜索，唯一与predict.gamlss的错误完全匹配的是一个russion网站，它显然会自动将SO问题翻译成俄语，在这个例子中是你的问题。有趣。谢谢@42-，我很想看看俄罗斯网站的进展。我还直接联系了GAMLSS，并将根据回复提供更新。

ZIPGAMM.dummy <- gamlss(Catch ~ cs(Datej, k=5) + cs(Year, k=5) + cs(Turbidity, k=5) + random(Site), family= ZIP(), data=mc.dummy, n.cyc=100)

MC.predicted <- predict(ZIPGAMM.dummy, newdata = MC.df.newdata, what="mu", type="response")    
new prediction 
Error in pred.s %*% rep(1, n.smooths) : 
requires numeric/complex matrix/vector arguments
In addition: Warning message:
contrasts dropped from factor random(Site)