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