gridded time series analysis
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.