bisect.beta.lim {QFASA} | R Documentation |
Find simultaneous and individual confidence intervals for diet proportions of a single prey species i.e. solve f(pio) = PVAL(pio) = alpha1 and f(pio) = PVAL(pio) = alpha2 using bisection.
Description
Assumptions: alpha1 < alpha2
Usage
bisect.beta.lim(alpha1, alpha2, par.list, R, p.mat, k)
Arguments
alpha1 |
simultaneous confidence level |
alpha2 |
individual confidence level |
par.list |
a list of R.p lists of I beta distribution parameters phi and theta that define diet proportion estimates for each of the prey species. Effectively R.p beta distibutions for each of the I prey species from which we bootstrap to calculate p-values. |
R |
number of bootstrap replicates to use in p-value estimation. |
p.mat |
predator diet estimates. |
k |
prey species index 1..I |
Details
Note: Tried to minimize number of times have to compute a pvalue since very slow.
Value
upper and lower limits for individual and simultaneous intervals respectively.
[Package QFASA version 1.2.0 Index]