<!DOCTYPE html>
<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"><meta name="description" content="Estimates of the initial population size, along with associated confidence intervals, are constructed from single or multiple census mark-recapture data using a variety of methods. For single census data, the initial population size (N) is estimated from the number of marked animals from a first sample (M), number of captured animals in a second sample (n), and the number of recaptured marked animals in the second sample (m) using either the &amp;#8216;naive&amp;#8217; Petersen method or Chapman, Ricker, or Bailey modifications of the Petersen method. Single census data can also be separated by group (e.g., size class) to estimate the initial population size by class and for the overall population size. For multiple census data, the initial population size is estimated from the number of captured animals (n), number of recaptured marked animals (m), the number of marked animals that are marked and returned to the population (R), or the number of extant marked animals prior to the sample (M) on each of several samples using either the Schnabel (1938) or Schumacher-Eschmeyer (1943) method."><title>Estimate initial population size for single or multiple census mark-recapture data. — mrClosed • FSA</title><!-- favicons --><link rel="icon" type="image/png" sizes="16x16" href="../favicon-16x16.png"><link rel="icon" type="image/png" sizes="32x32" href="../favicon-32x32.png"><link rel="apple-touch-icon" type="image/png" sizes="180x180" href="../apple-touch-icon.png"><link rel="apple-touch-icon" type="image/png" sizes="120x120" href="../apple-touch-icon-120x120.png"><link rel="apple-touch-icon" type="image/png" sizes="76x76" href="../apple-touch-icon-76x76.png"><link rel="apple-touch-icon" type="image/png" sizes="60x60" href="../apple-touch-icon-60x60.png"><script src="../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"><link href="../deps/bootstrap-5.2.2/bootstrap.min.css" rel="stylesheet"><script src="../deps/bootstrap-5.2.2/bootstrap.bundle.min.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- bootstrap-toc --><script src="https://cdn.jsdelivr.net/gh/afeld/bootstrap-toc@v1.0.1/dist/bootstrap-toc.min.js" integrity="sha256-4veVQbu7//Lk5TSmc7YV48MxtMy98e26cf5MrgZYnwo=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- search --><script src="https://cdnjs.cloudflare.com/ajax/libs/fuse.js/6.4.6/fuse.js" integrity="sha512-zv6Ywkjyktsohkbp9bb45V6tEMoWhzFzXis+LrMehmJZZSys19Yxf1dopHx7WzIKxr5tK2dVcYmaCk2uqdjF4A==" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/autocomplete.js/0.38.0/autocomplete.jquery.min.js" integrity="sha512-GU9ayf+66Xx2TmpxqJpliWbT5PiGYxpaG8rfnBEk1LL8l1KGkRShhngwdXK1UgqhAzWpZHSiYPc09/NwDQIGyg==" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mark.js/8.11.1/mark.min.js" integrity="sha512-5CYOlHXGh6QpOFA/TeTylKLWfB3ftPsde7AnmhuitiTX4K5SqCLBeKro6sPS8ilsz1Q4NRx3v8Ko2IBiszzdww==" crossorigin="anonymous"></script><!-- pkgdown --><script src="../pkgdown.js"></script><meta property="og:title" content="Estimate initial population size for single or multiple census mark-recapture data. — mrClosed"><meta property="og:description" content="Estimates of the initial population size, along with associated confidence intervals, are constructed from single or multiple census mark-recapture data using a variety of methods. For single census data, the initial population size (N) is estimated from the number of marked animals from a first sample (M), number of captured animals in a second sample (n), and the number of recaptured marked animals in the second sample (m) using either the &amp;#8216;naive&amp;#8217; Petersen method or Chapman, Ricker, or Bailey modifications of the Petersen method. Single census data can also be separated by group (e.g., size class) to estimate the initial population size by class and for the overall population size. For multiple census data, the initial population size is estimated from the number of captured animals (n), number of recaptured marked animals (m), the number of marked animals that are marked and returned to the population (R), or the number of extant marked animals prior to the sample (M) on each of several samples using either the Schnabel (1938) or Schumacher-Eschmeyer (1943) method."><meta property="og:image" content="https://fishr-core-team.github.io/FSA/logo.png"><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
<![endif]--></head><body>
    <a href="#main" class="visually-hidden-focusable">Skip to contents</a>
    

    <nav class="navbar fixed-top navbar-dark navbar-expand-lg bg-primary"><div class="container">
    
    <a class="navbar-brand me-2" href="../index.html">FSA</a>

    <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">0.9.5</small>

    
    <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation">
      <span class="navbar-toggler-icon"></span>
    </button>

    <div id="navbar" class="collapse navbar-collapse ms-3">
      <ul class="navbar-nav me-auto"><li class="active nav-item">
  <a class="nav-link" href="../reference/index.html">Reference</a>
</li>
      </ul><form class="form-inline my-2 my-lg-0" role="search">
        <input type="search" class="form-control me-sm-2" aria-label="Toggle navigation" name="search-input" data-search-index="../search.json" id="search-input" placeholder="Search for" autocomplete="off"></form>

      <ul class="navbar-nav"><li class="nav-item">
  <a class="nav-link" href="../news/index.html">
    <span class="fa fa-newspaper"></span>
     
  </a>
</li>
<li class="nav-item">
  <a class="nav-link" href="../index.html">
    <span class="fa fa-home fa"></span>
     
  </a>
</li>
<li class="nav-item">
  <a class="external-link nav-link" href="https://github.com/fishR-Core-Team/FSA">
    <span class="fa fa-github fa"></span>
     
  </a>
</li>
      </ul></div>

    
  </div>
</nav><div class="container template-reference-topic">
<div class="row">
  <main id="main" class="col-md-9"><div class="page-header">
      <img src="../logo.png" class="logo" alt=""><h1>Estimate initial population size for single or multiple census mark-recapture data.</h1>
      <small class="dont-index">Source: <a href="https://github.com/fishR-Core-Team/FSA/blob/HEAD/R/mrClosed.R" class="external-link"><code>R/mrClosed.R</code></a></small>
      <div class="d-none name"><code>mrClosed.Rd</code></div>
    </div>

    <div class="ref-description section level2">
    <p>Estimates of the initial population size, along with associated confidence intervals, are constructed from single or multiple census mark-recapture data using a variety of methods. For single census data, the initial population size (N) is estimated from the number of marked animals from a first sample (M), number of captured animals in a second sample (n), and the number of recaptured marked animals in the second sample (m) using either the ‘naive’ Petersen method or Chapman, Ricker, or Bailey modifications of the Petersen method. Single census data can also be separated by group (e.g., size class) to estimate the initial population size by class and for the overall population size. For multiple census data, the initial population size is estimated from the number of captured animals (n), number of recaptured marked animals (m), the number of marked animals that are marked and returned to the population (R), or the number of extant marked animals prior to the sample (M) on each of several samples using either the Schnabel (1938) or Schumacher-Eschmeyer (1943) method.</p>
    </div>

    <div class="section level2">
    <h2 id="ref-usage">Usage<a class="anchor" aria-label="anchor" href="#ref-usage"></a></h2>
    <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="fu">mrClosed</span><span class="op">(</span></span>
<span>  M <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  n <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  m <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  R <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  method <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"Petersen"</span>, <span class="st">"Chapman"</span>, <span class="st">"Ricker"</span>, <span class="st">"Bailey"</span>, <span class="st">"Schnabel"</span>,</span>
<span>    <span class="st">"SchumacherEschmeyer"</span><span class="op">)</span>,</span>
<span>  labels <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  chapman.mod <span class="op">=</span> <span class="cn">TRUE</span></span>
<span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for mrClosed1</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span></span>
<span>  <span class="va">object</span>,</span>
<span>  digits <span class="op">=</span> <span class="fl">0</span>,</span>
<span>  incl.SE <span class="op">=</span> <span class="cn">FALSE</span>,</span>
<span>  incl.all <span class="op">=</span> <span class="cn">TRUE</span>,</span>
<span>  verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span>
<span>  <span class="va">...</span></span>
<span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for mrClosed1</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span></span>
<span>  <span class="va">object</span>,</span>
<span>  parm <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  level <span class="op">=</span> <span class="va">conf.level</span>,</span>
<span>  conf.level <span class="op">=</span> <span class="fl">0.95</span>,</span>
<span>  digits <span class="op">=</span> <span class="fl">0</span>,</span>
<span>  type <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"suggested"</span>, <span class="st">"binomial"</span>, <span class="st">"hypergeometric"</span>, <span class="st">"normal"</span>, <span class="st">"Poisson"</span><span class="op">)</span>,</span>
<span>  bin.type <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"wilson"</span>, <span class="st">"exact"</span>, <span class="st">"asymptotic"</span><span class="op">)</span>,</span>
<span>  poi.type <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"exact"</span>, <span class="st">"daly"</span>, <span class="st">"byar"</span>, <span class="st">"asymptotic"</span><span class="op">)</span>,</span>
<span>  incl.all <span class="op">=</span> <span class="cn">TRUE</span>,</span>
<span>  verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span>
<span>  <span class="va">...</span></span>
<span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for mrClosed2</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">object</span>, digits <span class="op">=</span> <span class="fl">0</span>, verbose <span class="op">=</span> <span class="cn">FALSE</span>, <span class="va">...</span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for mrClosed2</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span></span>
<span>  <span class="va">object</span>,</span>
<span>  parm <span class="op">=</span> <span class="cn">NULL</span>,</span>
<span>  level <span class="op">=</span> <span class="va">conf.level</span>,</span>
<span>  conf.level <span class="op">=</span> <span class="fl">0.95</span>,</span>
<span>  digits <span class="op">=</span> <span class="fl">0</span>,</span>
<span>  type <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"suggested"</span>, <span class="st">"normal"</span>, <span class="st">"Poisson"</span><span class="op">)</span>,</span>
<span>  poi.type <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"exact"</span>, <span class="st">"daly"</span>, <span class="st">"byar"</span>, <span class="st">"asymptotic"</span><span class="op">)</span>,</span>
<span>  verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span>
<span>  <span class="va">...</span></span>
<span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for mrClosed2</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span></span>
<span>  <span class="va">x</span>,</span>
<span>  pch <span class="op">=</span> <span class="fl">19</span>,</span>
<span>  col.pt <span class="op">=</span> <span class="st">"black"</span>,</span>
<span>  xlab <span class="op">=</span> <span class="st">"Marked in Population"</span>,</span>
<span>  ylab <span class="op">=</span> <span class="st">"Prop. Recaptures in Sample"</span>,</span>
<span>  loess <span class="op">=</span> <span class="cn">FALSE</span>,</span>
<span>  lty.loess <span class="op">=</span> <span class="fl">2</span>,</span>
<span>  lwd.loess <span class="op">=</span> <span class="fl">1</span>,</span>
<span>  col.loess <span class="op">=</span> <span class="st">"gray20"</span>,</span>
<span>  trans.loess <span class="op">=</span> <span class="fl">10</span>,</span>
<span>  span <span class="op">=</span> <span class="fl">0.9</span>,</span>
<span>  <span class="va">...</span></span>
<span><span class="op">)</span></span></code></pre></div>
    </div>

    <div class="section level2">
    <h2 id="arguments">Arguments<a class="anchor" aria-label="anchor" href="#arguments"></a></h2>
    <dl><dt>M</dt>
<dd><p>A numeric representing the number of marked fish from the first sample (single-census), an object from <code><a href="capHistSum.html">capHistSum()</a></code> (single- or multiple-census), or numeric vector of marked fish prior to ith samples (multiple-census).</p></dd>


<dt>n</dt>
<dd><p>A numeric representing the number of captured fish in the second sample (single-census) or numeric vector of captured fish in ith sample (multiple-census).</p></dd>


<dt>m</dt>
<dd><p>A numeric representing the number of recaptured (marked) fish in the second sample (single-census) or numeric vector of recaptured (marked) fish in ith sample (multiple-census).</p></dd>


<dt>R</dt>
<dd><p>A numeric vector representing the number of marked fish returned to the population (multiple-census). Note that several references use the number of “new” marks returned to the population rather than the “total” number of marks returned to the population that is used here.</p></dd>


<dt>method</dt>
<dd><p>A single string that identifies the type of calculation method to use in the main function.</p></dd>


<dt>labels</dt>
<dd><p>A character or character vector used to label the rows of the resulting output matrix when using a single census method separated by groups. Must be the same length as <code>M</code>, <code>n</code>, and <code>m</code>. Defaults to upper-case letters if no values are given.</p></dd>


<dt>chapman.mod</dt>
<dd><p>A logical that represents whether the Chapman modification should be used (<code>=TRUE</code>, default) or not (<code>=FALSE</code>) when performing the Schnabel multiple census method.</p></dd>


<dt>object, x</dt>
<dd><p>An <code>mrClosed1</code> or <code>mrClosed2</code> object.</p></dd>


<dt>digits</dt>
<dd><p>The number of decimal digits to round the population estimates to. If <code>incl.SE=TRUE</code> then SE will be rounded to one more decimal place then given in <code>digits</code>.</p></dd>


<dt>incl.SE</dt>
<dd><p>A logical that indicates whether the results should include the calculated SE value. See details.</p></dd>


<dt>incl.all</dt>
<dd><p>A logical that indicates whether an overall population estimate should be computed when using a single census method that has been separated into sub-groups. See details.</p></dd>


<dt>verbose</dt>
<dd><p>A logical that indicates whether a reminder of the inputted values and what type of method was used should be printed with the summary and confidence interval results.</p></dd>


<dt>...</dt>
<dd><p>Additional arguments for methods.</p></dd>


<dt>parm</dt>
<dd><p>Not used here (included in <code>confint</code> generic).</p></dd>


<dt>level</dt>
<dd><p>Same as <code>conf.level</code> but used for compatibility with <code>confint</code> generic.</p></dd>


<dt>conf.level</dt>
<dd><p>A numeric representing the level of confidence to use for confidence intervals.</p></dd>


<dt>type</dt>
<dd><p>A single string that identifies the distribution to use when constructing confidence intervals in <code>confint</code>. See details.</p></dd>


<dt>bin.type</dt>
<dd><p>A string that identifies the method used to construct binomial confidence intervals (default is <code>"wilson"</code>). This is only used if <code>type="binomial"</code> in <code>confint</code>. See details of <code><a href="binCI.html">binCI</a></code>.</p></dd>


<dt>poi.type</dt>
<dd><p>A string that identifies the method used to construct Poisson confidence intervals (default is <code>"exact"</code>). This is only used if <code>type="Poisson"</code> in <code>confint</code>. See details of <code><a href="poiCI.html">poiCI</a></code>.</p></dd>


<dt>pch</dt>
<dd><p>A numeric used to indicate the type of plotting character.</p></dd>


<dt>col.pt</dt>
<dd><p>a string used to indicate the color of the plotted points.</p></dd>


<dt>xlab</dt>
<dd><p>A label for the x-axis.</p></dd>


<dt>ylab</dt>
<dd><p>A label for the y-axis.</p></dd>


<dt>loess</dt>
<dd><p>A logical that indicates if a loess smoother line (and approximate 95% confidence band) is fit to and shown on plot.</p></dd>


<dt>lty.loess</dt>
<dd><p>A single numeric used to indicate the type of line used for the loess line.</p></dd>


<dt>lwd.loess</dt>
<dd><p>A single numeric used to indicate the line width of the loess line.</p></dd>


<dt>col.loess</dt>
<dd><p>A single string used to indicate the color of the loess line.</p></dd>


<dt>trans.loess</dt>
<dd><p>A single numeric that indicates how transparent the loess band should be (larger numbers are more transparent).</p></dd>


<dt>span</dt>
<dd><p>A single numeric that controls the degree of smoothing. Values closer to 1 are more smooth.</p></dd>

</dl></div>
    <div class="section level2">
    <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2>
    

<p>A list with the following items</p>
<p></p>
<ul><li><p>M The number of marked fish from the first sample that was provided.</p></li>
<li><p>n The number of captured fish in the second sample that was provided.</p></li>
<li><p>m The number of recaptured (marked) fish in the second sample that was provided.</p></li>
<li><p>M1 The adjusted (depending on <code>type</code>) number of marked fish from the first sample.</p></li>
<li><p>n1 The adjusted (depending on <code>type</code>) number of captured fish in the second sample.</p></li>
<li><p>m1 The adjusted (depending on <code>type</code>) number of recaptured (marked) fish in the second sample.</p></li>
<li><p>cf A correction factor for the population estimate that depends on <code>type</code>.</p></li>
<li><p>method The type of method used (provided by the user).</p></li>
<li><p>methodLbl A label for the type of method used.</p></li>
<li><p>N The estimated initial population size.</p></li>
<li><p>labels Labels for the rows of summary matrix.</p></li>
</ul></div>
    <div class="section level2">
    <h2 id="details">Details<a class="anchor" aria-label="anchor" href="#details"></a></h2>
    <p>For single census data, the following methods can be used:</p><ul><li><p><code>method="Petersen"</code>. The ‘naive’ Petersen as computed using equation 2.1 from Krebs (1989).</p></li>
<li><p><code>method="Chapman"</code>. The Chapman (1951) modification of the Petersen method as computed using equation 2.2 from Krebs (1989).</p></li>
<li><p><code>method="Ricker"</code>. The Ricker (1975) modification of the Petersen as computed using equation 3.7 from Ricker (1975). This is basically the same <code>method="Chapman"</code> except that Ricker (1975) did NOT subtract a 1 from the answer in the final step. Thus, the estimate from <code>method="Chapman"</code> will always be one less than the estimate from <code>method="Ricker"</code>.</p></li>
<li><p><code>method="Bailey"</code>. The Bailey (1951, 1952) modification of the Petersen as computed using equation 2.3 from Krebs (1989).</p></li>
</ul><p>If <code>M</code> contains an object from <code><a href="capHistSum.html">capHistSum</a></code> and one of Petersen, Chapman, Ricker, or Bailey methods has been selected with <code>method=</code> then <code>n=</code> and <code>m=</code> can be left missing or will be ignored and the needed data will be extracted from the <code>sum</code> portion of the <code>CapHist</code> class object. If the data were not summarized with <code><a href="capHistSum.html">capHistSum</a></code> then all of <code>M=</code>, <code>n=</code>, and <code>m=</code> must be supplied by the user.</p>
<p>The population estimate (as computed with the formulas noted in the table above) is extracted with <code>summary</code>. In addition, the standard error of the population estimate (SE) can be extracted by including <code>incl.SE=TRUE</code>. The SE is from equation 3.6 (p. 78) in Ricker (1975) for the Petersen method, from p. 60 (near bottom) of Seber (2002) for the Chapman method, from p. 61 (middle) of Seber (2002) (and as noted on p. 79 of Ricker (1975)) for the Bailey method, and from equation 3.8 (p. 78) in Ricker (1975) for the Ricker method.</p>
<p>Confidence intervals for the initial population size from the single census methods can be constructed using four different distributions as chosen with <code>type=</code> in <code>confint</code>. If <code>type="suggested"</code> then the type of confidence interval suggested by the rules on p. 18 in Krebs (1989) are used. The general methods for constructing confidence intervals for N are described below</p>
<ul><li><p><code>type="hypergeometric"</code>. Uses <code><a href="hyperCI.html">hyperCI</a></code>. This is experimental at this point.</p></li>
<li><p><code>type="binomial"</code>. Use <code><a href="binCI.html">binCI</a></code> to construct a confidence interval for m/n (Petersen method) or (m+1)/(n+1) (Chapman, Bailey, Ricker methods), divides M or (M+1) by the CI endpoints, and subtract 1 (for the Chapman method).</p></li>
<li><p><code>type="Poisson"</code>. Use <code><a href="poiCI.html">poiCI</a></code> to construct a confidence interval for m (Petersen method) or (m+1) (Chapman, Bailey, Ricker methods), substitute the CI endpoints into the appropriate equation for estimating N, and subtract 1 (for the Chapman method).</p></li>
<li><p><code>type="normal"</code>. Used equation 2.4 (p.20) from Krebs (2002) for the Petersen method. For the other methods, used N+/- Z(0.975)*SE, where the SE was computed as noted above.</p></li>
</ul><p>If <code>incl.all=TRUE</code> in <code>summary</code> and population estimates have been constructed for multiple sub-groups then an overall population estimate is included by summing the population estimates for the multiple sub-groups. If <code>incl.SE=TRUE</code>, then an overall SE is computed by taking the square root of the summed VARIANCES for the multiple sub-groups.</p>
<p>For multiple census data, the following methods can be declared for use with the <code>method=</code> argument:</p><ul><li><p><code>method="Schnabel"</code>. The Schnabel (1938) method as computed with equation 3.15 from Ricker (1975).</p></li>
<li><p><code>method="SchumacherEschmeyer"</code>. The Schumacher and Eschmeyer (1943) method as computed with equation 3.12 from Ricker (1975) eqn 3.12.</p></li>
</ul><p>If <code>M</code> contains an object from <code><a href="capHistSum.html">capHistSum</a></code> and the Schnabel or Schumacher-Eschmeyer methods has been chosen then <code>n</code>, <code>m</code> and <code>R</code> can be left missing or will be ignored. In this case, the needed data is extracted from the <code>sum</code> portion of the <code>CapHist</code> class object. Otherwise, the user must supply vectors of results in <code>n</code>, <code>m</code>, and <code>R</code> or <code>M</code>.</p>
<p>The population estimate for each method is extracted with <code>summary</code>. Standard errors for the population estimate can NOT be computed for the Schnabel or Schumacher-Eschmeyer methods (a warning will be produced if <code>incl.SE=TRUE</code> is used).</p>
<p>Confidence intervals for the initial population size using multiple census methods can be constructed using the normal or Poisson distributions for the Schnabel method or the normal distribution for the Schumacher-Eschmeyer method as chosen with <code>type=</code>. If <code>type="suggested"</code> then the type of confidence interval suggested by the rule on p. 32 of Krebs (1989) is used (for the Schnabel method). If <code>type="Poisson"</code> for the Schnabel method then a confidence interval for the sum of m is computed with <code><a href="poiCI.html">poiCI</a></code> and the end points are substituted into the Schnabel equation to produce a CI for the population size. If <code>type="normal"</code> for the Schnabel method then the standard error for the <em>inverse</em> of the population estimate is computed as the square root of equation 2.11 from Krebs (1989) or equation 3.16 from Ricker (1975). The standard error for the Schumacher-Eschmeyer method is for the <em>inverse</em> of the population estimate and is computed with equation 2.14 from Krebs (1989) [Note that the divisor in Krebs (1989) is different than the divisor in equation 3.12 in Ricker (1975), but is consistent with equation 4.17 in Seber (2002).]  The confidence interval for the <em>inverse</em> population estimate is constructed from the inverse population estimate plus/minus a t critical value times the standard error for the inverse population estimate. The t critical value uses the number of samples minus 1 for the Schnabel method and the number of samples minus 2 when for the Schumacher-Eschmeyer method according to p. 32 of Krebs (1989) (note that this is different than what Ricker (1975) does). Finally, the confidence interval for the population estimate is obtained by inverting the confidence interval for the inverse population estimate. Note that confidence intervals for the population size when <code>type="normal"</code> may contain negative values (for the upper value) when the population estimate is relatively large and the number of samples is small (say, three) because the intervals are originally constructed on the inverted population estimate and they use the t-distribution.</p>
<p>The <code>plot</code> can be used to identify assumption violations in the Schnabel and Schumacher-Eschmeyer methods (an error will be returned if used with any of the other methods). If the assumptions ARE met then the plot of the proportion of marked fish in a sample versus the cumulative number of marked fish should look linear. A loess line (with approximate 95% confidence bands) can be added to aid interpretation with <code>loess=TRUE</code>. Note, however, that adding the loess line may return a number of warning or produce a non-informative if the number of samples is small (&lt;8).</p>
    </div>
    <div class="section level2">
    <h2 id="testing">Testing<a class="anchor" aria-label="anchor" href="#testing"></a></h2>
    <p>The results from the single census methods have had the following checks. The population estimates for all methods match reputable sources. The SE for the Chapman and Bailey methods match the results from <code><a href="https://rdrr.io/pkg/fishmethods/man/mrN.single.html" class="external-link">mrN.single</a></code> in <span class="pkg">fishmethods</span>,  The CI for the Petersen, Chapman, and Bailey methods partially match (are within 1<!-- % when they do not match) the CIs from reputable sources. --></p>
<p>The results for the multiple census methods have had the following checks. The population estimates for both methods match reputable sources. The intermediate calculations for both methods match those in Krebs (1989). The confidence interval for the Schnabel method using the Poisson distribution does NOT match Krebs (1989). This appears to be a difference in the use <code><a href="poiCI.html">poiCI</a></code> here versus distributional tables in Krebs (i.e., the difference appears to be completely in the critical values from the Poisson distribution). The confidence interval for the Schnabel method using the normal or the Poisson distribution do NOT match Ricker (1975), but there is not enough information in Ricker to determine why (it is likely due to numerical differences on the inverse scale). The confidence interval for the Schumacher-Eschmeyer method do match Krebs (1989) but not Ricker (1975). The Ricker result may be due to different df as noted above.</p>
    </div>
    <div class="section level2">
    <h2 id="ifar-chapter">IFAR Chapter<a class="anchor" aria-label="anchor" href="#ifar-chapter"></a></h2>
    <p>9-Abundance from Capture-Recapture Data.</p>
    </div>
    <div class="section level2">
    <h2 id="references">References<a class="anchor" aria-label="anchor" href="#references"></a></h2>
    <p>Ogle, D.H. 2016. <a href="https://fishr-core-team.github.io/fishR/pages/books.html#introductory-fisheries-analyses-with-r" class="external-link">Introductory Fisheries Analyses with R</a>. Chapman &amp; Hall/CRC, Boca Raton, FL.</p>
<p>Krebs, C.J. 1989. Ecological Methodology. Addison-Welsey Educational Publishing.</p>
<p>Ricker, W.E. 1975. Computation and interpretation of biological statistics of fish populations. Technical Report Bulletin 191, Bulletin of the Fisheries Research Board of Canada. [Was (is?) from http://www.dfo-mpo.gc.ca/Library/1485.pdf.]</p>
<p>Seber, G.A.F. 2002. The Estimation of Animal Abundance and Related Parameters. Edward Arnold, second edition.</p>
<p>Schnabel, Z.E. 1938. The estimation of the total fish population of a lake. American Mathematician Monthly, 45:348-352.</p>
<p>Schumacher, F.X. and R.W. Eschmeyer. 1943. The estimation of fish populations in lakes and ponds. Journal of the Tennessee Academy of Sciences, 18:228-249.</p>
    </div>
    <div class="section level2">
    <h2 id="see-also">See also<a class="anchor" aria-label="anchor" href="#see-also"></a></h2>
    <div class="dont-index"><p>See <code><a href="capHistSum.html">capHistSum</a></code> for generating input data from capture histories. See <code><a href="poiCI.html">poiCI</a></code>, <code><a href="binCI.html">binCI</a></code>, and <code><a href="hyperCI.html">hyperCI</a></code> for specifics on functions used in confidence interval construction. See <code><a href="mrOpen.html">mrOpen</a></code> for handling mark-recapture data in an open population. See <code><a href="https://fishr-core-team.github.io/FSAdata/reference/SunfishIN.html">SunfishIN</a></code> in <span class="pkg">FSAdata</span> for an example to test matching of results with Ricker (1975)'  See <code><a href="https://rdrr.io/pkg/fishmethods/man/mrN.single.html" class="external-link">mrN.single</a></code> and <code><a href="https://rdrr.io/pkg/fishmethods/man/schnabel.html" class="external-link">schnabel</a></code> in <span class="pkg">fishmethods</span> for similar functionality.</p></div>
    </div>
    <div class="section level2">
    <h2 id="author">Author<a class="anchor" aria-label="anchor" href="#author"></a></h2>
    <p>Derek H. Ogle, <a href="mailto:DerekOgle51@gmail.com">DerekOgle51@gmail.com</a></p>
    </div>

    <div class="section level2">
    <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-examples"></a></h2>
    <div class="sourceCode"><pre class="sourceCode r"><code><span class="r-in"><span><span class="co">### Single census with no sub-groups</span></span></span>
<span class="r-in"><span><span class="co">## Petersen estimate -- the default</span></span></span>
<span class="r-in"><span><span class="va">mr1</span> <span class="op">&lt;-</span> <span class="fu">mrClosed</span><span class="op">(</span><span class="fl">346</span>,<span class="fl">184</span>,<span class="fl">49</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr1</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1299</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr1</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used the 'naive' Petersen method with M=346, n=184, and m=49.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1299</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr1</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N  SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1299 159</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr1</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span>,digits<span class="op">=</span><span class="fl">1</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>           N     SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1299.3 158.99</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr1</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]    1034    1666</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr1</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> The binomial (wilson method) distribution was used.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]    1034    1666</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr1</span>,type<span class="op">=</span><span class="st">"hypergeometric"</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]    1049    1670</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">## Chapman modification of the Petersen estimate</span></span></span>
<span class="r-in"><span><span class="va">mr2</span> <span class="op">&lt;-</span> <span class="fu">mrClosed</span><span class="op">(</span><span class="fl">346</span>,<span class="fl">184</span>,<span class="fl">49</span>,method<span class="op">=</span><span class="st">"Chapman"</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr2</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N    SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1283 142.1</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr2</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used Chapman's modification of the Petersen method with M=346, n=184, and m=49.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N    SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1283 142.1</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">### Single census, using capHistSum() results</span></span></span>
<span class="r-in"><span><span class="co">## data in capture history format</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/utils/str.html" class="external-link">str</a></span><span class="op">(</span><span class="va">BluegillJL</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 'data.frame':	277 obs. of  2 variables:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>  $ first : int  1 0 1 0 1 1 1 1 1 1 ...</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>  $ second: int  0 1 0 1 0 0 0 0 0 0 ...</span>
<span class="r-in"><span><span class="va">ch1</span> <span class="op">&lt;-</span> <span class="fu"><a href="capHistSum.html">capHistSum</a></span><span class="op">(</span><span class="va">BluegillJL</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="va">mr3</span> <span class="op">&lt;-</span> <span class="fu">mrClosed</span><span class="op">(</span><span class="va">ch1</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr3</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used the 'naive' Petersen method with M=196, n=90, and m=9.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>         N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 1960</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr3</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> The Poisson (exact method) distribution was used.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]    1032    4286</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">### Single census with sub-groups</span></span></span>
<span class="r-in"><span><span class="va">marked</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">93</span>,<span class="fl">35</span>,<span class="fl">72</span>,<span class="fl">16</span>,<span class="fl">46</span>,<span class="fl">20</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="va">captured</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">103</span>,<span class="fl">30</span>,<span class="fl">73</span>,<span class="fl">17</span>,<span class="fl">39</span>,<span class="fl">18</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="va">recaps</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">20</span>,<span class="fl">23</span>,<span class="fl">52</span>,<span class="fl">15</span>,<span class="fl">35</span>,<span class="fl">16</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="va">lbls</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"YOY"</span>,<span class="st">"Juvenile"</span>,<span class="st">"Stock"</span>,<span class="st">"Quality"</span>,<span class="st">"Preferred"</span>,<span class="st">"Memorable"</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="va">mr4</span> <span class="op">&lt;-</span> <span class="fu">mrClosed</span><span class="op">(</span><span class="va">marked</span>,<span class="va">captured</span>,<span class="va">recaps</span>,method<span class="op">=</span><span class="st">"Ricker"</span>,labels<span class="op">=</span><span class="va">lbls</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr4</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>             N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY       466</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile   46</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock     102</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality    19</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred  52</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable  23</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> All       708</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr4</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>             N   SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY       466 88.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile   46  4.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock     102  7.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality    19  1.5</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred  52  2.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable  23  1.8</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> All       708 89.2</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr4</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used Ricker's modification of the Petersen method with observed inputs (by group) of:</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   YOY: M=93, n=103, and m=20</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Juvenile: M=35, n=30, and m=23</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Stock: M=72, n=73, and m=52</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Quality: M=16, n=17, and m=15</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Preferred: M=46, n=39, and m=35</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Memorable: M=20, n=18, and m=16</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>             N   SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY       466 88.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile   46  4.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock     102  7.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality    19  1.5</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred  52  2.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable  23  1.8</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> All       708 89.2</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr4</span>,incl.SE<span class="op">=</span><span class="cn">TRUE</span>,incl.all<span class="op">=</span><span class="cn">FALSE</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used Ricker's modification of the Petersen method with observed inputs (by group) of:</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   YOY: M=93, n=103, and m=20</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Juvenile: M=35, n=30, and m=23</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Stock: M=72, n=73, and m=52</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Quality: M=16, n=17, and m=15</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Preferred: M=46, n=39, and m=35</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span>   Memorable: M=20, n=18, and m=16</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>             N   SE</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY       466 88.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile   46  4.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock     102  7.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality    19  1.5</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred  52  2.7</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable  23  1.8</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr4</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>           95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY           325     691</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile       41      60</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock          91     121</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality        18      25</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred      49      61</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable      22      31</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> All           533     883</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr4</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> YOY - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Juvenile - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Stock - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Quality - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Preferred - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Memorable - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> All - The normal distribution was used.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>           95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY           325     691</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile       41      60</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock          91     121</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality        18      25</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred      49      61</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable      22      31</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> All           533     883</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr4</span>,incl.all<span class="op">=</span><span class="cn">FALSE</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> YOY - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Juvenile - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Stock - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Quality - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Preferred - The binomial (wilson method) distribution was used.</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Memorable - The binomial (wilson method) distribution was used.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>           95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> YOY           325     691</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Juvenile       41      60</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Stock          91     121</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Quality        18      25</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Preferred      49      61</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Memorable      22      31</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">### Multiple Census</span></span></span>
<span class="r-in"><span><span class="co">## Data in summarized form ... Schnabel method</span></span></span>
<span class="r-in"><span><span class="va">mr5</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/with.html" class="external-link">with</a></span><span class="op">(</span><span class="va">PikeNY</span>,<span class="fu">mrClosed</span><span class="op">(</span>n<span class="op">=</span><span class="va">n</span>,m<span class="op">=</span><span class="va">m</span>,R<span class="op">=</span><span class="va">R</span>,method<span class="op">=</span><span class="st">"Schnabel"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mr5</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mr5</span>,loess<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="mrClosed-1.png" alt="" width="700" height="433"></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr5</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>       N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 87</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr5</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Used the Schnabel method with Chapman modification.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>       N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 87</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr5</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]      71     113</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr5</span>,verbose<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> The normal distribution was used.</span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]      71     113</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">## Schumacher-Eschmeyer method</span></span></span>
<span class="r-in"><span><span class="va">mr6</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/with.html" class="external-link">with</a></span><span class="op">(</span><span class="va">PikeNY</span>,<span class="fu">mrClosed</span><span class="op">(</span>n<span class="op">=</span><span class="va">n</span>,m<span class="op">=</span><span class="va">m</span>,R<span class="op">=</span><span class="va">R</span>,method<span class="op">=</span><span class="st">"Schumacher"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr6</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>       N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 85</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr6</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]      76      96</span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">### Capture history data summarized by capHistSum()</span></span></span>
<span class="r-in"><span><span class="co"># ignore first column of ID numbers</span></span></span>
<span class="r-in"><span><span class="va">ch2</span> <span class="op">&lt;-</span> <span class="fu"><a href="capHistSum.html">capHistSum</a></span><span class="op">(</span><span class="va">PikeNYPartial1</span>,cols2ignore<span class="op">=</span><span class="st">"id"</span><span class="op">)</span></span></span>
<span class="r-in"><span></span></span>
<span class="r-in"><span><span class="co">## Schnabel method</span></span></span>
<span class="r-in"><span><span class="va">mr7</span> <span class="op">&lt;-</span> <span class="fu">mrClosed</span><span class="op">(</span><span class="va">ch2</span>,method<span class="op">=</span><span class="st">"Schnabel"</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mr7</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="mrClosed-2.png" alt="" width="700" height="433"></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">mr7</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>        N</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,] 128</span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/confint.html" class="external-link">confint</a></span><span class="op">(</span><span class="va">mr7</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span>      95% LCI 95% UCI</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> [1,]      75     238</span>
<span class="r-in"><span></span></span>
</code></pre></div>
    </div>
  </main><aside class="col-md-3"><nav id="toc"><h2>On this page</h2>
    </nav></aside></div>


    <footer><div class="pkgdown-footer-left">
  <p></p><p>Developed by <a href="http://derekogle.com" class="external-link">Derek H. Ogle</a>, Jason C. Doll, A. Powell Wheeler, Alexis Dinno.</p>
</div>

<div class="pkgdown-footer-right">
  <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.0.7.</p>
</div>

    </footer></div>

  

  

  </body></html>

