From 126658418f871ec06f1d17ecdd8adda764ab0f0e Mon Sep 17 00:00:00 2001
From: Maria Voigt <maria.voigt@idiv.de>
Date: Wed, 8 Feb 2017 16:24:02 +0100
Subject: [PATCH] replacing aerial index with value for ESW
---
src/model_fitting/abundance_model.R | 14 ++++++--------
1 file changed, 6 insertions(+), 8 deletions(-)
diff --git a/src/model_fitting/abundance_model.R b/src/model_fitting/abundance_model.R
index fbb54a1..27e94d2 100644
--- a/src/model_fitting/abundance_model.R
+++ b/src/model_fitting/abundance_model.R
@@ -30,10 +30,10 @@ option_list <- list (
make_option(c("-o", "--output-directory"), dest = "output_directory",
type = "character", help = "directory with output files",
metavar = "/path/to/output-dir"),
- make_option("--ESW-aerial-index",
- dest = "ESW_aerial_index",
+ make_option("--ESW-aerial",
+ dest = "ESW_aerial",
type = "double",
- help = "index aerial effective strip width"),
+ help = "aerial effective strip width"),
make_option("--exclude-year", dest = "exclude_year", type = "integer",
default = NA, help = "year to exclude", metavar = "2015"),
make_option("--include-aerial",
@@ -78,8 +78,8 @@ if(is_verbose){print(paste("indir", indir))}
outdir <- options$output_directory
if(is_verbose){print(paste("outdir", outdir))}
-ESW_aerial_index <- options$ESW_aerial_index
-if(is_verbose){print(paste("Aerial ESW index", ESW_aerial_index))}
+ESW_aerial <- options$ESW_aerial
+if(is_verbose){print(paste("Aerial ESW", ESW_aerial))}
do_stability <- options$stability
if(is_verbose){print(paste("stability", do_stability))}
@@ -109,9 +109,7 @@ source(file.path(indir_fun, "roger_functions/get_conf_set.r"))
#define offset ground
ESW <- 0.01595 #effective strip width in km
-
-ESW_aerial_options <- seq(from = 0.01, to = 0.15, by = 0.01)
-ESW_aerial <- ESW_aerial_options[ESW_aerial_index]
+# ESW_aerial already defined
print(paste("this is ESW aerial:", ESW_aerial))
NCS <- 1.12 #nest construction rate from Spehar et al. 2010
PNB <- 0.88 # proportion of nest builders from Spehar et al. 2010
--
GitLab