Code to reproduce the mean weight per species using the Zissou color palette (from #datacarpentry tweet)

zissouRodent.R

## Code to reproduce: https://twitter.com/FrancoisInvert/status/464493646151446528

library(wesanderson) # for the Zissou palette
library(RCurl) # needed to get data with HTTPS protocol

rodentURL <- getURL("https://raw.githubusercontent.com/NESCent/2014-05-08-datacarpentry/master/data/surveys.csv")
rodentData <- read.csv(textConnection(rodentURL))

meanPerSpecies <- aggregate(wgt ~ species, data=rodentData, mean, na.rm=TRUE)

zissouPal <- wes.palette(name = "Zissou", type="continuous")

png(file="zissouRodent.png", width=657, heigh=348)
par(mar=c(4, 2, 0, 1)) # if you want to see the legend
## par(mar=c(2, 2, 0, 1)) # for twitter (without legend)
barplot(sort(meanPerSpecies$wgt), horiz=TRUE, las=1,
        col=zissouPal(nlevels(meanPerSpecies$species[, drop=TRUE])),
        names.arg=meanPerSpecies$species, xlab="Mean weight (g)",
        cex.names=.7, cex.axis=.7)
dev.off()