## LOAD THE CUSTOM FUNCTIONS AND THE DATA
source("myFunctions.R")
load("Data/BirdS0.Rda") ; NDO <- read.csv("Data/NomsDoizos.csv")
Outliers <- c("ALAARV", "APUAPU") ; library(reshape)
for (i in levels(BirdS0$SPEC)[!levels(BirdS0$SPEC) %in% Outliers]){
    ## SHAPE THE DATA AND DROP OUTILERS SPECIES
    sbs0 <- subset(BirdS0, BirdS0$SPEC== i)
    sbs0 <- aggregate(sbs0[, 3: 8], by= list(sbs0$MAILLE), mean)
    names(sbs0)[ 1] <-  "MAILLE" ; yop <- melt(sbs0, id= "MAILLE")
    yop$time <- as.numeric(substr(yop$variable, 3, 4))
    ## LINEAR INTERPOLATION OF DECENIAL TRENDS
    yop$value <- ifelse(yop$value> 10, 10, yop$value) ; LinApprox(yop)
    ## MERGE WITH THE GEOGRAPHIC FILES
    yop <- merge(DAT, CRD[, 6: 8], by= "MAILLE")
    coordinates(yop) <- c("X1", "X2")
    ## CREATE THE SPATIO-TEMPORAL ANIMATION
    AnimHTML(yop,
             ttle= paste(NDO$Anglais[ NDO$SPEC== i],
                         " Index, Scenario S0", sep= ""), pal= "Spectral",
             brks= c(-1000, -250, -75, -15, 0, 15, 75, 250, 1000),
             t1=paste("Scenario S0, ",
                      NDO$Anglais[ NDO$SPEC== i], " Index: ", sep= ""),
             t2=", ref= 2003", dir= paste("./output/", i, "S0/", sep= ""),
             name= paste(i , "S0.html", sep= ""))
}