Performance and optimal design of N-mixture models for spatiotemporally replicated drone-based surveys
Repository for the manuscript OPTIMALLY DESIGNING DRONE SURVEYS FOR
WILDLIFE ABUNDANCE MODELING WITH N-MIXTURE MODELS (Brack et al. 2022;
Methods in Ecology and Evolution).
This repository contains the simulation results and R
code to explore,
in a very wide scan study, the performance (based on root mean squared
error) and optimal survey effort allocation for hierarchical N-mixture
models, focusing on their application for drone-based surveys. We also
investigate the use of a double-observer protocol in image reviewing to
decompose the detection process in availability and perception.
We provide an interactive and “ready-to-consult”
webpage hoping
to assist ecologists and conservationists in panning spatiotemporally
replicated drone-based surveys for abundance modeling. There we show in
detail simulation results and run some examples.
The simulation study is divided in three parts:
-
This first experiment assesses the performance of N-mixture models using double and single observer counts under different scenarios of population density and detection probability and address the optimal survey effort allocation in terms of spatial vs. temporal prioritization for each scenario.
-
In this experiment, we investigate how the use of double-observer protocol increases model performance and affects optimal survey effort allocation.
-
Here, we evaluate if the use of double-observer protocol can reduce the effort needed in fieldwork to achieve the same model performance as in a single observer approach.
data/
: simulation results filesraw_simul_resu/
: raw simulation results are not available in this repo (download info below)processed_simul_resu/
: extracted and cleaned simulation results
docs/
: files for the websitems/
: manuscript filesoutputs/
: figures, tables and intermediate outputsfigs/
: figures built with R scripts. RMSE curves for scenarios and RMSE relationshipstabs/
: tables built with R scripts. Optimal J tables and budget savings table
R/
: R scripts and functions. Script file names reference to each simulation experiment withPart
orP
ex.sumulation*.R
: script to run a single iteration simulating spatiotemporally replicated counts and analyzing with the respective N-mixture model. Examples with maximum likelihood estimation (*_likeli.R
) and Bayesian approach (*_bayes.R
).func_*.R
: functions to support R scriptssimulNmix_scenarios_parallel.R
: script to run several simulation iterations. This script was used to produce raw simulation results using maximum likelihood estimation.Part*/
: code for each simulation experimentC&E_.R*
: clean and extract raw simulation results and save in the processed foldercalcRMSE_*.R
: code to calculate the RMSE for each scenario (under a combination of local abundance and availability) and get the optimal number of visits for each scenario.fig*.R
: code to generate figures
download.file(url = "https://github.com/ismaelvbrack/designNmix4droneSurveys/archive/main.zip", destfile = "designNmix4droneSurveys.zip")
To reproduce the results from raw simulation results, download them from
zenodo (645 Mb).
Otherwise, reproduce starting from processed simulation results (cleaned
and extracted, available here) directly from this repo.
Manual download here and
paste in ~/data/raw_simul_resu
or
Download from RStudio:
library(zen4R)
download_zenodo(doi="10.5281/zenodo.5156592", path=here::here("data","raw_simul_resu"))