# # Written by: # -- # John L. Weatherwax 2009-04-21 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # # EPage 788 # #----- set.seed(0) epsilon = 1e-1 xs = seq( 0, 1, length.out=200 ) plot( xs, dbeta( xs, epsilon, epsilon ), main='original beta[eps,eps] distribution' ) plot( xs, dbeta( xs, epsilon+1, epsilon ), main='distribution with one positive example' ) plot( xs, dbeta( xs, epsilon+1, epsilon+1 ), main='with one positive and one negative example' )