添加 p 值和 R2 ggplot [后续]

Question

这是一个后续question。当我运行下面给出的代码时，我收到警告消息，我认为这是由于我的代码中没有方面的要求，而链接中提到的源代码包括方面。看看，请让我知道哪个部分需要修改。期待！

代码：

library(dplyr) 
library(ggplot2)
library(ggpmisc)

df <- diamonds %>%
  dplyr::filter(cut%in%c("Fair","Ideal")) %>%
  dplyr::filter(clarity%in%c("I1" ,  "SI2" , "SI1" , "VS2" , "VS1",  "VVS2")) %>%
  dplyr::mutate(new_price = ifelse(cut == "Fair", 
                                   price* 0.5, 
                                   price * 1.1))


p <- ggplot(df, aes(x,y, color=factor(cut))) 
p <- p + stat_smooth(method = "lm", formula = y ~ x-1, size = 1, level=0.95) 
p <- p + geom_point(alpha = 0.3) 
p <- p + stat_poly_eq(aes(label = paste(..rr.label..)),
                      label.x.npc = "right", label.y.npc = 0.15, formula = formula, 
                      parse = TRUE, size = 3) + 
          stat_fit_glance(method = 'lm', method.args = list(formula = formula),
                      geom = 'text', aes(label = paste("P-value = ", 
                      signif(..p.value.., digits = 4), sep = "")),label.x.npc = 'right',
                      label.y.npc = 0.35, size = 3)
print(p)

警告信息：

1：stat_poly_eq() 中的计算失败： “闭包”类型的对象不是子集

2：stat_fit_glance() 中的计算失败： “闭包”类型的对象不是子集

Answer 1

简答：你需要添加

formula <- y ~ x

（或任何您定义的公式）在您调用ggplot 之前（即在读取p <- ggplot(...) 的行之前。

---

A "closure" is a type of function in R. 所以警告消息“'closure' 类型的对象不是子集的”意味着您运行的任何代码都没有期望一个函数对象。

当我们仔细查看您的代码时，我们会在您对 stat_poly_eq 和 stat_fit_glance 的调用中看到 formula = formula。注意formula是R中的一个函数。如果你没有单独定义一个formula对象，R会认为你指的是formula函数。 stat_poly_eq() 和 stat_fit_glance() 抱怨是因为他们希望函数中的 formula 参数 是 formula 类对象，而不是函数。

更一般地说，您不应该将公式命名为“公式”，因为这会造成混淆。你可以使用例如改为“模型”。