knitr::opts_chunk$set(echo = TRUE)
knitr::opts_chunk$set(message = FALSE)
knitr::opts_chunk$set(warning = FALSE)
library(ProvidenciaChemo)
library(tidyverse)
colonies <- read_csv('extdata/Figure_1E_cfu.csv') %>%
  mutate(supe_total = n_sup * 10, # plated 1/10 volume of last wash after bleaching to get n cells left in supe
         cells_per_worm = ( (n_colonies * dilut) - supe_total ) / n_worms) # normalize based on # cells left over after bleaching

colonies %>% ggplot(aes(x = food, y = cells_per_worm)) +
  geom_bardots(fillvar = food) +
  scale_fill_plot(palette = "grey-blue", drop = TRUE) +
  labs(y = "cfu / worm")


mikeod38/ProvidenciaChemo documentation built on April 6, 2020, 11:57 p.m.