I recently gathered fish harvest data from the U.S. National Oceanic and Atmospheric Administarion (NOAA), which I downloaded from Infochimps. The data is fish harvest by weight and value, by species for 21 years, from 1985 to 2005. Here is a link to a google document of the data I used below: https://spreadsheets.google.com/ccc?key=0Aq6aW8n11tS_dFRySXQzYkppLXFaU2F5aC04d19ZS0E&hl=en

################# Fish harvest data ########################################

setwd("/Mac/R_stuff/Blog_etc/Infochimps/Fishharvest") # Set path

library(ggplot2)

library(googleVis)

library(Hmisc)

fish <- read.csv("fishharvest.csv") # read data

fish2 <- melt(fish,id=1:3,measure=4:24) # melt table

year <- rep(1985:2005, each = 117)

fish2 <- data.frame(fish2,year) # replace year with actual values

# Google visusalization API

fishdata <- data.frame(subset(fish2,fish2$var == "quantity_1000lbs",-4),value_1000dollars=subset(fish2,fish2$var == "value_1000dollars",-4)[,4])

names(fishdata)[4] <- "quantity_1000lbs"

fishharvest <- gvisMotionChart(fishdata, idvar="species", timevar="year")

plot(fishharvest)

Created by Pretty R at inside-R.org

Data: fishdata, Chart ID: MotionChart_2011-01-17-08-09-24

R version 2.12.1 (2010-12-16),

Google Terms of Use

fishdatagg2 <- ddply(fish2,.(species,var),summarise,

mean = mean(value),

se = sd(value)/sqrt(length(value))

)

fishdatagg2 <- subset(fishdatagg2,fishdatagg2$var %in% c("quantity_1000lbs","value_1000dollars"))

limit3 <- aes(ymax = mean + se, ymin = mean - se)

bysppfgrid <- ggplot(fishdatagg2,aes(x=reorder(species,rank(mean)),y=mean,colour=species)) + geom_point() + geom_errorbar(limit3) + facet_grid(. ~ var, scales="free") + opts(legend.position="none") + coord_flip() + scale_y_continuous(trans="log")

ggsave("bysppfgrid.jpeg")

Created by Pretty R at inside-R.org

*Related*

To

**leave a comment** for the author, please follow the link and comment on their blog:

** Recology**.

R-bloggers.com offers

**daily e-mail updates** about

R news and

tutorials on topics such as:

Data science,

Big Data, R jobs, visualization (

ggplot2,

Boxplots,

maps,

animation), programming (

RStudio,

Sweave,

LaTeX,

SQL,

Eclipse,

git,

hadoop,

Web Scraping) statistics (

regression,

PCA,

time series,

trading) and more...

If you got this far, why not

__subscribe for updates__ from the site? Choose your flavor:

e-mail,

twitter,

RSS, or

facebook...

**Tags:** datasets, Ecology, ggplot2, googleVis, R