R help es - Sep 2013 - usando RasterBricks

Necesitaría ayuda con una cuestión. Pongamos por ejemplo estos valores
aleatorios para generar un RasterBrick s de 3 capas.
library(raster)
r <- raster(ncol=10, nrow=10)
r[]=1:ncell(r)
s <- brick(r,r,r)
s <- s * 1:3

Habría alguna forma sin hacer un bucle for, que me devolviese un vector con
los valores medios de cada una de las 3 capas.

mean_s<-NULL


for (i in 1:3)
{
mean_s[i]<-mean(s[[i]][])
}

Un saludo y gracias

	[[alternative HTML version deleted]]

Prueba:

apply(values(s),2,mean)

Un saludo. Olivier
--   
____________________________________

Olivier G. Nuñez
Email: onunez@iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 09/09/2013, a las 9:35, Jorge Naya Macias escribió:
> Necesitaría ayuda con una cuestión. Pongamos por ejemplo estos valores
> aleatorios para generar un RasterBrick s de 3 capas.
> library(raster)
> r <- raster(ncol=10, nrow=10)
> r[]=1:ncell(r)
> s <- brick(r,r,r)
> s <- s * 1:3
>
> Habría alguna forma sin hacer un bucle for, que me devolviese un  
> vector con
> los valores medios de cada una de las 3 capas.
>
> mean_s<-NULL
>
>
> for (i in 1:3)
> {
> mean_s[i]<-mean(s[[i]][])
> }
>
> Un saludo y gracias
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

	[[alternative HTML version deleted]]