R中的回归残差和拟合值栅格输出

Question

我正在使用 10 个光栅文件计算回归。 我能够获得斜率和 Rsquared 的栅格输出，但在生成残差和拟合值时遇到了一些麻烦。

这是我的脚本：

##########################################
require(raster)
require(rgdal)

#Create list of files
rasters <- list.files(pattern='\\.tif$')


time <- 1:nlayers(rasters)
fun <- function(x) { if (is.na(x[1])){ NA } else { lm(x ~ time)$coefficients[2] }}
slope <- calc(rasters, fun)

time <- 1:nlayers(rasters)
fun <- function(x) { if (is.na(x[1])){ NA } else { lm(x ~ time)$;summary(m)$r.squared }}
Rsquared <- calc(rasters, fun)

time <- 1:nlayers(rasters)
fun <- function(x) { if (is.na(x[1])){ NA } else { lm(x ~ time)$;summary(m)$residuals }}
residuals <- calc(rasters, fun)

time <- 1:nlayers(rasters)
fun <- function(x) { if (is.na(x[1])){ NA } else { lm(x ~ time)$;summary(m)$fitted.values }}
fitted.values <- calc(rasters, fun)
##############################

请问，有人可以帮帮我吗？

非常感谢

吉安卡

Answer 1

你需要学习 lm 函数。例如，要获得拟合（预测）值，您可以这样做：

y = runif(10)
x = 1:10
m = lm(y~x)
predict(m)

但不是

summary(m)$fitted.values