Advances in peptide microarray technology have allowed for the creation of fast-paced and modular experiments within affinity ligand discovery. Previously, low density peptide arrays of 10,000 peptides were used to identify low affinity peptide ligands for a target protein; an approach that can be subsequently improved upon with a number of techniques. VDAP[a] offers more information about the relative affinity of protein-peptide interactions via signal intensity in contrast to high throughput screening (HTS) and display technologies which offer binary data. Now, high density peptide arrays with 130,000 to 330,000 peptides are available that allow screening across peptide libraries of greater diversity. With this increase in scale and diversity, faster analytical tools are needed to adequately characterize array data. Using the statistical power available in the R programming language, we have created a flexible analysis package that efficiently processes high density peptide array data from a variety of layouts, rank existing peptide hits, and utilize signal intensity data to generate new hits. This analysis provides a user-friendly method to efficiently analyze high density peptide array data, generate peptide leads for targeted therapeutic development, and further improve peptide array technologies.