회귀 잔차 및 적합 값 래스터 출력 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) 
############################## 

가, 누군가가 나를 도울 수하십시오 : 여기

내 스크립트입니다?

덕분에 많은

gianca

Answer 1

당신은 LM의 기능을 연구 할 필요가있다. 예를 장착하려면 (예측) 값의 경우, 당신은 할 수 있습니다 :

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

하지만

summary(m)$fitted.values