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?