forked from FEProgramGhana/FE-2024-Ghana-examples
-
Notifications
You must be signed in to change notification settings - Fork 0
/
spatial_plotter.Rmd
48 lines (40 loc) · 1.36 KB
/
spatial_plotter.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
---
title: "Plot Spatial Example Exercise"
author: "Tobias Holden"
date: "6/14/2023"
output: html_document
editor_options:
chunk_output_type: console
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(data.table)
```
```{r}
# Input path (yours may differ)
root <- paste('experiments','my_outputs','spatial',sep='/')
### Load Files
events <- fread(paste(root,"CountedEvents.csv", sep='/'))
malaria_channels <- fread(paste(root,"SpatialReportMalariaFiltered_all_ages.csv",sep='/'))
### Align Column Names
colnames(events) <- tolower(colnames(events))
colnames(malaria_channels) <- tolower(colnames(malaria_channels))
### Join datasets by time and node
df <- inner_join(events, malaria_channels)
### Plot Summary
df %>%
gather(key = "variable", value="value", -c(node, time, run_number, x_temporary_larval_habitat, experiment)) %>%
group_by(time, node, x_temporary_larval_habitat, variable) %>%
summarize(value = mean(value)) %>%
mutate(node = paste("Node",node),
variable=gsub("_","\n",variable)) %>%
ggplot(aes(x=time, y=value)) +
geom_path(aes(group=x_temporary_larval_habitat, color=x_temporary_larval_habitat)) +
facet_grid(variable~node, scales="free") +
scale_color_distiller(palette = "Spectral") +
ylab("Average Value") +
theme_minimal()
### Save plot
ggsave(paste(root,'SpatialSummary.png',sep = '/'))
```