1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
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"
|
