NB. beta binomial distribution NB. written by Giichiro Suzuki /Feb/2011 mean=:[:+/i.@#*]%+/ var=:[:+/([:*:i.@#-mean)*]%+/ bden=:3 :'(k!{:y)*(p^k)*(1-p=.{.y)^|.k=.i.>:{:y' NB. -------------------------------------------- S=:0 0 6 13 24 38 53 37 102 54 10 17 17 7 5 3 AA=:386*bden 0.5 15 A=:".3": AA NB. --------------------------------------------- beta=:([:*/!@<:"0)%[:!<:@+/ bbden=:3 :0 b=.(k+{.c),.({:c=.>{.y)+|.k=.i.>:n=.>{:y (k!n)*(beta"1 b)%beta c ) NB. --------------------------------------------- BB=:386*bbden(2#10.443);15 B=:".3": BB W=:".3":386*(0.3*bden 0.3 15)+(0.4*bden 0.5 15)+0.3*bden 0.7 15