images/comparison-comparators-16.gnuplot

   1 #!/usr/bin/gnuplot -persist
   2 #
   3 #
   4 #       G N U P L O T
   5 #       Version 4.4 patchlevel 0
   6 #       last modified March 2010
   7 #       System: Linux 2.6.32-5-amd64
   8 #
   9 #       Copyright (C) 1986-1993, 1998, 2004, 2007-2010
  10 #       Thomas Williams, Colin Kelley and many others
  11 #
  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 'comparison-comparators-16.pdf'
  18 unset clip points
  19 set clip one
  20 unset clip two
  21 set bar 1.000000 front
  22 set border 31 front linetype -1 linewidth 1.000
  23 set xdata
  24 set ydata
  25 set zdata
  26 set x2data
  27 set y2data
  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"
  34 set boxwidth
  35 set style fill  empty border
  36 set style rectangle back fc  lt -3 fillstyle   solid 1.00 border lt -1
  37 set dummy x,y
  38 set format x "% g"
  39 set format y "% g"
  40 set format x2 "% g"
  41 set format y2 "% g"
  42 set format z "% g"
  43 set format cb "% g"
  44 set angles radians
  45 unset grid
  46 set key title ""
  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
  49 unset label
  50 unset arrow
  51 set style increment default
  52 unset style line
  53 unset style arrow
  54 set style histogram clustered gap 2 title  offset character 0, 0, 0
  55 unset logscale
  56 set offsets 0, 0, 0, 0
  57 set pointsize 1
  58 set encoding default
  59 unset polar
  60 unset parametric
  61 unset decimalsign
  62 set view 60, 30, 1, 1
  63 set samples 100, 100
  64 set isosamples 10, 10
  65 set surface
  66 unset contour
  67 set clabel '%8.3g'
  68 set mapping cartesian
  69 set datafile separator whitespace
  70 unset hidden3d
  71 set cntrparam order 4
  72 set cntrparam linear
  73 set cntrparam levels auto 5
  74 set cntrparam points 5
  75 set size ratio 0 1,1
  76 set origin 0,0
  77 set style data points
  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
  84 set ticslevel 0.5
  85 set mxtics default
  86 set mytics default
  87 set mztics default
  88 set mx2tics default
  89 set my2tics default
  90 set mcbtics default
  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
  97 set nox2tics
  98 set noy2tics
  99 set cbtics border in scale 1,0.5 mirror norotate  offset character 0, 0, 0
 100 set cbtics autofreq  norangelimit
 101 set title ""
 102 set title  offset character 0, 0, 0 font "" norotate
 103 set timestamp bottom
 104 set timestamp ""
 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 [-10.0000:10.0000] )
 109 set vrange [ * : * ] noreverse nowriteback  # (currently [-10.0000:10.0000] )
 110 set xlabel "Komparatoren"
 111 set xlabel  offset character 0, 0, 0 font "" textcolor lt -1 norotate
 112 set x2label ""
 113 set x2label  offset character 0, 0, 0 font "" textcolor lt -1 norotate
 114 set xrange [ * : * ] noreverse nowriteback  # (currently [63.0000:124.000] )
 115 set x2range [ * : * ] noreverse nowriteback  # (currently [63.0000:124.000] )
 116 set ylabel "Prozent"
 117 set ylabel  offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
 118 set y2label ""
 119 set y2label  offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
 120 set yrange [ * : * ] noreverse nowriteback  # (currently [0.00000:12.0000] )
 121 set y2range [ * : * ] noreverse nowriteback  # (currently [0.00000:11.7344] )
 122 set zlabel ""
 123 set zlabel  offset character 0, 0, 0 font "" textcolor lt -1 norotate
 124 set zrange [ * : * ] noreverse nowriteback  # (currently [-10.0000:10.0000] )
 125 set cblabel ""
 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] )
 128 set zero 1e-08
 129 set lmargin  -1
 130 set bmargin  -1
 131 set rmargin  -1
 132 set tmargin  -1
 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
 139 set colorbox default
 140 set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front bdefault
 141 set loadpath
 142 set fontpath
 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)
 166 GNUTERM = "wxt"
 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))"
 169 n = 29
 170 p = 0.172562132555044
 171 GPFUN_bin = "bin(x) = binv(n,floor(x),p)"
 172 a1 = 0.398942280401433
 173 sigma = 1.88445909011715
 174 mu = 23.800757470569
 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
 180 n_erl = 54
 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)"
 186 offset = 19
 187 GPFUN_poisson_fit = "poisson_fit(x) = 100.0 * poisson(x - offset)"
 188 GPFUN_gauss_fit = "gauss_fit(x) = 100.0 * gauss(x)"
 189 s = 1.39871428165847
 190 t = 23.5320369849105
 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 = 17.9390580825481
 196 gamma_theta = 1.09113152299828
 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))"
 198 gamma_offset = 62
 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))"
 206 y = 0
 207 FIT_CONVERGED = 1
 208 a = 240000
 209 a0 = 2.0
 210 FIT_NDF = 52
 211 FIT_STDFIT = 0.15913245106829
 212 FIT_WSSR = 1.31680312311609
 213 lamdba = 14.0
 214 scale = 5000000
 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=63:124] 'markov-comparators-16-pct.data' title 'SN-Markov', 'evolution-comparators-16-bs-pct.data' title 'SN-Evolution (BS)' lt 3, 'evolution-comparators-16-oes-pct.data' title 'SN-Evolution (OES)' lt 4, gamma_fit(x) title "Gamma-Verteilung" lt 2
 218 plot [x=63:124] 'markov-comparators-16-pct.data' title 'SN-Markov', 'evolution-comparators-16-bs-pct.data' title 'SN-Evolution (BS)' lt 3, 'evolution-comparators-16-oes-pct.data' title 'SN-Evolution (OES)' lt 4
 219 #    EOF