Revision c1b1cb13
Added by Benoit Parmentier almost 10 years ago
| climate/research/oregon/interpolation/multi_timescale_paper_interpolation.R | ||
|---|---|---|
| 7 | 7 |
#Analyses, figures, tables and data for the paper are also produced in the script. |
| 8 | 8 |
#AUTHOR: Benoit Parmentier |
| 9 | 9 |
#CREATED ON: 10/31/2013 |
| 10 |
#MODIFIED ON: 10/23/2014
|
|
| 10 |
#MODIFIED ON: 01/06/2015
|
|
| 11 | 11 |
#Version: 6 |
| 12 | 12 |
#PROJECT: Environmental Layers project |
| 13 | 13 |
################################################################################################# |
| ... | ... | |
| 111 | 111 |
"gam_fss","kriging_fss","gwr_fss") |
| 112 | 112 |
|
| 113 | 113 |
y_var_name <- "dailyTmax" |
| 114 |
out_prefix<-"analyses_10062014"
|
|
| 114 |
out_prefix<-"analyses_01062015"
|
|
| 115 | 115 |
out_dir<-"/home/parmentier/Data/IPLANT_project/paper_multitime_scale__analyses_tables" |
| 116 | 116 |
create_out_dir_param = TRUE |
| 117 | 117 |
|
| ... | ... | |
| 430 | 430 |
height=500*layout_m[1],width=500*layout_m[2]) |
| 431 | 431 |
|
| 432 | 432 |
p<-xyplot(data ~ month | interp_method, groups=which,dd,type="b", |
| 433 |
par.settings = list(axis.text = list(font = 2, cex = 1.3),layout=layout_m, |
|
| 434 |
par.main.text=list(font=2,cex=2),strip.background=list(col="white")),par.strip.text=list(font=2,cex=1.5), |
|
| 433 |
par.settings = list(axis.text = list(font = 2, cex = 1.3), |
|
| 434 |
superpose.symbol = list( #allow changing symbol size, type and color |
|
| 435 |
pch=1:2,cex=1.5), |
|
| 436 |
layout=layout_m, |
|
| 437 |
par.main.text=list(font=2,cex=2),strip.background=list(col="white")), |
|
| 438 |
par.strip.text=list(font=2,cex=1.5), |
|
| 435 | 439 |
#strip=strip.custom(factor.levels=names_layers), |
| 436 | 440 |
xlab=list(label="Month", cex=2,font=2), |
| 437 | 441 |
ylab=list(label="\u0394RMSE from model with elevation", cex=2, font=2), |
| 438 |
auto.key=list(columns=2,space="top",cex=2.5,font=2), |
|
| 442 |
#key = list(points=list(pch=1:2)), |
|
| 443 |
|
|
| 444 |
auto.key=list(columns=2,space="top",cex=2.5,font=2,pch=1:2), |
|
| 439 | 445 |
|
| 440 | 446 |
#list(label="\u0394RMSE between mod1 and mod4",cex=1.5) |
| 441 | 447 |
as.table=TRUE) #as.table controls the order of the pannels!! |
| ... | ... | |
| 445 | 451 |
}) |
| 446 | 452 |
print(p) |
| 447 | 453 |
|
| 448 |
#p_dif <- xyplot(rmse_dif ~ 1:12, |
|
| 449 |
# col=c("black"),
|
|
| 450 |
# type="b", |
|
| 451 |
# ylab=list(label="\u0394RMSE between mod1 and mod4",cex=1.5), |
|
| 452 |
# xlab=list(label="Month",cex=1.5), |
|
| 453 |
# main=list(label="(b) Monthly \u0394RMSE betwen mod1 and mod4",cex=1.8), |
|
| 454 |
# par.settings = list(axis.text = list(font = 2, cex = 1.3), |
|
| 455 |
# par.main.text=list(font=2,cex=2),strip.background=list(col="white")),par.strip.text=list(font=2,cex=1.5)) |
|
| 456 |
#p_dif <- update(p_dif, panel = function(...) {
|
|
| 457 |
# panel.abline(h = 0, lty = 2, col = "gray") |
|
| 458 |
# panel.xyplot(...) |
|
| 459 |
# }) |
|
| 460 |
|
|
| 461 | 454 |
dev.off() |
| 462 | 455 |
|
| 463 | 456 |
####################################################### |
Also available in: Unified diff
paper proof multi-timescale Figure 11 modifications