[R-sig-Geo] gridded time series analysis
Martin
martin_brandt at gmx.net
Fri Nov 26 22:35:11 CET 2010
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.
More information about the R-sig-Geo
mailing list