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

# Reads an ".est" file produced by the predictor "pre".
# Plots true phase, estimated phase, and its estimated probability,
# 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,1
set output
set yrange [-1.1:3.1]
plot \
  "${datafile}" using 1:6 title "est. phase" with histeps lt 3, \
  "${datafile}" using 1:7 title "phase error" with histeps lt 4, \
  "${datafile}" using 1:8 title "Pr(phase)" with histeps lt 7
EOF

/bin/rm -f ${datafile}