blob: deb9c1ce131ef81809ae5d87beb6f4cfec88a777 (
plain)
1
2
3
4
5
6
7
8
9
10
11
|
#!/usr/bin/env Rscript
library(ggplot2)
data <- read.csv("tables/pa-table.csv")
for (i in c(2:10)) {
name <- names(data)[i]
cols <- append(c(15,17),c(19:28))
group <- data[data[i] == 1,cols]
freq <- 100*colSums(group,na.rm=TRUE)/colSums(!is.na(group))
plot <- ggplot(data.frame(freq),aes(x=names(data)[cols],y=freq)) + geom_bar(stat="identity") + ylab("% mutagenic") + xlab(name) + ylim(c(0,100))
ggsave(paste("figures/",name,".png",sep=""))
}
|