sanivult
Evaluating the impacts of the outbreak of a Bovine Spongiform Encephalopathy epidemic on the demographic and population dynamics of one of the world's largest colonies of the Eurasian Griffon vulture (Gyps fulvus).
Science Score: 67.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
✓DOI references
Found 23 DOI reference(s) in README -
✓Academic publication links
Links to: scholar.google, wiley.com, zenodo.org -
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (5.0%) to scientific vocabulary
Keywords
Repository
Evaluating the impacts of the outbreak of a Bovine Spongiform Encephalopathy epidemic on the demographic and population dynamics of one of the world's largest colonies of the Eurasian Griffon vulture (Gyps fulvus).
Basic Info
Statistics
- Stars: 1
- Watchers: 1
- Forks: 1
- Open Issues: 0
- Releases: 3
Topics
Metadata Files
README.html
README
The SaniVult project
Table of Contents
About
This is the GitHub hosting of the project SaniVult. The paper associated to the project is published in the journal Ecological Applications. See the
CITATIONfile for a BibTex entry to the article. This folder contains the files needed to reproduce all the results of the project, and compile the manuscript of the associated paper.Abstract
Scavenging is a key ecological process controlling energy flow in ecosystems and providing valuable ecosystem services worldwide. As long-lived species, the demographic dynamics of vultures can be disrupted by spatio-temporal fluctuations in food availability, with dramatic impacts on their population viability and the ecosystem services provided. In Europe, the outbreak of Bovine Spongiform Encephalopathy (BSE) in 2001 prompted a restrictive sanitary legislation banning the presence of livestock carcasses in the wild at a continental scale. In long-lived vertebrate species the buffering hypothesis predicts that the demographic traits with the largest contribution to population growth rate should be less temporally variable. The BSE outbreak provides a unique opportunity to test for the impact of demographic buffering in a keystone scavenger suffering abrupt but transient food shortages. We study the 42-year dynamics (1978-2020) of one of the worlds largest breeding colonies of Eurasian griffon vultures (Gyps fulvus). We fitted an inverse Bayesian state-space model with density-dependent demographic rates to the time-series of stage-structured abundances to investigate shifts in vital rates and population dynamics before, during and after the implementation of a restrictive sanitary regulation. Prior to the BSE outbreak the dynamics was mainly driven by adult survival: 83% of temporal variance in abundance was explained by variability in this rate. Moreover, during this period the regulation of population size operated through density-dependent fecundity and sub-adult survival. However, after the onset of the European ban, a one-month delay in average laying date, a drop in fecundity and a reduction in the number of fledglings induced a transient increase in the impact of fledgling and sub-adult recruitment on dynamics. Although adult survival rate remained constantly high, as predicted by the buffering hypothesis, its relative impact on the temporal variance in abundance dropped to 71% during the sanitary legislation and to 54% after the ban was lifted. A significant increase in the relative impact of environmental stochasticity on dynamics was modeled after the BSE outbreak. These results provide empirical evidence on how abrupt environmental deterioration may induce dramatic demographic and dynamic changes in the populations of keystone scavengers, with far-reaching impacts on ecosystem functioning worldwide.
Built With
This is a workflowr project bootstraped by a suite of open-source tools.
A suite of R packages were used in this project. I am grateful to all the people involved in the development of these open-source packages. Run the following R command from within the project for producing a reference list of the packages used:
grateful::cite_packages(out.format = "rmd", out.dir = file.path(getwd(), "analysis"))A list of these packages is placed at the end of this document.
Reproducible workflow
This section shows how to reproduce the results of the accompanying paper. The folder
./codehas the following structure:. ./code ./code/S4D3M_JAGS_Fitting.R ./code/S4D3M_JAGS_model.jags ./code/utilities.RIn this folder, the file
./code/utilities.Rcontains all the functions and utilities necessary to conduct the analyses. The file./code/S4D3M_JAGS_model.jagscontains the state-space stage-structured demographic density-dependent model (S4D3M) developed in the accompanying paper written in the JAGS language.The data folder has the following structure:
./data ./data/Breeding_output.csv ./data/BSE_cases.csv ./data/data.csvThe manuscript folder has the following structure:
ms appendix.pdf appendix.tex arxiv.sty biblio.bib DataS1.zip figs Fig1.pdf Fig2.pdf Fig3.pdf Fig4.pdf FigS1.pdf FigS3.pdf FigS4.pdf main_text.pdf main_text.tex MetadataS1.docx MetadataS1.pdfPrerequisites
Prior to reproducing the results, make sure to have installed all the necessary software. In particular, you need JAGS and R. The R libraries needed to reproduce the results (see below) will be automatically installed by the package pacman.
Workflow
Note that the S4D3M is fitted through Bayesian MCMC methods using Gibbs sampling, and runs in JAGS: even though JAGS is written in the C++ language, the code can take several hours to run depending on the architecture used. Note that there are relatively easy ways of parallelizing this code.
You can reproduce the results of the accompanying paper with three methods:
The first, easiest way to reproduce all the analyses in the project is to use the
Makefile. With simple GNU Make syntax, you can reproduce all the project, from statistical analyses to manuscript production. For example, in GNU/Linux based systems, you can point with the command shell to the project folder and run the following command:make allThis command will first conduct all the statistical analyses in the project, and produce all the figures. It then will assemble and compile the manuscript and associated supplementary materials with the necessary figures. Finally, it will open the files. Alternatively, note that you can run this command within RStudio from the Terminal tab.
From within
R, simplysourcethe file./code/S4D3M_JAGS_Fitting.R. This will perform all the analyses of the paper in the required order.The final method is to open the R Markdown file
./analysis/index.Rmdto interactively execute the workflow.
Contact
Pablo Almaraz - @palmaraz_Eco - palmaraz@gmail.com
Personal webpage: https://palmaraz.github.io/
Project Link: https://github.com/palmaraz/SaniVult
R packages used in this project
- base (R Core Team 2021)
- workflowr (Blischak, Carbonetto, and Stephens 2019)
- rmarkdown (Xie, Dervieux, and Riederer 2020)
- checkpoint (Ooi, de Vries, and Microsoft 2021)
- coda (Plummer et al.2006)
- data.table (Dowle and Srinivasan 2021)
- ggmcmc (Fernndez-i-Marn 2016)
- ggsci (Xiao 2018)
- grateful (Rodrguez-Snchez and Hutchins 2020)
- mvtnorm (Genz and Bretz 2009)
- pacman (Rinker and Kurkiewicz 2018)
- patchwork (Pedersen 2020)
- runjags (Denwood 2016)
- tidyverse (Wickham et al.2019)
- truncnorm (Mersmann et al.2018)
- viridis (Garnier et al.2021)
- xtable (Dahl et al.2019)
- dplyr (Wickham et al.2021)
- ggpubr (Kassambara 2020)
- plyr (Wickham 2011)
- readr (Wickham and Hester 2021)
- reshape2 (Wickham 2007)
- tibble (Mller and Wickham 2021)
References
Blischak, John D, Peter Carbonetto, and Matthew Stephens. 2019. Creating and Sharing Reproducible Research Code the Workflowr Way F1000Research 8 (1749). https://doi.org/10.12688/f1000research.20843.1.
Dahl, David B., David Scott, Charles Roosen, Arni Magnusson, and Jonathan Swinton. 2019. Xtable: Export Tables to LaTeX or HTML. https://CRAN.R-project.org/package=xtable.
Denwood, Matthew J. 2016. runjags: An R Package Providing Interface Utilities, Model Templates, Parallel Computing Methods and Additional Distributions for MCMC Models in JAGS. Journal of Statistical Software 71 (9): 125. https://doi.org/10.18637/jss.v071.i09.
Dowle, Matt, and Arun Srinivasan. 2021. Data.table: Extension of Data.frame. https://CRAN.R-project.org/package=data.table.
Fernndez-i-Marn, Xavier. 2016. ggmcmc: Analysis of MCMC Samples and Bayesian Inference. Journal of Statistical Software 70 (9): 120. https://doi.org/10.18637/jss.v070.i09.
Garnier, Simon, Ross, Noam, Rudis, Robert, Camargo, et al.2021. viridis - Colorblind-Friendly Color Maps for r. https://doi.org/10.5281/zenodo.4679424.
Genz, Alan, and Frank Bretz. 2009. Computation of Multivariate Normal and t Probabilities. Lecture Notes in Statistics. Heidelberg: Springer-Verlag.
Kassambara, Alboukadel. 2020. Ggpubr: Ggplot2 Based Publication Ready Plots. https://CRAN.R-project.org/package=ggpubr.
Mersmann, Olaf, Heike Trautmann, Detlef Steuer, and Bjrn Bornkamp. Truncnorm: Truncated Normal Distribution. https://CRAN.R-project.org/package=truncnorm.
Mller, Kirill, and Hadley Wickham. 2021. Tibble: Simple Data Frames. https://CRAN.R-project.org/package=tibble.
Ooi, Hong, Andrie de Vries, and Microsoft. 2021. Checkpoint: Install Packages from Snapshots on the Checkpoint Server for Reproducibility. https://CRAN.R-project.org/package=checkpoint.
Pedersen, Thomas Lin. 2020. Patchwork: The Composer of Plots. https://CRAN.R-project.org/package=patchwork.
Plummer, Martyn, Nicky Best, Kate Cowles, and Karen Vines. 2006. CODA: Convergence Diagnosis and Output Analysis for MCMC. R News 6 (1): 711. https://journal.r-project.org/archive/.
R Core Team. 2021. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
Rinker, Tyler W., and Dason Kurkiewicz. 2018. pacman: Package Management for R. Buffalo, New York. http://github.com/trinker/pacman.
Rodrguez-Snchez, Francisco, and Shaurita D. Hutchins. 2020. Grateful: Facilitate Citation of r Packages. https://github.com/Pakillo/grateful.
Wickham, Hadley. 2007. Reshaping Data with the reshape Package. Journal of Statistical Software 21 (12): 120. http://www.jstatsoft.org/v21/i12/.
. 2011. The Split-Apply-Combine Strategy for Data Analysis. Journal of Statistical Software 40 (1): 129. http://www.jstatsoft.org/v40/i01/.
Wickham, Hadley, Mara Averick, Jennifer Bryan, Winston Chang, Lucy DAgostino McGowan, Romain Franois, Garrett Grolemund, et al.2019. Welcome to the tidyverse. Journal of Open Source Software 4 (43): 1686. https://doi.org/10.21105/joss.01686.
Wickham, Hadley, Romain Franois, Lionel Henry, and Kirill Mller. 2021. Dplyr: A Grammar of Data Manipulation. https://CRAN.R-project.org/package=dplyr.
Wickham, Hadley, and Jim Hester. 2021. Readr: Read Rectangular Text Data. https://CRAN.R-project.org/package=readr.
Xiao, Nan. 2018. Ggsci: Scientific Journal and Sci-Fi Themed Color Palettes for Ggplot2. https://CRAN.R-project.org/package=ggsci.
Xie, Yihui, Christophe Dervieux, and Emily Riederer. 2020. R Markdown Cookbook. Boca Raton, Florida: Chapman; Hall/CRC. https://bookdown.org/yihui/rmarkdown-cookbook.
Owner
- Name: Pablo Almaraz
- Login: palmaraz
- Kind: user
- Location: Spain
- Company: CSIC, US
- Website: https://palmaraz.github.io/
- Twitter: palmaraz_Eco
- Repositories: 1
- Profile: https://github.com/palmaraz
A quantitative ecologist working in the mathematical and statistical analysis of biological data from real world applications.
Citation (CITATION)
@article{Almaraz2021,
author = {Almaraz, Pablo and Mart{\'{i}}nez, F{\'{e}}lix and Morales-reyes, Zebensui and S{\'{a}}nchez-Zapata, Jos{\'{e}} A.},
journal = {Ecological Applications},
title = {{Long-term demographic dynamics of a keystone scavenger disrupted by human-induced shifts in food availability}},
year = {2021}
}
GitHub Events
Total
- Push event: 9
- Fork event: 1
Last Year
- Push event: 9
- Fork event: 1