--------------------------------------------------------------
                Binomial, Bernoulli Trials  SW 3.8
          Albino ex SW Table 3.5, pp. 103-4
number albino children = sum two Bernoulli Trials with prob 1/4               
> dbinom(0:2, s=2, p=1/4)
[1] 0.5625 0.3750 0.0625
> k = 0:2
> list(k, dbinom(k, s=2, p=1/4))
[[1]]
[1] 0 1 2
[[2]]
[1] 0.5625 0.3750 0.0625
> pbinom(1, s=2, p=1/4)
[1] 0.9375

          Mutants SW Ex 3.45 pp106-7
> k = 0:5
> list(k, dbinom(k, s=5, p=.39))
[[1]]
[1] 0 1 2 3 4 5
[[2]]
[1] 0.08445963 0.26999390 0.34523810 0.22072600 0.07055995 0.00902242
    matches SW minitab output

          Blood Type SW Ex 3.47 pp.108-9
> k = 0:6
> list(k, dbinom(k, s=6, p=.85))
[[1]]
[1] 0 1 2 3 4 5 6
[[2]]
[1] 1.139063e-05 3.872813e-04 5.486484e-03 4.145344e-02 1.761771e-01 
                                           3.993348e-01 3.771495e-01
> 1 - dbinom(6, s=6, p=.85)
[1] 0.6228505
----------------------------------------------------------------------
Discrete Distributions:  Binomial-Poisson approx for np, p small
MTB > pdf;                             |MTB > pdf;                                       
SUBC> poisson 1.                       |SUBC> binomial 1000 .001.               
    POISSON WITH MEAN =   1.000        |    BINOMIAL WITH N =1000  P = 0.001000 
       K          P( X = K)            |       K          P( X = K)             
       0            0.3679             |       0            0.3677              
       1            0.3679             |       1            0.3681              
       2            0.1839             |       2            0.1840              
       3            0.0613             |       3            0.0613              
       4            0.0153             |       4            0.0153              
       5            0.0031             |       5            0.0030              
       6            0.0005             |       6            0.0005              
       7            0.0001             |       7            0.0001              
       8            0.0000             |       8            0.0000              
                                       |
                 Agresti Table 1.1 p.5
MTB > pdf;                         MTB > pdf ;
SUBC> poisson 2.                   SUBC> binomial 10 .2.
Poisson with mu = 2.00000          Binomial with n = 10 and p = 0.20 
                                                                     
         x      P( X = x)            x      P( X = x)         
         0        0.1353             0        0.1074          
         1        0.2707             1        0.2684          
         2        0.2707             2        0.3020          
         3        0.1804             3        0.2013          
         4        0.0902             4        0.0881          
         5        0.0361             5        0.0264          
         6        0.0120             6        0.0055          
         7        0.0034             7        0.0008          
         8        0.0009             8        0.0001          
         9        0.0002             9        0.0000          
        10        0.0000     
match less good.