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