[R-sig-Geo] gridded time series analysis

[Martin]

this is what I did to perform a regression between two bricks (each brick
represent a time series):

r <- raster(brick1)
for (i in 1:ncell(r)) { 
r[i] = lm(as.ts(cellValues(brick1, i)) ~ as.ts(cellValues(brick2,
i)))$coefficients[2]
}

The result will be a slope raster, but it really takes a lot of time, so
maybe there is a better solution..


-- 
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/gridded-time-series-analysis-tp5775651p5778472.html
Sent from the R-sig-geo mailing list archive at Nabble.com.