---
title: "02-Crab-qpcr"
output: github_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

```{r}
library("dplyr")
```

```{r}
qpcr <- read.csv("../data/hematqpcr_crabRNA.csv")

head(qpcr)
```
Want to get average sqallmeans for each FRP...

```{r}
avQpcr <- group_by(qpcr, FRP) %>%
summarize(avg_sq = mean(sq_all.runs_mean, na.rm=TRUE), trtmnt_tank, infection_status, maturity, sample_pos.neg) %>% 
distinct_all() %>%
filter(trtmnt_tank == "cold" | trtmnt_tank == "warm" | trtmnt_tank == "ambient")


avQpcr
```

```{r}
library(ggplot2)
```

```{r}
ggplot(avQpcr, aes(x = avg_sq)) +
  geom_histogram(bins = 100) +
  scale_x_log10() +
  facet_grid(~infection_status)
```


