Pessoal, Estou fazendo um loop para variáveis (tair, tsoil, biomass, lai) e dentro desse loop estou fazendo outro com anos (2005, 2006, 2007) etc. Ao final dos dois loops, gostaria de ter um data.frame assim:
d var date value 1 tair 2005-01-01 10 2 tair 2006-01-01 20 3 tair 2007-01-01 30 4 tair 2008-01-01 40 5 tsoil 2005-01-01 10
6 tsoil 2006-01-01 20 7 tsoil 2007-01-01 30 8 tsoil 2008-01-01 40 etc etc Qual é o meio mais eficiente para fazer isso? Talvez até exista uma forma mais rápida de fazer isso do que usar loops, então qualquer sugestão será muito bem-vinda. Segue o CMR (que não está eficiente, é só para fins de compreensão do problema): ------------------------------------------------------------------------------------------------------------- library(raster) # Create random rasterstack and set names r1 <- r2 <- r3 <- r4 <- raster(nrow = 50, ncol = 100) r1[] <- 10 r2[] <- 20 r3[] <- 30 r4[] <- 40 s <- stack(r1, r2, r3, r4) # Set years names(s) <- c('X2005.01.01', 'X2006.01.01', 'X2007.01.01', 'X2008.01.01') # Create random variable names vars <- c('tempsoi', 'tempair', 'biomass', 'lai') for( i in 1:length (vars) ) { # Loop over years (names), extract raster average values and put in the data.frame d <- NULL for ( j in 1:length (names(s)) ) { var <- vars[i] date <- format (as.Date (gsub (pattern='\\.', replacement='-', substring (names(s)[j],2)))) value <- mean (getValues (s[[j]]), na.rm = TRUE) d <- rbind (d, data.frame(var, date, value)) } # end loop over names } # end loop vars ------------------------------------------------------------------------------------------------------------- Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/inde... Phone: (612) 323 9898