Find root (i.e. solve f(pio) = PVAL(pio) = alpha1 using either uniroot or optimize.
root.beta(x1, x2, alpha, par.list, R, p.mat, k, low = TRUE)root which will be lower or upper CI limit
lower starting value
upper starting value
simultaneous confidence level
parameters needed to compute p-value by bootstrapping
number of bootstraps
predator diet estimates
prey species index 1...I
TRUE if computing lower interval.