tT<-read.csv(file.choose(),stringsAsFactors = F) plot(tT$logFC,-log10(tT$P.Value),pch=16,main="Volcano map",cex.lab=1.5, xlab="log2FoldChange",ylab="-log10Pvalue") points(tT[tT$logFC>1&tT$P.Value<0.05,]$logFC,pch=16, -log10(tT[tT$logFC>1&tT$P.Value<0.05,]$P.Value), col='red') points(tT[tT$logFC< -1&tT$P.Value<0.05,]$logFC,pch=16, -log10(tT[tT$logFC< -1&tT$P.Value<0.05,]$P.Value), col='green') abline(h=-log10(0.05),v=c(-1,1))