我正在运行的分析具有比例形式的因变量(即变量的范围为0到1)。我正在多级模型中运行此分析。
DV:个人在21天中每天报告五次焦虑(编码为1)或不焦虑(编码为0)。当我尝试评估个体在几天中的焦虑轨迹时,我平均每天要花费这5个时间点(即,使其成为比例变量)。
IV:我正在尝试运行增长曲线模型。因此,我将学习日(如上所述,有21天)作为我的IV。
为阐明我的模型,级别1在个人中处于日常级别,级别2在个体中处于日常级别。
我知道我必须使用LME4函数才能在多级别上运行GLM,但是我不确定我应该使用哪个发行系列。这不是Poisson分布(因为结果不是整数),并且由于没有二进制结果,因此我无法进行逻辑回归。
好,所以我对此进行了深入研究,发现零膨胀的Beta模型最适合我的数据,因为我的结果是按比例的形式,并且包含值0。下面是我运行的代码,但我以某种方式得到错误。 关于为什么我可能会收到错误的任何建议?
焦虑<-glmmTMB(imean.anxstate〜1 + studyday1,data = ema2,family = list(family =“ beta”,link =“ logit”),ziformula =〜1) 摘要(着急)
``家庭''中缺少的某些组件:下游方法可能无法指定``家庭'',因为不推荐使用普通列表nlminb中的错误(start = par,objective = fn,gradient = gr,control = control $ optCtrl): NA / NaN梯度评估 时间停止在:0.05 0 0.22
谢谢!