Revision a1a88eb3
Added by Benoit Parmentier over 8 years ago
| climate/research/oregon/interpolation/global_product_assessment_part1.R | ||
|---|---|---|
| 4 | 4 |
#Combining tables and figures for individual runs for years and tiles. |
| 5 | 5 |
#AUTHOR: Benoit Parmentier |
| 6 | 6 |
#CREATED ON: 05/15/2016 |
| 7 |
#MODIFIED ON: 09/04/2016
|
|
| 7 |
#MODIFIED ON: 09/05/2016
|
|
| 8 | 8 |
#Version: 1 |
| 9 | 9 |
#PROJECT: Environmental Layers project |
| 10 | 10 |
#COMMENTS: Initial commit, script based on part NASA biodiversity conferenc |
| ... | ... | |
| 19 | 19 |
# |
| 20 | 20 |
#setfacl -Rmd user:aguzman4:rwx /nobackupp8/bparmen1/output_run10_1500x4500_global_analyses_pred_1992_10052015 |
| 21 | 21 |
|
| 22 |
#COMMIT: testing plot for both training and testing and introducing a new function
|
|
| 22 |
#COMMIT: testing plot function for training residuals at at stations and fixing bugs, region 1
|
|
| 23 | 23 |
|
| 24 | 24 |
################################################################################################# |
| 25 | 25 |
|
| ... | ... | |
| 85 | 85 |
source(file.path(script_path,function_assessment_part3)) #source all functions used in this script |
| 86 | 86 |
|
| 87 | 87 |
#Product assessment |
| 88 |
function_product_assessment_part1_functions <- "global_product_assessment_part1_functions_09042016.R"
|
|
| 88 |
function_product_assessment_part1_functions <- "global_product_assessment_part1_functions_09052016.R"
|
|
| 89 | 89 |
source(file.path(script_path,function_product_assessment_part1_functions)) #source all functions used in this script |
| 90 | 90 |
|
| 91 | 91 |
############################### |
| ... | ... | |
| 300 | 300 |
list_param_plot_stations_val_by_date <- list(l_dates,df_combined_data_v_dates,countries_shp_tmp,CRS_locs_WGS84, |
| 301 | 301 |
station_type_name,model_name, |
| 302 | 302 |
var_selected,y_var_name,out_suffix,out_dir) |
| 303 |
names(list_param_plot_stations_val_by_date) <- c("l_dates","df_combined_data_points","countries_shp","proj_str","
|
|
| 304 |
station_type_name","model_name", |
|
| 303 |
names(list_param_plot_stations_val_by_date) <- c("l_dates","df_combined_data_points","countries_shp","proj_str",
|
|
| 304 |
"station_type_name","model_name",
|
|
| 305 | 305 |
"var_selected","y_var_name","out_suffix","out_dir") |
| 306 | 306 |
|
| 307 | 307 |
|
| 308 |
list_plot_obj <- mclapply(1:length(l_dates),FUN=plot_stations_val_by_date,list_param=list_param_plot_stations_val_by_date, |
|
| 308 |
list_plot_obj_data_v <- mclapply(1:length(l_dates),FUN=plot_stations_val_by_date,list_param=list_param_plot_stations_val_by_date,
|
|
| 309 | 309 |
mc.preschedule=FALSE, |
| 310 | 310 |
mc.cores = num_cores) |
| 311 | 311 |
|
| ... | ... | |
| 314 | 314 |
list_param_plot_stations_val_by_date <- list(l_dates,df_combined_data_s_dates,countries_shp_tmp,CRS_locs_WGS84, |
| 315 | 315 |
station_type_name,model_name, |
| 316 | 316 |
var_selected,y_var_name,out_suffix,out_dir) |
| 317 |
names(list_param_plot_stations_val_by_date) <- c("l_dates","df_combined_data_points","countries_shp","proj_str","
|
|
| 318 |
station_type_name","model_name", |
|
| 317 |
names(list_param_plot_stations_val_by_date) <- c("l_dates","df_combined_data_points","countries_shp","proj_str",
|
|
| 318 |
"station_type_name","model_name",
|
|
| 319 | 319 |
"var_selected","y_var_name","out_suffix","out_dir") |
| 320 | 320 |
#debug(plot_stations_val_by_date) |
| 321 |
#test <- plot_stations_val_by_date(1,list_param = list_param_plot_stations_val_by_date)
|
|
| 321 |
#test <- plot_stations_val_by_date(2,list_param = list_param_plot_stations_val_by_date)
|
|
| 322 | 322 |
|
| 323 | 323 |
list_plot_obj_data_s <- mclapply(1:length(l_dates),FUN=plot_stations_val_by_date,list_param=list_param_plot_stations_val_by_date, |
| 324 | 324 |
mc.preschedule=FALSE, |
Also available in: Unified diff
product assessment part1 testing plot function for training residuals at at stations and fixing bugs, region 1