diff --git a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-110-1.pdf b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-110-1.pdf index c38208f..6cdf3f3 100644 Binary files a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-110-1.pdf and b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-110-1.pdf differ diff --git a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-111-1.pdf b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-111-1.pdf index 480ce82..848963a 100644 Binary files a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-111-1.pdf and b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-111-1.pdf differ diff --git a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-131-1.pdf b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-131-1.pdf index ecacfa3..b367003 100644 Binary files a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-131-1.pdf and b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-131-1.pdf differ diff --git a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-132-1.pdf b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-132-1.pdf index 8beb34f..e827cc1 100644 Binary files a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-132-1.pdf and b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-132-1.pdf differ diff --git a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-169-1.pdf b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-169-1.pdf index d002ee8..62448cd 100644 Binary files a/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-169-1.pdf and b/_bookdown_files/banana-book_files/figure-html/unnamed-chunk-169-1.pdf differ diff --git a/dispersal.Rmd b/dispersal.Rmd index 67e7f56..e970e09 100644 --- a/dispersal.Rmd +++ b/dispersal.Rmd @@ -4,13 +4,11 @@ WORK IN PROGRESS. **I'm missing a section on gof w/ R2ucare and posterior predic ## Introduction -Blabla. +In this fifth chapter, you will learn about the Arnason-Schwarz model that allows estimating transitions between sites and states based on capture-recapture data. You will also see how to deal with uncertainty in the assignment of states to individuals. ## The Arnason-Schwarz (AS) model -+ Arnason 1972, 1973. Schwarz 1993. - -+ Nichols et al. 1992 (microtus and mass as state). Nichols et al. 1994 (microtus and reproductive costs). +In Chapter \@ref(survival), we got acquainted with Cormack-Jolly-Seber (CJS) model which accommodates transitions between the states alive and dead while accounting for imperfect detection. It is often the case that besides being alive, more detailed information is collected on the state of animals when they are detected. For example, if the study area is split into several discrete sites, you may record where an animal is detected, the state being now alive in this particular site. The Arnason-Schwarz (AS) model can be viewed as an extension to the CJS model in which we estimate movements between sites on top of survival. The AS model is named after the two statisticians -- Neil Arnason and Carl Schwarz. ## Multisite capture-recapture data @@ -1155,11 +1153,7 @@ Breeders are difficult to assigned to the correct state. Non-breeders are relati ## Suggested reading -+ Lebreton, J.-D., J. D. Nichols, R. J. Barker, R. Pradel and J. A. Spendelow (2009). [Modeling Individual Animal Histories with Multistate Capture–Recapture Models](https://multievent.sciencesconf.org/conference/multievent/pages/Lebretonetal2009AER.pdf). Advances in Ecological Research, 41:87-173. - -+ Seminal paper by Pradel (2005) [Multievent: An Extension of Multistate Capture–Recapture Models to Uncertain States](https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1541-0420.2005.00318.x). Biometrics, 61: 442-447. - -+ Dupuis (1995) had a similar idea for the Arnason-Schwarz model: Dupuis, J. (1995) [Bayesian estimation of movement and survival probabilities from capture-recapture data](https://academic.oup.com/biomet/article-abstract/82/4/761/252161). Biometrika. Vol. 82, pp 761-772. ++ Sites: Arnason 1972, 1973. Schwarz 1993. States: Nichols et al. 1992 (microtus and mass as state). Nichols et al. 1994 (microtus and reproductive costs). Lebreton, J.-D., J. D. Nichols, R. J. Barker, R. Pradel and J. A. Spendelow (2009). [Modeling Individual Animal Histories with Multistate Capture–Recapture Models](https://multievent.sciencesconf.org/conference/multievent/pages/Lebretonetal2009AER.pdf). Advances in Ecological Research, 41:87-173. -+ See also for a review Gimenez et al. (2012) [Estimating demographic parameters using hidden process dynamic models](https://oliviergimenez.github.io/pubs/Gimenezetal2012TPB.pdf). Theoretical Population Biology 82: 307-316. ++ Uncertainty: Seminal paper by Pradel (2005) [Multievent: An Extension of Multistate Capture–Recapture Models to Uncertain States](https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1541-0420.2005.00318.x). Biometrics, 61: 442-447. Dupuis (1995) had a similar idea for the Arnason-Schwarz model: Dupuis, J. (1995) [Bayesian estimation of movement and survival probabilities from capture-recapture data](https://academic.oup.com/biomet/article-abstract/82/4/761/252161). Biometrika. Vol. 82, pp 761-772. See also for a review Gimenez et al. (2012) [Estimating demographic parameters using hidden process dynamic models](https://oliviergimenez.github.io/pubs/Gimenezetal2012TPB.pdf). Theoretical Population Biology 82: 307-316. diff --git a/docs/404.html b/docs/404.html index 9153378..c13bb03 100644 --- a/docs/404.html +++ b/docs/404.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • diff --git a/docs/about-the-author.html b/docs/about-the-author.html index e9d7c18..141f938 100644 --- a/docs/about-the-author.html +++ b/docs/about-the-author.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • diff --git a/docs/covariateschapter.html b/docs/covariateschapter.html index 81c7a66..0c271e8 100644 --- a/docs/covariateschapter.html +++ b/docs/covariateschapter.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • diff --git a/docs/crashcourse.html b/docs/crashcourse.html index b1d44ba..0e853cd 100644 --- a/docs/crashcourse.html +++ b/docs/crashcourse.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • @@ -260,7 +260,7 @@

     sample_from_posterior <- rbeta(1000, 20, 39) # draw 1000 values from posterior survival beta(20,39)
     mean(sample_from_posterior) # compute mean with Monte Carlo integration
    -## [1] 0.3365
    +## [1] 0.3381

    You may check that the mean we have just calculated matches closely the expectation of a beta distribution10:

     20/(20+39) # expectation of beta(20,39)
    @@ -269,7 +269,7 @@ 

     quantile(sample_from_posterior, probs = c(2.5/100, 97.5/100))
     ##   2.5%  97.5% 
    -## 0.2266 0.4695
    +## 0.2231 0.4687

    diff --git a/docs/dispersal.html b/docs/dispersal.html index fd77578..dadf121 100644 --- a/docs/dispersal.html +++ b/docs/dispersal.html @@ -6,15 +6,15 @@ Chapter 5 Sites and states | Bayesian Analysis of Capture-Recapture Data with Hidden Markov Models - + - + - + @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • @@ -106,16 +106,13 @@

    5.1 Introduction

    -

    Blabla.

    +

    In this fifth chapter, you will learn about the Arnason-Schwarz model that allows estimating transitions between sites and states based on capture-recapture data. You will also see how to deal with uncertainty in the assignment of states to individuals.

    5.2 The Arnason-Schwarz (AS) model

    - +

    In Chapter 4, we got acquainted with Cormack-Jolly-Seber (CJS) model which accommodates transitions between the states alive and dead while accounting for imperfect detection. It is often the case that besides being alive, more detailed information is collected on the state of animals when they are detected. For example, if the study area is split into several discrete sites, you may record where an animal is detected, the state being now alive in this particular site. The Arnason-Schwarz (AS) model can be viewed as an extension to the CJS model in which we estimate movements between sites on top of survival. The AS model is named after the two statisticians – Neil Arnason and Carl Schwarz.

    @@ -35055,10 +35052,8 @@

    5.11 Suggested reading

    @@ -35066,7 +35061,7 @@

    - +
    -
    ## Time difference of 19.01 secs
    +
    ## Time difference of 19.63 secs

    We can have a look to numerical summaries:

     MCMCsummary(mcmc.output, round = 2)
    @@ -5627,7 +5627,7 @@ 

    ## |-------------------------------------------------------| end_time <- Sys.time() end_time - start_time -## Time difference of 17.84 secs

    +## Time difference of 17.86 secs

    The numerical summaries are similar to those we obtained with the complete likelihood, and effective samples sizes are larger denoting better mixing:

     MCMCsummary(mcmc.output, round = 2)
    @@ -5763,7 +5763,7 @@ 

    ## |-------------------------------------------------------| end_time <- Sys.time() end_time - start_time -## Time difference of 17.63 secs +## Time difference of 17.71 secs MCMCsummary(mcmc.output, round = 2) ## mean sd 2.5% 50% 97.5% Rhat n.eff ## p 0.61 0.06 0.49 0.61 0.72 1 1453 diff --git a/docs/index.html b/docs/index.html index 016b818..66e7816 100644 --- a/docs/index.html +++ b/docs/index.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • diff --git a/docs/introduction-4.html b/docs/introduction-4.html index 04f8e2f..866dd75 100644 --- a/docs/introduction-4.html +++ b/docs/introduction-4.html @@ -78,7 +78,7 @@

  • 3 Hidden Markov models
  • II. Transitions
  • Introduction
  • -
  • 4 Live and dead
  • +
  • 4 Alive and dead
  • 5 Sites and states
  • III. Case studies
  • Introduction
  • @@ -104,7 +104,7 @@

    Introduction

    -4 Live and dead +4 Alive and dead

    4.1 Introduction

    -

    In this fourth chapter, you will learn about the Cormack-Jolly-Seber model that allows estimating survival based on capture-recapture data. You will also see how to deal with covariates to try and explain temporal and/or individual variation in survival. This chapter will also be the opportunity to illustrate how to incorporate prior information to improve model inference.

    +

    In this fourth chapter, you will learn about the Cormack-Jolly-Seber model that allows estimating survival based on capture-recapture data. You will also see how to deal with covariates to try and explain temporal and/or individual variation in survival. This chapter will also be the opportunity to introduce tools to compare models and assess their quality of fit to data.

    4.2 The Cormack-Jolly-Seber (CJS) model

    -

    In chapter 3, we introduced a capture-recapture model with constant survival and detection probabilities which we formulated as a HMM and fitted to data in NIMBLE. Historically, however, it was a slightly more complicated model that was first proposed – the so-called Cormack-Jolly-Seber (CJS) model – in which survival and recapture probabilities are time-varying. This feature of the CJS model is useful to account for variation due to environmental conditions in survival or to sampling effort in detection. Schematically the CJS model can be represented this way:

    +

    In Chapter 3, we introduced a capture-recapture model with constant survival and detection probabilities which we formulated as a HMM and fitted to data in NIMBLE. Historically, however, it was a slightly more complicated model that was first proposed – the so-called Cormack-Jolly-Seber (CJS) model – in which survival and recapture probabilities are time-varying. This feature of the CJS model is useful to account for variation due to environmental conditions in survival or to sampling effort in detection. Schematically the CJS model can be represented this way:

    Note that the states (in gray) and the observations (in white) do not change. We still have \(z = 1\) for alive, \(z = 2\) for dead, \(y = 1\) for non-detected, and \(y = 2\) for detected.

    Parameters are now indexed by time. The survival probability is defined as the probability of staying alive (or “ah, ha, ha, ha, stayin’ alive” like the Bee Gees would say) over the interval between \(t\) and \(t+1\), that is \(\phi_t = \Pr(z_{t+1} = 1 | z_t = 1)\). The detection probability is defined as the probability of being observed at \(t\) given you’re alive at \(t\), that is \(p_t = \Pr(y_{t} = 1 | z_t = 1)\). It is important to bear in mind that survival operates over an interval while detection occurs at a specific time (see Section 4.8).

    -

    The CJS model is named after three statisticians who each published independently a paper introducing more or less the same approach, a year apart ! In fact, Richard Cormack and George Jolly were working in the same corridor in Scotland back in the 1960’s. They would meet every day at coffee and would play some game together, but never mention work and were not aware of each other’s work.

    +

    The CJS model is named after the three statisticians – Richard Cormack, George Jolly and George Seber – who each published independently a paper introducing more or less the same approach, a year apart ! In fact, Richard Cormack and George Jolly were working in the same corridor in Scotland back in the 1960’s. They would meet every day at coffee and would play some game together, but never mention work and were not aware of each other’s work.

    @@ -9896,7 +9896,7 @@