/
IC_inc_data.Rmd
50 lines (43 loc) · 1.06 KB
/
IC_inc_data.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
49
---
title: "IC_inc_data"
author: "Tristan Caro"
date: "7/4/2021"
output: html_document
---
```{r}
# Load Libraries
library(tidyverse)
library(readxl)
library(ggsci)
library(ggdist)
library(magrittr)
```
```{r}
IC_tbl <- read_xlsx("ML_incubation_IC_data_EE.xlsx")
# Data represented in µM
# Cleanup the .xlsx
IC_tbl <- IC_tbl %>% na_if("n.a.")
IC_tbl$Ammonium %<>% as.double()
IC_tbl$Nitrate %<>% as.double()
IC_tbl$Formate %<>% as.double()
IC_tbl$`Time (hrs)` %<>% as.integer()
```
```{r}
IC_tbl %>% filter(N_source == "+ Nitrate") %>%
ggplot(aes(x = `Time (hrs)`, y = Nitrate)) +
geom_point(aes(color = ID)) +
geom_smooth(method = "lm") +
facet_wrap(vars(light)) +
ggtitle("Nitrate (µM)")
IC_tbl %>% filter(N_source == "+ Ammonium") %>%
ggplot(aes(x = `Time (hrs)`, y = Ammonium)) +
geom_point(aes(color = ID)) +
geom_smooth(method = "lm") +
facet_wrap(vars(light)) +
ggtitle("Ammonium (µM)")
IC_tbl %>% ggplot(aes(x = `Time (hrs)`, y = Formate)) +
geom_point() +
geom_smooth(method = "lm") +
facet_wrap(vars(light)) +
ggtitle("Formate (µM)")
```