library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.2     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.2     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.1     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
qpcrd <- read.csv("../data/qPCR-data-2-21.csv", header = TRUE)
ggplot(data = qpcrd, mapping = aes(x = Subgroup, y = CqMean)) +   geom_boxplot() +
  facet_wrap('Target') +
  coord_cartesian(ylim = c(10, 40))

qpcrd %>%
  filter(Target == "Cg_HSP70(598/9)") %>%
ggplot(mapping = aes(x = Subgroup, y = CqMean)) + 
  geom_boxplot()  +
  geom_jitter(aes(color = History), width = 0.1, alpha = 0.7)

qpcrd %>%
  filter(Target == "Cg_HSP90(1532/3)") %>%
ggplot(mapping = aes(x = Subgroup, y = CqMean)) + 
  geom_boxplot()  +
  geom_jitter(aes(color = History), width = 0.1, alpha = 0.7)

qpcrd %>%
  filter(Target == "Cg_VIPERIN(1828/9)") %>%
ggplot(mapping = aes(x = Subgroup, y = CqMean)) + 
  geom_boxplot()  +
  geom_jitter(aes(color = History), width = 0.1, alpha = 0.7)