Построение линии регрессии и r2 на диаграмме рассеяния между двумя растрами в R? - PullRequest
Новая
       75

Построение линии регрессии и r2 на диаграмме рассеяния между двумя растрами в R?

1 голос
/ 21 марта 2020

У меня есть два rasters in r

> lpjENLF
class      : RasterLayer 
dimensions : 2803, 5303, 14864309  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : memory
names      : VegC 
values     : 0, 17.99169  (min, max)

> geocarbon2
class      : RasterLayer 
dimensions : 2803, 5303, 14864309  (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : memory
names      : layer 
values     : 0, 26.53035  (min, max)

Я делаю диаграмму рассеяния между ними следующим образом: 

plot(lpjENLF, geocarbon2,maxpixels=900000)

Давая выход выше scatterplot

Я узнаю корреляцию по: 

> cor(values(geocarbon2), values(lpjENLF), use="complete.obs", method = 'pearson')
[1] 0.6883869

Как я могу построить линию регрессии и коэффициент r2 Пирсона на диаграмме рассеяния?

...