interactions ggplot

interactions.R

npg <- 20
mu <- c(200,400,400,800,800,1600)
cond <- rep(rep(c("A","B"),each=npg),3)
geno <- rep(c("X","Y","Z"),each=2*npg)
table(cond, geno)
counts <- rnbinom(6*npg, mu=rep(mu,each=npg), size=1/.01)
d <- data.frame(counts, cond, geno)
library(ggplot2)
plotit <- function(d) {
  ggplot(d, aes(x=cond, y=counts, group=geno)) + 
    geom_jitter(position = position_jitter(width=.1)) +
    facet_wrap(~ geno) + 
    scale_y_log10() + 
    stat_summary(fun.y=mean, geom="line", colour="red", size=1)
}
plotit(d)
lm(log2(counts) ~ cond + geno + geno:cond)

mu[6] <- 400
counts <- rnbinom(6*npg, mu=rep(mu,each=npg), size=1/.01)
d <- data.frame(counts, cond, geno)
plotit(d)
lm(log2(counts) ~ cond + geno + geno:cond)