| mr.closed1.sim {FSA} | R Documentation |
This function is used to simulate multiple results for a single census mark-recapture data. The user can control many aspects of the simulation including simulating the loss of marks, differential survival rates for tagged and untagged fish, the addition of recruits, and differential probabilities of capture for tagged and untagged fish. In addition, Petersen, Chapman, and Bailey estimates of population size can be made.
mr.closed1.sim(type=c("Petersen","P","Chapman","C","Ricker","R","Bailey","B"),
N=1000, rsmpls=500, incl.final=TRUE)
type |
a string that identifies the type of calculation method to use (see details). |
N |
a number that represents the ``known'' size of the simulated population just prior to the first sample. |
rsmpls |
a number indicating the number of simulations to run. |
incl.final |
a logical indicating whether the mean final population size and histogram of percent error from final population size should be shown. |
The user can use slider bars to choose values for the expected number of fish to capture in the first (i.e., the marking or tagging) sample, the expected number of fish to be captured in the second or final (i.e., the recapture) sample, a probablity that a tagged fish loses the tag, the probability of survival for tagged fish, the probability of survival for untagged fish, a proportion of the initial population size to recruit to the population between the first and final samples, and the ratio of the probability of capture of tagged fish to the probability of capture of untagged fish in the final sample. The user selects the “true” initial population size, the method of population estimator, and the number of resamples to use with arguments to the function. The simulation can be re-run without changing any of the slider buttons by pressing the ‘Re-Randomize’ button.
In general, the simulation follows these steps: (1) a population of N fish is created; (2) randomly select M fish in the first sample to be tagged such that the average of all M values should be equal to the user-supplied ‘Tagged (M)’ value; (3) randomly select m fish to lose tags according to the user-supplied probability of tag loss (‘PR(Tag Loss)’); (4) randomly select fish to die according to the user-supplied probabilities of survival for tagged (‘PR(Surv Tagged)’) and untagged fish (‘PR(Surv UNTagged)’); (5) add recruits to the population in proportion to N and in accordance to the user-supplied proportion of recruits (‘Proportion Recruit’) to add; (6) identify the actual population size just before taking the final sample (N1); (7) randomly select n fish in the second sample such that the average of all n values should be equal to the user-supplied ‘Captured (n)’ value (see note below about how differential probabilties of capture are incorporated into the model); (8) count the number of tagged fish in the second sample; and (9) compute the population estimate using the chosen method (see below).
The methods for estimating the population size are
type="P" | naive Petersen. |
type="C" | Chapman(1951) modification of the Petersen. |
type="CR" | Ricker(1975) modification of the Chapman modification. |
type="B" | Bailey(1951,1952) modification of the Petersen. |
The effect of violating the assumption of tag loss can be simulated by changing the probability of tag loss slider to a value greater than 0 but less than 1. For example, setting ‘PR(Tag Loss)’ to 0.1 is used to simulate a 10 percent probability of losing the tag.
The effect of violating the assumption of no mortality from the first to second sample can be simulated by changing one or both of the probabilities of survival for tagged and untagged fish to values less than 1 (but greater than 0). For example, setting ‘PR(Surv Tagged)’ AND ‘PR(Surv UNTagged)’ to 0.8 will simulate mortality between the first and final sample but NOT differential mortality between the tagged and untagged fish (i.e., the mortalities are the same for both groups of fish). The effect of differential mortalities between tagged and untagged fish can be simulated by using different survival probabilities for tagged and untagged fish.
The effect of violating the assumption of no recruitment from the first to final sample can be simulated by changing the Proportion Recruit slider to a value greater than 0. For example, setting ‘Proportion Recruit’ to 0.1 will simulate 10 percent of N recruiting to the population just before the final sample.
The effect of violating the assumption of equal catchabilities in the final sample for tagged and untagged fish can be simulated by changing the ‘Ratio PR(Capture)’ slider to a value different than 1. Values greater than 1 indicate that the catchability of tagged fish is greater than the catchability of untagged fish. Values less than 1 indicate that the catchability of tagged fish is less than that of untagged fish. For example, setting ‘Ratio PR(Capture)’ to 0.8 will simulate a situation where the capture probability of tagged fish is 80 percent of the capture probablity of untagged fish (i.e., simulates the situation where marked fish are less likely to be captured).
The probability of capture in the final sample is equal to the expected number of fish to be collected in the final sample (set with ‘Captured (n)’) divided by the actual population size just prior to the final sample (N1) as long as ‘Ratio PR(Capture)’ is 1. The probabilities of capture for the tagged and untagged fish must be carefully adjusted if the ‘Ratio PR(Capture)’ value is different than 1. Because of the different numbers of tagged and untagged individuals in the population, the probability of capture for tagged and untagged individuals must be computed by adjusting the overall probability of capture to assure that, on average, the user-provided expected number captured in the final sample is met. This modification is found by solving the following system of equations for the probabilities of capture for the tagged (PM) and untagged fish (PU), respectively,
PM*M + PU*(N-M) = P*N
frac{PM}{PU} = k
where M is the number of tagged animals, and N-M is the number of untagged animals, k is the ‘Ratio PR(Capture)’ value, and P is the overall probability of capture if there was no difference in catchability between tagged and untagged animals. The solutions to this system are
PU = frac{PN}{M+N}
PM = PU*k
These solutions are used in the simulation model.
None. However, a dynamic graphic is produced that is controlled by slider bars as described in the details. The dynamic graphic is a histogram is the population estimate from all resamples with a red vertical line at the initial population size (provided by the user), a blue vertical line at the population size just prior to the final sample (N1; if incl.final=TRUE), and a green vertical line at the mean population estimate. The vertical blue and red lines may not be visible under some scenarios because of overprinting.
Derek H. Ogle, dogle@northland.edu
## Not run: ## This example cannot be run by examples() but should work in an interactive R session ## Not run: mr.closed1.sim() ## Not run: mr.closed1.sim(type="C") # use Chapman modification