# # 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. # #----- Ms = c(5,10,20) eps = c(0.1,0.2,0.4) for( m in Ms ){ for( e in eps ){ print( sprintf( "%5d %8.2f %10.6f", m, e, 1 - pbinom( floor( m/2 ), m, e ) ) ) } }