=================================================================== RCS file: /home/cvs/OpenXM_contrib/gnuplot/demo/Attic/prob.dem,v retrieving revision 1.1.1.1 retrieving revision 1.1.1.3 diff -u -p -r1.1.1.1 -r1.1.1.3 --- OpenXM_contrib/gnuplot/demo/Attic/prob.dem 2000/01/09 17:01:01 1.1.1.1 +++ OpenXM_contrib/gnuplot/demo/Attic/prob.dem 2003/09/15 07:09:29 1.1.1.3 @@ -1,15 +1,12 @@ # -# $Id: prob.dem,v 1.1.1.1 2000/01/09 17:01:01 maekawa Exp $ +# $Id: prob.dem,v 1.1.1.3 2003/09/15 07:09:29 ohara Exp $ # # Demo Statistical Functions version 2.3 # -# Permission granted to distribute freely for non-commercial purposes only -# # Copyright (c) 1991, 1992 Jos van der Woude, jvdwoude@hut.nl print " Statistical Library Demo, version 2.3" print "\n Copyright (c) 1991, 1992, Jos van de Woude, jvdwoude@hut.nl" -print "Permission granted to distribute freely for non-commercial purposes only" print "\n\n\n\n\n\n\n" print "NOTE: contains 54 plots and consequently takes a lot of time to run" print " Press Ctrl-C to exit right now" @@ -51,7 +48,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 +65,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 +107,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 +135,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 +147,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 +167,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 +178,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 +196,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 +207,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 +226,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 +237,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 +255,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 +266,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 +454,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 +465,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 +537,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 +548,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 +648,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 +659,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 +780,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