Looking at number of cell, and expression data.

gast <- read.delim("../data/gastrula_geneBYcellsexpressed.txt", sep = "\t")

bla <- read.delim("../data/CPbla_geneBYcellsexpressed.txt", sep = "\t")

lc <- full_join(bla, gast, by = "id")

lc %>%
  filter(num_cells_expressed.x == 0, num_cells_expressed.y > 5)

lc %>%
  filter(num_cells_expressed.y == 0, num_cells_expressed.x > 5)

Using the given filters there are 1637 genes not expressed in blastula stage and expresses >5 cells in gastrula library.

And 260 genes expressed in greater that 5 cells in blastula library and not expressed in gastrula library.

Also have an expression comparison but not quite sure how to come up with a statistic.