# $CVSHeader: pal8e/dat/stat/plot-QSynCurve,v 1.2 2005/02/15 01:10:24 indraneel Rel $

## Run using GNUPlot

set data style lines
set mxtics 5
show mxtics
set mytics 5
show mytics
#set key below box

#a == amplitude
#s == standard deviation
#c == center

#fs(x) = as * exp (-0.5 * (((x-cs)/ss)**2))

#fa1(x) = aa1 * exp (-0.5 * (((x-ca1)/sa)**2))
#fa2(x) = aa2 * exp (-0.5 * (((x-ca2)/sa)**2))
#fa(x) = fa1(x) + fa2(x)

#as = 0.170
#cs = 4.810
#ss = 0.227

#aa1 = 0.073
#aa2 = 0.077
#ca1 = 4.463
#ca2 = 5.245
#sa = 0.261

#fit [3.2:7.2] fs(x) "QdistSyn" using 1:2 via as, ss, cs
#fit [3.2:7.2] fa(x) "QdistAnti" using 1:2 via aa1, aa2, ca1, ca2, sa

set xrange [0:120]
#set yrange [0:0.2]
set xlabel "Quadruplet score"
set ylabel "Normalised Count"
set title "True and False Parallel Quadruplet Scores from DSSP"

set arrow from 45,graph 0 to 45,graph 1 nohead lt 1 lw 4
set label "c1" at 40,graph 0.96
set arrow from 46,graph 0 to 46,graph 1 nohead lt 1 lw 4
set label "c2" at 47,graph 0.96

set terminal x11
plot "QgoodSyn" with lines lt 2 lw 4 title "True quadruplets", \
	"QbadSyn" with lines lt 3 lw 4 title "False quadruplets"

#set terminal postscript enhanced color blacktext solid "Times-Roman" 12
set terminal postscript eps enhanced "Times-Roman" 22
set output "QSynCurve.eps"
replot


