source("biossay.S") DD_input(LL=1000) PP _ grid.value(DD$biossay,DD$alpha,DD$beta) draws_sampling(M=1000,PP,DD$alpha,DD$beta) plot.joint.post(w=0,DD$bioassay,draws)