指定对数进行对数变换以绘制反变换趋势

问题描述 投票:0回答:1

我已经在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

r transformation nlme emmeans
1个回答
0
投票

您不更新模型对象,而是更新参考网格:

© www.soinside.com 2019 - 2024. All rights reserved.