diff options
| author | Christoph Helma <helma@in-silico.ch> | 2017-12-20 16:16:40 +0100 |
|---|---|---|
| committer | Christoph Helma <helma@in-silico.ch> | 2017-12-20 16:16:40 +0100 |
| commit | 0983a03eaf0df05b464a4f537947692825560908 (patch) | |
| tree | e66e24d11f6636673f6ff909c1a5441f9a9fb949 /scripts | |
| parent | a5a9144dd7eb4cb9455c5674325ce6e0cc17af61 (diff) | |
database names in figures
Diffstat (limited to 'scripts')
| -rwxr-xr-x | scripts/dataset-variability.R | 3 | ||||
| -rwxr-xr-x | scripts/functional-groups4R.rb | 5 | ||||
| -rwxr-xr-x | scripts/median-correlation-plot.R | 2 |
3 files changed, 5 insertions, 5 deletions
diff --git a/scripts/dataset-variability.R b/scripts/dataset-variability.R index 68271d9..2fa7327 100755 --- a/scripts/dataset-variability.R +++ b/scripts/dataset-variability.R @@ -3,10 +3,9 @@ library(ggplot2) library(grid) library(gridExtra) -data <- read.csv("data/test_log10.csv",header=T) +data <- read.csv("data/test_log10_database_fix.csv",header=T) 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")) ggsave(file='figures/dataset-variability.pdf', plot=img, width=12,height=8) diff --git a/scripts/functional-groups4R.rb b/scripts/functional-groups4R.rb index 0b14b7a..3203b0f 100755 --- a/scripts/functional-groups4R.rb +++ b/scripts/functional-groups4R.rb @@ -1,3 +1,4 @@ +#!/usr/bin/env ruby require 'csv' csv = [] exclude = [ @@ -20,8 +21,8 @@ CSV.foreach("data/functional-groups.csv") do |row| keep = false if row[0].match(patt) end if keep and [row[1].to_i,row[2].to_i].max >= 25 - csv << [row[0].gsub('_',' '),row[1].to_i,"Mazzatorta"] - csv << [row[0].gsub('_',' '),row[2].to_i,"Swiss Federal Office"] + csv << [row[0].gsub('_',' '),row[1].to_i,"Nestle"] + csv << [row[0].gsub('_',' '),row[2].to_i,"FSVO"] else p row end diff --git a/scripts/median-correlation-plot.R b/scripts/median-correlation-plot.R index f4b28c2..3ea5ddd 100755 --- a/scripts/median-correlation-plot.R +++ b/scripts/median-correlation-plot.R @@ -3,6 +3,6 @@ library(ggplot2) experimental <- read.csv("data/median-correlation.csv",header=T) -img = qplot(mazzatorta,swiss,data=experimental,xlab="-log10(LOAEL Mazzatorta median)",ylab="-log10(LOAEL Swiss Federal Office median)") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) +img = qplot(mazzatorta,swiss,data=experimental,xlab="-log10(LOAEL Nestle median)",ylab="-log10(LOAEL FSVO median)") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) ggsave(file='figures/median-correlation.pdf', plot=img,width=12, height=8) |