diff options
Diffstat (limited to 'report/figures/results/plot.gnu')
| -rw-r--r-- | report/figures/results/plot.gnu | 104 |
1 files changed, 104 insertions, 0 deletions
diff --git a/report/figures/results/plot.gnu b/report/figures/results/plot.gnu new file mode 100644 index 0000000..53517b3 --- /dev/null +++ b/report/figures/results/plot.gnu @@ -0,0 +1,104 @@ +#set term x11 persist +set term cairolatex pdf size 6in,6in +set output "size_pts.tex" +set linetype 1 lc rgb '#E41A1C' # red +set linetype 2 lc rgb '#377EB8' # blue +set linetype 3 lc rgb '#4DAF4A' # green + +set key samplen 2 +set grid +set xlabel "x" +set ylabel "p(x)" + +set style fill transparent solid 0.5 border -1 + +set multiplot layout 3,2 + +set title "$N$=1750, $L$=128, $S$=$-$0.169, $K$=0.139" +plot 'N1750_L128.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +set title "$N$=500, $L$=128, $S$=0.227, $K$=0.753" +plot 'N500_L128.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +set title "$N$=1750, $L$=2048, $S$=$-$0.044, $K$=0.065" +plot 'N1750_L2048.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +set title "$N$=500, $L$=2048, $S$=0.187, $K$=0.204" +plot 'N500_L2048.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +set title "$N$=1750, $L$=4096, $S$=$-$0.045, $K$=0.072" +plot 'N1750_L4096.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +set title "$N$=500, $L$=4096, $S$=0.140, $K$=0.218" +plot 'N500_L4096.txt' using 1:2 with boxes lt 2 notitle, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +unset multiplot + +set terminal cairolatex pdf size 6in,3in +set output "single.tex" +unset title +set lmargin at screen 0.15 +set rmargin at screen 0.85 + +plot 'N1750_L4096.txt' using 1:2 with boxes title '$N=1750,\, L=4096$' lt 2, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines title 'TW\textsubscript{1}' lt 1, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:3 with lines title 'N(0,1)' dt 2 lt 3 + + +set terminal cairolatex pdf size 6in, 6in +set output "maxcomp.tex" +#set term x11 persist +set multiplot layout 3,2 + +unset lmargin +unset rmargin + +set title "$N$=500" +plot '500_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle +set title "$N$=750" +plot '750_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle +set title "$N$=1000" +plot '1000_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle +set title "$N$=1250" +plot '1250_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle +set title "$N$=1500" +plot '1500_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle +set title "$N$=1750" +plot '1750_max.txt' using 1:2 with boxes notitle lt 3, \ + '../tracy-widom-approx/tracy_widom_normal.txt' using 1:2 with lines lt 1 lw 5 notitle + +unset multiplot + +#set term x11 persist +set terminal cairolatex pdf size 6in,3in +set output "parameters.tex" +unset title +set lmargin at screen 0.15 +set rmargin at screen 0.85 + +set arrow from 400,0.2935 to 1800,0.2935 nohead lt 1 dt 3 lw 8 +set arrow from 400,0.1292 to 1800,0.1292 nohead lt 2 dt 3 lw 8 + +set xlabel "$N$" +set ylabel "Skewness/Excess Kurtosis" + +set y2tics +set ytics nomirror +set y2label "Error [\\%]" +set key outside top center horizontal + +plot 'props.txt' using 1:2 with linespoints lt 1 pt 7 ps 0.4 title "Skewness" axes x1y1, \ + 'props.txt' using 1:(($2-0.2935)/(0.2935)*100) with lines lt 1 dt 2 axes x1y2 title "Skewness Error", \ + 'props.txt' using 1:3 with linespoints lt 2 pt 7 ps 0.4 title "Excess Kurtosis" axes x1y1, \ + 'props.txt' using 1:(($3-0.1292)/(0.1292)*100) with lines lt 2 dt 2 axes x1y2 title "Excess Kurtosis Error" |