version 1.1.1.1, 2000/01/09 17:01:01 |
version 1.1.1.3, 2003/09/15 07:09:29 |
|
|
# |
# |
# Demo Statistical Functions version 2.3 |
# 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 |
# Copyright (c) 1991, 1992 Jos van der Woude, jvdwoude@hut.nl |
|
|
print " Statistical Library Demo, version 2.3" |
print " Statistical Library Demo, version 2.3" |
print "\n Copyright (c) 1991, 1992, Jos van de Woude, jvdwoude@hut.nl" |
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 "\n\n\n\n\n\n\n" |
print "NOTE: contains 54 plots and consequently takes a lot of time to run" |
print "NOTE: contains 54 plots and consequently takes a lot of time to run" |
print " Press Ctrl-C to exit right now" |
print " Press Ctrl-C to exit right now" |
Line 51 pause -1 "Hit return to continue" |
|
Line 48 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 65 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 107 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 135 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 147 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 167 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 178 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 196 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 207 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 226 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 237 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 255 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 266 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 454 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 465 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 537 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 548 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 648 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 659 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 780 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 |