Dear all,
First, I apologize for this long-winded post.
I am more seeking to confirm that I am doing things correct for a simulation than anything. I am aware of the assumptions I am making for this closed population model, but interested in making sure my methodology makes sense.
The scenario:
Basically I have used a Huggins CC model to assess detection probability of a known population (lizards that were surveyed in enclosures within their natural habitat), and it is really low (<0.1) for both males and females. I know these low detection probabilities are of limited use, but my goal is to simulate encounter histories for all individuals in populations of size "x" for two and three occasions. Then I will be using those encounter histories to generate hypothetical population estimates for these "known" population sizes using a simple lincoln-peterson calculation to show how much we typically underestimate population sizes of this species in this area. I am perfectly fine with MARK simply outputting the N(hat) for a number of iterations of these simulations.
What I have done:
I am simulating a Huggins CC Model, simply p(g)=c(g) where g=sex. I set up 2 groups, and 2 encounter occasions. I am inputting the respective detection probabilities for all parameters relating to males and females, and checking the boxes for "Derived Estimates" and "Input Data in Output", and running 100 simulations for 100 Males and 100 Females. I am getting Deriv1 and Deriv2 in my output which are reading numbers that I would generally expect (~25 males, ~17 females).
Does this make sense? Or am I missing something?
Thank you to all who take the time to look at this- I greatly appreciate your input.
Best,
Mike