diff options
author | Christoph Helma <helma@in-silico.ch> | 2016-03-04 15:27:20 +0100 |
---|---|---|
committer | Christoph Helma <helma@in-silico.ch> | 2016-03-04 15:27:20 +0100 |
commit | 7ad7c10c1e708f6b5a3473de24dbeab03d0b74a3 (patch) | |
tree | 14e7a8c37343f3d878e8116873978861b226b5bc /dataset-variability.R | |
parent | d3071896a7116670756199f0df7c2a618de2aea3 (diff) |
prediction interval, discussion
Diffstat (limited to 'dataset-variability.R')
-rw-r--r-- | dataset-variability.R | 15 |
1 files changed, 12 insertions, 3 deletions
diff --git a/dataset-variability.R b/dataset-variability.R index b0e3c76..775fd03 100644 --- a/dataset-variability.R +++ b/dataset-variability.R @@ -18,7 +18,16 @@ s.dup$SMILES <- reorder(s.dup$SMILES,s.dup$LOAEL) p1 <- ggplot(m.dup, aes(SMILES,LOAEL),ymin = min(LOAEL), ymax=max(LOAEL)) + ylab('-log(LOAEL mg/kg_bw/day)') + xlab('Compound') + theme(axis.text.x = element_blank()) + geom_point() + ggtitle("Mazzatorta") + ylim(-1,4) p2 <- ggplot(s.dup, aes(SMILES,LOAEL),ymin = min(LOAEL), ymax=max(LOAEL)) + ylab('-log(LOAEL mg/kg_bw/day)') + xlab('Compound') + theme(axis.text.x = element_blank()) + geom_point() + ggtitle("Swiss Federal Office") + ylim(-1,4) -pdf('figure/dataset-variability.pdf') -grid.arrange(p1,p2,ncol=1) -dev.off() +#pdf('figure/dataset-variability.pdf') +#grid.arrange(p1,p2,ncol=1) +#dev.off() +data <- read.csv("data/test.csv",header=T) +data$LOAEL = -log(data$LOAEL) +data$SMILES <- reorder(data$SMILES,data$LOAEL) +img = ggplot(data,aes(SMILES,LOAEL,ymin = min(LOAEL), ymax=max(LOAEL),color=Dataset)) + geom_point() +img <- img + ylab('-log(LOAEL mg/kg_bw/day)') + xlab('Compound') + theme(axis.text.x = element_blank()) + theme(legend.title=element_blank()) +img = img + scale_fill_discrete(breaks=c("Mazzatorta", "Both", "Swiss Federal Office")) +img = img + +ggsave(file='figure/dataset-variability.pdf', plot=img, width=12,height=8) |