sim_cdf.R
In cdfinv: Confidence Interval Estimation via CDF Inversion

## -----------------------------------------------------------------------------
# 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)