Revision 893b8668
Added by Benoit Parmentier over 9 years ago
| climate/research/oregon/interpolation/presentation_global_run_scalingup_assessment.Rmd | ||
|---|---|---|
| 1 | 1 |
--- |
| 2 |
title : Summary for year 1992, region 4 -Africa
|
|
| 2 |
title : Summary for year 1982, region 4 -Africa
|
|
| 3 | 3 |
subtitle : Global assessment interpolation |
| 4 | 4 |
author : Benoit Parmentier and Alberto Guzman |
| 5 | 5 |
job : NEX job world |
| ... | ... | |
| 15 | 15 |
|
| 16 | 16 |
#First let's set up the input parameters for the script in R: |
| 17 | 17 |
|
| 18 |
in_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1992_09012015"
|
|
| 19 |
|
|
| 18 |
#in_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pr#ed_1982_09152015"
|
|
| 19 |
in_dir <- "/home/bparmentier/Dropbox/Data/NCEAS/env_layers_documents/output_run10_1500x4500_global_analyses_pred_1982_09152015" |
|
| 20 | 20 |
#region_name <- "reg4" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3 |
| 21 | 21 |
#out_suffix <- paste(region_name,"_","pred_1992_09032015",sep="") |
| 22 |
out_suffix <- "pred_1992_09032015"
|
|
| 22 |
out_suffix <- "pred_1982_09152015"
|
|
| 23 | 23 |
#day_to_mosaic <- c("20100831",
|
| 24 | 24 |
# "20100901") #PARAM7 |
| 25 | 25 |
|
| ... | ... | |
| 31 | 31 |
There were 43 tiles for region 4 (Africa) |
| 32 | 32 |
|
| 33 | 33 |
1. Tiles were run for 365 days. |
| 34 |
2. Year tested was 1992.
|
|
| 34 |
2. Year tested was 1982.
|
|
| 35 | 35 |
3. Assessments for the number of days predicted and accuracy were carried out. |
| 36 | 36 |
|
| 37 | 37 |
|
| ... | ... | |
| 51 | 51 |
```{r, echo=FALSE}
|
| 52 | 52 |
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015") |
| 53 | 53 |
|
| 54 |
fname <- paste("Figure1_tile_processed_region_world_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 54 |
fname <- paste("Figure1_tile_processed_region_world_run10_1500x4500_global_analyses_",out_suffix)
|
|
| 55 | 55 |
|
| 56 | 56 |
image_file1 <- file.path(in_dir,fname) |
| 57 | 57 |
|
| ... | ... | |
| 69 | 69 |
|
| 70 | 70 |
```{r, echo=FALSE}
|
| 71 | 71 |
|
| 72 |
fname <- paste("Figure2a_boxplot_with_oultiers_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 72 |
fname <- paste("Figure2a_boxplot_with_oultiers_by_tiles_mod1_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 73 | 73 |
image_file2 <- file.path(in_dir,fname) |
| 74 | 74 |
|
| 75 |
fname <- paste("Figure2b_boxplot_scaling_by_tiles_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 75 |
fname <- paste("Figure2b_boxplot_scaling_by_tiles_mod1_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 76 | 76 |
image_file3 <- file.path(in_dir,fname) |
| 77 | 77 |
|
| 78 | 78 |
``` |
| ... | ... | |
| 96 | 96 |
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png") |
| 97 | 97 |
|
| 98 | 98 |
|
| 99 |
fname <- paste("Figure3a_boxplot_overall_region_with_oultiers_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 99 |
fname <- paste("Figure3a_boxplot_overall_region_with_oultiers_mod_kr_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 100 | 100 |
image_file4 <- file.path(in_dir,fname) |
| 101 | 101 |
|
| 102 |
fname <- paste("Figure3b_boxplot_overall_region_scaling_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 102 |
fname <- paste("Figure3b_boxplot_overall_region_scaling_mod_kr_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 103 | 103 |
image_file5 <- file.path(in_dir,fname) |
| 104 | 104 |
|
| 105 | 105 |
``` |
| ... | ... | |
| 126 | 126 |
#image_file2 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_edge_weighted_20100831_world_mosaic_07092015.png") |
| 127 | 127 |
|
| 128 | 128 |
|
| 129 |
fname <- paste("Figure5_ac_metrics_ranked_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 129 |
fname <- paste("Figure5_ac_metrics_ranked_mod1_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 130 | 130 |
image_file6 <- file.path(in_dir,fname) |
| 131 | 131 |
|
| 132 |
fname <- paste("Figure5_ac_metrics_ranked_mod_kr_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 132 |
fname <- paste("Figure5_ac_metrics_ranked_mod_kr_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 133 | 133 |
image_file7 <- file.path(in_dir,fname) |
| 134 | 134 |
|
| 135 | 135 |
``` |
| ... | ... | |
| 155 | 155 |
```{r, echo=FALSE}
|
| 156 | 156 |
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015") |
| 157 | 157 |
|
| 158 |
fname <- paste("Figure6_ac_metrics_map_centroids_tile_mod1_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 158 |
fname <- paste("Figure6_ac_metrics_map_centroids_tile_mod1_run10_1500x4500_global_analyses_",out_suffix, sep="")
|
|
| 159 | 159 |
|
| 160 | 160 |
image_file8 <- file.path(in_dir,fname) |
| 161 | 161 |
|
| ... | ... | |
| 174 | 174 |
|
| 175 | 175 |
--- |
| 176 | 176 |
|
| 177 |
# Predicted tiles for 1992
|
|
| 177 |
# Predicted tiles for 1982
|
|
| 178 | 178 |
|
| 179 | 179 |
```{r, echo=FALSE}
|
| 180 | 180 |
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015") |
| 181 | 181 |
|
| 182 |
fname <- paste("Figure7a_ac_metrics_map_centroids_tile_mod1_missing_day_367_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 182 |
fname <- paste("Figure7a_ac_metrics_map_centroids_tile_mod1_missing_day_367_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 183 | 183 |
|
| 184 | 184 |
image_file9 <- file.path(in_dir,fname) |
| 185 | 185 |
|
| ... | ... | |
| 203 | 203 |
```{r, echo=FALSE}
|
| 204 | 204 |
#image_file1 <- file.path(in_dir,"Figure9_clim_mosaics_day_test_unweighted_20100831_world_mosaic_07092015") |
| 205 | 205 |
|
| 206 |
fname <- paste("Figure7a_ac_metrics_map_centroids_tile_mod1_missing_day_365_run10_1500x4500_global_analyses_pred_1992_09012015.png")
|
|
| 206 |
fname <- paste("Figure7a_ac_metrics_map_centroids_tile_mod1_missing_day_365_run10_1500x4500_global_analyses_",out_suffix,sep="")
|
|
| 207 | 207 |
|
| 208 | 208 |
image_file10 <- file.path(in_dir,fname) |
| 209 | 209 |
|
| ... | ... | |
| 226 | 226 |
|
| 227 | 227 |
```{r echo = F}
|
| 228 | 228 |
|
| 229 |
#tb <- read.table(file=file.path(out_dir,paste("tb_diagnostic_v_NA","_",out_suffix,".txt",sep="")),sep=",")
|
|
| 230 |
tb <- read.table(file.path("/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1992_09012015","tb_diagnostic_v_NA_run10_1500x4500_global_analyses_pred_1992_09012015.txt"),sep=",")
|
|
| 229 |
#tb <- read.table(file.path(in_dir,"tb_diagnostic_v_NA_run10_1500x4500_global_analyses_pred#_1992_09012015.txt"),sep=",") |
|
| 230 |
|
|
| 231 |
tb_s <- read.table(file=file.path(in_dir,paste("tb_diagnostic_v_NA","_",out_suffix,".txt",sep="")),sep=",")
|
|
| 232 |
|
|
| 231 | 233 |
|
| 232 | 234 |
#table(tb$pred_mod) |
| 233 | 235 |
#table(tb$index_d) |
Also available in: Unified diff
presentation slidify for 1982 first modifications