1

我必须创建大量（100 多个）线性模型的 ggplots。我想将 p 值（可能还有 R2）添加到每个图中。我知道可以使用ggpmisc. 在这里，我使用stat_fit_glance添加 p 值。我的“问题”是这两个都需要我先运行lm才能插入为公式= my_lm。

由于我必须创建大量绘图，我想知道是否有办法避免先创建 lm 对象，而在生成 ggplot 时简单地计算它？我可以使用 对箱线图进行 t 检验stat_compare_means，并且真的希望能找到一种使用 lm 的方法。

我的代码在下面。我希望能够跳过第一行代码：

my_lm <- lm(y ~ x)


ggplot(data = complete, aes(x= x, y = y))+  
geom_point()+
theme_classic()+
geom_smooth(method = "lm")+
labs(x="Ellenberg F", y = "Species richness")+
stat_fit_glance(method = 'lm',
              method.args = list(data = complete, formula = my_lm),
              geom = 'text',
              aes(label = paste("p-value = ", signif(..p.value.., digits = 4), sep = "")),
              label.x = 8.5, label.y = 25, size = 3)

我试过简单地把公式 = y ~ x 没有运气。

4

1 回答 1

2

ggpmisc::stat_fit_glance在:的帮助下method.args = list(formula = y ~ x)。
这意味着您不需要先运行lm。
您只能指定线性模型的公式。

library(ggpmisc)
set.seed(1)
n <- 100
x <- 8+rnorm(n)
y <- 11+x+2*rnorm(n)
complete <- data.frame(x, y)

summary(lm(y~x))
ggplot(data = complete, aes(x= x, y = y))+  
geom_point()+
theme_classic()+
geom_smooth(method = "lm")+
labs(x="Ellenberg F", y = "Species richness")+
stat_fit_glance(method = 'lm',
       method.args = list(formula = y ~ x),  geom = 'text', 
       aes(label = paste("p-value=", signif(..p.value.., digits = 4), 
                      "   R-squared=", signif(..r.squared.., digits = 3), sep = "")),
       label.x = 8.5, label.y = 25, size = 5)

在此处输入图像描述

于 2019-02-07T21:45:21.733 回答