Misc sketches

The gaol of integration

n_cells <- 600
expression <- data.frame(y=rpois(n_cells, 2), cell_type=rep(c("B cells", "T cells", "NK cells"), n_cells/3), Dataset=rep(c("A", "B", "C"), each=n_cells/3))
expression$y[expression$cell_type == "T cells"] <- expression$y[expression$cell_type == "T cells"] + 2
expression$y[expression$cell_type == "NK cells"] <- expression$y[expression$cell_type == "NK cells"] + 8

expression_norm <- expression
expression$y[expression$Dataset == "A"] <- expression$y[expression$Dataset == "A"] + 20
expression$y[expression$Dataset == "B"] <- expression$y[expression$Dataset == "B"] + 30


ggplot(expression, aes(x=cell_type, y=y, colour=Dataset)) +
  xlab("Cell type") + ylab("Expression") +
  labs(subtitle = "Gene X") +
  geom_jitter() +
  theme_classic() +
  theme(axis.text.y = element_blank(), axis.ticks = element_blank())

ggplot(expression_norm, aes(x=cell_type, y=y, colour=Dataset)) +
  xlab("Cell type") + ylab("Expression") +
  labs(subtitle = "Gene X") +
  geom_jitter() +
  theme_classic() +
  theme(axis.text.y = element_blank(), axis.ticks = element_blank())