## 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", "SYLMEL") ; library(reshape)
for (i in levels(BirdS2$SPEC)[ !levels(BirdS2$SPEC) %in% Outliers]){
    ## SHAPE THE DATA AND DROP OUTILERS SPECIES
    sbs2 <- subset(BirdS2, BirdS2$SPEC== i)
    sbs2 <- aggregate(sbs2[, 3: 8], by= list(sbs2$MAILLE), mean)
    names(sbs2)[ 1] <-  "MAILLE" ; yop <- melt(sbs2, 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 S2", sep= ""), pal= "Spectral",
             brks= c(-1000, -250, -75, -15, 0, 15, 75, 250, 1000),
             t1=paste("Scenario S2, ",
                      NDO$Anglais[ NDO$SPEC== i], " Index: ", sep= ""),
             t2=", ref= 2003", dir= paste("./output/", i, "S2", sep= ""),
             name= paste(i , "S2.html", sep= ""))
}