#################################################
set term png
set out 'self-mi.png'
# set term postscript eps enhanced color 20
# set out 'self-mi.eps'
set key right top
set title "Self MI"
set xlabel "MI (bits)"
set ylabel "Frequency"
set logscale y
set format y "10^{%L}"
# set yrange [1.0e-6:1]
nbins = 40.0
width = 40
ncls3 = 530
normcla3 = nbins/(width*ncls3)
ncls4 = 536
normcla4 = nbins/(width*ncls4)
nwrds3 = 9490
normw3 = nbins/(width*nwrds3)
nwrds4 = 9486
normw4 = nbins/(width*nwrds4)
plot "self-mi-n3-classes.dat" using 2:($3*normcla3) with lines lw 2 title "classes noise=3", \
"self-mi-n4-classes.dat" using 2:($3*normcla4) with lines lw 2 title "classes noise=4", \
"self-mi-n3-words.dat" using 2:($3*normw3) with lines lw 2 title "words noise=3", \
"self-mi-n4-words.dat" using 2:($3*normw4) with lines lw 2 title "words noise=4"
#################################################
set term png
set out 'class-mi.png'
# set term postscript eps enhanced color 20
# set out 'class-mi.eps'
set key right top
set title "Class-Class MI"
set xlabel "MI (bits)"
set ylabel "Frequency"
set logscale y
set format y "10^{%L}"
set yrange [2.0e-4:0.2]
nbins = 100.0
width = 40
ncls3 = 10154
normcla3 = nbins/(width*ncls3)
ncls4 = 9988
normcla4 = nbins/(width*ncls4)
# Normal distribution
g(x,mu,sigma) = exp(-(x-mu)**2/(2*sigma**2))/(sigma * sqrt(2*3.14159))
plot "class-mi-n3.dat" using 2:($3*normcla3) with lines lw 2 title "classes noise=3", \
"class-mi-n4.dat" using 2:($3*normcla4) with lines lw 2 title "classes noise=4", \
"class-mi-n3.dat" using 2:(g($2,-3,3.5)) with lines lw 2 title "G(-3,3.5)"
#################################################
set term png
set out 'word-mi.png'
# set term postscript eps enhanced color 20
# set out 'word-mi.eps'
set key left top
set title "Word-Word MI"
set xlabel "MI (bits)"
set ylabel "Frequency"
set logscale y
set format y "10^{%L}"
set yrange [2.0e-5:0.2]
nbins = 100.0
width = 30
nwrds3 = 117218
normw3 = nbins/(width*nwrds3)
nwrds4 = 117023
normw4 = nbins/(width*nwrds4)
plot "word-mi-n3.dat" using 2:($3*normw3) with lines lw 2 title "words noise=3", \
"word-mi-n4.dat" using 2:($3*normw4) with lines lw 2 title "words noise=4", \
"word-mi-n3.dat" using 2:(g($2,-1.5,3.1)) with lines lw 2 title "G(-1.5,3.1)"
#################################################
set term png
set out 'clawrd-mi.png'
# set term postscript eps enhanced color 20
# set out 'clawrd-mi.eps'
set key right top
set title "Class-Word MI"
set xlabel "MI (bits)"
set ylabel "Frequency"
set logscale y
set format y "10^{%L}"
set yrange [4.0e-5:0.2]
nbins = 100.0
width = 30
nclwr3 = 58555
normclw3 = nbins/(width*nclwr3)
nclwr4 = 55194
normclw4 = nbins/(width*nclwr4)
plot "clawrd-mi-n3.dat" using 2:($3*normclw3) with lines lw 2 title "noise=3", \
"clawrd-mi-n4.dat" using 2:($3*normclw4) with lines lw 2 title "noise=4", \
"clawrd-mi-n3.dat" using 2:(g($2,-1.5,3.1)) with lines lw 2 title "G(-1.5,3.1)"
#################################################
set term png
#set out 'all-mi.png'
set term postscript eps enhanced color 20
set out 'all-mi.eps'
set key right top
set title "All-in-One MI"
set xlabel "MI (bits)"
set ylabel "Frequency"
set logscale y
set format y "10^{%L}"
set yrange [4.0e-5:0.2]
plot "class-mi-n3.dat" using 2:($3*normcla3) with lines lw 2 title "classes n=3", \
"class-mi-n4.dat" using 2:($3*normcla4) with lines lw 2 title "classes n=4", \
"word-mi-n3.dat" using 2:($3*normw3) with lines lw 2 title "words n=3", \
"word-mi-n4.dat" using 2:($3*normw4) with lines lw 2 title "words n=4", \
"clawrd-mi-n3.dat" using 2:($3*normclw3) with lines lw 2 title "class-word n=3", \
"clawrd-mi-n4.dat" using 2:($3*normclw4) with lines lw 2 title "class-word n=4", \
#################################################