![](https://i0.wp.com/quantixed.org/wp-content/uploads/2023/06/volcanoPlot_igor_q.png?resize=364%2C512&ssl=1)
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?