#! /usr/bin/env gnuplot

set terminal postscript eps color
set size 5.5/10,4.0/7
set key outside samplen 1
set xrange [-0.5:*]

set output "histrmsd.eps"
set title "Backbone RMSD"
set xlabel "RMSD"
set ylabel "Fraction" 1.0
plot \
    "< cut -f2 verifymatch.out.rmsd.avg.100  | histogram.py -r 0:5:0.5 -N 1" \
        using 2:1 title "Functional" with linespoints lw 2 \
    , "< cut -f2 verifymatch.R1000.out.rmsd.avg.100  | histogram.py -r 0:5:0.5 -N 1" \
        using 2:1 title "Random" with linespoints lw 2

system "convert -density 120 histrmsd.eps histrmsd.png"

set output "histfeatured.eps"
set title "Physicochemical Distance"
set xlabel "FEATURE distance"
set ylabel "Fraction" 1.0
set xtics rotate
plot \
    "< cut -f2 verifymatch.out.featured.avg.100  | histogram.py -r 0:2000:200 -N 1" \
        using 2:1 title "Functional" with linespoints lw 2 \
    , "< cut -f2 verifymatch.R1000.out.featured.avg.100  | histogram.py -r 0:2000:200 -N 1" \
        using 2:1 title "Random" with linespoints lw 2

system "convert -density 120 histfeatured.eps histfeatured.png"

set output "histroc.eps"
set size square 6.0/10,4.0/7
set title "Performance"
set xlabel "Specificity" ,.5
set ylabel "Sensitivity" 1.0
set xtics norotate
set key outside samplen 1
set xrange [0:1]
set yrange [0:1]
plot \
    "< echo 0 0 0 0; cat rmsd.roc.txt" using 3:1 title "RMSD" with linespoints lw 2 \
    , "< echo 0 0 0 0; cat feature.roc.txt" using 3:1 title "Physicochemical" with linespoints lw 2 \
    , "< echo 0 0 0 0; echo 1 1 1 1" using 3:1 title "random" with linespoints lw 1

system "convert -density 120 histroc.eps histroc.png"