[BACK]Return to prob.dem CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gnuplot / demo

Diff for /OpenXM_contrib/gnuplot/demo/Attic/prob.dem between version 1.1.1.1 and 1.1.1.2

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

Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.1.1.2

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>