Question

我是R的新手。现在，我想绘制数据（两个变量）并显示包括boxplot的回归线。我可以显示那些数据，除了r平方值和等式图。

下面是我显示图形的脚本

library (car)
scatterplot(FIRST_S2A_NDVI, MEAN_DRONE_NDVI, 
            main = "NDVI Value from Sentinel and Drone", 
            xlab = "NDVI Value from Sentinel", 
            ylab = "NDVI Value from Drone", 
            pch = 15, col = "black", 
            regLine = list(col="green"), smooth = FALSE)

图是这样的。 Comparison NDVI Value

现在，最后的任务是在我的图形上添加方程式和r平方值。我需要写什么脚本。我从Add regression line equation and R^2 on graph开始尝试了此脚本，但仍然不知道如何显示它们。

感谢您的阅读，希望对我有帮助。

p.s。

我的数据内容

OBJECTID SAMPLE_GRID FIRST_S2A_NDVI MEAN_DRONE_NDVI
1         1           1      0.6411405       0.8676092
2         2           2      0.4335293       0.5697814
3         3           3      0.7350439       0.7321858
4         4           4      0.7268013       0.8271566
5         5           5      0.3638939       0.5682631
6         6           6      0.1953890       0.3168246
7         7           7      0.4841993       0.7380627
8         8           8      0.4137447       0.3239288
9         9           9      0.8219178       0.8676065
10       10          10      0.2647872       0.2296441
11       11          11      0.8126657       0.8519964
12       12          12      0.2648504       0.2465738
13       13          13      0.5992035       0.8016030
14       14          14      0.2420299       0.3933670
15       15          15      0.5059137       0.7593807
16       16          16      0.7713419       0.8026068
17       17          17      0.3762540       0.5941540
18       18          18      0.5876435       0.7763927
19       19          19      0.2491609       0.5095306
20       20          20      0.3213648       0.4456958
21       21          21      0.2101466       0.1960858
22       22          22      0.3749034       0.4956361
23       23          23      0.5712630       0.7350484
24       24          24      0.8444895       0.8577550
25       25          25      0.3331450       0.4390229
26       26          26      0.1851611       0.4573663
27       27          27      0.4914998       0.2750837
28       28          28      0.7121390       0.7780228

Answer 1

用于将方程式和R平方值添加到当前图。您可以简单地使用y和x变量创建模型，格式化方程式，然后使用mtext函数粘贴到绘图上。

m <- lm(MEAN_DRONE_NDVI~FIRST_S2A_NDVI)

eq <- paste0("y = ",round(coef(m)[2],3),"x ",
             ifelse(coef(m)[1]<0,round(coef(m)[1],3),
                    paste("+",round(coef(m)[1],3))))         

mtext(eq, 3,-1)
mtext(paste0("R^2 = ",round(as.numeric(summary(m)[8]),3)), 3, -3)

您可以更改模型中的变量，还可以使用mtext函数中的第二个和第三个参数来更改文本的位置

散点图-添加方程式和r平方值

1 个答案: