Hi all,
Is there a way to get longitude and latitude information related to each BOLD accesion? the function bold_seq gives back four entries; id, name, gene, and sequence but not lon and lat.
Thanks in advance!
Hi all,
Is there a way to get longitude and latitude information related to each BOLD accesion? the function bold_seq gives back four entries; id, name, gene, and sequence but not lon and lat.
Thanks in advance!
Thanks for your question about the bold package @nereal
Checkout the function bold_seqspec()
, which gets sequence data and specimen data, e.g.,
res <- bold_seqspec(taxon = "Coelioxys")
tibble::as_tibble(res[, c("lat", "lon", "nucleotides")])
#> # A tibble: 315 x 3
#> lat lon nucleotides
#> <dbl> <dbl> <chr>
#> 1 47.4 16.3 ---------------------------------------------------------------------------------------------------------------------------------------------AATAATGATATAATTTATAACā¦
#> 2 52.9 -118. TATAATATATATAATTTTTGCAATATGATCAGGAATAATTGGATCTTCTATAAGTATAATTATTCGAATAGAATTAAGAATCCCAGGATCATGAATTAATAATGATCAAATTTATAACTCTTTTATTACAGCACATGCATTTTTAATAATTTTTTTTTTAGTā¦
#> 3 NA NA -------------------------------GGTATAATTGGATCATCTTTAAGAATAATTATTCGCATAGAATTAAGAATCCCGGGTTCTTGAATTAACAATGATCAAATTTATAATTCTTTTATTACAGCTCATGCCTTTTTAATAATTTTTTTCCTAGTā¦
#> 4 NA NA CATCCTATATATAATTTTTGCCATATGGTCAGGAATAATTGGATCTTCATTAAGAATAATTATTCGTATAGAATTAAGAATCCCAGGCTCTTGGATTAGTAATGACCAAATTTATAATTCTTTTATTACTGCTCATGCATTTTTAATAATTTTTTTTTTAGTā¦
#> 5 51.5 -0.967 ATAAAGATATTGGTATAATATATATAATTTTTGCTATATGATCAGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGGTCTTGAATTAATAATGACCAAATTTATAATTCTTTCATTACAGCCCATGCATTTTTAATAAā¦
#> 6 NA NA ----------------------------------ATAATTGGATCTTCTTTAAGAATAATCATTCGAATAGAATTAAGAACTCCAGGATCTTGAATTAATAATGATCAAATCTATAACTCATTTATTACCGCCCATGCATTTCTAATAATTTTTTTTTTAGTā¦
#> 7 NA NA AATATTATATATAATTTTTGCTATATGATCTGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGGTCATGAATTAATAATGATCAAATTTATAATTCTTTCATTACAGCTCATGCATTTTTAATAATTTTTTTTTTAGTā¦
#> 8 18.8 -96.9 -ATATTATATATAATTTTTGCTATATGATCTGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGATCATGAATCAATAATGATCAAATTTATAATTCTTTTATTACAGCTCATGCATTTTTAATAATTTTTTTTTTAGTā¦
#> 9 NA NA ---------------------------------AATAATTGGATCATCAATTAGAATATTAATTCGAATAGAATTAAGAATTCCTGGAACATGAATTAATAATGATCAAATCTATAATTCATTTATTACCGCTCATGCATTTTTAATAATTTTTTTTTTAGTā¦
#> 10 52.2 13.5 AATAATATATATAATTTTTGCTATATGATCAGGAATAATTGGTTCTTCATTAAGAATAATTATTCGTATAGAATTAAGAATTCCTGGAGCTTGAATTAATAATGATCAAATTTATAATTCATTTATTACAGCTCACGCATTTATAATAATTTTTTTTTTAGTā¦
#> # ⦠with 305 more rows
Does that get what you want?
Yes, fantastic! sorry i missed that funtion out
Thank you Sckott