data <- read.table("gos.fn.GS020.rarefaction")
attach(data)
plot(numsampled, unique, type="b", xlab="# sequences sampled", ylab="Observed OTUs", col="blue")
lines(numsampled, X0.03, type="b", col="green")
lines(numsampled, X0.05, type="b", col="red")
lines(numsampled, X0.10, type="b", col="magenta")
legend("topleft",c('Unique','97%','95%', '90%'),
   pch=c(16,16,16,16), lty=c(1,1,1,1),
  col=c('blue','green','red','magenta')
  ,ncol=1)

http://www.plob.org/2014/05/11/7264.html

http://docs.ggplot2.org/current/index.html

http://www.tuicool.com/articles/nUf2ae