# volcano plot:
plot(M,Bstat,pch=".")
abline(h=c(0,2,4),lty=c(1,2,3))
abline(v=c(-1,1),lty=3)
#
bigM <- (abs(M) > 1)
bigB <- (Bstat > 0)
points(M[bigM&!bigB],Bstat[bigM&!bigB],pch=16,col="green")
points(M[bigB&!bigM],Bstat[bigB&!bigM],pch=16,col="magenta")
points(M[bigM&bigB],Bstat[bigM&bigB],pch=16,col="blue")
#
#
# It looks like this: