我们正在拟合位置、比例和形状的广义加性模型,以访问月份(连续(和年份(分类(对响应变量 Y 的影响。 为了研究Y在各年中的季节变化,纳入了年和月之间的交互作用项。使用的 R 代码如下所示:
M1<-gamlss(Y~cs(Month(+pvc(Month, by=Year_N(,family=NBI,data = data(
函数 pvc(( 是可变系数函数。这种交互作用取解释变量x的线性系数根据另一个解释变量r平滑变化。
当我运行此代码时,发生了以下错误:"有 50 个或更多警告(使用 warnings(( 查看前 50 个("。
当包"gamlss"上可用的以下数据集和代码时,也会出现此错误 (2020 年 2 月 3 日(用于:
租金数据
首先有一个因素
数据(出租(
现在有了因素
r2 <- gamlss(R~pvc(Fl, by=B(, data=rent(
有人可以解释一下会发生什么吗?
乔纳塔斯·普拉多
警告不是错误,因此模型可能没有失败,但在达到收敛之前,可能性优化在最大迭代阈值处停止。增加迭代拟合的参数值可能有助于避免出现警告消息。例如"方法 = RS(100(">