1 #!/usr/bin/gnuplot -persist
5 # Version 4.4 patchlevel 0
6 # last modified March 2010
7 # System: Linux 2.6.32-5-amd64
9 # Copyright (C) 1986-1993, 1998, 2004, 2007-2010
10 # Thomas Williams, Colin Kelley and many others
12 # gnuplot home: http://www.gnuplot.info
13 # faq, bugs, etc: type "help seeking-assistance"
14 # immediate help: type "help"
15 # plot window: hit 'h'
16 set terminal pdfcairo size 15.00cm, 9.27cm
17 set output 'markov-comparators-18-pct.pdf'
21 set bar 1.000000 front
22 set border 31 front linetype -1 linewidth 1.000
28 set timefmt x "%d/%m/%y,%H:%M"
29 set timefmt y "%d/%m/%y,%H:%M"
30 set timefmt z "%d/%m/%y,%H:%M"
31 set timefmt x2 "%d/%m/%y,%H:%M"
32 set timefmt y2 "%d/%m/%y,%H:%M"
33 set timefmt cb "%d/%m/%y,%H:%M"
35 set style fill empty border
36 set style rectangle back fc lt -3 fillstyle solid 1.00 border lt -1
47 set key outside center bottom horizontal Right noreverse enhanced autotitles box linetype -1 linewidth 1.000
48 set key noinvert samplen 4 spacing 1 width 0 height 0
51 set style increment default
54 set style histogram clustered gap 2 title offset character 0, 0, 0
56 set offsets 0, 0, 0, 0
69 set datafile separator whitespace
73 set cntrparam levels auto 5
74 set cntrparam points 5
78 set style function lines
79 set xzeroaxis linetype -2 linewidth 1.000
80 set yzeroaxis linetype -2 linewidth 1.000
81 set zzeroaxis linetype -2 linewidth 1.000
82 set x2zeroaxis linetype -2 linewidth 1.000
83 set y2zeroaxis linetype -2 linewidth 1.000
91 set xtics border in scale 1,0.5 mirror norotate offset character 0, 0, 0
92 set xtics autofreq norangelimit
93 set ytics border in scale 1,0.5 mirror norotate offset character 0, 0, 0
94 set ytics autofreq norangelimit
95 set ztics border in scale 1,0.5 nomirror norotate offset character 0, 0, 0
96 set ztics autofreq norangelimit
99 set cbtics border in scale 1,0.5 mirror norotate offset character 0, 0, 0
100 set cbtics autofreq norangelimit
102 set title offset character 0, 0, 0 font "" norotate
105 set timestamp offset character 0, 0, 0 font "" norotate
106 set rrange [ * : * ] noreverse nowriteback # (currently [8.98847e+307:-8.98847e+307] )
107 set trange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] )
108 set urange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] )
109 set vrange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] )
110 set xlabel "Komparatoren"
111 set xlabel offset character 0, 0, 0 font "" textcolor lt -1 norotate
113 set x2label offset character 0, 0, 0 font "" textcolor lt -1 norotate
114 set xrange [ * : * ] noreverse nowriteback # (currently [82.0000:139.000] )
115 set x2range [ * : * ] noreverse nowriteback # (currently [82.0000:139.000] )
117 set ylabel offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
119 set y2label offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
120 set yrange [ * : * ] noreverse nowriteback # (currently [0.00000:10.0000] )
121 set y2range [ * : * ] noreverse nowriteback # (currently [0.00000:9.74040] )
123 set zlabel offset character 0, 0, 0 font "" textcolor lt -1 norotate
124 set zrange [ * : * ] noreverse nowriteback # (currently [8.98847e+307:-8.98847e+307] )
126 set cblabel offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
127 set cbrange [ * : * ] noreverse nowriteback # (currently [8.98847e+307:-8.98847e+307] )
133 set locale "en_US.UTF-8"
134 set pm3d explicit at s
135 set pm3d scansautomatic
136 set pm3d interpolate 1,1 flush begin noftriangles nohidden3d corners2color mean
137 set palette positive nops_allcF maxcolors 0 gamma 1.5 color model RGB
138 set palette rgbformulae 7, 5, 15
140 set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front bdefault
143 set fit noerrorvariables
144 binom(n,k) = (n!) / ((k!) * ((n-k)!))
145 binv(n,k,p) = binom(n,k) * (p**k) * ((1-p)**(n-k))
146 bin(x) = binv(n,floor(x),p)
147 gauss(x) = (a1 / sigma) * exp((-0.5) * ((x - mu) / sigma)**2.0)
148 poisson_int(x) = exp((-1.0) * lambda) * (lambda**x) / (x!)
149 poisson(x) = poisson_int(floor(x + 0.5))
150 erl(x) = gamma * exp(-1.0 * gamma * x) * (gamma * x)**(n_erl - 1) / ((n_erl - 1)!)
151 binomial_int(x) = binomial_coeff(n, x) * p**x * (1-p)**(n-x)
152 binomial_coeff(n,k) = (n!) / ((k!) * ((n-k)!))
153 binomial(x) = binomial_int(floor(x + 0.5))
154 binomial_fit(x) = 100.0 * binomial(x-19.0)
155 poisson_fit(x) = 100.0 * poisson(x - offset)
156 gauss_fit(x) = 100.0 * gauss(x)
157 cauchy(x) = (1.0 / pi) * s / (s**2.0 + (x - t)**2)
158 cauchy_fit(x) = 100.0 * cauchy(x)
159 erlang_fit(x) = 100.0 * erl(x - erl_offset)
160 gamma_dist(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma(gamma_k))
161 gamma_fit(x) = 100.0 * dgamma(x - gamma_offset, gamma_k, gamma_theta)
162 dgamma(x, shape, rate) = (x<0)? 0 : (x==0)? ((shape<1)? 1/0 : (shape==1)? rate : 0) : (rate==0)? 0 : exp(_ln_dgamma(x, shape, rate))
163 foo_fit(x) = foo_a0 * exp(foo_a1 * (x - 62))
164 _ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x
165 pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)
167 GPFUN_binom = "binom(n,k) = (n!) / ((k!) * ((n-k)!))"
168 GPFUN_binv = "binv(n,k,p) = binom(n,k) * (p**k) * ((1-p)**(n-k))"
170 p = 0.172562132555044
171 GPFUN_bin = "bin(x) = binv(n,floor(x),p)"
172 a1 = 0.398942280401433
173 sigma = 1.88445909011715
175 GPFUN_gauss = "gauss(x) = (a1 / sigma) * exp((-0.5) * ((x - mu) / sigma)**2.0)"
176 lambda = 5.10554056666183
177 GPFUN_poisson_int = "poisson_int(x) = exp((-1.0) * lambda) * (lambda**x) / (x!)"
178 GPFUN_poisson = "poisson(x) = poisson_int(floor(x + 0.5))"
179 gamma = 1.90073273091251
181 GPFUN_erl = "erl(x) = gamma * exp(-1.0 * gamma * x) * (gamma * x)**(n_erl - 1) / ((n_erl - 1)!)"
182 GPFUN_binomial_int = "binomial_int(x) = binomial_coeff(n, x) * p**x * (1-p)**(n-x)"
183 GPFUN_binomial_coeff = "binomial_coeff(n,k) = (n!) / ((k!) * ((n-k)!))"
184 GPFUN_binomial = "binomial(x) = binomial_int(floor(x + 0.5))"
185 GPFUN_binomial_fit = "binomial_fit(x) = 100.0 * binomial(x-19.0)"
187 GPFUN_poisson_fit = "poisson_fit(x) = 100.0 * poisson(x - offset)"
188 GPFUN_gauss_fit = "gauss_fit(x) = 100.0 * gauss(x)"
191 GPFUN_cauchy = "cauchy(x) = (1.0 / pi) * s / (s**2.0 + (x - t)**2)"
192 GPFUN_cauchy_fit = "cauchy_fit(x) = 100.0 * cauchy(x)"
193 erl_offset = 49.8425987940659
194 GPFUN_erlang_fit = "erlang_fit(x) = 100.0 * erl(x - erl_offset)"
195 gamma_k = 10.7234521744736
196 gamma_theta = 0.76572583751246
197 GPFUN_gamma_dist = "gamma_dist(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma_gamma(gamma_k))"
199 GPFUN_gamma_fit = "gamma_fit(x) = 100.0 * dgamma(x - gamma_offset, gamma_k, gamma_theta)"
200 foo_a0 = 0.000540217815483465
201 foo_a1 = 0.838098058527587
202 GPFUN_foo_fit = "foo_fit(x) = foo_a0 * exp(foo_a1 * (x - 62))"
203 GPFUN__ln_dgamma = "_ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x"
204 GPFUN_pgamma = "pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)"
205 GPFUN_dgamma = "dgamma(x, shape, rate) = (x<0)? 0 : (x==0)? ((shape<1)? 1/0 : (shape==1)? rate : 0) : (rate==0)? 0 : exp(_ln_dgamma(x, shape, rate))"
211 FIT_STDFIT = 0.0927264844441511
212 FIT_WSSR = 0.481499251372798
215 GPFUN_gamma_gamma = "gamma_gamma(x) = (floor(x - 0.5))!"
216 GPFUN_gamma = "gamma(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma_gamma(k))"
217 plot [x=82:139] 'markov-comparators-18-pct.data' title 'Gemessene Daten', gamma_fit(x) title "Gamma-Verteilung"
218 ## fit gamma_fit(x) 'markov-comparators-18-pct.data' via gamma_k, gamma_theta