Mostly finished results and lots of plots

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"