version 1.1.1.1, 2000/01/09 17:01:02 |
version 1.1.1.3, 2003/09/15 07:09:29 |
|
|
# |
# |
# Demo Statistical Approximations version 1.1 |
# Demo Statistical Approximations version 1.1 |
# |
# |
# Permission granted to distribute freely for non-commercial purposes only |
|
# |
|
# Copyright (c) 1991, Jos van der Woude, jvdwoude@hut.nl |
# Copyright (c) 1991, Jos van der Woude, jvdwoude@hut.nl |
|
|
pause 0 "" |
pause 0 "" |
|
|
pause 0 " Copyright (c) 1991, 1992, Jos van de Woude, jvdwoude@hut.nl" |
pause 0 " Copyright (c) 1991, 1992, Jos van de Woude, jvdwoude@hut.nl" |
pause 0 "" |
pause 0 "" |
pause 0 "" |
pause 0 "" |
pause 0 " Permission granted to distribute freely for non-commercial purposes only" |
|
pause 0 "" |
pause 0 "" |
pause 0 "" |
pause 0 "" |
pause 0 "" |
pause 0 "" |
Line 45 xmax = ceil(mu + 4.0 * sigma) |
|
Line 42 xmax = ceil(mu + 4.0 * sigma) |
|
ymax = 1.1 * binom(mu) #mode of binomial PDF used |
ymax = 1.1 * binom(mu) #mode of binomial PDF used |
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 77 xmax = ceil(mu + 4.0 * sigma) |
|
Line 74 xmax = ceil(mu + 4.0 * sigma) |
|
ymax = 1.1 * binom(mu) #mode of binomial PDF used |
ymax = 1.1 * binom(mu) #mode of binomial PDF used |
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 111 xmax = ceil(mu + 4.0 * sigma) |
|
Line 108 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * p |
ymax = 1.1 * p |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 143 xmax = ceil(mu + 4.0 * sigma) |
|
Line 140 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * p |
ymax = 1.1 * p |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 176 xmax = ceil(mu + 4.0 * sigma) |
|
Line 173 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * hypgeo(mu) #mode of binom PDF used |
ymax = 1.1 * hypgeo(mu) #mode of binom PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 209 xmax = ceil(mu + 4.0 * sigma) |
|
Line 206 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * hypgeo(mu) #mode of binom PDF used |
ymax = 1.1 * hypgeo(mu) #mode of binom PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 243 xmax = int(mu + 4.0 * sigma) |
|
Line 240 xmax = int(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * g((rho - 1) / lambda) #mode of gamma PDF used |
ymax = 1.1 * g((rho - 1) / lambda) #mode of gamma PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 275 xmax = ceil(mu + 4.0 * sigma) |
|
Line 272 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * negbin(mu - 1/p) #mode of gamma PDF used |
ymax = 1.1 * negbin(mu - 1/p) #mode of gamma PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 304 lambda = pi / (sqrt(3.0) * sigma) |
|
Line 301 lambda = pi / (sqrt(3.0) * sigma) |
|
xmin = mu - 4.0 * sigma |
xmin = mu - 4.0 * sigma |
xmax = mu + 4.0 * sigma |
xmax = mu + 4.0 * sigma |
ymax = 1.1 * logistic(mu) #mode of logistic PDF used |
ymax = 1.1 * logistic(mu) #mode of logistic PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin: xmax] |
set xrange [xmin: xmax] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 335 xmax = ceil(mu + 4.0 * sigma) |
|
Line 332 xmax = ceil(mu + 4.0 * sigma) |
|
xinc = ceil((xmax - xmin) / 10) |
xinc = ceil((xmax - xmin) / 10) |
xinc = xinc > 1 ? xinc : 1 |
xinc = xinc > 1 ? xinc : 1 |
ymax = 1.1 * poisson(mu) #mode of poisson PDF used |
ymax = 1.1 * poisson(mu) #mode of poisson PDF used |
set key |
set key box |
set nozeroaxis |
set nozeroaxis |
set xrange [xmin - 1 : xmax + 1] |
set xrange [xmin - 1 : xmax + 1] |
set yrange [0 : ymax] |
set yrange [0 : ymax] |
Line 353 set arrow from mu + 0.5, normal(mu + sigma) \ |
|
Line 350 set arrow from mu + 0.5, normal(mu + sigma) \ |
|
set label "mu" at mu + 1, ymax / 10 |
set label "mu" at mu + 1, ymax / 10 |
set label "sigma" at mu + 1 + sigma, normal(mu + sigma) |
set label "sigma" at mu + 1 + sigma, normal(mu + sigma) |
plot poisson(x), normal(x - 0.5) |
plot poisson(x), normal(x - 0.5) |
set noarrow |
reset |
set nolabel |
|