# $CVSHeader: pal8e/dat/stat/plot-Hrmsdbrok,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

set xrange [7:60]
#set yrange [0:0.2]
set xlabel "Helix length"
set ylabel "RMSD"
set title "Broken Helix RMSD vs length from PSSD"

f(x) = m*x + c

m = 0.0118
c = 2.1756

set terminal x11
plot "Hrmsd-brok" with points title "RMSD", \
	'Hrmsd-brok' smooth unique linewidth 2 title "Average RMSD", \
	f(x) linewidth 2 title "Approximated RMSD of linear helix"

set terminal postscript "Times-Roman" 12
set output "Hrmsd-brok.ps"
replot


