bisect.beta.lim: Find simultaneous and individual confidence intervals for...

In justinkamerman/QFASA: Quantitative Fatty Acid Signature Analysis

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.

justinkamerman/QFASA documentation built on Nov. 17, 2023, 12:33 p.m.