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