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))])
