=================================================================== RCS file: /home/cvs/OpenXM_contrib/gnuplot/demo/Attic/prob.dem,v retrieving revision 1.1.1.1 retrieving revision 1.1.1.2 diff -u -p -r1.1.1.1 -r1.1.1.2 --- OpenXM_contrib/gnuplot/demo/Attic/prob.dem 2000/01/09 17:01:01 1.1.1.1 +++ OpenXM_contrib/gnuplot/demo/Attic/prob.dem 2000/01/22 14:16:08 1.1.1.2 @@ -1,5 +1,5 @@ # -# $Id: prob.dem,v 1.1.1.1 2000/01/09 17:01:01 maekawa Exp $ +# $Id: prob.dem,v 1.1.1.2 2000/01/22 14:16:08 maekawa Exp $ # # Demo Statistical Functions version 2.3 # @@ -51,7 +51,7 @@ pause -1 "Hit return to continue" #xmax = 1.0 #Mode of beta PDF used #ymax = (p < 1.0 || q < 1.0) ? 2.0 : 1.1 * beta((p - 1.0)/(p + q - 2.0)) -set key +set key right box set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -68,6 +68,7 @@ plot [0:1] [0:5] p = 0.5, q = 0.7, beta(x) title "p = p = 0.5, q = 2.5, beta(x) title "p = 0.5, q = 2.5" pause -1 "Hit return to continue" set title "incomplete beta CDF" +set key left box plot [0:1] [0:1.1] p = 0.5, q = 0.7, cbeta(x) title "p = 0.5, q = 0.7", \ p = 5.0, q = 3.0, cbeta(x) title "p = 5.0, q = 3.0", \ p = 0.5, q = 2.5, cbeta(x) title "p = 0.5, q = 2.5" @@ -109,7 +110,7 @@ pause -1 "Hit return to continue" #xmin = a - 4.0 * b #xmax = a + 4.0 * b #ymax = 1.1 * cauchy(a) #Mode of cauchy PDF used -set key +set key left box set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -137,7 +138,7 @@ pause -1 "Hit return to continue" #xmin = xmin < 0 ? 0 : xmin #xmax = mu + 4.0 * sigma #ymax = 1.1 * (df1 > 2.0 ? chi(df1 - 2.0) : 1.0) #Mode of chi PDF used -set key +set key right box set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -149,10 +150,12 @@ set format x "%.1f" set format y "%.2f" set sample 100 set title "chi-square PDF" +set key right box plot [0:15] [0:0.2] df1 = 4, chi(x) title "df = 4", \ df1 = 6, chi(x) title "df = 6", \ df1 = 8, chi(x) title "df = 8" pause -1 "Hit return to continue" +set key left box set title "chi-square CDF" plot [0:15] [0:1.1] df1 = 4, cchi(x) title "df = 4", \ df1 = 6, cchi(x) title "df = 6", \ @@ -167,7 +170,6 @@ pause -1 "Hit return to continue" #xmin = xmin < 0 ? 0 : xmin #xmax = mu + 4.0 * sigma #ymax = n < 2.0 ? 1.0 : 1.1 * erlang((n - 1.0) / lambda) #Mode of erlang PDF used -set key set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -179,10 +181,12 @@ set format x "%.1f" set format y "%.1f" set sample 100 set title "erlang PDF" +set key right box plot [0:10] [0:1] lambda = 1, n = 2, erlang(x) title "lambda = 1, n = 2", \ lambda = 2, n = 2, erlang(x) title "lambda = 2, n = 2" pause -1 "Hit return to continue" set title "erlang CDF" +set key left box plot [0:10] [0:1.1] lambda = 1, n = 2, cerlang(x) title "lambda = 1, n = 2", \ lambda = 2, n = 2, cerlang(x) title "lambda = 2, n = 2" pause -1 "Hit return to continue" @@ -195,7 +199,6 @@ pause -1 "Hit return to continue" #xmin = mu - 4.0 * sigma #xmax = mu + 4.0 * sigma #ymax = 1.1 * extreme(u) #Mode of extreme PDF used -set key set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -207,10 +210,12 @@ set format x "%.1f" set format y "%.2f" set sample 100 set title "extreme PDF" +set key left box plot [-10:10] [0:0.4] alpha = 0.5, u = 1.0, extreme(x) title "alpha = 0.5, u = 1.0", \ alpha = 1.0, u = 0.0, extreme(x) title "alpha = 1.0, u = 0.0" pause -1 "Hit return to continue" set title "extreme CDF" +set key left box plot [-10:10] [0:1.1] alpha = 0.5, u = 1.0, cextreme(x) title "alpha = 0.5, u = 1.0", \ alpha = 1.0, u = 0.0, cextreme(x) title "alpha = 1.0, u = 0.0" pause -1 "Hit return to continue" @@ -224,7 +229,6 @@ pause -1 "Hit return to continue" #xmax = mu + 4.0 * sigma #Mode of F PDF used #ymax = df1 < 3.0 ? 1.0 : 1.1 * f((df1 / 2.0 - 1.0) / (df1 / 2.0 + df1 / df2)) -set key set zeroaxis #set xrange [xmin : xmax] #set yrange [0 : ymax] @@ -236,10 +240,12 @@ set format x "%.1f" set format y "%.2f" set sample 100 set title "F PDF" +set key right box plot [0:4] [0:0.8] df1 = 5.0, df2 = 9.0, f(x) title "df1 = 5, df2 = 9", \ df1 = 7.0, df2 = 6.0, f(x) title "df1 = 7, df2 = 6" pause -1 "Hit return to continue" set title "F CDF" +set key left box plot [0:4] [0:1.1] df1 = 5.0, df2 = 9.0, cf(x) title "df1 = 5, df2 = 9", \ df1 = 7.0, df2 = 6.0, cf(x) title "df1 = 7, df2 = 6" pause -1 "Hit return to continue" @@ -252,7 +258,6 @@ pause -1 "Hit return to continue" #xmin = xmin < 0 ? 0 : xmin #xmax = mu + 4.0 * sigma #ymax = rho < 1.0 ? 2.0 : 1.1 * g((rho - 1.0) / lambda) #Mode of gamma pdf used -set key set zeroaxis #set xrange [xmin: xmax] #set yrange [0: ymax] @@ -264,11 +269,13 @@ set format x "%.1f" set format y "%.1f" set sample 100 set title "gamma PDF" +set key right plot [0:5] [0:1.5] rho = 0.5, lambda = 1.0, g(x) title "rho = 0.5, lambda = 1.0", \ rho = 1.0, lambda = 1.0, g(x) title "rho = 1.0, lambda = 1.0", \ rho = 2.0, lambda = 2.0, g(x) title "rho = 2.0, lambda = 2.0" pause -1 "Hit return to continue" set title "incomplete gamma CDF (lambda == 1.0)" +set key right bottom plot [0:5] [0:1.1] rho = 0.5, cgamma(x) title "rho = 0.5", \ rho = 1.0, cgamma(x) title "rho = 1.0", \ rho = 2.0, cgamma(x) title "rho = 2.0" @@ -450,7 +457,6 @@ pause -1 "Hit return to continue" #xmin = xmin < 0 ? 0 : xmin #xmax = mu + 4.0 * sigma #ymax = 1.1 * maxwell(1.0 / a) #Mode of maxwell PDF used -set key set zeroaxis #set xrange[xmin: xmax] #set yrange[0: ymax] @@ -462,11 +468,13 @@ set format x "%.1f" set format y "%.1f" set sample 100 set title "maxwell PDF" +set key right top box plot [0:6] [0:1.4] a = 1.5, maxwell(x) title "a = 1.5", \ a = 1.0, maxwell(x) title "a = 1.0", \ a = 0.5, maxwell(x) title "a = 0.5" pause -1 "Hit return to continue" set title "maxwell CDF" +set key right bottom box plot [0:6] [0:1.1] a = 1.5, cmaxwell(x) title "a = 1.5", \ a = 1.0, cmaxwell(x) title "a = 1.0", \ a = 0.5, cmaxwell(x) title "a = 0.5" @@ -532,7 +540,6 @@ pause -1 "Hit return to continue" #xmin = mu - 4.0 * sigma #xmax = mu + 4.0 * sigma #ymax = 1.1 * normal(mu) #Mode of normal PDF used -set key set zeroaxis #set xrange [xmin: xmax] #set yrange [0: ymax] @@ -544,11 +551,13 @@ set format x "%.1f" set format y "%.1f" set sample 100 set title "normal (also called gauss or bell-curved) PDF" +set key left top box plot [-4:4] [0:1] mu = 0, sigma = 1.0, normal(x) title "mu = 0, sigma = 1.0", \ mu = 2, sigma = 0.5, normal(x) title "mu = 2, sigma = 0.5", \ mu = 1, sigma = 2.0, normal(x) title "mu = 1, sigma = 2.0" pause -1 "Hit return to continue" set title "normal (also called gauss or bell-curved) CDF" +set key left top box plot [-4:4] [0:1.1] mu = 0, sigma = 1.0, cnormal(x) title "mu = 0, sigma = 1.0", \ mu = 2, sigma = 0.5, cnormal(x) title "mu = 2, sigma = 0.5", \ mu = 1, sigma = 2.0, cnormal(x) title "mu = 1, sigma = 2.0" @@ -642,7 +651,6 @@ pause -1 "Hit return to continue" #xmin = 0.0 #xmax = a #ymax = 1.1 * 2.0 / a #Mode of sine PDF used -set key set zeroaxis #set xrange [xmin: xmax] #set yrange [0: ymax] @@ -654,10 +662,12 @@ set format x "%.2f" set format y "%.1f" set sample 100 set title "sine PDF" +set key bottom outside plot [0:2] [0:1.1] a = 2.0, n = 1, sine(x) title "a = 2.0, n = 1", \ a = 2.0, n = 3, sine(x) title "a = 2.0, n = 3" pause -1 "Hit return to continue" set title "sine CDF" +set key top left plot [0:2] [0:1.1] a = 2.0, n = 1, csine(x) title "a = 2.0, n = 1", \ a = 2.0, n = 3, csine(x) title "a = 2.0, n = 3" pause -1 "Hit return to continue" @@ -773,3 +783,4 @@ plot [0:3] [0:1.2] lambda = 1, n = 0.5, cweibull(x) ti lambda = 1, n = 1.0, cweibull(x) title "lambda = 1, n = 1.0", \ lambda = 1, n = 2.0, cweibull(x) title "lambda = 1, n = 2.0", \ lambda = 3, n = 2.0, cweibull(x) title "lambda = 3, n = 2.0" +reset