diff options
Diffstat (limited to 'scripts/test-correlation-plot.R')
| -rwxr-xr-x | scripts/test-correlation-plot.R | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/scripts/test-correlation-plot.R b/scripts/test-correlation-plot.R index 74a2739..ef69058 100755 --- a/scripts/test-correlation-plot.R +++ b/scripts/test-correlation-plot.R @@ -1,15 +1,17 @@ +#!/usr/bin/Rscript + library(ggplot2) library(grid) library(gridExtra) experimental <- read.csv("data/median-correlation.csv",header=T) -p1 = qplot(-log10(mazzatorta),-log10(swiss),data=experimental,xlab="-log10(LOAEL Mazzatorta median)",ylab="-log10(LOAEL Swiss Federal Office median)",main="Experimental data") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) +p1 = qplot(mazzatorta,swiss,data=experimental,xlab="-log10(LOAEL Mazzatorta median)",ylab="-log10(LOAEL Swiss Federal Office median)",main="Experimental data") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) training = read.csv("data/training-test-predictions.csv",header=T) -p2 = qplot(-log10(LOAEL_predicted),-log10(LOAEL_measured_median),data=training,xlab="-log10(LOAEL predicted)",ylab="-log10(LOAEL measured median)",main="Combined") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) +p2 = qplot(LOAEL_predicted,LOAEL_measured_median,data=training,xlab="-log10(LOAEL predicted)",ylab="-log10(LOAEL measured median)",main="Combined") + geom_point() + geom_abline(intercept=0.0) + xlim(-1,4) + ylim(-1,4) -pdf('figure/test-correlation.pdf') +pdf('figures/test-correlation.pdf') grid.arrange(p1,p2,ncol=1,respect=T) dev.off() |