Guest User

Untitled

a guest
Jan 23rd, 2019
102
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.80 KB | None | 0 0
  1. require(move)
  2. loginStored = movebankLogin(username = "Diego_Ellis", password = "")
  3. all_studies <- getMovebank(entity_type = "study", login=loginStored)
  4. dim(all_studies)
  5. head(all_studies$name)
  6.  
  7. library(ggplot2)
  8. world <- map_data("world")
  9.  
  10. ggplot() +
  11. geom_polygon(data = world,
  12. fill = "grey38",
  13. aes(x = long, y = lat, group = group)) +
  14. geom_point(data = all_studies,
  15. aes(x = main_location_long, y = main_location_lat, color = as.factor(there_are_data_which_i_cannot_see))) +
  16. scale_color_discrete(guide = guide_legend(title = "Data Restrictions"))
  17.  
  18. getMovebank("tag_type", loginStored)
  19. getMovebank("study", loginStored)
  20.  
  21. searchMovebankStudies(x="coyote", login=loginStored)
  22.  
  23. head(searchMovebankStudies(x="oose", loginStored))
  24.  
  25. head(searchMovebankStudies(x="bird", login = loginStored))
  26.  
  27.  
  28. ID <- getMovebankID(study = "Black-backed jackal, Etosha National Park, Namibia", login=loginStored)
  29. head(getMovebankAnimals(study = ID, login = loginStored))
  30.  
  31. # be patient
  32. jackals <- getMovebankData(study = "Black-backed jackal, Etosha National Park, Namibia", login = loginStored, removeDuplicatedTimestamps = TRUE)
  33.  
  34. idData <- jackals@idData
  35. jackal_movement <- jackals@data
  36. jackal_id_movement <- jackals@trackId
  37. jackal_movement$local_identifier <- jackal_id_movement
  38. str(jackal_movement)
  39.  
  40. jackal_movement <- read.table(file = "black-backed-jackal-Namibia.csv", header=T, sep=",")
  41.  
  42. library(tidyverse)
  43. library(ggmap)
  44.  
  45. jackal_box <- make_bbox(lon = location_long, lat = location_lat, data=jackal_movement, f= 0.1)
  46. jackal_map <- get_map(location = jackal_box, source = 'google', maptype = 'terrain')
  47. # ggmap(jackal_map) + geom_point(data=jackal_movement, aes(x=location_long, y=location_lat))
  48.  
  49. ggmap(jackal_map) + geom_point(data=jackal_movement, aes(x=location_long, y=location_lat, color=local_identifier), alpha = 0.5)
Add Comment
Please, Sign In to add comment