#Analyses, figures, tables and data for the  paper are also produced in the script.

#AUTHOR: Benoit Parmentier 

#CREATED ON: 10/31/2013  

#Version: 6

#PROJECT: Environmental Layers project                                     

#################################################################################################

                              par.main.text=list(font=2,cex=2),strip.background=list(col="white")),par.strip.text=list(font=2,cex=1.5),

          #strip=strip.custom(factor.levels=names_layers),

          xlab=list(label="Month", cex=2,font=2),

          auto.key=list(columns=2,space="top",cex=2.5,font=2),

          #list(label="\u0394RMSE between mod1 and mod4",cex=1.5)

p1 <- list_plots_spt[[1]]

p2 <- list_plots_spt[[2]]

#p3 <- list_plots_spt[[3]]

#grid.arrange(p1,p2,p3,ncol=1)

grid.arrange(p1,p2,ncol=1)

################################################

#### Figure 5: Spatial lag profiles  

#for Jan 1 and Sept 1 in 2010 for all models (1 to 7) and methods

index <- 1 #index corresponding to Jan 1 #For now create Moran's I for only one date...

nb_lag <- 10

list_filters<-lapply(1:nb_lag,FUN=autocor_filter_fun,f_type="queen") #generate lag 10 filters

list_param_stat_moran_gam_daily <- list(filter=list_filters[[10]],lf_list=list_lf_gam_daily)

#tt <- stat_moran_std_raster_fun(1,list_param=list_param_stat_moran)

#tt <- mclapply(1:11, list_param=list_param_stat_moran_CAI, FUN=stat_moran_std_raster_fun,mc.preschedule=FALSE,mc.cores = 11) #This is the end bracket from mclapply(...) statement

png(filename=paste("Figure_annex1_paper_Comparison_daily_monthly_mean_lst",out_prefix,".png",sep=""),width=960,height=480)

par(mfrow=c(1,2))

plot(lst_md)

title("Mean temperature (°C) for January 1")

plot(lst_mm_01)

title("Mean temperature (°C) for month of January")

dev.off()