data-raw/01b_data-process_taxonomy.R

library("seabirdPrey")
library(dplyr)

## read-in-raw-prey-taxonomy ----
seabirddiet <- readr::read_csv(here::here("data-raw", "intermediate", "seabirddiet.csv"),
                               guess_max = Inf, na = c("", "NA")) 

## process-worms_validate-prey-taxonomy ----
prey_taxonomy <-  worrms_validate(seabirddiet$prey_species)
pred_taxonomy <-  worrms_validate(seabirddiet$species)

## prepare classification table ----
classification <- classify_validname(unique(prey_taxonomy$valid_name,
                                            pred_taxonomy$valid_name), 
                                     ranks = c("phylum", "subphylum", "class", "order", 
                                               "family", "genus", "species"))  
classification %>%
    dplyr::arrange(phylum, subphylum, class, order, family, genus, species) %>%
    readr::write_csv(here::here("data-raw", "metadata", "classification.csv"))
annakrystalli/seabirddietDB documentation built on Nov. 2, 2019, 1:54 p.m.