In WendyLiuLab/Hsieh2016PCR: PCR data analysis for Hsieh 2006

knitr::opts_chunk$set(echo = FALSE)

library(dplyr)
library(ggplot2)
library(readr)
library(reshape2)

devtools::load_all(".")

cqs = load_all_cqs()
melts = load_all_melt_curves()

for(target in unique(cqs$Target)) {
  cat(paste0("\n\n", target, " --------------------------"))
  target.cqs = cqs %>% filter(Target == target)
  target.melts = melts %>% filter(Target == target)

  g = ggplot(target.melts, aes(Temperature, dF, color=as.factor(Col))) +
    facet_wrap(~Condition) +
    geom_line(alpha=0.7) +
    theme_bw()
  print(g)

  g = ggplot(target.cqs, aes(Condition, Cq, color=Mouse)) +
    geom_point() +
    facet_grid(Timepoint~.) +
    theme_bw()
  print(g)

  g = target.cqs %>%
    fold_vs_gapdh %>%
      ggplot(aes(Condition, FoldVsGapdh, color=Mouse)) +
      geom_point() +
      facet_grid(Timepoint~., scales="free_y") +
      theme_bw()
  print(g)
}