We have an analysis routine for proteomics data written for IgorPro. One output is a volcano plot . These plots show the fold change in one sample compared to another and plot that against a p-value to estimate how reproducible any changes observed are. This post is not about that software, but on the topic of how we can recreate this plot in R . What steps need to be considered?