First, thank you so much for your response!
We used program MARK after capturing 14 individuals in a particular trapping session to estimate the actual number of individuals at 23 with 95% confidence interval of 18-36. We need to run this on multiple trapping sessions with different numbers of individuals captured and different CIs for each one.
The probability of that number being 23 is unknown and when I used binom functions, I just guessed to see if it would give us anything of use. We need the probability of every number between 18 and 36 so that we can use this in an energetics model sampling from this distribution each time we run 10,000 MCMC replicates.
I have tried:
q <- seq(18:36)
pbinom(q, 23, 0.6)
and got probabilities that make no sense because 36 is more likely than 23.
When I try
rbinom(100, 23, 0.6)
I get a bunch of numbers lower than 23 and certainly lower than 18 so it's not giving me the distribution of numbers I am looking for.
I have tried using the package: sn but my numbers are out of the feasible range for the cp2dp command.
Does that help clarify things? Again, thank you so much for your help!