PLOS Biology 2010 Citation Counts

# code for figure 1: density plots for citation counts for
# PLOS Biology articles published in 2010

# load May 20, 2013 ALM report
alm <- read.csv("../data/alm_report_plos_biology_2013-05-20.csv", stringsAsFactors=FALSE)

# only look at research articles
alm <- subset(alm, alm$article_type == "Research Article")

# only look at papers published in 2010
alm$publication_date <- as.Date(alm$publication_date)
alm <- subset(alm, alm$publication_date > "2010-01-01" & alm$publication_date <= "2010-12-31")

# labels
colnames <- dimnames(alm)[[2]]
plos.color <- "#1ebd21"
plos.source <- "scopus"

plos.xlab <- "Scopus Citations"
plos.ylab <- "Probability"

quantile <- quantile(alm[,plos.source], c(.10,.50,.90), na.rm=TRUE)

# plot the chart
opar <- par(mai=c(0.5,.75,0.5,0.5), omi=c(0.25,0.1,0.25,0.1), mgp=c(3,.5,.5), fg="black", cex.main=2, cex.lab=1.5, col=plos.color, col.main=plos.color, col.lab=plos.color, xaxs="i", yaxs="i")

d <- density(alm[,plos.source], from=0, to=100)
d$x <- append(d$x, 0)
d$y <- append(d$y, 0)
plot(d, type="n", main=NA, xlab=NA, ylab=NA, xlim=c(0,100), frame.plot=FALSE)
polygon(d, col=plos.color, border=NA)
mtext(plos.xlab, side=1, col=plos.color, cex=1.25, outer=TRUE, adj=1, at=1)
mtext(plos.ylab, side=2, col=plos.color, cex=1.25, outer=TRUE, adj=0, at=1, las=1)

par(opar)