## ----------------------------------------------------------------------------- # give the function a name that will not conflict with others in the namespace beta.stat <- function(x) { -sum(log(x)) } ## ----------------------------------------------------------------------------- require(cdfinv) cdfinv.sim("beta","shape1",10.5,bound="lower",lpb=0,nsamp=6,stat.func=beta.stat,shape2=2.8)