# NOT RUN {
##
## Example: Distribution Function
## Creates Table I in Swed and Eisenhart (1943), p. 70,
## with n1 = 2 and n1 <= n2 <= 20
##
m <- NULL
for (i in 2:20){
m <- rbind(m, pruns(2:5,2,i))
}
rownames(m)=2:20
colnames(m)=2:5
#
# 2 3 4 5
# 2 0.333333333 0.6666667 1.0000000 1
# 3 0.200000000 0.5000000 0.9000000 1
# 4 0.133333333 0.4000000 0.8000000 1
# 5 0.095238095 0.3333333 0.7142857 1
# 6 0.071428571 0.2857143 0.6428571 1
# 7 0.055555556 0.2500000 0.5833333 1
# 8 0.044444444 0.2222222 0.5333333 1
# 9 0.036363636 0.2000000 0.4909091 1
# 10 0.030303030 0.1818182 0.4545455 1
# 11 0.025641026 0.1666667 0.4230769 1
# 12 0.021978022 0.1538462 0.3956044 1
# 13 0.019047619 0.1428571 0.3714286 1
# 14 0.016666667 0.1333333 0.3500000 1
# 15 0.014705882 0.1250000 0.3308824 1
# 16 0.013071895 0.1176471 0.3137255 1
# 17 0.011695906 0.1111111 0.2982456 1
# 18 0.010526316 0.1052632 0.2842105 1
# 19 0.009523810 0.1000000 0.2714286 1
# 20 0.008658009 0.0952381 0.2597403 1
#
# }
Run the code above in your browser using DataLab