按组可视化混合模型中随机效应的方差

Question

我已经使用lmer在R中运行了线性混合模型。我正在尝试可视化随机效果结构。为了生成图，我使用了print(dotplot(ranef(RT.model.4, condVar=T))[['part_no']])，其中part_no是混合模型的随机效应。它会创建如下内容：

太好了。但是，我希望能够在图表中直观地看出两组参与者之间的差异（正在讨论随机效应）。我有A组和B组。在我的数据集中，有一列用于参与者类型，每行给出一个A或B值。

我想对图形进行颜色编码以显示来自A组和B组的参与者。或者更好的方法是创建两个单独的面板，每个组一个。

非常感谢您提出任何有关此操作的建议。

我已经使用lmer在R中运行了线性混合模型。我正在尝试可视化随机效果结构。要生成图，我使用了print（dotplot（ranef（RT.model.4，condVar = T））[['part_no']]）...

Answer 1

您应该可以使用groups=中的dotplot()选项。假设您的数据位于名为df的数据帧中，并且组变量位于group中，则可以使用

Answer 2

这是使用ggplot而不是lattice的一种方法（只是因为我比较熟悉），使用?dotplot.ranef.mer中示例的代码。您需要将数据中的治疗组