data(BivalvePBDB)
BivalveMatrix<-CoordList_PBDB(BivalvePBDB)
BivalveGeo<-GeoRarefaction_MultiTaxa(nLocCut=50,OccMatrix=BivalveMatrix,TaxaStart=3,iter=20)
PEE_SingleTaxon(GeoRare=BivalveGeo,TName=names(BivalveGeo)[3])
Run the code above in your browser using DataLab