我已经在R中拟合了lme模型,并具有logit转换后的响应。我无法找到执行logit转换的直接命令,因此我已经手动完成了。
logitr<-log(r/1-r)
然后,我将其用作我的lme模型中的两个因素和一个数值变量之间相互作用的响应。
model<-lme(logitr<-factor1*factor2*numeric,random=1|random)
现在,R显然不知道该模型是对数变换的。如何将其指定为R?
我没有运气尝试过:
update(model, tran="logit")
之所以要指定模型是对数变换的原因是因为我想使用emmeans软件包中的函数emmip绘制逆变换的结果,以显示变量之间相互作用的趋势。
通常(如果我只有因素的话,我只会使用:
update_refgrid_model<-ref_grid(update(model, tran="logit"))
但是当我想使用emmip绘制数值变量和因子之间相互作用的趋势时,这种方法不起作用。如果我指定:
emmip(update_refgrid_model, factor1~numeric|factor2, cov.reduce = range, type = "response")
然后我没有绘制任何趋势,仅显示了数字变量的平均水平的估计值。
所以,我如何指定对数变换并绘制具有与数值变量相互作用的因子的lme模型的反变换趋势?
我已经在R中拟合了lme模型,并具有logit转换后的响应。我无法找到执行logit转换的直接命令,因此我已经手动完成了。 logitr
您不更新模型对象,而是更新参考网格: