X-Git-Url: https://git.octo.it/?a=blobdiff_plain;f=images%2Fmarkov-comparators-12-pct.gnuplot;h=bde07ade9b3320d9d05e88b009c82acdea4522cf;hb=fc5785676841c1186ea1240150e9e27622fbabd2;hp=718b7a008adc1a3df84e540ac9981f28de6228bc;hpb=f0c46ce4d21d333772c7a967b3c8965379fd8699;p=diplomarbeit.git diff --git a/images/markov-comparators-12-pct.gnuplot b/images/markov-comparators-12-pct.gnuplot index 718b7a0..bde07ad 100644 --- 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