如何控制lmerTest输出结果的绘图布局？

Question

我使用lme4和lmerTest运行一个混合模型，然后对我的模型使用反向变量消除(step)。这似乎很管用。在lmerTest中运行'step‘函数之后，我绘制了最终模型。“绘图”结果与ggplot2输出类似。

我想改变情节的布局。显而易见的答案是自己手动创建一个使用ggplot2的原始地块。如果可能的话，我只想简单地更改输出的布局，以便每个图(即在最终模型中绘制的因变量)都在它们自己的行中。

请参阅下面的代码和图表，以查看我的结果。注: plot有三列，我想要三行。此外，我没有提供样本数据(如果我也需要的话，请告诉我！)

library(lme4)
library(lmerTest)

# Full model
Female.Survival.model.1 <- lmer(Survival.Female ~ Location + Substrate + Location:Substrate + (1|Replicate), data = Transplant.Survival, REML = TRUE)

# lmerTest - backward stepwise elimination of dependent variables
Female.Survival.model.ST <- step(Female.Survival.model.1, reduce.fixed = TRUE, reduce.random = FALSE, ddf = "Kenward-Roger" )
Female.Survival.model.ST
plot(Female.Survival.model.ST)

​

​

Answer 1

创建这些图的函数称为plotLSMEANS。您可以通过lmerTest:::plotLSMEANS查看函数的代码。查看代码的原因是: 1)验证这些情节确实是基于ggplot2代码；2)查看是否可以确定需要更改哪些内容才能得到所需的内容。

在这种情况下，听起来您希望facet_wrap有一列而不是三列。我使用**lmerTest*函数step帮助页面中的示例进行了测试，看起来您可以简单地向图中添加一个新的facet_wrap层。

library(ggplot2)
plot(Female.Survival.model.ST) + 
    facet_wrap(~namesforplots, scales = "free", ncol = 1)

Answer 2

试一试: plot(difflsmeans(Female.Survival.model.ST$model，test.effs = "Location“)