I am trying to show the limestone treatment amount of Ca on treatment sites (2,000 kg/ha) and I would like to remove the line from the Untreated site facet.
The code I have been using is:
soillimebudget %>%
mutate(`Stand Age` = factor(`Stand Age`, levels = c("NL", "14", "19", "23", "28", "32", "33", "37"))) %>%
mutate(`Treatment` = factor(`Treatment`, levels = c("Pre-Treatment Conditions", "Limestone, Stable", "Limestone, Eroded"))) %>% mutate(Treatment = recode(Treatment, "Pre-Treatment Conditions" = "Untreated")) %>%
mutate(Pool = factor(Pool, levels = c("ABG", "LFH", "Mineral Soil"))) %>%
mutate(Pool = recode(Pool, "ABG" = "AGB")) %>%
group_by(Treatment, `Stand Age`, Pool) %>%
ggplot(aes(x = Mg, y = `Stand Age`, fill = `Pool`, position = "stack")) +
labs(title = "Mg Budget") +
labs(x = "Stand Age", y = "Mg (kg/ha)") +
geom_bar(aes(`Stand Age`, `Mg`, fill = Pool),
position = "stack", stat = "identity") +
#scale_fill_grey(start = 0.3, end = 0.7) +
theme_bw() +
geom_hline(linetype = "dashed", yintercept = 1000) +
theme(legend.position = "top") +
theme(panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.background = element_blank(),
panel.border = element_rect(colour = "Black", fill = NA)) +
theme(text = element_text(size = 30)) +
facet_grid(~`Treatment`, scales = "free_x", space = "free_x")
Currently my figure is how I would like, barring removing the one H line from the Facet.
Thanks