add shape files

OJWatson · Feb 2, 2024 · fb89b07 · fb89b07
1 parent 0fc8e26
commit fb89b07
Show file tree

Hide file tree

Showing 12 changed files with 14 additions and 0 deletions.
diff --git a/analysis/11_simulating_spread.R b/analysis/11_simulating_spread.R
@@ -382,3 +382,11 @@ write.csv(central_df, "analysis/data-out/central_times_prospective.csv", row.nam
 write.csv(worst_df, "analysis/data-out/pessimistic_times_prospective.csv", row.names = FALSE)
 write.csv(best_df, "analysis/data-out/optimistic_times_prospective.csv", row.names = FALSE)
 
+# and out put the shape files
+dir.create("analysis/data-out/shape_files")
+afr_map_full <- readRDS(here::here("analysis/data-derived/map_for_CAR.rds"))
+saveRDS(afr_map_full %>% select(iso:country_level), "analysis/data-out/shape_files/admin1.rds")
+sf::write_sf(afr_map_full %>% select(iso:country_level), "analysis/data-out/shape_files/admin1.shp")
+saveRDS(map_0, "analysis/data-out/shape_files/admin0.rds")
+sf::write_sf(map_0, "analysis/data-out/shape_files/admin0.shp")
+
diff --git a/analysis/data-out/README.md b/analysis/data-out/README.md
@@ -26,3 +26,7 @@ In all the outputs above, the **central** times provide the estimated times base
 The **optimistic** times assume the upper or lower estimate value for each parameter (depending on the direction of its effect on selection) such that the selection of ArtR will increase at its slowest. Conversely, the **pessimistic** times assume the upper or lower value for each parameter (depending on the direction of its effect on selection) such that the selection of ArtR will increase at its fastest.
 
 Lastly, the uncertainty provided in each of the files, reflects the uncertainty that arises from our use of a stochastic, individual based model for the selection of ArtR. 
+
+### Shape Files
+
+The linked shape files for these outputs are in `shape_files/`. The admin level 0 and admin level 1 shape files are provided, both as a saved R object (`.rds`) and as a shape file (`.shp`). To link these with the risk outputs, these can be matched using `id_1`.
diff --git a/analysis/data-out/shape_files/admin0.dbf b/analysis/data-out/shape_files/admin0.dbf
diff --git a/analysis/data-out/shape_files/admin0.prj b/analysis/data-out/shape_files/admin0.prj
@@ -0,0 +1 @@
+GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
diff --git a/analysis/data-out/shape_files/admin0.rds b/analysis/data-out/shape_files/admin0.rds
diff --git a/analysis/data-out/shape_files/admin0.shp b/analysis/data-out/shape_files/admin0.shp
diff --git a/analysis/data-out/shape_files/admin0.shx b/analysis/data-out/shape_files/admin0.shx
diff --git a/analysis/data-out/shape_files/admin1.dbf b/analysis/data-out/shape_files/admin1.dbf
diff --git a/analysis/data-out/shape_files/admin1.prj b/analysis/data-out/shape_files/admin1.prj
@@ -0,0 +1 @@
+GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
diff --git a/analysis/data-out/shape_files/admin1.rds b/analysis/data-out/shape_files/admin1.rds
diff --git a/analysis/data-out/shape_files/admin1.shp b/analysis/data-out/shape_files/admin1.shp
diff --git a/analysis/data-out/shape_files/admin1.shx b/analysis/data-out/shape_files/admin1.shx
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]