Learn R Programming

bayess (version 1.4)

thresh: Bound for the accept-reject algorithm in Chapter 5

Description

This function is used in ardipper to determine the bound for the accept-reject algorithm simulating the non-standard conditional distribution of $r_1$.

Usage

thresh(k, n1, c2, c3, r2, q1)

Arguments

k
current proposal for the number of individuals vanishing between the first and second experiments
n1
first capture population size
c2
number of individuals recaptured during the second experiment
c3
number of individuals recaptured during the third experiment
r2
number of individuals vanishing between the second and third experiments
q1
probability of disappearing from the population

Value

  • numerical value of the upper bound, to be compared with the uniform random draw

Details

This upper bound is equal to $$\frac{{n_1-c_2 \choose k} {n_1-k \choose c_3+r_2}}{{\bar r\choose k}}$$

See Also

ardipper

Examples

Run this code
if (runif(1) < thresh(y,n1,c2,c3,r2,q1))

Run the code above in your browser using DataLab