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