GAMLSS中pvc ()的相互作用项

Question

我们对位置、尺度和形状的广义加法模型进行拟合，以了解月份(连续)和年份(范畴)对响应变量Y的影响，以研究Y在不同年份的季节变化，包括年与月之间的相互作用项。所使用的R代码如下：

M1<-gamlss(Y~cs(月份)+pvc(月份，by=Year_N)，family=NBI，data = data)

函数pvc()是变系数函数。这种相互作用取的是解释变量x的线性系数，根据另一个解释变量r，x是平稳变化的。

当我运行此代码时，出现以下错误：“有50个或更多警告(使用警告()查看前50个)”。

当使用包“gamlss”(2020年2月3日)上可用的下列数据集和代码时，也会出现此错误：

租金数据

先考虑一个因素

数据(租金)

现在考虑到这个因素

r2 <- gamlss(R~pvc(F1，by=B)，data=rent)

有人能解释一下可能发生的事吗？

乔纳塔斯普拉多

Answer 1

警告不是错误，因此模型很可能没有失败，但在达到收敛之前，可能在最大迭代阈值处停止了优化。增加迭代拟合的参数值有助于避免警告信息。例如，"method = RS(100)“