![enter image description here][2]res_lumbar <-read.csv("res_lumbar_results_new.csv",header=TRUE)

sig_lumbar<- subset(res_lumbar, padj < 0.1 & abs(log2FoldChange) >= 1)

up_lumbar <- subset(sig_lumbar, log2FoldChange >= 1)

down_lumbar<- subset(sig_lumbar, log2FoldChange <= -1)

cat("Lumbar Upregulated:", nrow(up_lumbar), "\n")

cat("Lumbar Downregulated:", nrow(down_lumbar), "\n")

EnhancedVolcano(res_lumbar, lab = NA, x = "log2FoldChange", y = "padj", title = "Volcano Plot: Thoracic CT vs IT", xlab = "Log2 Fold Change", ylab = "-log10(p-value)", pCutoff = 0.05, FCcutoff = 1, pointSize = 1, colAlpha = 0.75, legendPosition = "right", legendLabSize = 12, legendIconSize = 4, drawConnectors = FALSE, ylim=c(0,130) ) This is the code i have used. While plotting the Volcano plots, due to few outliers, my whole plot is getting skewed, can anyone please suggest me how to include nreaks in the y-axis. Any help is appreciated