#! /bin/bash 
# Last edited on 2005-06-04 00:29:08 by stolfi
 
usage="plot-estimates < INFILE.est > CPLOT.eps"

# Reads an ".est" file produced by the predictor "pre".
# Plots the true coding/noncoding indicator Is(C)
# and the estimated probability Pr(C) for each base.
 
tmp="/tmp/$$"
 
datafile="${tmp}.dat"

# Extract fields from ".est" file:
cat \
  | plot-estimates-aux \
  > ${datafile}
             
# plot
gnuplot <<EOF
set terminal postscript eps color lw 2 "TimesRoman" 18
set size 1.5,0.5
set output
set yrange [-0.05:1.05]
plot \
  "${datafile}" using 1:3 title "Is(C)" with histeps lt 4, \
  "${datafile}" using 1:4 title "Pr(C)" with histeps lt 3
EOF

/bin/rm -f ${datafile}