## LOAD THE CUSTOM FUNCTIONS source("myFunctions.R") ## LOAD THE DATA load("Data/BirdS2.Rda") ; NDO <- read.csv("Data/NomsDoizos.csv") Outliers <- c("ALAARV", "APUAPU", "MOTFLA") ; library(reshape) ## RESHAPE THE DATA AND DROP SOME OUTLIERS SPECIES SBS2 <- subset(BirdS2, !BirdS2$SPEC %in% Outliers) SBS2 <- aggregate(SBS2[, 3: 8], by= list(SBS2$MAILLE), mean) names(SBS2)[ 1] <- "MAILLE" ; PAGS2 <- melt(SBS2, id= "MAILLE") PAGS2$time <- as.numeric(substr(PAGS2$variable, 3, 4)) ## LINEAR INTERPOLATION OF DECENIAL TRENDS PAGS2$value <- ifelse(PAGS2$value> 10, 10, PAGS2$value); LinApprox(PAGS2) ## MERGE WITH THE GEOGRAPHIC FILES MAGS2 <- merge(DAT, CRD[, 6: 8], by= "MAILLE") coordinates(MAGS2) <- c("X1", "X2") ## CREATE THE SPATIO-TEMPORAL ANIMATION AnimHTML(MAGS2, ttle= "Aggregate Bird Index, Scenario S2", pal= "Spectral", brks= c(-671, -85, -50, -25, -5, 1, 5, 40, 150), t1= "Scenario S2, Aggregate Bird Index: ", t2= ", ref= 2003", dir= "./output/AGGINDS2/", name= "AGGINDS2.html")