确定R中两个光栅之间r2的正确方法



我正试图在两个光栅之间找到r2。例如,我有两个光栅r1和r2,哪个代码计算r2?

cor(values(r1), values(r2), use="complete.obs", method = 'pearson')

堆叠光栅并将其转换为数据帧并找到r2。

stack1 <- stack(r1, r2)
df1 <- data.frame(na.omit(values(stack1)))
head(df1)
valueofdf1 = lm(gc ~ vc, data=df1)
summary(valueofdf1)$r.squared 

两个过程给出不同的r2。这就是我困惑的原因。我不想找到相关性,但我想找到两个光栅之间的r2值。

首先,cor不返回R^2,而是返回R。正如注释中所述,您得到的是预期的输出,R2=R*R。

因此,只要进行必要的调整,这两种方法都是可互换的(出于所有实际原因(。

最新更新