我目前正在研究混合效应逻辑中一些固定效应因子的边际效应。为此,我使用了非常有用的ggpredict
包中的ggeffects
函数。然后,我还使用了关联的tab_model
包的sjPlot
函数来生成包含奇数比率的表。但是,我对每个输出的输出感到惊讶:
1)现在,我看到我的因子预测变量的所有级别都包含在输出中(这与R常用的one level of each factor serves as a reference用于对比的伪编码相反)。是否可以在ggpredict
输出中保留参考水平?我希望将其用于i)对照人工计算,ii)将其与根据随机效应有条件地计算得出的glmer模型系数进行比较。
2)tab_model
提供的几率与我对原始glmer
模型(per the IDRE example procedure)提供的系数取幂后获得的几率相同。但是,我给人的印象是,计算出的OR是由边际系数得出的,该边际系数并未考虑模型中随机效应的影响(请参见以“许多人更喜欢” here,概率和图形”部分here,以及答案here的顶部,以获取更多信息)。反过来,这是否意味着tab_model
提供的固定效应变量的OR不能同样考虑随机效应的影响?如果是这样,是否有论据或其他手段可以做到这一点?
谢谢!