#! /usr/bin/env gnuplot

set terminal postscript eps color
set output "plotroc_combine.eps"
set size square 6.0/10,4.0/7
set xlabel "False Positive Rate" ,.5
set ylabel "True Positive Rate" 1.0
set title "Serine endopeptidase"
set key outside samplen 1

plot \
    "pfam-astral40.trypsin.hmm.out.max-sprot42.ec_3_4_21.act_site_ser.astral40.id-roc" \
        using (1-$2):(1-$1) title "Pfam" with linespoints lw 2 \
    , "seqfeature-astral40.scores.max_sprot42.ec_3_4_21.act_site_ser.astral40.id_roc" \
        using (1-$2):(1-$1) title "SeqFEATURE" with linespoints lw 2 \
    , "<perl -e 'print 22/45,\" \",1/(5733-45),\"\n\"'" \
        using ($2):($1) title "PROSITE" with points lw 2 pt 5

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