FECRT <-
structure(list(Farm = structure(c(1L, 1L, 1L, 2L, 2L, 3L, 3L,
3L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 7L, 7L, 7L), .Label = c("5",
"11", "15", "20", "24", "36", "47"), class = "factor"), Drug =
structure(c(3L,
2L, 1L, 3L, 2L, 2L, 4L, 1L, 3L, 4L, 2L, 4L, 1L, 3L, 4L, 1L, 3L,
4L, 1L), .Label = c("Closantel", "Febendazole", "Ivermectin",
"Moxidectin"), class = "factor"), Reduction = structure(c(10L,
2L, 1L, 15L, 12L, 11L, 15L, 9L, 5L, 16L, 8L, 14L, 4L, 7L, 16L,
3L, 6L, 16L, 13L), .Label = c("-64", "-6", "7", "17", "46", "49",
"51", "53", "55", "64", "67", "82", "89", "90", "99", "100"), class =
"factor")), .Names = c("Farm",
"Drug", "Reduction"), row.names = c(NA, -19L), class = "data.frame")

barplot(FECRT$"Reduction", name = Farmgroup, main="% EPG Reduction after treatment - FECRT Farms", xlab = "Farm ID", ylab = "EPG Percentual Reduction")