Species abundance curves and histograms

# Code for these Graphs # David LeBauer dlebauer@uci.edu / @gmail.com setwd ( "/media/sda3/Users/David/Research/NCEAS_DGS/Meeting" ) sp.i.abun <- read.csv ( 'spp.i.abun.csv', header = TRUE) attach ( sp.i.abun ) # Picture of Curves #pdf( png ( file = "sppabun.png" ) #, height = 11 ,width= 8.5); # par ( omi = c(1,1,1,1)) #NB remember to turn dev.off() when done writing figure # X11( ) <- use X11 for real figures # par( mfrow = (c(2,2))) plot( rank , 100*sum_abun, xlab = 'rank abundance of sp_i', ylab = 'cumulative % abundance' ) plot( rank, 100*rel.abun, xlab = 'rank abundance of sp_i', ylab = 'absolute % abundance' ) hist ( number_individuals , breaks= 25, ylab = 'number of spp. in bin' ) hist (log10( number_individuals ), breaks= 25, ylab = 'number of spp. in bin' ) dev.off()

Click to view full-size image…