GAMLSS 中通过 pvc 的交互项 ()



我们正在拟合位置、比例和形状的广义加性模型,以访问月份(连续(和年份(分类(对响应变量 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(">

最新更新