Here is what my dataframe looks like:
> head(full_malaria_data)
  X index     lonx      laty District population malaria_cases river_length distance_to_river distance_to_coast
1 0     1 6.470243 0.2406650    Lemba         10             0     3098.054          136.9634         210.53670
2 1     2 6.474831 0.2397604    Lemba        395            23    15498.375          240.7952         214.72492
3 2     3 6.460882 0.2677222    Lemba       1862             8    13198.230          583.5622          65.33937
4 3     4 6.471704 0.2500610    Lemba        302             0    13198.230          231.2028         523.73073
I am trying to find the total number of malaria cases per district (for Lemba, Canta Calo, Principe, Agua Grande).
So far my code does one district at a time:
full_malaria_data %>% 
  summarise(out = sum(malaria_cases[District == "Lobata"])) %>%
  pull(out)
How can I augment my code so that it finds the sum of malaria cases for all districts at once?
