JAMP
JAMP copied to clipboard
Update BOLD_web_hack.R
- New option to use output from Bold identification with Species Level Barcode Records database by setting
db = "sp"
. With default options (db="all"
) it should work just as original code does. Withdb="sp"
, the script expects input in format of species level database.
*see attachment for input and output ofBold_web_hack("BOLDtest.txt",db="sp")
*if db is messed up, you either get an error or table full of NAs
Other minor optimizing changes are:
- one
write.csv(temp_tab,...)
is after theif
, instead of two inside theif
-
OTU <- sub(...)
now gives "OTU_1234"; previously it gave "Query: OTU_1234" - object
whois_num
removed, because it wasn't used anywhere in the code and it was identical withwhois