End Result


Code

library(scales) #to use the alpha function
a=rnorm(250 , mean=0 , sd=4)
b=rnorm(250 , mean=0 , sd=4)
c=rnorm(250 , mean=1 , sd=8)
par(bg="black")
par(mar=c(0,0,0,0))
# Hex color #003366 is midnight blue
#cex maginfies the shape diameter
#pch sets the shape 16=circle, 14=triangle,
#xaxt, yaxt, xlab, ylab, xlim, ylim all control the appearnce and scaling of the axis
plot(a ,b , cex=abs(c) ,col=alpha("#003366",0.25), pch=16 , xaxt="n" , yaxt="n" , xlab="" , 
     ylab="" , xlim=c(-3,10) , ylim=c(-3,10))