# -*- gnuplot -*-
#
# LikeLor.gpt -- plots the log likelihood of quantized Lorenz data vs number
#                of states.  Derived from EM1.gpt
#

# Copyright (c) 2005, 2007, 2008 Andrew Fraser
#
# This file is part of HMM_DS_Code.
#
# HMM_DS_Code is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# HMM_DS_Code is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
# 02111-1307, USA.

set terminal epslatex monochrome dashed "normal" 9
set output "`echo -n $GP_OUTPUT`"
set size 0.98,0.98

# Not this time!  There are ytics and y2tics on different scales for
# this plot.
#
### set border 3

set logscale x
set xrange [16:4.0e6] 
set xtics 100
set xtics ('$10^2$' 100, '$10^4$' 10000, '$10^6$' 1e+06)
set xlabel '\Nstates'

set yrange [0:0.8]
set ytics .2
set ytics nomirror
set y2tics ('$0$'   0,		\
	    '$0.2$'  .2*log(2),	\
            '$0.4$'  .4*log(2),	\
            '$0.6$'  .6*log(2),	\
            '$0.8$'  .8*log(2))
set ylabel '\h'
set y2label '\bits'

# 'data/LikeLor'
plot "`echo -n $GP_INPUT1`" index 0 using ($3):($5) title 'HMM' with lines, \
     0.1359 title 'Bound' with lines

set output