server_test2/BuckleyAFB.R at main · jessiereese/server_test2 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
library(tidyverse)

data <- read_csv("2021_Buckley.csv")

# list of species
data %>%
  filter(substring(BirdCode,1,2)!='UN') %>%
  filter(BirdCode != "NOBI") %>%
  distinct(Species)

# species
data %>%
  filter(substring(BirdCode,1,2)!='UN') %>%
  filter(BirdCode != "NOBI") %>%
  summarise(n_species = n_distinct(Species))

# n individuals
data %>%
  filter(substring(BirdCode,1,2)!='UN') %>%
  filter(BirdCode != "NOBI") %>%
  summarise(n_indiv = n())

# of individuals per species
n.species <- data %>%
  filter(substring(BirdCode,1,2)!='UN') %>%
  filter(BirdCode != "NOBI") %>%
  group_by(Species) %>%
  summarise(n_indiv = n())

# start date
data %>%
  summarise(min(date))

# end data
data %>%
  summarise(max(date))

species <- read_csv("CPW_SpeciesStatus_Designations.csv")

data <- data %>%
  left_join(species)

data %>%
  select(Species, CPW_Designation) %>%
  filter(CPW_Designation != "NA") %>%
  distinct()

n.species <- n.species %>%
  left_join(species)

write.csv(n.species, file = "n.species.csv" )