def chip_ps_ref(ps,mean_frag_length,cells=10000): """Do a chip seq experiment given the distribution ps""" G = len(ps) return concat(chip_ps(rfd_xs(ps),mean_frag_length) for cell in verbose_gen(xrange(cells)))
def make_chip_dataset(num_cells): return concat([chip(genome,rfd_xs(ps),MEAN_FRAGMENT_LENGTH) for i in verbose_gen(xrange(num_cells))])