我目前正在研究混合效应逻辑中一些固定效应因子的边际效应。为此,我使用ggpredict
了非常有用的ggeffects
软件包的功能。然后,我还使用tab_model
关联sjPlot
包的功能来生成包含奇数比的表格。然而,我对每个的输出有点惊讶:
1)我现在看到我的因子预测变量的所有级别都包含在输出中(与 R 通常的虚拟编码相反,其中每个因子的一个级别用作对比的参考)。是否可以在ggpredict
输出中保留参考电平?我希望用它来 i) 检查手动计算和 ii) 将其与 glmer 模型系数进行比较,这些系数不是根据随机效应有条件地计算的。
2) 提供的优势比tab_model
与我通过对原始glmer
模型提供的系数求幂获得的优势比相同(根据 IDRE 示例程序)。但是,我的印象是,计算的 OR 是从边际系数得出的,这些系数没有考虑到我的模型中随机效应的影响(请参阅此处以“许多人喜欢”开头的段落,“预测概率和图形”此处的部分,以及此处的最佳答案以获取更多信息)。反过来,这是否意味着固定效应变量的 OR 由tab_model
同样不考虑随机效应的影响?如果是这种情况,是否有论据或其他方法可以这样做?
谢谢!