SN-Markov-Graphen ausgebaut.

[diplomarbeit.git] / images / markov-comparators-12-pct.gnuplot

diff --git a/images/markov-comparators-12-pct.gnuplot b/images/markov-comparators-12-pct.gnuplot
index 718b7a0..bde07ad 100644 (file)
--- a/images/markov-comparators-12-pct.gnuplot
+++ b/images/markov-comparators-12-pct.gnuplot
@@ -13,8 +13,8 @@
 #      faq, bugs, etc:   type "help seeking-assistance"
 #      immediate help:   type "help"
 #      plot window:      hit 'h'
-# set terminal wxt 0
-# set output
+set terminal pdfcairo  size 15.00cm, 9.27cm 
+set output 'markov-comparators-12-pct.pdf'
 unset clip points
 set clip one
 unset clip two
@@ -44,7 +44,7 @@ set format cb "% g"
 set angles radians
 unset grid
 set key title ""
-set key inside right top vertical Right noreverse enhanced autotitles nobox
+set key outside center bottom horizontal Right noreverse enhanced autotitles box linetype -1 linewidth 1.000
 set key noinvert samplen 4 spacing 1 width 0 height 0 
 unset label
 unset arrow
@@ -105,8 +105,8 @@ set timestamp ""
 set timestamp  offset character 0, 0, 0 font "" norotate
 set rrange [ * : * ] noreverse nowriteback  # (currently [8.98847e+307:-8.98847e+307] )
 set trange [ * : * ] noreverse nowriteback  # (currently [-5.00000:5.00000] )
-set urange [ * : * ] noreverse nowriteback  # (currently [-5.00000:5.00000] )
-set vrange [ * : * ] noreverse nowriteback  # (currently [-5.00000:5.00000] )
+set urange [ * : * ] noreverse nowriteback  # (currently [-10.0000:10.0000] )
+set vrange [ * : * ] noreverse nowriteback  # (currently [-10.0000:10.0000] )
 set xlabel "Komparatoren" 
 set xlabel  offset character 0, 0, 0 font "" textcolor lt -1 norotate
 set x2label "" 
@@ -121,7 +121,7 @@ set yrange [ * : * ] noreverse nowriteback  # (currently [0.00000:16.0000] )
 set y2range [ * : * ] noreverse nowriteback  # (currently [0.00000:14.0879] )
 set zlabel "" 
 set zlabel  offset character 0, 0, 0 font "" textcolor lt -1 norotate
-set zrange [ * : * ] noreverse nowriteback  # (currently [8.98847e+307:-8.98847e+307] )
+set zrange [ * : * ] noreverse nowriteback  # (currently [-10.0000:10.0000] )
 set cblabel "" 
 set cblabel  offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
 set cbrange [ * : * ] noreverse nowriteback  # (currently [8.98847e+307:-8.98847e+307] )
@@ -144,9 +144,9 @@ set fit noerrorvariables
 gamma_gamma(x) = (floor(x - 0.5))!
 gamma_dist(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma(gamma_k))
 gamma_fit(x) = 100.0 * dgamma(x - gamma_offset, gamma_k, gamma_theta)
+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))
 _ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x
 pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)
-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))
 GNUTERM = "wxt"
 GPFUN_gamma_gamma = "gamma_gamma(x) = (floor(x - 0.5))!"
 gamma_k = 8.26729758801638
@@ -154,13 +154,12 @@ gamma_theta = 0.961930249498844
 GPFUN_gamma_dist = "gamma_dist(x) = x**(gamma_k - 1) * exp((-1.0) * x / gamma_theta) / (gamma_theta**gamma_k * gamma(gamma_k))"
 gamma_offset = 40
 GPFUN_gamma_fit = "gamma_fit(x) = 100.0 * dgamma(x - gamma_offset, gamma_k, gamma_theta)"
+GPFUN__ln_dgamma = "_ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x"
+GPFUN_pgamma = "pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)"
+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))"
 FIT_CONVERGED = 1
 FIT_NDF = 39
 FIT_STDFIT = 0.222936770551606
 FIT_WSSR = 1.9383313428952
-GPFUN__ln_dgamma = "_ln_dgamma(x, a, b) = a*log(b) - lgamma(a) + (a-1)*log(x) - b*x"
-GPFUN_pgamma = "pgamma(x, shape, rate) = (x<0)? 0 : igamma(shape, x*rate)"
-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))"
-plot [x=41:81] './images/markov-comparators-12-pct.data' title 'Gemessene Daten', gamma_fit(x) title 'Gamma-Verteilung'
-## fit gamma_fit(x) 'images/markov-comparators-12-pct.data' via gamma_k, gamma_theta
+plot [x=41:81] 'markov-comparators-12-pct.data' title 'Gemessene Daten', gamma_fit(x) title 'Gamma-Verteilung'
 #    EOF