Figure 1.
An iterative peak calling procedure yields more peaks from the complex mix of worm cell types. (A) The core peak calling procedure is to model the data using latent Dirichlet allocation (LDA), cluster the cells, and call peaks based on the clusters. (B) A flow chart represents the overall peak-calling strategy. First, bulk peaks are called, followed by two iterations of clustering and peak calling based on an LDA model. Then, we group the cells by tissue and repeat the two steps of clustering and peak calling. The number of cells included and the number of peaks called at each step are given in the inset table.
