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. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo