Ultradense peptide binding arrays that can probe millions of linear peptides comprising the entire proteomes or immunomes of human or mouse, or numerous microbes, are powerful tools for studying the abundance of different antibody repertoire in serum samples to understand adaptive immune responses. There are few statistical analysis tools for exploring high-dimensional, significant and reproducible antibody targets for ultradense peptide binding arrays at the linear peptide, epitope (grouping of adjacent peptides), and protein level across multiple samples/subjects (I.e. epitope spread or immunogenic regions within each protein) for understanding the heterogeneity of immune responses. We developed HERON (Hierarchical antibody bindingEpitopes and pROteins from liNear peptides), an R package, which allows users to identify immunogenic epitopes using meta-analyses and spatial clustering techniques to explore antibody targets at various resolution and confidence levels, that can be found consistently across a specified number of samples through the entire proteome to study antibody responses for diagnostics or treatment. Our approach estimates significance values at the linear peptide (probe), epitope, and protein level to identify top candidates for validation. We test the performance of predictions on all three levels using correlation between technical replicates and comparison of epitope calls on 2 datasets, which shows HERON’s competitiveness in estimating false discovery rates and finding general and sample-level regions of interest for antibody binding. The code is available as an R package downloadable fromhttp://github.com/Ong-Research/HERON.