是否有R-ggplot函数可用于打印每个子批次的R2和图例？

我正在处理的数据如下：

      "low_rd"

id  source      x   y
-----------------------
id1 source a    1   1
id1 source a    2   2
id1 source a    3   3
id1 source x    1   4
id1 source x    2   5
id1 source x    3   6
id2 source x    1   1
id2 source x    2   2
id2 source a    3   3
id2 source x    4   7
id2 source y    5   8
id2 source z    6   9
id3 source y    1   1
id3 source y    2   2
id3 source y    3   3
id3 source y    4   7
id3 source z    5   8
id3 source z    6   9
id4 source a    1   1
id4 source x    2   2
id4 source x    3   3
id4 source a    4   7
id4 source y    5   8
id4 source y    6   9
id5 source y    1   1
id5 source z    2   2
id5 source x    3   3
id5 source y    4   7
id5 source z    5   8
id5 source y    6   9

我在运行下面的代码时，在上面的数据上得到的曲线图

ggplot(low_rd,aes(x = x, y = y)) +
  geom_point(aes(colour = factor(source)), show.legend = FALSE) + 
  geom_smooth(method = 'lm', se = FALSE) + 
  geom_smooth(method = 'lm', aes(colour = factor(source)), se = FALSE, show.legend = FALSE) +
  facet_wrap(~id, scales = "free", shrink = FALSE)

详情如下:

      "low_rd"

id  source      x   y
-----------------------
id1 source a    1   1
id1 source a    2   2
id1 source a    3   3
id1 source x    1   4
id1 source x    2   5
id1 source x    3   6
id2 source x    1   1
id2 source x    2   2
id2 source a    3   3
id2 source x    4   7
id2 source y    5   8
id2 source z    6   9
id3 source y    1   1
id3 source y    2   2
id3 source y    3   3
id3 source y    4   7
id3 source z    5   8
id3 source z    6   9
id4 source a    1   1
id4 source x    2   2
id4 source x    3   3
id4 source a    4   7
id4 source y    5   8
id4 source y    6   9
id5 source y    1   1
id5 source z    2   2
id5 source x    3   3
id5 source y    4   7
id5 source z    5   8
id5 source y    6   9

每个图中的深蓝色线是我通过对每个“id”的所有数据点进行回归得到的总体回归线。但是，对于每个id，每个数据点都与特定的“源”关联，我绘制了与每个源关联的每个数据点的斜率，如图所示

我现在要做的是显示一个图例（最好是右上角），特定于每个子地块，它显示与特定id相关的所有源。此外，我还尝试显示每个ID的R2和总坡度值

我将非常感谢任何关于如何在我的代码中实现上述方面的输入

---

干杯

也许这可以帮助您：[）对于R2，您可以从

ggpubr

包中查看

stat\u cor（）

函数。

code至今+stat\u cor（aes（label=…rr.label..）