that's cool, I'm also interested in a similar problem but just with one
brick ending up with a slope raster as the output. It may be possible with
stackApply(). have a look. or maybe robert will chime in



On Fri, Nov 26, 2010 at 1:35 PM, Martin <martin_brandt@gmx.net> wrote:

>
> 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
>

	[[alternative HTML version deleted]]

