# HG changeset patch
# User nicolas
# Date 1477420864 14400
# Node ID d961f726b619a65e415378a17450ebd9541bc447
# Parent  2e66da6efc4183fdab007feb22dd6b5566b91684
Uploaded

diff -r 2e66da6efc41 -r d961f726b619 plotPrediction.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/plotPrediction.R	Tue Oct 25 14:41:04 2016 -0400
@@ -0,0 +1,30 @@
+########################################################
+#
+# creation date : 07/06/16
+# last modification : 07/06/16
+# author : Dr Nicolas Beaume
+# owner : IRRI
+#
+########################################################
+
+library("miscTools")
+# scatterplot of the prediction vs target
+r2.plot <- function(true, predicted) {
+  # the scatterplot
+  plot(true, predicted, xlab="trait value", ylab="predicted value", main="", pch=16, 
+       ylim=c(min(min(true), min(predicted)), max(max(true), max(predicted))))
+  # add a red lines with ideal case
+  lines(true, true, col="red")
+}
+
+############################ main #############################
+# load argument
+cmd <- commandArgs(trailingOnly = T)
+source(cmd[1])
+# load prediction and target
+phenotype <- read.table(phenotype, sep="\t", h=T)[,1] 
+predicted <- read.table(predicted, sep = "\t", h=T)[,2]
+# plot in a pdf that will be available in galaxy history panel
+pdf(out)
+r2.plot(phenotype, predicted = predicted)
+dev.off()
\ No newline at end of file