如何在LME4中获取收敛错误消息或最大梯度 笔记

编程问答 2022-08-11

问题描述

用lme4做随机效应模型后，如何获得最大梯度值？默认情况下会显示它是否大于0.002，但是我想制作一个原始函数，如下所示，无论该值是否大于0.002，该函数都需要返回最大梯度值。

  1. 制作原始模型
MM<- lmer(Y ~ 1 + X + (1 |cluster_ID ),data=data)
summary(MM)
  1. 如果上述原始模型的最大坡度大于0.002，请执行以下过程。
MM_par <- getME(MM,c("theta","fixef"))
RS_MM <- update(MM,start=MM_par,control=lmerControl(optCtrl = list(maxeval=1e4)))
summary(RS_MM)

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

mixedrr