#original code used to obtain this dataset on March 21st, 2015
# using version 1.2 of the Paleobiology Database API
# (sorry, URLs removed as they lead to the PBDB test server...)
save(graptOccPBDB,graptTaxaPBDB,file="graptPBDB.rdata")
# load archived example data
data(graptPBDB)
# let's visualize who entered the majority of the occurrence data
pie(sort(table(graptOccPBDB$enterer)))
# and now who authorized it
pie(sort(table(graptOccPBDB$authorizer)))
# I apologize for using pie charts.
# Let's look at age resolution of these occurrences
hist(graptOccPBDB$early_age-graptOccPBDB$late_age,
main="Age Resolution of Occurrences", xlab="Ma")
#distribution of taxa among taxonomic ranks
table(graptTaxaPBDB$taxon_rank)
barplot(table(graptTaxaPBDB$taxon_rank))
Run the code above in your browser using DataLab