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